KMP--君住长江头，我住长江尾，日日思君不见君，共饮长江水

POJ 3461: Oulipo

题意：

求出第一个串在第二个串中的出现次数...

分析：

KMP板子题...

代码：

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 //by NeighThorn
 using namespace std;
 //眉眼如初，岁月如故

 const int maxn=+;

 int cas,lens,lenp,nxt[maxn];

 char s[maxn],p[maxn];

 inline void getnxt(void){
     nxt[]=nxt[]=;int k;
     for(int i=;i<lenp;i++){
         k=nxt[i];
         while(k&&p[k+]!=p[i+])
             k=nxt[k];
         if(p[k+]==p[i+])
             nxt[i+]=k+;
         else
             nxt[i+]=;
     }
 }

 inline void kmp(void){
     int ans=,posp=,poss=;
     while(poss<=lens){
         if(p[posp]==s[poss])
             posp++,poss++;
         else if(posp==)
             poss++;
         else
             posp=nxt[posp-]+;
         if(posp==lenp+)
             ans++,posp=nxt[posp-]+;
     }
     printf("%d\n",ans);
 }

 signed main(void){
     scanf("%d",&cas);
     while(cas--){
         memset(nxt,,sizeof(nxt));
         scanf("%s%s",p+,s+);
         lenp=strlen(p+);lens=strlen(s+);
         getnxt();kmp();
     }
     return ;
 }//Cap ou pas cap. Pas cap.

POJ 2406: Power Strings

题意：

求出每个串的最短循环节出现次数...

分析：

此题两种解法...然而本质上是一样的...

No.1 Hash

对于一个串其前缀A和后缀B相同并且有重叠，并且len%strlen(str-B)==0，那么str-B一定是循环节...所以我们可以用hash水过...

No.2 next数组...

然而更机智的做法是next数组...

代码：

No.1 Hash

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #define int long long
 using namespace std;
 const int maxn=+,MOD=;
 char str[maxn];
 int len,ans;
 unsigned int hash[maxn],pow[maxn];
 signed main(void){
     pow[]=(long long);
     for(int i=;i<maxn;i++)
         pow[i]=(pow[i-]*)%MOD;
     while(scanf("%s",str+)&&str[]!='.'){
         len=strlen(str+),hash[]=,ans=;
         for(int i=;i<=len;i++)
             hash[i]=(hash[i-]*%MOD+str[i]-'a'+)%MOD;
         for(int i=;i<len;i++){
             if(len%(len-i)==){
                 if(hash[i]%MOD==(hash[len]-hash[len-i]*pow[i]%MOD+MOD)%MOD){
                     ans=len/(len-i);
                 }
             }
         }
         cout<<ans<<endl;
     }
     return ;
 }

No.2 next数组

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 //by NeighThorn
 using namespace std;
 //眉眼如初，岁月如故

 const int maxn=+;

 int len,nxt[maxn];

 char str[maxn];

 inline void getnxt(void){
     nxt[]=nxt[]=;int k;
     for(int i=;i<len;i++){
         k=nxt[i];
         while(k&&str[k+]!=str[i+])
             k=nxt[k];
         if(str[k+]==str[i+])
             nxt[i+]=k+;
         else
             nxt[i+]=;
     }
 }

 signed main(void){
     while(scanf("%s",str+)&&str[]!='.'){
         len=strlen(str+);getnxt();
         if(len%(len-nxt[len])==)
             printf("%d\n",len/(len-nxt[len]));
         else
             puts("");
     }
     return ;
 }//Cap ou pas cap. Pas cap.

POJ 2752: Seek the Name, Seek the Fame

题意：

输出所有s的合法前缀，合法前缀的定义为前缀等于后缀...

分析：

不断nxt...

代码：

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 //by NeighThorn
 using namespace std;
 //眉眼如初，岁月如故 

 const int maxn=+;

 int len,tail,nxt[maxn],stk[maxn];

 char str[maxn];

 inline void getnxt(void){
     nxt[]=nxt[]=;int k;
     for(int i=;i<len;i++){
         k=nxt[i];
         while(k&&str[k+]!=str[i+])
             k=nxt[k];
         if(str[k+]==str[i+])
             nxt[i+]=k+;
         else
             nxt[i+]=;
     }
 }

 signed main(void){
     while(scanf("%s",str+)!=EOF){
         len=strlen(str+);getnxt();tail=;stk[++tail]=len;
         while(nxt[len])
             stk[++tail]=nxt[len],len=nxt[len];
         for(int i=tail;i>=;i--)
             printf("%d ",stk[i]);
         puts("");
     }
     return ;
 }//Cap ou pas cap. Pas cap.

POJ 3167: Cow Patterns

题意：

给出两个串，如果a和b匹配当且仅当ab中的每个对应位置的数字rank1和rank2相同，rank1代表当前元素前面小于他的元素个数，rank2代表当前元素前面等于他的元素个数，求出匹配串和模式串的匹配位置（开头位置）

分析：

KMP…裸的KMP…就是计算一下rank就好… 
因为数字不会超过25，所以直接暴力求rank就好…

代码：

 #include<algorithm>
 #include<iostream>
 #include<cstring>
 #include<cstdio>
 //by NeighThorn
 using namespace std;
 const int maxn=+,maxk=+,maxs=+;
 int n,k,s,a[][maxn],nxt[maxk],num[][maxn][maxs],vis[maxk];
 inline int rank1(int len,int en,int id){
     int st=en-len,ans=;
     for(int i=;i<a[id][en];i++)
         ans+=num[id][en][i]-num[id][st][i];
     return ans;
 }
 inline int rank2(int len,int en,int id){
     return num[id][en][a[id][en]]-num[id][en-len][a[id][en]];
 }
 inline void getnext(void){
     nxt[]=nxt[]=;int p;
     for(int i=;i<k;i++){
         p=nxt[i];int a=rank1(p+,p+,),b=rank1(p+,i+,),c=rank2(p+,p+,),d=rank2(p+,i+,);
         while(p&&(a!=b||c!=d))
             p=nxt[p],a=rank1(p+,p+,),b=rank1(p+,i+,),c=rank2(p+,p+,),d=rank2(p+,i+,);
         if(rank1(p+,p+,)==rank1(p+,i+,)&&rank2(p+,p+,)==rank2(p+,i+,))
             nxt[i+]=p+;
         else
             nxt[i+]=;
     }
 }
 inline void kmp(void){
     int posa=,posb=,ans=,stk[maxn];
     while(posa<=n){
         if(rank1(posb,posa,)==rank1(posb,posb,)&&rank2(posb,posa,)==rank2(posb,posb,))
             posa++,posb++;
         else if(posb==)
             posa++;
         else
             posb=nxt[posb-]+;
         if(posb==k+)
             stk[++ans]=posa-k,posb=nxt[posb-]+;
     }
     cout<<ans<<endl;
     for(int i=;i<=ans;i++)
         cout<<stk[i]<<endl;
 }
 signed main(void){
     scanf("%d%d%d",&n,&k,&s);memset(num,,sizeof(num));
     for(int i=;i<=n;i++)
         scanf("%d",&a[][i]),num[][i][a[][i]]=;
     memset(vis,,sizeof(vis));
     for(int i=;i<=n;i++)
         for(int j=;j<=s;j++)
             num[][i][j]+=num[][i-][j];
     for(int i=;i<=k;i++)
         scanf("%d",&a[][i]),num[][i][a[][i]]=;
     for(int i=;i<=k;i++)
         for(int j=;j<=s;j++)
             num[][i][j]+=num[][i-][j];
     getnext();kmp();
     return ;
 }

---

By NeighThorn