hdu 3689 Infinite monkey theorem dp

dp[i][j] 表示长度为i的时候匹配到单词的第j位，并且没有出现过该单词的概率。

f[i][j] 表示匹配到第i位时加入‘a'+j时会匹配到哪一位。

先暴力出f[i][j]。然后转移就非常简答了。

#include <iostream>

#include <cstdio>
#include <cstring>
using namespace std;
int n,m;
double a[1001];
char t[101];
double dp[1011][14];
int f[101][101];


int check(int r,int s)
{
    for(int i=s,j=0;i>=1;i--,j++)
    {
//        printf("%d %d\n",i,r-j);
        if(t[i]!=t[r-j])
        return(0);
    }
    return(1);
}


int work(int r,int s)
{
    for(int j=r+1;j>=1;j--)
    if(t[j]==s+'a')
    {
        if(check(r,j-1))
        return(j);
    }
    return(0);
}


int main()
{
//    freopen("in.txt","r",stdin);
    while(scanf("%d %d",&n,&m),n||m)
    {
        memset(a,0,sizeof(a));
        for(int i=1;i<=n;i++)
        {
            getchar();
            char tmp;
            scanf("%c",&tmp);
            scanf("%lf",&a[tmp]);
        }
        scanf("%s",&t[1]);


        int lon=strlen(&t[1]);


//        printf("%s",&t[1]);


//          f[2]['o'-'a']=work(2,'o'-'a');


        for(int i=0;i<lon;i++)
        for(int j=0;j<26;j++)
        {
            f[i][j]=work(i,j);
        }


//        printf("%d\n",f[2]['o'-'a']);


        memset(dp,0,sizeof(dp));


        dp[0][0]=1;


        for(int k=0;k<m;k++)
        for(int i=0;i<lon;i++)
        for(int j=0;j<26;j++)
        {
            dp[k+1][f[i][j]]+=dp[k][i]*a[j+'a'];
        }


//        for(int k=0;k<m;k++)
//        for(int i=0;i<lon;i++)
//        printf("%d %d %lf\n",k,i,dp[k][i]);


        double ans=0;


        for(int i=0;i<lon;i++)
        ans+=dp[m][i];


        printf("%.2f",100-ans*100);
        printf("%c\n",'%');
    }
    return 0;
}