首页 诗词 字典 板报 句子 名言 友答 励志 学校 网站地图
当前位置: 首页 > 教程频道 > 开发语言 > 编程 >

zoj 3587 蛮美的KMP题

2012-08-25 
zoj 3587 蛮好的KMP题我以前做过的KMP http://www.cnblogs.com/wuyiqi/tag/KMP/题意:从S串中拿出两个子串

zoj 3587 蛮好的KMP题

我以前做过的KMP http://www.cnblogs.com/wuyiqi/tag/KMP/

题意:从S串中拿出两个子串能拼成T串的方案有多少种,两个子串可以重复

这道题其实要解决的问题就是求T串的某个前缀在S串中出现了几次(cnt[])    以及T串的某个后缀在S串中出现了几次(cnt2[]),那么答案就是cnt[i]*cnt2[len-i]的和

具体求的过程就很简单了,cnt[next[i]]+=cnt[i]即可,最后的cnt[i]就表示了i前缀出现的次数,后缀的话将字符串反转再KMP即可

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;typedef long long lld;const int maxn = 100001;int n,m;char a[maxn],b[maxn];int p[maxn],cnt[maxn],cnt2[maxn];void getp(){    p[1]=0;    int i,j=0;    for(i=2;i<=m;i++){        while(j>0&&b[j+1]!=b[i]) j=p[j];        if(b[j+1]==b[i]) j+=1;        p[i]=j;    }}void kmp(int cnt[]){    int i,j=0;    for(i=1;i<=n;i++){        while(j>0&&b[j+1]!=a[i]) j=p[j];        if(b[j+1]==a[i]) j+=1,cnt[j]++;        if(j==m)    j=p[j];    }for(i=m;i>=1;i--){if(p[i]!=0) cnt[p[i]]+=cnt[i];}}void RE(char s[]){int len=strlen(s+1);for(int i=1;i<=len/2;i++){swap(s[i],s[len-i+1]);}}lld solve(){lld ans=0;for(int i=1;i<=m;i++)ans+=(lld)cnt[i]*cnt2[m-i];return ans;}int main(){int t;scanf("%d",&t);while(t--)    {scanf("%s",a+1);scanf("%s",b+1);        m=strlen(b+1);        n=strlen(a+1);memset(cnt,0,sizeof(cnt));        getp();        kmp(cnt);RE(a);RE(b);memset(cnt2,0,sizeof(cnt2));getp();kmp(cnt2);printf("%lld\n",solve());    }}


热点排行