首先我们应该是枚举 L个位置上的每个字符来得到最终概率

然后AC自动机的作用就是为了判断你枚举的地方是否对应了单词节点,如果对应了,就肯定要不得

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int max_node = ;
const int N =;
const int node_num=;
int idx[N],k,n;
char str[][];
double prob[]; //这个地方一开始开小了,题目也没说n为多少。坑啊。。一直RE没找到原因
int vis[max_node][];
double d[max_node][];
struct ACtrie
{
int ch[max_node][];
int sz;
int val[max_node];
int fail[max_node];
//int last[max_node];
void init(){
memset(ch,,sizeof ch);
memset(val,,sizeof val);
memset(fail,,sizeof fail);
//memset(last,0,sizeof last);
sz=;
}
void insert(char* s,int v){
int u=;
for (int i=;s[i];i++){
int c=idx[s[i]];
if (!ch[u][c]){
memset(ch[sz],,sizeof ch[sz]);
val[sz]=;
ch[u][c]=sz++;
}
u=ch[u][c];
}
val[u]=;
}
void getfail(){
queue<int> q;
fail[]=;
for (int i=;i<node_num;i++)
if (ch[][i]){
fail[ch[][i]]=;
q.push(ch[][i]);
} while (!q.empty()){
int u=q.front();
q.pop();
for (int c=;c<node_num;c++){
int v=ch[u][c];
if (!v) {ch[u][c] = ch[fail[u]][c];continue;} //这个个普通AC自动机不一样,因为要一视同仁所有字符,当前字符在该子树如果找不到就可以转到其失败节点去,因为失败节点有的话,同样是满足条件的
q.push(v);
int p=fail[u];
while (p && !ch[p][c]) p=fail[p];
fail[v]=ch[p][c];
//last[v]=val[fail[p]]?fail[p] : last[fail[p]];
val[v]|=val[fail[v]];
}
}
}
} T;
double getprob(int u,int L)
{
if (L==) return 1.0;
if (vis[u][L]) return d[u][L];
vis[u][L]=;
d[u][L]=0.0;
for (int i=;i<n;i++){
if (!T.val[T.ch[u][i]]) d[u][L]+=prob[i]*getprob(T.ch[u][i],L-);
}
return d[u][L];
}
int main()
{
int t,L;
char cc[];
scanf("%d",&t);
int kase=;
while (t--){
T.init();
scanf("%d",&k);
for (int i=;i<k;i++)
scanf("%s",str[i]);
scanf("%d",&n);
for (int i=;i<n;i++){
scanf("%s%lf",&cc,&prob[i]);
idx[cc[]]=i;
}
for (int i=;i<k;i++){
T.insert(str[i],);
}
T.getfail();
memset(vis,,sizeof vis);
scanf("%d",&L);
double ans=getprob(,L);
printf("Case #%d: %.6lf\n",++kase,ans);
}
return ;
}

