UVA11468-一个AC自动机+概率dp的题目;
题意:给出一些字符和对应的选择几率,随机选择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