學習kmp算法我最後是看的數據結構書上的一本教材學會的。。我認爲kmp相對於普通的BF算法就是避免了很是多沒必要要的匹配。而kmp算法的精髓天然就在於next數組的運用。。。而next數組簡而言之就是存儲的就是模式串中第j個字符與主串中對應字符「失配」時,在模式串中需要又一次和主串中失配的字符相比較的位置。php
。ios
。我認爲這句歸納挺好的。算法
。數組
。數據結構
題1:post
hdu 1711 number sequence 學習
題目連接:spa
http://acm.hdu.edu.cn/showproblem.php?pid=1711
code
#include<cstdio> #include<cstring> const int maxn=1000000+10; int s[maxn],p[maxn]; int next[maxn]; int n,m; /*算法複雜度爲o(m+n)*/ void Getnext() { int j,k; next[0]=-1; j=0; k=-1; /*n=strlen(p)*/ while(j<m) { if(k==-1||p[j]==p[k]) //匹配的狀況下,p[j]==p[k] { j++; k++; next[j]=k; } else //p[j]!=p[k] k=next[k]; } } int KMP() { int i,j; i=0,j=0; Getnext(); while(i<n) { if(j==-1||s[i]==p[j]) { i++; j++; } else j=next[j]; if(j==m) return i-m+1; } return -1; } int main() { int t; scanf("%d",&t); while(t--) { scanf("%d%d",&n,&m); for(int i=0;i<n;i++) scanf("%d",&s[i]); for(int i=0;i<m;i++) scanf("%d",&p[i]); int ans=KMP(); printf("%d\n",ans); } return 0; }
#include<cstdio> #include<cstring> #include<iostream> using namespace std; const int maxn=1000000+10; char s[maxn],p[maxn]; int next[maxn]; int n,m; /*算法複雜度爲o(m+n)*/ void Getnext() { int j,k; next[1]=0; j=1; k=0; while(j<=m) { if(k==0||p[j]==p[k]) //匹配的狀況下,p[j]==p[k] { j++; k++; next[j]=k; } else //p[j]!=p[k] k=next[k]; } } int KMP() { int i,j; i=1,j=1; Getnext(); while(i<n) { if(j==0||s[i]==p[j]) { i++; j++; } else j=next[j]; if(j==m+1) return i-m; } return -1; } int KMP_Count() { int ans=0; int i,j=1; if(n==1&&m==1) { if(s[1]==p[1]) return 1; else return 0; } Getnext(); for(i=1;i<=n;i++) { while(j>1&&s[i]!=p[j]) j=next[j]; if(s[i]==p[j]) j++; if(j==m+1) { ans++; j=next[j]; } } return ans; } int main() { int t; scanf("%d",&t); while(t--) { scanf("%s",p+1); scanf("%s",s+1); m=strlen(p+1); n=strlen(s+1); int ans2=KMP_Count(); printf("%d\n",ans2); } return 0; } /* ababcabcacbab abcac abababc ab */
。求週期。。。
。
。
。
<h2><pre name="code" class="cpp">#include<cstdio> #include<cstring> const int maxn=1000000+10; char s[maxn],p[maxn]; int next[maxn]; int m; /*算法複雜度爲o(m+n)*/ void Getnext() { int j,k; next[1]=0; j=1; k=0; while(j<=m) { if(k==0||p[j]==p[k]) //匹配的狀況下,p[j]==p[k] { j++; k++; next[j]=k; } else //p[j]!=p[k] k=next[k]; } } int KMP() { int i,j; i=0,j=0; Getnext(); while(i<strlen(s)) { if(j==-1||s[i]==p[j]) { i++; j++; } else j=next[j]; if(j==strlen(p)) return i-strlen(p); } return -1; } int main() { int cas=1; while(~scanf("%d",&m),m) { scanf("%s",p+1); Getnext(); printf("Test case #%d\n",cas++); for(int i=2;i<=m+1;i++) { if(next[i]>1&&(i-1)%(i-next[i])==0) printf("%d %d\n",i-1,(i-1)/(i-next[i])); } printf("\n"); } return 0; } /* 7 abcabcc */
。
因此跟上題同樣。。僅僅需要推斷一下末尾就可以了。
。。
即n%(n+1-next[n+1])==0成立與否。。
。
#include<cstdio> #include<cstring> const int maxn=1000000+10; char p[maxn]; int next[maxn],n; void Getnext() { int j,k; next[1]=0; k=0; j=1; while(j<=n) { if(k==0||p[k]==p[j]) { k++; j++; next[j]=k; } else k=next[k]; } } int main() { int ans; while(~scanf("%s",p+1)) { if(strcmp(p+1,".")==0) return 0; ans=-1; n=strlen(p+1); Getnext(); if(n%(n+1-next[n+1])==0) ans=n/(n+1-next[n+1]); if(ans==-1) ans=1; printf("%d\n",ans); } return 0; }
。之後作的題我也會補充。。。