Vijos1982 NOIP2015Day2T2 子串 substring 动态规划

子串 (substring.cpp/c/pas)　　　题目链接

【问题描述】
有两个仅包含小写英文字母的字符串 A 和 B。现在要从字符串 A 中取出 k 个 互不重叠 的非空子串，然后把这 k 个子串按照其在字符串 A 中出现的顺序依次连接起来得到一
个新的字符串，请问有多少种方案可以使得这个新串与字符串 B 相等？注意：子串取出的位置不同也认为是不同的方案 。
【输入格式】
输入文件名为 substring.in。
第一行是三个正整数 n，m，k，分别表示字符串 A 的长度，字符串 B 的长度，以及问题描述中所提到的 k，每两个整数之间用一个空格隔开。
第二行包含一个长度为 n 的字符串，表示字符串 A。
第三行包含一个长度为 m 的字符串，表示字符串 B。
【输出格式】
输出文件名为 substring.out。
输出共一行，包含一个整数，表示所求方案数。 由于答案可能很大，所以这里要求对输出答案对 1,000,000,007 取模 的结果。
【输入输出样例 1】
substring.in
6 3 1
aabaab
aab
substring.out
2
【输入输出样例 2】
substring.in
6 3 2
aabaab
aab
substring.out
7
【输入输出样例 3】
substring.in
6 3 3
aabaab
aab
substring.out
7

【题解】
NOIP2015Day2T2
一道好好的DP题
我们用dp[i][j][k]表示在B串中匹配i个，在A串中匹配到的位置为j，共使用k个子串的方案总数，则dp[i][j][k]=Σdp[i-1][j'][k-1] +dp[i-1][j-1][k]
那么，对于Σ可以用前缀和优化，这样的时间就可以卡进去了，但是空间还是要炸，所以我们采用滚动数组来优化空间即可。详见代码。

#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
using namespace std;
const int N=+,M=+;
const int mod=1e9+;
char s1[N],s2[M];
int n,m,K;
int dp[][N][M],sum[][N][M];
//dp[i][j][k]=Σdp[i-1][j'][k-1](1<=j'<=j-1) +dp[i-1][j-1][k]
int main(){
    scanf("%d%d%d%s%s",&n,&m,&K,s1+,s2+);
    memset(dp,,sizeof dp);
    memset(sum,,sizeof sum);
    int I=,J=;
    for (int i=;i<=n;i++){
        if (s2[]==s1[i])
            dp[][i][]=;
        sum[][i][]=sum[][i-][]+dp[][i][];
    }
    for (int i=;i<=m;i++,I^=,J^=){
        memset(dp[J],,sizeof dp[J]);
        memset(sum[J],,sizeof sum[J]);
        for (int j=;j<=n;j++){
            if (s2[i]!=s1[j])
                continue;
            for (int k=;k<=K;k++)
                if (j>=)
                    dp[J][j][k]=(sum[I][j-][k-]+dp[I][j-][k])%mod;
                else
                    dp[J][j][k]=dp[I][j-][k];
        }
        for (int k=;k<=K;k++)
            for (int j=;j<=n;j++)
                sum[J][j][k]=(sum[J][j-][k]+dp[J][j][k])%mod;
    }
    printf("%d",sum[I][n][K]);
    return ;
}