bzoj 4556: [Tjoi2016&Heoi2016]字符串

Posted 2020-08-27 SD_le

      二分ans，二分区间长度，st表查，最后主席树判断。

      这题最大收获学到了一个nb的卡常技巧，主席树元素个数为0直接返回，不敢相信快了一倍。

    

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<algorithm>
  4 #include<cstring>
  5 #define N 200005
  6 using namespace std;
  7 inline int read()
  8 {
  9     int p=0;char c=getchar();
 10     while(c<‘0‘||c>‘9‘)c=getchar();
 11     while(c>=‘0‘&&c<=‘9‘)p=p*10+c-‘0‘,c=getchar();
 12     return p;
 13 }
 14 int wb[N],rank[N],sa[N],sum[N];
 15 int mnn(int x,int y)
 16 {
 17     if(x<y)return x;return y;
 18 }
 19 char s[N];
 20 void da(int n,int m)
 21 {
 22     int *x=rank,*y=wb;
 23     for(int i=0;i<m;i++)sum[i]=0;
 24     for(int i=0;i<n;i++)sum[x[i]=s[i]]++;
 25     for(int i=1;i<m;i++)sum[i]+=sum[i-1];
 26     for(int i=n-1;i>=0;i--)sa[--sum[x[i]]]=i;
 27     int p=1;
 28     for(int j=1;p<n;j<<=1,m=p)
 29     {
 30         p=0;
 31         for(int i=n-j;i<n;i++)y[p++]=i;
 32         for(int i=0;i<n;i++)if(sa[i]>=j)y[p++]=sa[i]-j;
 33         for(int i=0;i<m;i++)sum[i]=0;
 34         for(int i=0;i<n;i++)sum[x[i]]++;
 35         for(int i=1;i<m;i++)sum[i]+=sum[i-1];
 36         for(int i=n-1;i>=0;i--)sa[--sum[x[y[i]]]]=y[i];
 37         swap(x,y);x[sa[0]]=0;p=1;
 38         for(int i=1;i<n;i++)
 39          x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i-1]+j]==y[sa[i]+j]?p-1:p++;
 40     }
 41     return ;
 42 }
 43 int h[N];
 44 void cal(int n)
 45 {
 46     int k=0;
 47     for(int i=1;i<=n;i++)rank[sa[i]]=i;
 48     for(int i=0;i<n;i++)
 49     {
 50         if(k)k--;
 51         int j=sa[rank[i]-1];
 52         while(s[j+k]==s[i+k])k++;
 53         h[rank[i]]=k;
 54     }
 55 }
 56 int n,m;
 57 int mn[N][20];
 58 int lg[N];
 59 int root[N];
 60 void st()
 61 {
 62     int now=0;
 63     for(int i=1;i<=n;i<<=1)
 64     {
 65         for(int j=i;j<i<<1;j++)
 66         {
 67             lg[j]=now;
 68         }
 69         now++;
 70     }
 71     for(int i=1;i<=n;i++)mn[i][0]=h[i];
 72     for(int i=1;i<=18;i++)
 73     {
 74         int k=(1<<(i-1));
 75         for(int j=1;j<=n;j++)
 76         {
 77             if(j+k<=n)mn[j][i]=mnn(mn[j][i-1],mn[j+k][i-1]);
 78             else mn[j][i]=mn[j][i-1];
 79         }
 80     }return ;
 81 }
 82 struct node
 83 {
 84     int l,r,sum;
 85 }a[N*20];int cnt;
 86 void merge(int x,int y,int l,int r,int pos)
 87 {
 88     if(l==r)
 89     {
 90         a[x].sum=1;
 91         return ;
 92     }
 93     int mid=(l+r)>>1;
 94     if(pos<=mid)
 95     {
 96         a[x].l=++cnt;
 97         a[x].r=a[y].r;
 98         merge(a[x].l,a[y].l,l,mid,pos);
 99     }
100     else
101     {
102         a[x].r=++cnt;
103         a[x].l=a[y].l;
104         merge(a[x].r,a[y].r,mid+1,r,pos);
105     }
106     a[x].sum=a[a[x].l].sum+a[a[x].r].sum;
107 }
108 bool qur(int x,int y,int l,int r,int ll,int rr)
109 {
110     if(a[x].sum==a[y].sum)return 0;
111     if(l>=ll&&r<=rr)return a[x].sum!=a[y].sum;
112     int mid=(l+r)>>1;
113     if(ll<=mid)if(qur(a[x].l,a[y].l,l,mid,ll,rr))return 1;
114     if(rr>mid)if(qur(a[x].r,a[y].r,mid+1,r,ll,rr))return 1;
115     return 0;
116 }
117 int t1,t2,t3,t4;
118 inline int qr(int l,int r)
119 {
120     if(l>r)return N;
121     int k=lg[r-l+1];
122     return mnn(mn[l][k],mn[r-(1<<k)+1][k]);
123 }
124 bool pan(int x)
125 {
126      int pos=rank[t3-1];
127      int ll,rr;
128      int l=pos,r=n,mid;
129      while(l<=r)
130      {
131          mid=(l+r)>>1;
132          if(qr(pos+1,mid)<x)r=mid-1;
133          else l=mid+1;
134      }
135      rr=r;
136      l=1,r=pos;
137      while(l<=r)
138      {
139          mid=(l+r)>>1;
140          if(qr(mid+1,pos)<x)l=mid+1;
141          else r=mid-1;
142      }
143      ll=l;
144      return qur(root[rr],root[ll-1],1,n,t1,t2-x+1);
145 }
146 
147 int main()
148 {
149     scanf("%d%d",&n,&m);
150     scanf("%s",s);
151     da(n+1,256);cal(n);
152     st();
153     for(int i=1;i<=n;i++)
154     {
155         root[i]=++cnt;
156         merge(root[i],root[i-1],1,n,sa[i]+1);
157     }
158     for(int i=1;i<=m;i++)
159     {
160         t1=read();t2=read();t3=read();t4=read();
161         int l=0;int r=min(t2-t1+1,t4-t3+1);
162         while(l<=r)
163         {
164             int mid=(l+r)>>1;
165             if(pan(mid))l=mid+1;
166             else r=mid-1;
167         }
168         printf("%d\n",r);
169     }
170     return 0;
171 }
172 /*
173 5 2
174 ababaa
175 1 4 2 5
176 1 4 2 3
177 */