子串 [NOIP2015]

Description

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

Input

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

Output

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

Sample Input

样例输入1：
6 3 1
aabaab
aab

样例输入2：
6 3 2
aabaab
aab

Sample Output

样例输出1：
2

样例输出2：
7

Hint

数据范围：
对于第 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。

Solution

这道题一看tm就是道DP题，而且第一感觉会很复杂

设状态

首先先来观察，我们设状态的话，要存储的应该有当前位置i，被匹配到的位置j，用了几个子串k（dp[i][j][k]）

但我们发现还有一个棘手的问题，如何判断当前的是一个独立的子串还是连接到前面的子串

那么我们就要知道某一个位置是否被匹配，再加上一维[0/1]代表是否被匹配

状态转移

对于dp[i][j][k][0]，他没被匹配，那匹配位置没变，子串个数没变，从i-1转移就是 dp[i-1][j][k][0]+dp[i-1][j][k][1]

对于dp[i][j][k][1]，前提条件就是A[i]==B[j]，满足条件后我们发现，它可以分别从 (1)前一位没被匹配 (2)前一位被匹配，这个字符连接到前面一个子串 (3)前一位没被匹配，这一个字符单独成为一个子串，那么转移方程为dp[i-1][j-1][k-1][0]+dp[i-1][j-1][k-1][1]+dp[i-1][j-1][k][1]

细节&优化

边界条件：dp[i][1][1][0]=sigma(j=1->i-1)(A[j]==B[1])（显然可以累计，不需要每次求）　　dp[i][1][1][1]=(A[i]==B[1])

dp[1000][200][200][2]的空间是显然爆炸的，由于只要用到前一维，可以滚一滚第一维

而每次计算完dp[i]后，要清空dp[i-1]，否则会对dp[i+1]产生干扰(因为是滚动的)

Code

 #include<set>
 #include<map>
 #include<stack>
 #include<queue>
 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define RG register int
 #define rep(i,a,b)    for(RG i=a;i<=b;i++)
 #define per(i,a,b)    for(RG i=a;i>=b;i--)
 #define ll long long
 #define inf (1<<30)
 #define maxm 205
 using namespace std;
 const ll mo=;
 int n,m,t,p,q;
 ll sum;
 ll dp[][maxm][maxm][];
 char A[],B[maxm];
 inline int read()
 {
     int x=,f=;char c=getchar();
     while(c<''||c>''){if(c=='-')f=-;c=getchar();}
     while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
     return x*f;
 }

 int main()
 {
     n=read(),m=read(),t=read();
     scanf("%s",A+);scanf("%s",B+);
     p=,q=;
     rep(i,,n)
     {
         swap(p,q);
         dp[p][][][]=sum;
         if(A[i]==B[])    dp[p][][][]=,sum++;
         rep(j,,m)
             rep(k,,t)
             {
                 if(A[i]==B[j])    (dp[p][j][k][]=(dp[q][j-][k-][]+dp[q][j-][k][])%mo+dp[q][j-][k-][])%=mo;
                 dp[p][j][k][]=(dp[q][j][k][]+dp[q][j][k][])%mo;
             }
        rep(j,,m)    rep(k,,t)    dp[q][j][k][]=dp[q][j][k][]=;
     }
     cout<<(dp[p][m][t][]+dp[p][m][t][])%mo;
     return ;
 }