noip2015 洛谷 P2679 子串
阿新 • • 發佈:2018-12-23
題目:子串
思路:
大佬的思路
令f[i][j][kk][0/1]表示a的前i個,b的前j個,分為kk段,第i個取不取的方案數。
轉移方程:
if(a[i]!=b[j]) {
f[val][j][kk][0]=(f[val^1][j][kk][0]+f[val^1][j][kk][1])%md;
f[val][j][kk][1]=0;
} else {
f[val][j][kk][0]=(f[val^1][j][kk][0]+f[val^1][j][kk][1])%md;
f[val][j][kk][1]=((f[val^1][j-1][kk-1][0]+f[val^1][j-1][kk-1][1]) %md
+f[val^1][j-1][kk][1])%md;
}
程式碼:
#include<bits/stdc++.h>
using namespace std;
#define read(x) scanf("%d",&x)
#define maxn 1000
#define maxm 200
#define md 1000000007
int n,m,k;
char a[maxn+5],b[maxm+5];
int f[5][maxm+5][maxm+5][2];
int main() {
read(n),read(m),read(k);
scanf("%s",a+1);
scanf("%s",b+1);
f[1][0][0][0]=f[0][0][0][0]=1;
int val=1;
for(int i=1;i<=n;i++) {
for(int j=1;j<=m;j++) {
for(int kk=1;kk<=k;kk++) {
if(a[i]!=b[j]) {
f[val][j][kk][0]=(f[val^1][j][kk][0]+f[val^1][j][kk][1])%md;
f[val][j][kk][1]=0;
} else {
f[val][j][kk][0]=(f[val^1] [j][kk][0]+f[val^1][j][kk][1])%md;
f[val][j][kk][1]=((f[val^1][j-1][kk-1][0]+f[val^1][j-1][kk-1][1])%md+f[val^1][j-1][kk][1])%md;
}
}
}
val^=1;
}
printf("%d",(f[n&1][m][k][1]+f[n&1][m][k][0])%md);
return 0;
}