洛谷 P2679 子串

题目背景

无

题目描述

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

输入输出格式

输入格式：

输入文件名为 substring.in。

第一行是三个正整数 n,m,k,分别表示字符串 A 的长度,字符串 B 的长度,以及问

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

输出格式：

输出文件名为 substring.out。 输出共一行,包含一个整数,表示所求方案数。由于答案可能很大,所以这里要求[b]输出答案对 1,000,000,007 取模的结果。[/b]

输入输出样例

输入样例#1： 

6 3 1 
aabaab 
aab

输出样例#1： 

2

输入样例#2：

6 3 2 
aabaab 
aab

输出样例#2：

7

输入样例#3： 

6 3 3 
aabaab 
aab

输出样例#3： 

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 ]为A用到了 i ，B用到了 j ，已经用了 k 个子串， 并且一定用了当前字符(A[i])时的方案数。

设f[ i ][ j ][ k ]为A用到了 i ，B用到了 j ，已经用了 k 个子串， 无论用不用当前字符(A[i])时的方案数总和。

接下来这个转移可就有蛮难想了。

一个一个来，

先分析一下 s 的转移。

能转移的前提自然是 A[ i ] == B [ j ]啦。

既然 A[i] 一定要用，那么依旧是两种情况：独自成一串 或 与前面的成一串。

独自成一串，方案数为：f[ i-1 ][ j-1 ][ k-1]

与前方共成一串，方案数为：dp[ i-1 ][ j-1 ][ k ]，因为前一个字符串(A[i-1])也一定要用！

所以合并一下： dp[ i ][ j ][ k ] = f[ i-1 ][ j-1 ][ k-1 ] + dp[ i-1 ][ j-1 ][ k ]；

接着分析 f 的转移。

f[ i ][ j ][ k ] 的来源也有两种： 使用当前字符 或 不使用当前字符

对于使用当前字符，方案数算法如上，答案即：dp[ i ][ j ][ k ]；

对于不使用当前字符，则从f[ i-1 ]转来，即：f[ i -1 ][ j ][ k ]；

合并一下： f[ i ][ j ][ k ] = f[ i-1 ][ j ][ k ] + dp[ i ][ j ][ k ]；

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<algorithm>
 5 using namespace std;
 6 const int MAXN=1001;
 7 const int INF=0x7fffff;
 8 const int mod=1e9+7;
 9 inline int read()
10 {
11     char c=getchar();int flag=1,x=0;
12     while(c<'0'||c>'9')    {if(c=='-')    flag=-1;c=getchar();}
13     while(c>='0'&&c<='9')    x=x*10+c-48,c=getchar();return x*flag;
14 }
15 int dp[3][MAXN][MAXN];//一定要用 
16 int f[3][MAXN][MAXN];//可以不用 
17 char s1[MAXN],s2[MAXN];
18 int n,m,t;
19 int main()
20 {
21     n=read();m=read();t=read();
22     scanf("%s",s1+1);scanf("%s",s2+1);
23     int now=1,past=0;
24     f[0][0][0]=1;
25     for(int i=1;i<=n;i++)
26     {
27         f[now][0][0]=1;
28         for(int j=1;j<=m;j++)
29         {
30             for(int k=1;k<=t;k++)
31             {
32                 if(s1[i]==s2[j])    dp[now][j][k]=(dp[past][j-1][k]+f[past][j-1][k-1])%mod;
33                 else dp[now][j][k]=0;
34                 f[now][j][k]=(f[past][j][k]+dp[now][j][k])%mod;
35             }
36         }
37         swap(now,past);
38     }
39     printf("%d",f[past][m][t]);
40     return 0;
41 }