NOIP 2015 D2 T2 子串 substring

来源:互联网 发布:dream it possible 编辑:程序博客网 时间:2024/05/16 00:46


:
ABAKkA使B

f(i,j,k)Aibjk

f(i,j,k)=f(i1,j1,k),A[i]=B[j]A[i1]B[j1]

f(i,j,k)=f(i1,j1,k)+f(i,j,k1),A[i]=B[j]A[i1]=B[j1]

#include <stdio.h>#include <string.h>#include <iostream>#include <algorithm>using namespace std;#define mod 1000000007llchar a[1000+5],b[200+5];int n,m,K;long long f[2][1001][201];long long tmp[2][1001][201];int main(){    cin>>n>>m>>K;    scanf("%s%s",a+1,b+1);    long long ans = 0;    f[0][0][0]=tmp[0][0][0]=1;    for(int i=0;i<=n;++i)tmp[0][i][0]=1;    for(int k=1;k<=K;++k)    {        register int o= k&1,op = (k+1)&1;        memset(f[o],0,sizeof f[o]);        memset(tmp[o],0,sizeof tmp[o]);        for(register int i=1;i<=n;++i)            for(register int j=1;j<=m;++j)            {                if(a[i]==b[j])                {                    f[o][i][j]=tmp[op][i-1][j-1];                     if(a[i-1]==b[j-1])                         (f[o][i][j]+=f[o][i-1][j-1])%=mod;                }                (tmp[o][i][j] +=  f[o][i][j] + tmp[o][i-1][j])%=mod;            }    }    for(int i=1;i<=n;++i)        (ans += f[K&1][i][m])%=mod;    cout<<ans<<endl;}

这里写图片描述

1 0
原创粉丝点击