CF的每道DP题目都能让我回味良久,太经典了!!!!!!!!
题意:给两个字符串A,B,每次操作能把A分为两部分C,D,组成DC,比如A是abcdefg,abc | defg ----->defg | abc,然后指定操作次数K,问在K步之内将A转化为B的方法数。
链接:http://codeforces.com/problemset/problem/176/B
解析:先暴力求出循环点个数(循环点:在这个点切断之后可以生成B),然后dp【k】【2】,dp【i】【0】表示现在和B一样,dp【i】【1】表示现在和B不一样,转移方程看代码。
#include<stdio.h>
#include<string.h>
#define LL __int64
#define MOD 1000000007LL
const int maxn = 1005;
char p[maxn],s[maxn];
LL dp[maxn*100][2];
int main()
{
int i,j,cnt;
int k,lenp,lens;
scanf("%s%s",p,s);
scanf("%d",&k);
lenp = strlen(p);
lens = strlen(s);
if(lenp - lens){
puts("0");
return 0;
}
cnt = 0;
for(i = 0;i < lenp;i ++){
p[i + lenp] = p[i];
if(!strncmp(p + i,s,lenp)) cnt ++;
}
if(!strncmp(s,p,lenp)) dp[0][0] = 1,dp[0][1] = 0;
else dp[0][0] = 0,dp[0][1] = 1;
for(i = 1;i <= k;i ++){
dp[i][0] = ( dp[i - 1][0] * (cnt - 1) % MOD + dp[i - 1][1] * cnt % MOD) % MOD;
dp[i][1] = (dp[i - 1][0] *(lenp - cnt) % MOD + dp[i - 1][1] * (lenp - cnt - 1) % MOD) % MOD;
}
printf("%I64d\n",dp[k][0]%MOD);
return 0;
}