题面

题目描述

小阳阳发明了一个有趣的游戏:有n个玩家,每一个玩家均有一个长度为 l 的字母序列,任何两个玩家的字母序列不同。共有m种不同的字母,所有的字母序列都由这m种字母构成,为了方便,我们取大写字母的前。个字母。例如。m =3 , l = 4 , ABAA , CBCA 为两个合法的字母序列。现在由小阳阳来操控一台神奇的机器,每个时刻机器会随机产生一个字母,其中第i种字母随机出来的概率为pi/qi,显然 sum(pi/qi)=1 。这样T个时刻后机器会产生一个长度为 T 的字母序列。如果某个时刻某个玩家发现自己的字母序列在机器产生的字母序列中出现了,“出现”的定义是玩家的字母序列是机器产生的字母序列中连续的一段,那么我们称这个玩家获胜,游戏结束。现在小阳阳感兴趣的一个问题是,每个玩家分别有多大的概率能获得这场游戏的胜利呢?

输入格式

第一行有三个正整数n,l,m表示有n个人,每个人的字母序列长度为l,共有m个字母。

接下来m行,每行有两个非负整数p,q,表示随机到第i个字母的概率为p/q(0<=p<=q<=10,(p,q)=1)。数据保证m个字母的随机概率之和为1。

接下来n行,每行有一个长度为l的字母序列,表示第i个人的字母序列。数据保证所有的字母为大写字母的前m个且没有两个字母序列完全相同。

输出格式

包含n行,每行包含一个实数,表示第i个人获胜的概率,输出结果四舍五入到两位小数。

样例输入1

3 2 2
1 2
1 2
AB
BA
AA

样例输出1

0.25
0.50
0.25

样例输入2

3 4 2
1 2
1 2
AABA
ABAA
BAAA

样例输出2

0.31
0.33
0.37

样例说明1

两种字母 A 和 B ,概率均为 1/2。若前两个字母为 AB , BA 或AA,均有一个人获胜,获胜概率为 1 / 4 ;若前两个字母为 BB ,那么之后随机到 BBA , BBBA , BBBB 入都一定是 BA 获胜。因此 BA 的获胜概率为 1/4 + 1/4 = 1/2 。

样例说明 2

三个人的获胜概率分别为 4/13 , 17/52 , 19/52 ,注意输出结果四舍五入到两位小数。、

数据范围

100%的数据保证, n , l, m≤ 10.

题解

涉及概率 / 期望的题, 无非就是概率转期望, 期望转概率, DP一下就好了.

比如说这一题, 我们只需要建立trie图, 求出每个节点的经过的期望次数即可.

注意, 当一个节点的儿子为NULL时, 应该把这个概率加到根节点上.

