bzoj2553 [BeiJing2011]禁忌（ACAM+概率dp+矩阵快速幂）

先把所有禁忌串建成ACAM，搞出Trie图。然后我们首先对于一个字符串，他的伤害是多少，就可以转化成选最多的线段两两不相交。这个问题可以贪心地按右端点排序后能取就取。体现在ACAM上就是走到一个结束节点就返回根，++ans。f[i][j]表示长度为i的串匹配到ACAM的节点j的伤害期望值。然后长度超级大！我们需要矩阵倍增！但是我们还要记录ans，怎么办呢！再多开一个点累计答案！如果一个点再走到他的儿子就结束了，就向根连一条边，向新点连一条边。新点向自己连一条权值为1的边，用来累计答案。嗯，假装我理解了。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ld long double
#define inf 0x3f3f3f3f
#define N 80
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,K,son[N][26],fail[N],rt,k;
ld tmp;char s[N];
bool ed[N];
struct Matrix{
    ld a[N][N];
    Matrix(){}
    Matrix(bool t){memset(a,0,sizeof(a));if(t) 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;
    }friend Matrix operator^(Matrix a,int k){
        Matrix res(1);
        for(;k;k>>=1,a=a*a) if(k&1) res=res*a;return res;
    }
}trans;
inline void ins(){
    scanf("%s",s+1);int m=strlen(s+1),p=rt;
    for(int i=1;i<=m;++i){
        int &y=son[p][s[i]-'a'];if(!y) y=++n;p=y;
    }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;i<k;++i){
            int &y=son[x][i];
            if(!y){y=son[fail[x]][i];continue;}q.push(y);
            fail[y]=son[fail[x]][i];ed[y]|=ed[fail[y]];
        }
    }
}
int main(){
//  freopen("a.in","r",stdin);
    int m=read();K=read();k=read();
    rt=++n;for(int i=0;i<k;++i) son[0][i]=rt;
    while(m--) ins();buildAC();tmp=1.0/k;
    for(int i=1;i<=n;++i){
        if(ed[i]) continue;
        for(int j=0;j<k;++j)
            if(ed[son[i][j]]) trans.a[i][rt]+=tmp,trans.a[i][n+1]+=tmp;
            else trans.a[i][son[i][j]]+=tmp;
    }++n;trans.a[n][n]=1;
    trans=trans^K;
    printf("%.10Lf\n",trans.a[1][n]);
    return 0;
}