noip2015day2 子串

对于dp蒟蒻的我，这道题把我难倒了。翻看了网上的各种题解，现自己总结一份如下：
我们用sum[i][j][k]为a串前i个中选出k个字串拼成b串前j个的方案数，f[i][j][k]是a串第i个必须用的方案数。

当且仅当a[i] == b[j]，f[i][j][k]=sum[i-1][j-1][k-1]+f[i-1][j-1][k]

f[i-1][j-1][k]代表第i个用于接上上一个字串
sum[i-1][j-1][k-1]代表第i个用于新开一个字串
显然f[i][j][k]=g[i][j][k]+f[i-1][j][k]，因为此时有两个决策，要么用第i个，要么不用第i个

f 等于 a 中选到第 i-1 个和 b 中选到第 j-1 个用 k-1 段的总和新开一段加上 a 中选到第 i-1 个和 b 中选到第 j-1 个接上第i个。
若不等于则第i个一定不用，所以此时f[i][j][k]=0

#include<cstdio>
#include<iostream>
using namespace std;
int f[2][201][201],sum[2][201][201]; 
char a[1001],b[201];
const int mod=1000000007;
int main()
{
    int lena,lenb,K;
    cin>>lena>>lenb>>K;
    for(int i=1;i<=lena;i++)
        cin>>a[i];
    for(int i=1;i<=lenb;i++)
        cin>>b[i];
        int l,s;
    int now=0,pre=1;
    sum[0][0][0]=1;
    for(int i=1;i<=lena;i++)
    {
        swap(now,pre);
        sum[now][0][0]=1;
        for(int j=1;j<=lenb;j++)
        {
            for(int k=1;k<=K;k++)
            {
                if(a[i]==b[j])
                    f[now][j][k]=(sum[pre][j-1][k-1]+f[pre][j-1][k])%mod;

                else f[now][j][k]=0;
                sum[now][j][k]=(sum[pre][j][k]+f[now][j][k])%mod;
            }
        }
    }
    printf("%d",sum[now][lenb][K]);
}