#include <cstdio>
#include <cstring>
#include <deque>
#include <algorithm> const int N = 10, M = 10, L = 10;
int n, l, m;
double p[M];
struct matrix
{
int n;
double a[N * L][N * L + 1];
inline matrix()
{
memset(a, 0, sizeof(a));
}
inline void gauss()
{
for(int i = 0; i < n; ++ i)
{
int p;
for(p = i; p < n && a[p][i] == 0; ++ p);
if(p ^ i)
for(int j = 0; j <= n; ++ j)
std::swap(a[i][j], a[p][j]);
for(int j = 0; j < n; ++ j)
if(j ^ i)
{
double tmp = a[j][i] / a[i][i];
for(int k = 0; k <= n; ++ k)
a[j][k] -= a[i][k] * tmp;
}
}
}
}A;
struct ACautomaton
{
int cnt;
struct node
{
node *suc[10], *fl;
int ed, id;
inline node(int _id)
{
for(int i = 0; i < 26; ++ i)
suc[i] = NULL;
ed = 0;
id = _id;
}
}*rt;
inline ACautomaton()
{
cnt = 0;
rt = new node(cnt ++);
rt->fl = rt;
}
inline int insert(char *str, int len, int id)
{
node *u = rt;
for(int i = 0; i < len; u = u->suc[str[i] - 'A'], ++ i)
if(u->suc[str[i] - 'A'] == NULL)
u->suc[str[i] - 'A'] = new node(cnt ++);
u->ed = 1;
return u->id;
}
inline void build()
{
std::deque<node*> que;
que.clear();
for(int i = 0; i < 26; ++ i)
if(rt->suc[i] != NULL)
rt->suc[i]->fl = rt, que.push_back(rt->suc[i]), A.a[rt->suc[i]->id][0] = p[i];
else if(i < m)
A.a[0][0] += p[i];
for(; ! que.empty(); que.pop_front())
{
node *u = que.front();
for(int i = 0; i < 26; ++ i)
if(u->suc[i] != NULL)
{
if(! u->ed)
A.a[u->suc[i]->id][u->id] = p[i];
node *p = u->fl;
for(; p != rt && p->suc[i] == NULL; p = p->fl);
u->suc[i]->fl = p->suc[i] == NULL ? p : p->suc[i];
que.push_back(u->suc[i]);
}
else
{
u->suc[i] = u->fl->suc[i];
if(u->suc[i] != NULL && ! u->ed)
A.a[u->suc[i]->id][u->id] = p[i];
else if(i < m && ! u->ed)
A.a[0][u->id] += p[i];
}
}
A.n = cnt;
A.a[0][cnt] -= 1;
for(int i = 0; i < cnt; ++ i)
A.a[i][i] += -1;
}
}ACA;
int main()
{ #ifndef ONLINE_JUDGE
freopen("BZOJ1444.in", "r", stdin);
freopen("BZOJ1444.out", "w", stdout);
#endif scanf("%d%d%d", &n, &l, &m);
for(int i = 0; i < m; ++ i)
{
int x, y;
scanf("%d%d\n", &x, &y);
p[i] = (double)x / y;
}
static int ed[N];
for(int i = 0; i < n; ++ i)
{
static char str[L];
scanf("%s", str);
ed[i] = ACA.insert(str, l, i);
}
ACA.build();
A.gauss();
for(int i = 0; i < n; ++ i)
printf("%.2lf\n", A.a[ed[i]][A.n] / A.a[ed[i]][ed[i]] == 0 ? 0 : A.a[ed[i]][A.n] / A.a[ed[i]][ed[i]]);
}

