bzoj1444 [Jsoi2009]有趣的游戏（ACAM+矩阵+概率dp）

给定n个长度为l的模式串，现在要用前m个大写字母生成一个随机串，每个字符有自己的出现几率，第一次出现的字符串获胜，求最终每个字符串的获胜几率

建出ACAM，搞出Trie图，弄出转移矩阵：
如果某个节点是模式串结尾那么这个节点只向自己连一条概率为1的出边（因为只要一走到结尾节点游戏就停止了）否则的话向每个转移到的节点，加上转移的概率。

然后把转移矩阵自乘50遍即可。

吐槽：这题假的假的【再见】，我写矩阵倍增10000次就wa了，连乘50次就A了，也不知道是谁挂精度了【再见】

#include 
#include 
#include 
#include 
using namespace std;
#define ll long long
#define N 110
#define inf 0x3f3f3f3f
#define ld long double
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n=0,nn,l,m,fail[N],son[N][10],pos[11],rt,ed[N];
ld pro[11];
char s[N];
struct Matrix{
    ld a[N][N];
    Matrix(){}
    Matrix(bool flag){
        memset(a,0,sizeof(a));
        if(flag) for(int i=1;i<=n;++i) a[i][i]=1;
    }
    friend Matrix operator*(Matrix a,Matrix b){
        Matrix res(0);
        for(int i=1;i<=n;++i)
            for(int j=1;j<=n;++j)
                for(int k=1;k<=n;++k)
                    res.a[i][j]+=a.a[i][k]*b.a[k][j];return res;
    }
}a;
inline void ins(char *s,int id){
    int len=strlen(s+1),p=rt;
    for(int i=1;i<=len;++i){
        int &y=son[p][s[i]-'A'];if(!y) y=++n;p=y;
    }pos[id]=p;ed[p]=1;
}
inline void buildAC(){
    queue<int>q;q.push(rt);
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=0;iint &y=son[x][i];
            if(!y){y=son[fail[x]][i];continue;}q.push(y);
            fail[y]=son[fail[x]][i];
        }
    }
}
int main(){
//  freopen("a.in","r",stdin);
    nn=read();l=read();m=read();rt=++n;
    for(int i=0;i0][i]=rt;
    for(int i=0;iint x=read(),y=read();pro[i]=x*1.0/y;}
    for(int i=1;i<=nn;++i) scanf("%s",s+1),ins(s,i);buildAC();
    for(int i=1;i<=n;++i){
        if(ed[i]){a.a[i][i]=1;continue;}
        for(int j=0;jfor(int i=1;i<=50;++i) a=a*a;
    for(int i=1;i<=nn;++i) printf("%.2Lf\n",a.a[1][pos[i]]);
    return 0;
}