codechef [snackdown2017 Onsite Final] AND Graph

传送门

题解给出了一个很强势的dp：

i<K

$$dp[i][len]*Fib[len+2-(t[i]==1)] -> dp[i+1][len]$$

$$dp[i][len]*Fib[len+1-(t[i]==1)] -> dp[i+1][len+1]$$

i>=K

$$dp[i][len]*Fib[len+2-(t[i]==1)-(s[i]==1)] -> dp[i+1][len]$$

其中K是s的前导0个数

意思是单独考虑每一位，相邻两个不能同时为1，方案数是斐波那契数，然后就直接搞了。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define MN 5001
using namespace std;

const int MOD=;
int T,n,dp[MN][MN],k,F[MN],MMH;
char s[MN],t[MN];
int main(){
    F[]=;F[]=;
    for (int i=;i<=;i++) F[i]=F[i-]+F[i-],F[i]-=F[i]>=MOD?MOD:;
    scanf("%d",&T);
    while (T--){
        scanf("%s%s",s,t);MMH=;
        n=strlen(s);k=;
        while (s[k]=='') k++;
        for (int i=;i<n;i++)
        for (int j=;j<=n;j++) dp[j][i]=;
        dp[][]=;
        for (int i=;i<n;i++){
            for (int j=;j<k;j++){
                dp[j+][i]=(1LL*dp[j][i]*F[i+-(t[j]=='')]+dp[j+][i])%MOD;
                dp[j+][i+]=(1LL*dp[j][i]*F[i+-(t[j]=='')]+dp[j+][i+])%MOD;
            }
            for (int j=k;j<n;j++) dp[j+][i]=(1LL*dp[j][i]*F[i+-(t[j]=='')-(s[j]=='')]+dp[j+][i])%MOD;
            (MMH+=dp[n][i])%=MOD;
        }
        printf("%d\n",MMH);
    }
}