JSOI 2009 BZOJ 1444 有趣的游戏的更多相关文章

  1. BZOJ 1444 有趣的游戏(AC自动机+矩阵快速幂)

    真的是很有趣的游戏... 对每个单词构建好AC自动机后,由于单词都是相同长度的且不同,所以不会出现互相为子串的形式. 那么我们对AC自动机上的节点构建转移矩阵.对于每个单词末尾的节点.该节点的出边仅仅 ...

  2. BZOJ 1444:[JSOI2009]有趣的游戏

    BZOJ 1444:[JSOI2009]有趣的游戏 题目链接 首先我们建出Trie图,然后高斯消元. 我们设\(f_i\)表示经过第\(i\)个点的期望次数: \[ f_x=\sum i\cdot p ...

  3. BZOJ:4820: [Sdoi2017]硬币游戏&&BZOJ:1444: [Jsoi2009]有趣的游戏(高斯消元求概率)

    1444: [Jsoi2009]有趣的游戏 4820: [Sdoi2017]硬币游戏 这两道题都是关于不断随机生成字符后求出现给定字符串的概率的问题. 第一题数据范围较小,将串建成AC自动机以后,以A ...

  4. BZOJ 1444: [Jsoi2009]有趣的游戏 [AC自动机 高斯消元]

    1444: [Jsoi2009]有趣的游戏 题意:每种字母出现概率\(p_i\),有一些长度len的字符串,求他们出现的概率 套路DP的话,\(f[i][j]\) i个字符走到节点j的概率,建出转移矩 ...

  5. BZOJ 1444 [Jsoi2009]有趣的游戏 (AC自动机 + 概率DP + Gauss)

    1444: [Jsoi2009]有趣的游戏 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 1382  Solved: 498[Submit][Statu ...

  6. 1444: [Jsoi2009]有趣的游戏

    1444: [Jsoi2009]有趣的游戏 链接 分析: 如果一个点回到0号点,那么会使0号点的概率增加,而0号点的概率本来是1,不能增加,所以这题用期望做. 设$x_i$表示经过i的期望次数,然后初 ...

  7. BZOJ.4820.[SDOI2017]硬币游戏(思路 高斯消元 哈希/AC自动机/KMP)

    BZOJ 洛谷 建出AC自动机,每个点向两个儿子连边,可以得到一张有向图.参照 [SDOI2012]走迷宫 可以得到一个\(Tarjan\)+高斯消元的\(O((nm)^3)\)的做法.(理论有\(6 ...

  8. 【BZOJ1444】[JSOI2009]有趣的游戏(高斯消元,AC自动机)

    [BZOJ1444][JSOI2009]有趣的游戏(高斯消元,AC自动机) 题面 BZOJ 题解 先把\(AC\)自动机构建出来,最好构成\(Trie\)图.然后这样子显然是在一个有向图中有一堆概率的 ...

  9. 【bzoj1444】[Jsoi2009]有趣的游戏

    1444: [Jsoi2009]有趣的游戏 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 1007  Solved: 334[Submit][Statu ...

随机推荐

  1. Android 自定义debug.keystore

    场景分析: 有时候,我们要使用第三方的服务,需要提供自己的包名以及keystore的sha1值,比如微信支付,百度地图,都需要包名和keystore的sha1值作为唯一标识.这时候我们测试的时候,如果 ...

  2. “帮你APP”团队冲刺2

    1.整个项目预期的任务量 (任务量 = 所有工作的预期时间)和 目前已经花的时间 (所有记录的 ‘已经花费的时间’),还剩余的时间(所有工作的 ‘剩余时间’) : 所有工作的预期时间:88h 目前已经 ...

  3. Linux之匿名FTP服务器搭建

    FTP(File Transfer Protocol)是在服务器与客户端进行文件传输的一种传输协议.本次介绍的是vsftpd的软件体验ftp服务. FTP服务器默认情况下依据用户登录情况分为三种不同的 ...

  4. Python 连接数据库失败

    什么是 PyMySQL? PyMySQL 是在 Python3.x 版本中用于连接 MySQL 服务器的一个库,Python2中则使用mysqldb. PyMySQL 遵循 Python 数据库 AP ...

  5. 10 Java 对象的内存布局

    Java 创建对象的方式 1:new 语句和反射机制创建.该方式会调用类的构造器,同时满足诸多约束.如果一个类没有构造器的话,Java 编译器会自动添加一个无参数的构造器.子类的构造器需要调用父类的构 ...

  6. [python][django学习篇][5]选择数据库版本(默认SQLite3) 与操作数据库

    推荐学习博客:http://zmrenwu.com/post/6/ 选择数据库版本(SQLite3) 如果想选择MySQL等版本数据库,请先安装MySQL并且安装python mysql驱动,这里不做 ...

  7. 【转载】zookeeper使用和原理探究(一)

    最近开始看到一些公司在使用zookeeper,本身对此了解的很少,这里看到一篇非常好的文章,因此转载 原贴地址:http://www.blogjava.net/BucketLi/archive/201 ...

  8. Lua 语法要点

    table 默认键值都是从1开始 table array = { "A", "B" } array2 = array array[] = "D&quo ...

  9. sass和postcss

    sass是css预处理器 需要安装node-sass支持 核心是c++编写 集成 sass-loader 把scss装换成css css-loader 找出@import和url()导入的语法,告诉w ...

  10. Threadlocal_笔记

    参考:https://www.jianshu.com/p/377bb840802f https://www.cnblogs.com/dreamroute/p/5034726.html ThreadLo ...