NOIP2015 子串 【动态规划】

题目描述

有两个仅包含小写英文字母的字符串 A 和 B 。现在要从字符串 A 中取出 k 个互不重叠的非空子串，然后把这 k 个子串按照其在字符串 A 中出现的顺序依次连接起来得到一个新的字符串，请问有多少种方案可以使得这个新串与字符串 B 相等？注意：子串取出的位置不同也认为是不同的方案。

输入格式

第一行是三个正整数 n，m，k，分别表示字符串 A 的长度，字符串 B 的长度，以及问题描述中所提到的 k ，每两个整数之间用一个空格隔开。
第二行包含一个长度为 n 的字符串，表示字符串 A 。
第三行包含一个长度为 m 的字符串，表示字符串 B 。

输出格式

输出共一行，包含一个整数，表示所求方案数。由于答案可能很大，所以这里要求输出答案对 1,000,000,007 取模的结果。

样例数据 1

输入　

6 3 1
aabaab
aab

输出

2

样例数据 2

输入　

6 3 2
aabaab
aab

输出

7

样例数据 3

输入　

6 3 3
aabaab
aab

输出

7

【数据范围】

对于第1组数据：1≤n≤500，1≤m≤50，k=1；
对于第2组至第3组数据：1≤n≤500，1≤m≤50，k=2；
对于第4组至第5组数据：1≤n≤500，1≤m≤50，k=m；
对于第1组至第7组数据：1≤n≤500，1≤m≤50，1≤k≤m；
对于第1组至第9组数据：1≤n≤1000，1≤m≤100，1≤k≤m；
对于所有10组数据：1≤n≤1000，1≤m≤200，1≤k≤m。

解题思路：

dp[i][j][k][0/1]表示A串的前i个字符中取出 k 个互不重叠的非空子串匹配到了B串的前j个字符，0/1表示A[i]是否和B[j]匹配，则：

如果A[i]=B[j]：dp[i][j][k][1]=dp[i-1][j-1][k][1]+dp[i-1][j-1][k-1][1]+dp[i-1][j][k][0]
如果A[i]$\ne$B[j]：dp[i][j][k][1]=0;
dp[i][j][k][0]=dp[i-1][j][k][0]+dp[i-1][j][k][1]。

要开滚动数组。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c<'0'||c>'9')&&c!='-';c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=1005,M=205,p=1e9+7;
int n,m,k,dp[2][M][M][2];
char s[N],t[N];

int main()
{
    //freopen("lx.in","r",stdin);
    n=getint(),m=getint(),k=getint();
    scanf("%s%s",s+1,t+1);
    dp[0][0][0][0]=dp[1][0][0][0]=1;
    int now=1;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
            for(int l=1;l<=k;l++)
            {
                if(s[i]==t[j])dp[now][j][l][1]=((dp[now^1][j-1][l][1]+dp[now^1][j-1][l-1][1])%p+dp[now^1][j-1][l-1][0])%p;
                else dp[now][j][l][1]=0;
                dp[now][j][l][0]=(dp[now^1][j][l][0]+dp[now^1][j][l][1])%p;
            }
        now^=1;
    }
    now^=1;
    int ans=(dp[now][m][k][0]+dp[now][m][k][1])%p;
    cout<<ans;
    return 0;
}