NOIP2015Day2T2 子串
题目描述
有两个仅包含小写英文字母的字符串A和B。
现在要从字符串A中取出k个互不重叠的非空子串，然后把这 k 个子串按照其在字符串 A 中出现的顺序依次连接起来得到一个新的字符串，请问有多少种方案可以使得这个新串与字符串 B 相等？
注意：子串取出的位置不同也认为是不同的方案。
输入格式
第一行是三个正整数 n，m，k，分别表示字符串 A 的长度，字符串 B 的长度，以及问题描述中所提到的 k，每两个整数之间用一个空格隔开。
第二行包含一个长度为 n 的字符串，表示字符串A。
第三行包含一个长度为 m 的字符串，表示字符串B。
输出格式
输出共一行，包含一个整数，表示所求方案数。
由于答案可能很大，所以这里要求输出答案对 1,000,000,007 取模的结果。
数据范围
输入样例：
6 3 1
aabaab
aab
输出样例：
2
解题报告
题意理解
就是让你从字符串A中,拿出\(k\)个子串,子串之间不得有重复,要求\(k\)个子串,按照拿出来的顺序,可以拼接成为字符串B.
要你求出符合条件的方案总数.
所有合法方案如下：（加下划线的部分表示取出的子串）
样例一：aab aab / aab aab
样例二：a ab aab / a aba ab / a a ba ab / aab a ab / aa b aab / aa baa b / aab aa b
样例三：a a b aab / a a baa b / a ab a a b / a aba a b / a a b a a b / a a ba a b / aab a a b
思路解析
算法确定
一般来说,当我们看到由于答案可能很大，所以这里要求输出答案对 1,000,000,007 取模的结果。这句话的时候,一般来说动态规划算法是可以基本确定的了.
每一个字符串,其实我们都可以看作是字符区间.
我们接着来分析题目的内容.
我们发现,对于题目而言,一个非常重要的点就是,他是一个区间.
区间有关的动态规划,最基本的有哪些呢?
- 区间DP
- 线性DP
- 数位DP
我们发现线性DP似乎是一个不错的选择,因为区间DP,这道题目似乎并木有给出明显的区间划分,反而都是线性的思想比较明显.
至于数位DP跟本题没有任何关系.
所以现在我们完全可以确定,本题的算法就是线性动态规划
状态设置
我们观察这个题目,题目的大致意思,就是让我们选取一些字符串出来,去构造字符串B.
那么我们状态中,必须有以下这些信息.
- 关于字符串A的信息.
- 关于字符串B的信息
- 我们当前用了多少个字符串,因为我们最多只能选取k个字符串.
这就是我们可以确定的必备信息条件,那么接下来我们如何设置状态表示呢?
根据上面的信息表,我们不难想到我们的基本状态如何设置.
举个例子,让我们更好的理解一下.
但是我们发现我们的状态转移方程不太好想出来,为什么呢?
我们发现对于K的转移,我们很难划分.
比如说我们当前选择了第\(j\)个字符.
但是我们不知道我们之前的\(j-1\)这个字符是否选择了.
如果选择了,那么我们的k不需要变化.
没有选择,我们的k就需要+1.
综上所述,我们不得不处理以下这个选取问题.
我们可以增加一个维度,去存储当前字符是否选择了.
转移方程
1.当前\(a[i]==b[j]\)
这是什么意思?就是我们当前这个字符,可以匹配B串了.
f[i][j][k][1]=(f[i][j-1][k][1]+f[i][j-1][k-1][0])%Mod;
f[i][j][k][0]=(f[i][j][k][0]+f[i][j][k][1])%Mod;
我们来慢慢分析每一句话.
f[i][j][k][1]=(f[i][j-1][k][1]+f[i][j-1][k-1][0])%Mod;
对于这句话而言.
假如说,我们当前选择第\(i\)个字符.
接着我们分析第二句话.
f[i][j][k][0]=(f[i-1][j][k][0]+f[i-1][j][k][1])%Mod;
假如说我们当前不选择第\(i\)个字符.
就是这一位取的+这一位没取的（也就是上一段的0）
Hint:空间太大,所以我们要状态压缩.
我们发现第一维\(i\)其实并没有什么作用.
所以我们使用背包问题的优化思想,让循环反过来,然后这道题目就解决了.
#include <bits/stdc++.h>
using namespace std;
const int Mod=1000000007;
#define read() ios::sync_with_stdio(false),cin.tie(0),cout.tie(0) //读入优化
string a,b;
int n,m,k2,f[210][210][2];
int main()
{
read();
cin>>n>>m>>k2>>a>>b;
a=" "+a;//字符串有用部分从1开始
b=" "+b;//字符串有用部分从1开始
f[0][0][0]=1;
for(int i=1; i<=n; i++)
for(int j=m; j; j--)//循环反过来
if (a[i]==b[j])//相等判定
{
for(int k=min(j,k2); k; k--)//k必然小于j,而且不能超过范围
{
f[j][k][1]=(f[j-1][k][1]+f[j-1][k-1][0])%Mod;
f[j][k][0]=(f[j][k][0]+f[j][k][1])%Mod;
}
}
else
for(int k=1; k<=m; k++)//初始化把
f[j][k][1]=0;
cout<<(f[m][k2][0])%Mod;
return 0;
}
