BZOJ2099: [Usaco2010 Dec]Letter 恐吓信

给两个长度不超过50000的串，A串可每次截连续一段复制出来，求最少复制几次能得到B串。

方法一：SAM。不会。

嗯好会了。

 #include<stdio.h>
 #include<string.h>
 #include<algorithm>
 #include<stdlib.h>
 //#include<iostream>
 //#include<time.h>
 using namespace std;

 int n,m;
 #define maxn 200011*2
 char s[maxn],p[maxn];
 struct SAM
 {
     int ch[maxn][],size,last,pre[maxn],pos[maxn];
     SAM()
     {
         size=;
         memset(ch[],,sizeof(ch[]));
         pre[]=-;
     }
     int idx(char c) {return c-'A';}
     void insert(char c,int p)
     {
         int id=idx(c),x=++size;
         memset(ch[x],,sizeof(ch[x])); pos[x]=p;
         int y=last;
         for (;y!=- && !ch[y][id];y=pre[y]) ch[y][id]=x;
         last=x;
         if (y==-) pre[x]=;
         else
         {
             if (pos[ch[y][id]]==pos[y]+) pre[x]=ch[y][id];
             else
             {
                 int z=ch[y][id],w=++size;
                 memcpy(ch[w],ch[z],sizeof(ch[z]));
                 pos[w]=pos[y]+; pre[w]=pre[z];
                 pre[z]=pre[x]=w;
                 for (;y!=- && ch[y][id]==z;y=pre[y]) ch[y][id]=w;
             }
         }
     }
 }sam;

 int main()
 {
     scanf("%d%d",&n,&m);
     for (int i=;i<=n;i++)
     {
         char c;while (!((c=getchar())>='A' && c<='Z'));
         sam.insert(c,i);
     }
     for (int i=;i<=m;i++) while (!((p[i]=getchar())>='A' && p[i]<='Z'));
     int pos=,ans=;
     while (pos<=m)
     {
         int now=;
         while (pos<=m && sam.ch[now][p[pos]-'A']) now=sam.ch[now][p[pos]-'A'],pos++;
 //        cout<<endl;
         ans++;
     }
     printf("%d\n",ans);
     return ;
 }

方法二：其实就是拿B去A上面匹配若干次。多次匹配的话可以把A先建个后缀数组，然后每次二分到B的位置搞匹配。匹配一次怕太慢的话，用hash吧！二分匹配长度找到第一个不可匹配的位置，就可以算匹配长度啦！

 #include<stdio.h>
 #include<string.h>
 #include<stdlib.h>
 #include<algorithm>
 //#include<iostream>
 using namespace std;

 int n,m;
 #define maxn 50011
 char t[maxn],p[maxn],c;
 int ht[maxn],hp[maxn];
 int sa[maxn],cnt[maxn],rank[maxn],y[maxn];
 void makesa(char* s)
 {
     int n=strlen(s),m=;
     for (int i=;i<m;i++) cnt[i]=;
     for (int i=;i<n;i++) cnt[rank[i]=s[i]-'A']++;
     for (int i=;i<m;i++) cnt[i]+=cnt[i-];
     for (int i=n-;i>=;i--) sa[--cnt[rank[i]]]=i;
     for (int k=;k<=n;k<<=)
     {
         int p=;
         for (int i=n-k;i<n;i++) y[p++]=i;
         for (int i=;i<n;i++) if (sa[i]>=k) y[p++]=sa[i]-k;
         for (int i=;i<m;i++) cnt[i]=;
         for (int i=;i<n;i++) cnt[rank[y[i]]]++;
         for (int i=;i<m;i++) cnt[i]+=cnt[i-];
         for (int i=n-;i>=;i--) sa[--cnt[rank[y[i]]]]=y[i];
         memcpy(y,rank,sizeof(y));
         p=;rank[sa[]]=;
         for (int i=;i<n;i++)
             rank[sa[i]]=y[sa[i]]==y[sa[i-]] && ((sa[i]+k>=n && sa[i-]+k>=n)
             || (sa[i]+k<n && sa[i-]+k<n && y[sa[i]+k]==y[sa[i-]+k]))?p:++p;
         if (p==n-) break;
         m=p+;
     }
 }
 void makeh(char* s,int* h)
 {
     int n=strlen(s);
     h[]=s[]-'A';
     for (int i=;i<n;i++)
         h[i]=h[i-]*+s[i]-'A';
 }
 int pw[maxn];
 void makepw(int n)
 {
     pw[]=;
     for (int i=;i<=n;i++)
         pw[i]=pw[i-]*;
 }
 bool ok(int x,int y)
 {
     int L=,R=min(n-x,m-y);
     while (L<R)
     {
         if (ht[x+mid-]-ht[x-]*pw[mid]==hp[y+mid-]-hp[y-]*pw[mid]) L=mid;
         else R=mid-;
     }
     return p[y+L]<=t[x+L];
 }
 int lcp(int x,int y)
 {
     int now=;
     while (x<n && y<m && t[x]==p[y]) now++,x++,y++;
     return now;
 }
 void play(int &x)
 {
     int L=,R=n;
     while (L<R)
     {
         int mid=(L+R)>>;
         if (ok(sa[mid],x)) R=mid;
         else L=mid+;
     }
     if (!L) x+=lcp(sa[L],x);
     else if (L<n) x+=max(lcp(sa[L-],x),lcp(sa[L],x));
     else x+=lcp(sa[n-],x);
 }
 bool isalpha(char c) {return c>='A' && c<='Z';}
 int main()
 {
     scanf("%d%d",&n,&m);
     t[n]='\0';p[m]='\0';
     for (int i=;i<n;i++)
     {
         while (!isalpha(c=getchar()));
         t[i]=c;
     }
     for (int i=;i<m;i++)
     {
         while (!isalpha(c=getchar()));
         p[i]=c;
     }
     makesa(t);
     makeh(t,ht);makeh(p,hp);makepw(max(n,m));
     int now=,ans=;
     while (now<m) play(now),ans++;
     printf("%d\n",ans);
     return ;
 }