#include<bits/stdc++.h> using namespace std; const int mod=1e9+7; const int N = 1000005; int Next[N]; long long ex[N],tong[N]; //即extand[] char p[N],t[N]; int T; long long ans; void pre() // next[i]: 以第i位置开始的子串 与 T的公共前缀 { int lp=strlen(p); Next[0]=lp; int j=0,k=1; while(j+1<lp && p[j]==p[j+1]) j++; Next[1]=j; for(int i=2; i<lp; i++) { int P=Next[k]+k-1; int L=Next[i-k]; if(i+L<P+1) Next[i]=L; else { j=max(0,P-i+1); while(i+j<lp && p[i+j]==p[j]) j++; // 枚举(p+1,length) 与(p-k+1,length) 区间比较 Next[i]=j; k=i; } } } void exkmp() { int lp=strlen(p),lt=strlen(t); pre(); //next数组初始化 int j=0,k=0; while(j<lt && j<lp && p[j]==t[j]) j++; ex[0]=j; for(int i=1; i<lt; i++) { int P=ex[k]+k-1; int L=Next[i-k]; if(i+L<P+1) ex[i]=L; else { j=max(0,P-i+1); while(i+j<lt && j<lp && t[i+j]==p[j]) j++; ex[i]=j; k=i; } } } int main() { scanf("%d",&T); while(T--) { memset(tong,0,sizeof(tong)); scanf("%s%s",&t,&p); int lt=strlen(t); int lp=strlen(p); reverse(p,p+lp); reverse(t,t+lt); exkmp(); ans=0; for(int i=0; i<lt; i++) tong[ex[i]]++; for(int i=lp;i;i--) { tong[i]=(tong[i]+tong[i+1])%mod; ans=(ans+tong[i]*i%mod)%mod; } printf("%lld\n",ans); } }
原文:https://www.cnblogs.com/SDUTNING/p/10388792.html