【bzoj1444】[Jsoi2009]有趣的游戏

1444: [Jsoi2009]有趣的游戏

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 1007  Solved: 334
[Submit][Status][Discuss]

Description

Input

注意 是0<=P

Output

Sample Input

Sample Output

 

 

 

 

【题解】

 

AC自动机+矩阵乘法

 

首先把模式串建成AC自动机，构建出转移矩阵。

 

构造方法：a[i][j]表示从第i个结点转移到第j个结点的概率，则如果j被标记过，f[i][j]=1，否则f[i][j]=possble[ch[j]]

 

具体见代码：

 

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<ctime>
 #include<cmath>
 #include<algorithm>
 using namespace std;
 #define MAXN 510
 struct node{double p[MAXN][MAXN]; node(){memset(p,,sizeof(p));}}a;
 int n,l,m,cnt,id,fail[MAXN],end[MAXN],pos[MAXN],q[MAXN],tr[MAXN][];
 double chty[MAXN];
 char ch[MAXN];
 inline int read()
 {
     int x=,f=;  char ch=getchar();
     while(!isdigit(ch))  {if(ch=='-')  f=-;  ch=getchar();}
     while(isdigit(ch))  {x=x*+ch-'';  ch=getchar();}
     return x*f;
 }
 void insert()
 {
     int now=;
     for(int i=;i<=l;i++)
     {
         if(!tr[now][ch[i]-'A'])  tr[now][ch[i]-'A']=++cnt;
         now=tr[now][ch[i]-'A'];
     }
     end[now]=;  pos[++id]=now;
 }
 void build()
 {
     int head=,tail=;
     for(int i=;i<m;i++)  if(tr[][i])  q[++tail]=tr[][i];
     while(++head<=tail)
     {
         int x=q[head];  end[x]|=end[fail[x]];
         for(int i=;i<m;i++)
         {
             if(!tr[x][i])  tr[x][i]=tr[fail[x]][i];
             else {fail[tr[x][i]]=tr[fail[x]][i];  q[++tail]=tr[x][i];}
         }
     }
 }
 void get()
 {
     for(int i=;i<=cnt;i++)
     {
         if(end[i]) a.p[i][i]=;
         else for(int x=;x<m;x++) a.p[i][tr[i][x]]+=chty[x];
     }
 }
 inline node operator *(node &x,node &y)
 {
     node z;
     for(int i=;i<=cnt;i++)
         for(int j=;j<=cnt;j++)
             for(int k=;k<=cnt;k++)
                 z.p[i][j]+=x.p[i][k]*y.p[k][j];
     return z;
 }
 int main()
 {
     //freopen("cin.in","r",stdin);
     //freopen("cout.out","w",stdout);
     n=read();  l=read();  m=read();
     for(int i=;i<m;i++)  chty[i]=(double)read()/read();
     for(int i=;i<=n;i++)  {scanf("%s",ch+);  insert();}
     build();
     get();
     for(int i=;i<=;i++)  a=a*a;
     for(int i=;i<=n;i++)  printf("%.2lf\n",(double)a.p[][pos[i]]);
     return ;
 }