UVA 11468 AC 自动机的更多相关文章

  1. Uva 11468 AC自动机或运算

    AC自动机 UVa 11468 题意:给一些字符和各自出现的概率,在其中随机选择L次,形成长度为L的字符串S,给定K个模板串,求S不包含任意一个串的概率. 首先介绍改良版的AC自动机: 传统的AC自动 ...

  2. UVa 11468 (AC自动机 概率DP) Substring

    将K个模板串构成一个AC自动机,那些能匹配到的单词节点都称之为禁止节点. 然后问题就变成了在Tire树上走L步且不经过禁止节点的概率. 根据全概率公式用记忆化搜索求解. #include <cs ...

  3. UVA - 11468 (AC自动机+动态规划)

    建立AC自动机,把AC自动机当做一张图,在上面跑L个节点就行了. 参考了刘汝佳的代码,发现可能有一个潜在的Bug--如果模式串中出现了没有指定的字符,AC自动机可能会建立出错. 提供一组关于这个BUG ...

  4. Substring UVA - 11468 AC自动机+概率DP

    题意: 给出一些字符和各自对应的选择概率,随机选择L次后得到一个长度为L的随机字符串S. 给出K个模板串,计算S不包含任何一个模板串的概率 dp[i][j]表示走到AC自动机 i 这个节点 还需要走 ...

  5. UVA 11468 AC自动机入门题 记忆化概率dp+ac自动机

    /** 链接:https://vjudge.net/problem/UVA-11468 详见lrj训练指南P218 我的是反向求存在模板串的概率. dp[i][j]表示当前i位置选择字符,前面i-1个 ...

  6. uva 11468 AC自动机+概率DP

    #include<cstdio> #include<cstring> #include<queue> #include<cstdio> #include ...

  7. UVa 11019 (AC自动机 二维模式串匹配) Matrix Matcher

    就向书上说得那样,如果模式串P的第i行出现在文本串T的第r行第c列,则cnt[r-i][c]++; 还有个很棘手的问题就是模式串中可能会有相同的串,所以用repr[i]来记录第i个模式串P[i]第一次 ...

  8. AC自动机——Uva 11468 子串

    题目链接:http://vjudge.net/contest/142513#problem/A 题意:给出一些字符和各自对应的选择概率,随机选择L次后将得到一个长度为L的随机字符串S.给出K个模版串, ...

  9. 沉迷AC自动机无法自拔之:[UVA 11468] Substring

    图片加载可能有点慢,请跳过题面先看题解,谢谢 这个鬼题目,上一波套路好了 先用题目给的模板串建\(AC\)自动机,把单词结尾标记为 \(val=1\),然后在建好的\(AC\)自动机上跑 \(dp\) ...

随机推荐

  1. 10 Class文件结构

  2. Educational Codeforces Round 63 选做

    D. Beautiful Array 题意 给你一个长度为 \(n\) 的序列.你可以选择至多一个子段,将该子段所有数乘上给定常数 \(x\) .求操作后最大的最大子段和. 题解 考虑最大子段和的子段 ...

  3. python进阶 廖雪峰(慕课网)

    1.函数式编程 变量名可以指向函数,那么函数就可以通过一个变量传递给另一个函数或者变量. map()函数:接收一个函数 f 和一个 list,并通过把函数 f 依次作用在 list 的每个元素上,得到 ...

  4. 实验吧-杂项-pilot-logic、ROT-13变身了

    1.pilot-logic 题上说password藏在文件里,直接丢到Winhex里,搜索pass就拿到flag了. 有的大佬提供了另一种方法,题上说是一个磁盘文件,有一个处理磁盘文件的软件autop ...

  5. python-处理json、处理时间

    1.处理json:import json #json串就是字符串. d = { 'car1':{'color':'red','price':100,'count':50}, 'car2':{'colo ...

  6. Laravel框架的学习

    用xampp环境 1.Composer的安装 http://www.phpcomposer.com/ 下载Composer的安装exe php.ini中extension_dir.browscap路径 ...

  7. bzoj 2451 Uyuw's Concert

    裸的半平面交.感觉这些东西,纯属在考代码能力啊.. #include<cstdio> #include<algorithm> #include<cmath> #de ...

  8. 《方方格子》(WPS版) _v3.6.6.0

    <方方格子>(WPS版) 下载地址(b2a7) SHA1:35AE4D99B77613D9E2BAF912477DC74C5C2B8389 版本信息    发行版本 3.6.6.0    ...

  9. reduce()、filter()、map()、some()、every()、...展开属性

    reduce().filter().map().some().every()....展开属性   这些概念属于es5.es6中的语法,跟react+redux并没有什么联系,我们直接在https:// ...

  10. 121-PHP类成员函数(二)

    <?php class ren{ //定义人类 const name='ren'; public function classname(){ echo '这个类的名称:'.self::name; ...