[BZOJ2553][BeiJing2011]禁忌 dp+AC自动机+矩阵快速幂
2553: [BeiJing2011]禁忌
Time Limit: 20 Sec Memory Limit: 128 MBSec Special Judge
Submit: 1206 Solved: 607
[Submit][Status][Discuss]
Description
Magic Land上的人们总是提起那个传说:他们的祖先John在那个东方岛屿帮助Koishi与其姐姐Satori最终战平。而后,Koishi恢复了读心的能力……
如今,在John已经成为传说的时代,再次造访那座岛屿的人们却发现Koishi遇到了新麻烦。
这次她遇到了Flandre Scarlet——她拥有可以使用禁忌魔法而不会受到伤害的能力。
为了说明什么是禁忌魔法及其伤害,引入以下概念:
1.字母集A上的每个非空字符串对应了一个魔法。
其中A是包含了前alphabet个小写字母的集合。
2.有一个集合T,包含了N个字母集A上的字符串
T中的每一串称为一个禁忌串(Taboo string)
3.一个魔法,或等价地,其对应的串s因为包含禁忌而对使用者造成的伤害按以下方式确定:
把s分割成若干段,考虑其中是禁忌串的段的数目,不同的分割可能会有不同的数目,其最大值就是这个伤害。
由于拥有了读心的能力,Koishi总是随机地使用Flandre Scarlet的魔法,可以确定的是,她的魔法正好对应字母集A上所有长度为len的串。
但是,Flandre Scarlet所使用的一些魔法是带有禁忌的,由于其自身特性,她可以使用禁忌魔法而不受到伤害,而Koishi就不同了。可怜的Koishi每一次使用对方的魔法都面临着受到禁忌伤害的威胁。
你现在需要计算的是如果Koishi使用对方的每一个魔法的概率是均等的,那么每一次随机使用魔法所受到的禁忌伤害的期望值是多少。
Input
第一行包含三个正整数N、len、alphabet。
接下来N行,每行包含一个串Ti,表示禁忌串。
Output
一个非负实数,表示所受到禁忌伤害的期望值。
Sample Input
aa
abb
Sample Output
【样例1解释】
一共有2^4 = 16种不同的魔法。
需要注意的是“aabb”的禁忌伤害是1而不是2。
HINT
100%的数据中N ≤ 5,len ≤109,1 ≤ alphabet ≤ 26。
在所有数据中,有不少于40%的数据中:N = 1。
数据保证每个串Ti的长度不超过15,并且不是空串。
数据保证每个Ti均仅含有前alphabet个小写字母。
数据保证集合T中没有相同的元素,即对任意不同的i和j,有Ti≠Tj。
【评分方法】
对于每一组数据,如果没有得到正确的输出(TLE、MLE、RTE、输出格式错误等)得0分。
否则:设你的输出是YourAns,标准输出是StdAns:
记MaxEPS = max(1.0 , StdAns)×10-6
如果|YourAns – StdAns| ≤ MaxEPS则得10分,否则得0分。
即:你的答案需要保证相对误差或绝对误差不超过10-6。
Source
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define maxn 200
using namespace std;
int q[maxn],al,n,m;
struct mat {
long double f[maxn][maxn];
mat() {for(int i=;i<maxn;i++) for(int j=;j<maxn;j++) f[i][j]=;}
}Ans;
struct AC {
int ch[maxn][],fail[maxn],cnt,val[maxn];
void build(char *s) {
int now=,len=strlen(s+);
for(int i=;i<=len;i++) {
int to=s[i]-'a';
if(!ch[now][to]) ch[now][to]=++cnt;
now=ch[now][to];
}
val[now]=;
}
void getfail() {
int h=,t=;
for(int i=;i<al;i++) if(ch[][i]) q[t++]=ch[][i];
while(h!=t) {
int now=q[h++];if(h==) h=;
for(int i=;i<al;i++) {
int to=ch[now][i];
if(!to) {ch[now][i]=ch[fail[now]][i];continue;}
int tmp=fail[now];
while(tmp&&!ch[tmp][i]) tmp=fail[tmp];
fail[to]=ch[tmp][i];
val[to]|=val[fail[to]];
q[t++]=to;if(t==) t=;
}
}
}
}a;
mat mul(mat &x,mat &y) {
mat ans;
for(int i=;i<=a.cnt+;i++)
for(int j=;j<=a.cnt+;j++)
for(int k=;k<=a.cnt+;k++) ans.f[i][j]+=x.f[i][k]*y.f[k][j];
return ans;
}
mat power(mat x,int c) {
mat ans;
for(int i=;i<=a.cnt+;i++) ans.f[i][i]=;
while(c) {
if(c&) ans=mul(ans,x);
x=mul(x,x);c>>=;
}
return ans;
}
void build() {
long double hh=/(long double)al;
Ans.f[a.cnt+][a.cnt+]=;
for(int i=;i<=a.cnt;i++) {
if(a.val[i]) continue;
for(int j=;j<al;j++) {
int to=a.ch[i][j];
if(a.val[to]) {
Ans.f[a.cnt+][i]+=hh;
Ans.f[][i]+=hh;
}
else Ans.f[to][i]+=hh;
}
}
}
int main() {
scanf("%d%d%d",&n,&m,&al);
for(int i=;i<=n;i++) {
char s[];scanf("%s",s+);
a.build(s);
}
a.getfail();
build();
Ans=power(Ans,m);
/* for(int i=0;i<=a.cnt+1;i++) {
for(int j=0;j<=a.cnt+1;j++) cout<<Ans.f[i][j]<<' ';
cout<<endl;
}*/ printf("%.7Lf\n",Ans.f[a.cnt+][]);
}
[BZOJ2553][BeiJing2011]禁忌 dp+AC自动机+矩阵快速幂的更多相关文章
- BZOJ2553 Beijing2011禁忌(AC自动机+动态规划+矩阵快速幂+概率期望)
考虑对一个串如何分割能取得最大值.那么这是一个经典的线段覆盖问题,显然每次取右端点尽量靠前的串.于是可以把串放在AC自动机上跑,找到一个合法串后就记录并跳到根. 然后考虑dp.设f[i][j]表示前i ...
- HDU 2243考研路茫茫——单词情结 (AC自动机+矩阵快速幂)
背单词,始终是复习英语的重要环节.在荒废了3年大学生涯后,Lele也终于要开始背单词了. 一天,Lele在某本单词书上看到了一个根据词根来背单词的方法.比如"ab",放在单词前一般 ...
- HDU 2243 考研路茫茫――单词情结 ——(AC自动机+矩阵快速幂)
和前几天做的AC自动机类似. 思路简单但是代码200余行.. 假设solve_sub(i)表示长度为i的不含危险单词的总数. 最终答案为用总数(26^1+26^2+...+26^n)减去(solve_ ...
- POJ 2778 DNA Sequence(AC自动机 + 矩阵快速幂)题解
题意:给出m个模式串,要求你构造长度为n(n <= 2000000000)的主串,主串不包含模式串,问这样的主串有几个 思路:因为要不包含模式串,显然又是ac自动机.因为n很大,所以用dp不太好 ...
- POJ2778 DNA Sequence(AC自动机+矩阵快速幂)
题目给m个病毒串,问不包含病毒串的长度n的DNA片段有几个. 感觉这题好神,看了好久的题解. 所有病毒串构造一个AC自动机,这个AC自动机可以看作一张有向图,图上的每个顶点就是Trie树上的结点,每个 ...
- poj2778DNA Sequence (AC自动机+矩阵快速幂)
转载请注明出处: http://www.cnblogs.com/fraud/ ——by fraud DNA Sequence Time Limit: 1000MS Memory ...
- POJ2778(SummerTrainingDay10-B AC自动机+矩阵快速幂)
DNA Sequence Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 17160 Accepted: 6616 Des ...
- poj2778 ac自动机+矩阵快速幂
给m个子串,求长度为n的不包含子串的母串数,最直接的应该是暴搜,肯定tle,考虑用ac自动机 将子串建成字典树,通过next表来构造矩阵,然后用矩阵快速幂求长度为n的数量 邻接矩阵https://we ...
- POJ - 2778 ~ HDU - 2243 AC自动机+矩阵快速幂
这两题属于AC自动机的第二种套路通过矩阵快速幂求方案数. 题意:给m个病毒字符串,问长度为n的DNA片段有多少种没有包含病毒串的. 根据AC自动机的tire图,我们可以获得一个可达矩阵. 关于这题的t ...
随机推荐
- oracle ocp 052考试学习
1.数据字典存储在SYSTEM表空间中. 2.SYSAUX可以offline: SQL>alter tablespace sysaux offline; 3.SYSTEM和SYSAUX都是永久表 ...
- windows服务那些事
前一段时间由于项目需求,写了一个windows服务.下面总结如下: windows服务其实就是一些后台程序,和其他程序的主要区别是它运行于系统后台.微软公司为了方便我们自己定制我们的服务,提供了很多借 ...
- TCP的挥手协议和握手协议
三次握手协议:三次握手协议的主要过程是交互彼此之间的初始序列号,如果没有确认的ACK帧可以么?肯定是可以的 client A -------> server B client A 发送了自己的初 ...
- php 文件上传失败
使用OSX系统,在使用MAMP Pro作为虚拟服务器,并使用PHP作为后端语言进行文件上传,从临时文件夹拷贝文件的方法为 move_uploaded_file 代码如下: if($_FILES['fi ...
- Java 如何正确停止一个线程
自己在做实验性小项目的时候,发现自己遇到一个问题:如何控制线程的"死亡"? 首先,如何开启一个线程呢? 最简单的代码: public class Main { public sta ...
- 【题解】SDOI2014旅行
洛谷P3313 大概是一道树链剖分的裸题.可以看出如果不是查询相同宗教的这一点,就和普通的树链剖分毫无两样了.所以针对每一个宗教都单独开一棵线段树,变成单点修改+区间查询.只不过宗教数目很多,空间消耗 ...
- ZOJ 2314 Reactor Cooling | 无源汇可行流
题目: 无源汇可行流例题 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=1314 题解: 证明什么的就算了,下面给出一种建图方式 ...
- 【ZJ选讲·画山】
给出一张纸(N × M),你要在上面画山,但不能画出界(N,M<=100) Like this: 起点为(0,0),终点为(N,0) 给出w种线段画法(x,y),表示用了这种画法后,笔迹末 ...
- 【NOIP模拟赛】Evensgn 的债务 乱搞
biubiu~~~ 我们发现按照这道题的题意我们把一个个人的前后(欠钱,被欠钱)都缩一下,那么他对其他人没有影响,那么我们就可以依次缩完每个人,而且每个人最后的状态都是要买欠要么被欠,那么我们可以知道 ...
- Codeforces Round #328 (Div. 2) A
A. PawnChess time limit per test 1 second memory limit per test 256 megabytes input standard input o ...