UVA 11468-Substring(AC自动机+概率dp)

UVA11468-一个AC自动机+概率dp的题目；

 

题目链接：http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=26&page=show_problem&problem=2463

 

题意：给出一些字符和对应的选择几率，随机选择L次将得到一个长度为L的随机字符串S(每次都是独立且随机)，给出K个模板串，计算S不包含任何一个串的概率(即任何一个模板串都不是S的连续子串)；

 

思路：构造出AC自动机之后，没随机生成一个字母，相当于在自动机内部随即走了一步，将所有的单词结点标记为“禁止”，本题就是求在AC自动机中从0处走L步，不进入任何“禁止结点”的概率。设dp[i][j]为当前在结点i，还要走j步，不碰到任何禁止结点的概率，由全概率公式可以得出下面的记忆化搜索过程。

 

/****************************
* author:crazy_石头
* date：2014/01/14
* algorithm:Aho-Corasick自动机+概率dp
* Pro：UVA11468 - Substring
*****************************/
#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
        
          using namespace std; #define INF 1<<29 #define eps 1e-8 #define A system("pause") #define rep(i,h,n) for(int i=(h);i<(n);i++) #define ms(a,b) memset((a),(b),sizeof(a)) const int maxm=150+10; const int maxn=1000000+10; const int maxlen=2000000+5; const int MAX_NODE=444; const int SIGMA_SIZE=62; double p[SIGMA_SIZE]; struct ACAutomation { int sz; int fail[MAX_NODE]; int ch[MAX_NODE][SIGMA_SIZE]; int val[MAX_NODE]; double dp[444][200]; inline int idx(char c) { if(c>='0'&&c<='9')return c-'0'; if(c>='A'&&c<='Z')return c-'A'+10; if(c>='a'&&c<='z')return c-'a'+36; } inline void init() { sz=1; memset(ch[0],0,sizeof(ch[0])); val[0]=0; memset(dp,-1,sizeof(dp)); } inline void insert(const char *s) { int u=0; for(int i=0; s[i]!='\0'; i++) { int v=idx(s[i]); if(!ch[u][v]) { memset(ch[sz],0,sizeof(ch[sz])); val[sz]=0; ch[u][v]=sz++; } u=ch[u][v]; } val[u]=1; } inline void construct() { fail[0]=0; queue 
         
           q; for(int c=0; c 
          
            =0)return dp[u][l]; double &ans=dp[u][l]; ans=0.0; for(int c=0;c