[BJWC2011]禁忌

IX.[BJWC2011]禁忌

老套路，建出AC自动机，并求出每个节点的$ok$。

之后呢？

考虑dp。

令$f_{i,j}$表示在一个长度为$i$的串，当前在点$j$的可能性。

则有：

对于任意$k\in [0,alphabet)$，

如果$t[t[j].ch[k]].ok$为$true$，那么$f_{i+1,k}+=f_{i,j}/alphabet$；

否则，在统计答案的变量中加上$f_{i,j}/alphabet$，同时$f_{i+1,1}+=f_{i,j}/alphabet$（匹配成功，默认跳到根节点）。

但是！！！这里的$len$巨大，我们不能这么跑，否则复杂度是$O(size*len*alphabet)$的。

考虑一下，这每次的转移的步骤，都与$i$和具体的$f_{i,j}$的值无关。因此可以考虑矩阵快速幂。

构建转移矩阵$base$。

for(int i=1;i<=cnt;i++){
	for(int j=0;j<m;j++){
		if(t[t[i].ch[j]].ok)base.num[i][t[i].ch[j]]+=1;
		else base.num[i][1]+=1,base.num[i][cnt+1]+=1;
	}
	for(int i=1;i<=base.N;i++)for(int j=1;j<=base.N;j++)base.num[i][j]/=m;
	base.num[base.N][base.N]=1.0;
}

（$m$即为$alphabet$。）

该矩阵的最后一列为统计答案列。但是，该矩阵的右下角处必需强制赋成$1$。为什么呢？

看到矩阵乘法$x*y$的过程中什么时候要用到右下角：

or(int i=1;i<=x.N;i++)for(int j=1;j<=y.N;j++)for(int k=1;k<=y.N;k++)z.num[i][j]+=x.num[i][k]*y.num[k][j];

则在做最后一行和最后一列的乘法时，它调用了右下角。

最后一行全为$0$（因为原本应该是一个$size*size$的方阵，但是我们强制补了最后一列），可以忽略。

但是，计算最后一列时，这一点是必须的。

然后我们就可以暴力矩阵快速幂了。

最后答案为矩阵的右上角。

因为最右一列是答案列，而最上一行为根节点行。

代码：

#include<bits/stdc++.h>
using namespace std;
int n,len,m,S,cnt=1;
struct AC_Automaton{
	int ch[26],fail;
	bool ok;
}t[80];
char s[20];
void ins(){
	int x=1;
	for(int i=0;i<S;i++){
		if(!t[x].ch[s[i]-'a'])t[x].ch[s[i]-'a']=++cnt,t[cnt].ok=true;
		x=t[x].ch[s[i]-'a'];
	}
	t[x].ok=false;
}
queue<int>q;
void build(){
	for(int i=0;i<m;i++){
		if(t[1].ch[i])t[t[1].ch[i]].fail=1,q.push(t[1].ch[i]);
		else t[1].ch[i]=1;
	}
	while(!q.empty()){
		int x=q.front();q.pop();
		for(int i=0;i<m;i++){
			if(t[x].ch[i])t[t[x].ch[i]].fail=t[t[x].fail].ch[i],q.push(t[x].ch[i]);
			else t[x].ch[i]=t[t[x].fail].ch[i];
		}
		t[x].ok&=t[t[x].fail].ok;
	}
}
struct Matrix{
	long double num[80][80];
	int N;
	Matrix(){
		N=0;
		memset(num,0,sizeof(num));
	}
	friend Matrix operator *(const Matrix &x,const Matrix &y){
		Matrix z;
		z.N=x.N;
		for(int i=1;i<=x.N;i++)for(int j=1;j<=y.N;j++)for(int k=1;k<=y.N;k++)z.num[i][j]+=x.num[i][k]*y.num[k][j];
		return z;
	}
}base,I;
Matrix pov(Matrix x,int y){
	if(!y)return I;
	Matrix tmp=pov(x,y>>1);
	return tmp*tmp*(y&1?x:I);
}
int main(){
	scanf("%d%d%d",&n,&len,&m),t[1].ok=true;
	for(int i=1;i<=n;i++)scanf("%s",s),S=strlen(s),ins();
	build();
	for(int i=1;i<=cnt;i++){
		for(int j=0;j<m;j++){
			if(t[t[i].ch[j]].ok)base.num[i][t[i].ch[j]]+=1;
			else base.num[i][1]+=1,base.num[i][cnt+1]+=1;
		}
	}
	base.N=I.N=cnt+1;
	for(int i=1;i<=base.N;i++)for(int j=1;j<=base.N;j++)base.num[i][j]/=m;
	base.num[base.N][base.N]=1.0;
	for(int i=1;i<=I.N;i++)I.num[i][i]=1;
	base=pov(base,len);
	printf("%.8Lf\n",base.num[1][base.N]);
	return 0;
}