NOIP2015 D2T2子串

题目描述
有两个仅包含小写英文字母的字符串 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 取模的结果。

这个题还是比较好看出是DP的，因为子串是按A中出现顺序连接的。但是怎么设状态和转移呢？用dp[i][j]表示前i个位置用了j个串？可是i个位置是a串的还是b串的？dp[i][j][k]表示a串第i个位置，b串第j个位置，用了k个子串。如果i位置的字母等于j位置的字母时，就匹配上了。这时有两种情况。这个字母是跟着上一个串的，这个字母是自己串的开头。
所以：dp[i][j][k]=dp[i-1][j-1][k]+dp[i-1][j-1][k-1] (a[i]==b[j])
如果当前位置字母不一样呢？dp【i】【j】【k】=dp【i-1】【j】【k】；
那么问题来了，即使一样我们也不一定要选这个字母呀。。所以上面的状态转移方程，是错的。。
然后我就不会了。
然后我就思考了一下（看了下题解）。
我们可以用再用一个状态转移方程。
f【i】【j】【k】表示当前位置下字母一定要用。
g【i】【j】【k】表示当前位置下字母可以用可以不用。
f【i】【j】【k】=f【i-1】【j-1】【k】+g【i-1】【j-1】【k-1】；两种情况
g【i】【j】【k】=f【i】【j】【k】+g【i-1】【j】【k】；两种情况
如果i j不匹配则为0 不更新。
而初始状态就为没选时等于1：f【N】【0】【0】=g【N】【0】【0】=1；
这道题空间开大会MLE，因此需要滚动数组优化。

但这道题。。我们可以倒序循环来压维。把i的一维直接压掉。其他代码一样。下面是代码，短的惊人。（借用某题解的一句话，这真是我见过思维难度和代码难度相差最大的一道题）

#include<bits/stdc++.h>
using namespace std;

const int MAXN=1e3+5;
const int MINN=2e2+5;
const int MOD=1e9+7;

char a[MAXN],b[MINN];
int n,m,k,dp1[MINN][MINN],dp2[MINN][MINN]; 

int main(){
    scanf("%d%d%d",&n,&m,&k);
    scanf("%s%s",a,b);
    dp1[0][0]=1;
    dp2[0][0]=1;
    for(int i=1;i<=n;i++)
        for(int j=m;j>=1;j--)
            for(int k1=1;k1<=k;k1++){
                if(a[i-1]!=b[j-1]){
                    dp1[j][k1]=0;continue;
                }       
                dp1[j][k1]=(dp1[j-1][k1]%MOD+dp2[j-1][k1-1]%MOD)%MOD;
                dp2[j][k1]=(dp1[j][k1]%MOD+dp2[j][k1]%MOD)%MOD;
            }
    printf("%d\n",dp2[m][k]%MOD);
    return 0;
}