原题链接:https://loj.ac/problem/2424
子串
题目描述
有两个仅包含小写英文字母的字符串 A A A 和 B B B。现在要从字符串 A A A中取出 k k k个互不重叠的非空子串,然后把这 k k k个子串按照其在字符串 A A A中出现的顺序依次连接起来得到一个新的字符串,请问有多少种方案可以使得这个新串与字符串 B B B相等?
注意:子串取出的位置不同也认为是不同的方案。
输入格式
第一行是三个正整数
n
,
m
,
k
n,m,k
n,m,k,分别表示字符串
A
A
A的长度,字符串
B
B
B的长度,以及问题描述中所提到的
k
k
k,每两个整数之间用一个空格隔开。
第二行包含一个长度为
n
n
n的字符串,表示字符串
A
A
A。
第三行包含一个长度为
m
m
m的字符串,表示字符串
B
B
B。
输出格式
输出共一行,包含一个整数,表示所求方案数。由于答案可能很大,所以这里要求输出答案对 1 , 000 , 000 , 007 1,000,000,007 1,000,000,007取模的结果。
样例
样例输入 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
1组数据:
1
≤
n
≤
500
,
1
≤
m
≤
50
,
k
=
1
1 \leq n \leq 500,1 \leq m \leq 50,k = 1
1≤n≤500,1≤m≤50,k=1;
对于第
2
2
2组至第
3
3
3组数据:
1
≤
n
≤
500
,
1
≤
m
≤
50
,
k
=
2
1 \leq n \leq 500,1 \leq m \leq 50,k = 2
1≤n≤500,1≤m≤50,k=2;
对于第
4
4
4组至第
5
5
5组数据:
1
≤
n
≤
500
,
1
≤
m
≤
50
,
k
=
m
1 \leq n \leq 500,1 \leq m \leq 50,k = m
1≤n≤500,1≤m≤50,k=m;
对于第
1
1
1组至第
7
7
7组数据:
1
≤
n
≤
500
,
1
≤
m
≤
50
,
1
≤
k
≤
m
1 \leq n \leq 500,1 \leq m \leq 50,1 \leq k \leq m
1≤n≤500,1≤m≤50,1≤k≤m;
对于第
1
1
1组至第
9
9
9组数据:
1
≤
n
≤
1000
,
1
≤
m
≤
100
,
1
≤
k
≤
m
1 \leq n \leq 1000,1 \leq m \leq 100,1 \leq k \leq m
1≤n≤1000,1≤m≤100,1≤k≤m;
对于所有
10
10
10组数据:
1
≤
n
≤
1000
,
1
≤
m
≤
200
,
1
≤
k
≤
m
1 \leq n \leq 1000,1 \leq m \leq 200,1 \leq k \leq m
1≤n≤1000,1≤m≤200,1≤k≤m。
题解
d p [ i ] [ j ] [ k ] dp[i][j][k] dp[i][j][k]表示考虑到 A A A串第 i i i位, B B B串的第 j j j位,拆成 k k k份的方案数, d i v i [ i ] [ j ] [ k ] divi[i][j][k] divi[i][j][k]表示 A [ i ] = B [ j ] A[i]=B[j] A[i]=B[j]且分配方案中包含 A [ i ] , B [ j ] A[i],B[j] A[i],B[j],分成 k k k份的方案数,就有转移如下:
if(a[i]==b[j])divi[i][j][k]=divi[i-1][j-1][k]+dp[i-1][j-1][k-1];
else divi[i][j][k]=0;
dp[i][j][k]=divi[i][j][k]+dp[i-1][j][k];
发现只用得到 i i i与 i − 1 i-1 i−1,于是就可以滚动数组,又发现只会用到 j j j和 j − 1 j-1 j−1,所以我们倒过来更新,省掉一维。
代码
#include<cstdio>
const int N=1005,M=205,mod=1e9+7;
long long dp[M][M]={1},divi[M][M];
int n,m,p;
char a[N],b[M];
void in(){scanf("%d%d%d%s%s",&n,&m,&p,a,b);}
void ac()
{
for(int i=1;i<=n;++i)for(int j=m;j;--j)for(int k=p;k;--k)
(dp[j][k]+=(divi[j][k]=a[i-1]==b[j-1]?divi[j-1][k]+dp[j-1][k-1]:0))%=mod;
printf("%lld",dp[m][p]);
}
int main(){in(),ac();}