刷题记录-luoguP2679 子串
2017-09-21 00:43
246 查看
f[i][j][k]表示用A前i个字符且必选A[i]来匹配B前j个字符,分成k段
状态转移方程:
f[i][j][k]=f[i-1][j-1][k]+∑(p=j-1~i-1)f[p][j-1][k-1]
发现f[p][j-1][k-1]可以边转移边计算,于是时间复杂度O(n*m*k)
开过滚动之后空间也是可以的。
边界不大好想,转移的细节也有不少,真的要注意了
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define MAXN 1005
#define MAXM 205
#define MOD 1000000007
using namespace std;
int f[MAXN][MAXM];
int g[MAXN][MAXM];
char A[MAXN],B[MAXM];
int n,m,K;
int main()
{
// freopen("T1.in","r",stdin);
scanf("%d%d%d",&n,&m,&K);
scanf("%s%s",A+1,B+1);
for(int i=0;i<=n;i++){
f[i][0]=1;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
if(A[i]==B[j]){
f[i][j]=f[i-1][j-1];
}
}
}
for(int i=0;i<=n;i++){
f[i][0]=0;
}
for(int k=2;k<=K;k++){
for(int j=1;j<=m;j++){
int sum=0;
for(int i=j;i<=n;i++){
sum=(sum+f[i-1][j-1])%MOD;
if(A[i]==B[j]){
g[i][j]=(g[i-1][j-1]+sum)%MOD;
}
}
}
memcpy(f,g,sizeof(f));
memset(g,0,sizeof(g));
}
// for(int k=1;k<=K;k++){
// for(int i=1;i<=n;i++){
// for(int j=1;j<=m;j++){
// printf("%d ",f[i][j][k]);
// }
// printf("\n");
// }
// printf("\n\n");
// }
int ans=0;
for(int i=1;i<=n;i++){
ans=(ans+f[i][m])%MOD;
}
printf("%d\n",ans);
return 0;
}
状态转移方程:
f[i][j][k]=f[i-1][j-1][k]+∑(p=j-1~i-1)f[p][j-1][k-1]
发现f[p][j-1][k-1]可以边转移边计算,于是时间复杂度O(n*m*k)
开过滚动之后空间也是可以的。
边界不大好想,转移的细节也有不少,真的要注意了
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define MAXN 1005
#define MAXM 205
#define MOD 1000000007
using namespace std;
int f[MAXN][MAXM];
int g[MAXN][MAXM];
char A[MAXN],B[MAXM];
int n,m,K;
int main()
{
// freopen("T1.in","r",stdin);
scanf("%d%d%d",&n,&m,&K);
scanf("%s%s",A+1,B+1);
for(int i=0;i<=n;i++){
f[i][0]=1;
}
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
if(A[i]==B[j]){
f[i][j]=f[i-1][j-1];
}
}
}
for(int i=0;i<=n;i++){
f[i][0]=0;
}
for(int k=2;k<=K;k++){
for(int j=1;j<=m;j++){
int sum=0;
for(int i=j;i<=n;i++){
sum=(sum+f[i-1][j-1])%MOD;
if(A[i]==B[j]){
g[i][j]=(g[i-1][j-1]+sum)%MOD;
}
}
}
memcpy(f,g,sizeof(f));
memset(g,0,sizeof(g));
}
// for(int k=1;k<=K;k++){
// for(int i=1;i<=n;i++){
// for(int j=1;j<=m;j++){
// printf("%d ",f[i][j][k]);
// }
// printf("\n");
// }
// printf("\n\n");
// }
int ans=0;
for(int i=1;i<=n;i++){
ans=(ans+f[i][m])%MOD;
}
printf("%d\n",ans);
return 0;
}
相关文章推荐
- [luoguP2679] 子串(DP)
- LuoguP2679[NOIP2015] 子串 解题报告【多维DP】
- 刷题记录luoguP1972 [SDOI2009]HH的项链
- 【NOIP2015-Day2-T2】洛谷2679:子串 题解
- [luogu-2679]noip2015day2-T2 子串 题解
- 刷题记录-luoguP1525 关押罪犯
- 刷题记录luoguP1983 车站分级
- 刷题记录-luoguP1265 公路修建
- 【NOIP2015】洛谷2679 子串
- <记录>删除字符串中的子串
- 【p1032-子串变换】解题记录
- 刷题记录-luoguP2149 [SDOI2009]Elaxia的路线
- java 记录一个子串在整串中出现的次数实例
- Manarcher 求 字符串 的最长回文子串 【记录】
- 刷题记录-luoguP1991 无线通讯网
- Luogu 2679 子串 (动态规划)
- noip2015 子串 luogu2679
- java-记录一个子串在整串中出现的次数
- 【LeetCode题目记录-3】字符串中最长的没有重复字符的子串
- 刷题记录-luoguP2055 [ZJOI2009]假期的宿舍