BZOJ 1444: [Jsoi2009]有趣的游戏

1444: [Jsoi2009]有趣的游戏

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 1114  Solved: 386
[Submit][Status][Discuss]

Description

Input

注意 是0<=P

Output

Sample Input

Sample Output

HINT

 

30%的数据保证, n ≤ 2. 50%的数据保证, n ≤ 5. 100%的数据保证, n , l, m≤ 10.

Source

 

[Submit][Status][Discuss]

先做AC自动机，然后本来想gauss，但是听说做50次自乘就行，那岂不是贼简单。

#include<cstdio>
#include<cstring>
#define siz 125
char s[15]; double d[15],f[siz][siz],p[siz][siz];
int n,m,l,q[siz],hd,tl,fl[siz],ch[siz][15],ed[siz],bp[15],tot=1;
main(){
    scanf("%d%d%d",&n,&l,&m);
    for(int j=0,x,y;j<m;++j)
        scanf("%d%d",&x,&y),d[j]=1.0*x/y;
    for(int j=1,t;j<=n;++j){
        scanf("%s",s);t=1;
        for(int i=0;i<l;++i){
            if(!ch[t][s[i]-'A'])
                ch[t][s[i]-'A']=++tot;
            t=ch[t][s[i]-'A'];
        }
        ed[t]=j,bp[j]=t;
    }
    for(int j=0;j<m;++j)
        if(!ch[1][j])ch[1][j]=1;
        else fl[q[tl++]=ch[1][j]]=1;
    while(hd!=tl){
        int t=q[hd++];
        for(int j=0;j<m;++j)
            if(!ch[t][j])ch[t][j]=ch[fl[t]][j];
            else fl[q[tl++]=ch[t][j]]=ch[fl[t]][j];
    }
    for(int i=1;i<=tot;++i)
        if(!ed[i])for(int j=0;j<m;++j)
            f[i][ch[i][j]]+=d[j];
        else f[i][i]=1;
    for(int t=1;t<=50;++t){
        memset(p,0,sizeof(p));
        for(int k=1;k<=tot;++k)
            for(int i=1;i<=tot;++i)
                for(int j=1;j<=tot;++j)
                    p[i][j]+=f[i][k]*f[k][j];
        memcpy(f,p,sizeof(f));
    }
    for(int i=1;i<=n;++i)printf("%.2lf
",f[1][bp[i]]);
}

@Author: YouSiki