(本人本题完成于2016-7-22)
题目:子串-题目
做法:设f[i][j][k]为在A串前i个字符中分k份取出B串前j个字符的方案数,再设一个数组g[i][j][k]为在A串前i个字符中分k份取出B串前j个字符,其中必取A串的第i个字符的方案数,由此得到状态转移方程:
g[i][j][k]=f[i-1][j-1][k-1]+g[i-1][j-1][k] (A[i]=B[j])
g[i][j][k]=0 (A[i]≠B[j])
f[i][j][k]=f[i-1][j][k]+g[i][j][k]
由于数据较大,所以我们可以以i划分阶段,用滚动数组存储,可以大大节省空间。
以下是本人代码:
#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define mod 1000000007
using namespace std;
long n,m,K,f[2][210][210]={0},g[2][210][210]={0},now,past; //使用滚动数组,节省空间
char a[1010],b[210]; //存储A串和B串
int main()
{
scanf("%ld %ld %ld\n",&n,&m,&K);
for(int i=1;i<=n;i++) scanf("%c",&a[i]);
scanf("\n");
for(int i=1;i<=m;i++) scanf("%c",&b[i]);
f[0][0][0]=1; //初始化
now=1;past=0;
for(int i=1;i<=n;i++)
{
f[now][0][0]=1;
for(int j=1;j<=m;j++)
for(int k=1;k<=K;k++)
{
if (a[i]==b[j]) g[now][j][k]=(f[past][j-1][k-1]+g[past][j-1][k])%mod;
else g[now][j][k]=0;
f[now][j][k]=(f[past][j][k]+g[now][j][k])%mod;
}
int t=now;now=past;past=t;
}
printf("%ld",f[past][m][K]%mod);
return 0;
}