BZOJ 1444 [Jsoi2009] 有趣的游戏

Description

Input

注意 是0<=P

Output

Sample Input

Sample Output

HINT

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

Source

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

AC自动机+矩阵乘法~

把所有串在AC自动机上建出来，记录终止点位置，然后建出方程。这里应该是要用高斯消元的，但是数据小可以用50次矩阵乘法水过~

#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;

int n,l,m,x,y,cnt,c[125][11],last[125],ed[125],bp[125],kkzv;
double gl[11],a[2][125][125];
char s[11];

int read()
{
	int totnum=0,f=1;char ch=getchar();
	while(ch<'0' || ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0' && ch<='9') {totnum=(totnum<<1)+(totnum<<3)+ch-'0';ch=getchar();}
	return totnum*f;
}

struct ac{
	void add(int u)
	{
		scanf("%s",s+1);
		int x=1;
		for(int i=1;i<=l;i++)
		{
			if(!c[x][s[i]-'A']) c[x][s[i]-'A']=++cnt;
			x=c[x][s[i]-'A'];
		}
		ed[x]=u;bp[u]=x;
	}
	void getlast()
	{
		queue<int> q;
		for(int i=0;i<m;i++)
		  if(!c[1][i]) c[1][i]=1;
		  else last[c[1][i]]=1,q.push(c[1][i]);
		while(!q.empty())
		{
			int x=q.front();q.pop();
			for(int i=0;i<m;i++)
			  if(!c[x][i]) c[x][i]=c[last[x]][i];
			  else last[c[x][i]]=c[last[x]][i],q.push(c[x][i]);
		}
	}
}ac;

int main()
{
	n=read();l=read();m=read();cnt=1;
	for(int i=0;i<m;i++) x=read(),y=read(),gl[i]=(double)x/y;
	for(int i=1;i<=n;i++) ac.add(i);ac.getlast();
	for(int i=1;i<=cnt;i++)
	  if(!ed[i]) for(int j=0;j<m;j++) a[kkzv][i][c[i][j]]+=gl[j];
	  else a[kkzv][i][i]=1;
	for(int kkz=1;kkz<=50;kkz++)
	{
		kkzv^=1;
		memset(a[kkzv],0,sizeof(a[kkzv]));
		for(int i=1;i<=cnt;i++)
		  for(int j=1;j<=cnt;j++)
		    for(int k=1;k<=cnt;k++) a[kkzv][i][j]+=a[kkzv^1][i][k]*a[kkzv^1][k][j];
	}
	for(int i=1;i<=n;i++) printf("%.2lf\n",a[kkzv][1][bp[i]]);
	return 0;
}