[BZOJ 4820] [SDOI2017] 硬币游戏(高斯消元+概率论+字符串hash)

题面

扔很多次硬币后,用H表示正面朝上,用T表示反面朝上,会得到一个硬币序列。比如HTT表示第一次正面朝上,后两次反面朝上。

选出n个同学,每个同学猜一个长度为m的序列,当某一个同学猜的序列在硬币序列中出现时(匹配时的序列必须连续),就不再扔硬币了,并且这个同学胜利。猜的n个序列两两不同。

假设硬币正反面朝上的概率相同,求每个同学胜利的概率。

\(n \leq 300\)

分析

(注意,本题中不区分序列和串,因为只讨论连续的情况)

设与所有猜的序列都不匹配的串S期望长度是H(这个定义比较玄学,虽然是无限抛下去的,但是理解的时候当成有限长度抛出这种不匹配的情况比较好理解)

引理: 结尾包含一个长度为L的确定的串但没有包含任意一个给定序列的概率为\(\frac{1}{2^L}\).

感性理解一下就好了,证明略.

那么在这样的S串后面加上第i个人猜的串\(A_i\),那么\(SA_i\)就有可能成为一个胜利的状态。根据引理,抛出\(SA_i\)的概率为\(\frac{H}{2^m}\)。但是,有可能提前抛出其他的串。如\(S=\mathrm{HTH},A_1=\mathrm{HTT},A_2=\mathrm{THT}\),那么当串为\(\mathrm{HTHT}\)时,在抛出\(A_1\)之前就会抛出\(A_2\).同样,如果S本来就包含\(A_i\)的一部分,也可能提前抛出\(A_i\)导致提前胜利。把提前胜利的情况从\(\frac{H}{2^m}\)里减掉,就可以得到第i个人在抛出\(SA_i\)时胜利的概率。

注意到如果提前抛出串\(A_j\),那么\(A_i\)的某个前缀和\(A_j\)的某个后缀一定相等(可能不止1个),比如HTT和THT的公共前后缀为TH。

设第\(i\)个人胜利的概率为\(p_i\),\(pre(S,i)\)表示串S长度为i的前缀,\(suf(S,i)\)表示串S长度为i的后缀。那么有

\[p_i=\frac{H}{2^m}-\sum_{j=1}^n \sum_{k=1}^{m-[i=j]} [pre(A_i,k)=suf(A_j,k)]p_j \frac{1}{2^{m-k}}
\]

解释一下\(\sum_{k=1}^{m-[i=j]} [pre(A_i,k)=suf(A_j,k)]\frac{1}{2^{m-k}}\)

我们枚举公共前后缀的长度k,当k确定时,\(pre(A_i,k)\)确定,那么后面长度为\(m-k\)的后缀也确定了。根据引理,这样的概率为\(\frac{1}{2^{m-k}}\)。

\(m-[i=j]\)的意义是, 当\(i=j\)时长度<=m-1时才是提前结束,如果长度为m,就是胜利的状态了。当\(k=m\)时,\([pre(A_i,k)=suf(A_i,k)]p_i \frac{1}{2^{m-k}}=p_i\),因此可以把左边的\(p_i\)去掉,然后把循环上界改成\(m\)

\[\frac{H}{2^m}-\sum_{j=1}^n \sum_{k=1}^{m} [pre(A_i,k)=suf(A_j,k)]p_j \frac{1}{2^{m-k}}=0
\]

移项,

\[\sum_{j=1}^n \sum_{k=1}^m [pre(A_i,k)=suf(A_j,k)]p_j\frac{1}{2^{m-k}}-\frac{H}{2^m}=0
\]

对于每个i,我们都可以得到这样一个方程。现在我们有\(p_1,p_2 \dots p_n,H\)共n+1个未知数,有n个方程.又因为\(\sum_{i=1}^n p_i=1\),我们就有n+1个方程了,直接高斯消元即可.suf和pre可以用哈希求。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define eps 1e-10
#define maxn 300
#define maxm 300
using namespace std;
typedef long long ll;
typedef double db;
int n,m;
char a[maxn+5][maxm+5]; const ll seed=2;
const ll mod=998244353;
ll hs[maxn+5][maxm+5];
ll pows[maxn+5]; db pow2[maxn+5];
void ini_hash(){
pow2[0]=1;
for(int i=1;i<=m;i++) pow2[i]=pow2[i-1]*0.5;
pows[0]=1;
for(int i=1;i<=m;i++) pows[i]=pows[i-1]*seed%mod;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++) hs[i][j]=(hs[i][j-1]*seed+(a[i][j]=='T'))%mod;
}
}
db calc(int x,int y){
db ans=0;
for(int i=1;i<=m;i++){
if(hs[x][i]==(hs[y][m]-hs[y][m-i]*pows[i]%mod+mod)%mod) ans+=pow2[m-i];
}
return ans;
} db mat[maxn+5][maxm+5];
void gauss(int n,int m){
for(int i=1;i<=n;i++){
int id=i;
for(int j=i+1;j<=n;j++){
if(mat[j][i]>mat[id][i]) id=j;//把系数最大的行j交换到第i行
}
for(int k=1;k<=m;k++) swap(mat[i][k],mat[id][k]);
for(int j=1;j<=n;j++){
if(i==j) continue;
db r=mat[j][i]/mat[i][i];//把第j行第i个未知数的系数消成0
for(int k=i;k<=m;k++) mat[j][k]-=mat[i][k]*r;
}
}
for(int i=1;i<=n;i++){
mat[i][m]/=mat[i][i];
}
} int main(){
//#ifdef LOCAL
// freopen("game5.in","r",stdin);
//#endif
scanf("%d %d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%s",a[i]+1);
}
ini_hash();
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
mat[i][j]=calc(i,j); //1/2^(m-a)
}
mat[i][n+1]=-pow2[m];//H也是未知数
mat[i][n+2]=0;
}
for(int i=1;i<=n;i++) mat[n+1][i]=1;//sum(p[i])=1
mat[n+1][n+2]=1;
//#ifdef DEBUG
for(int i=1;i<=n+1;i++){
for(int j=1;j<=n+2;j++) printf("%.4f ",mat[i][j]);
printf("\n");
}
//#endif
gauss(n+1,n+2);
for(int i=1;i<=n;i++){
printf("%.10lf\n",mat[i][n+2]);
}
}

[BZOJ 4820] [SDOI2017] 硬币游戏(高斯消元+概率论+字符串hash)的更多相关文章

  1. [Sdoi2017]硬币游戏 [高斯消元 KMP]

    [Sdoi2017]硬币游戏 题意:硬币序列,H T等概率出现,\(n \le 300\)个人猜了一个长为$ m \le 300$的字符串,出现即获胜游戏结束.求每个人获胜概率 考场用了[1444: ...

  2. [BZOJ4820][SDOI2017]硬币游戏(高斯消元+KMP)

    比较神的一道题,正解比较难以理解. 首先不难得出一个(nm)^3的算法,对所有串建AC自动机,将在每个点停止的概率作为未知数做高斯消元即可. 可以证明,AC自动机上所有不是模式串终止节点的点可以看成一 ...

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

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

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

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

  5. BZOJ 4820 [Sdoi2017]硬币游戏 ——期望DP 高斯消元

    做法太神了,理解不了. 自己想到的是建出AC自动机然后建出矩阵然后求逆计算,感觉可以过$40%$ 用一个状态$N$表示任意一个位置没有匹配成功的概率和. 每种匹配不成功的情况都是等价的. 然后我们强制 ...

  6. bzoj 4820: [Sdoi2017]硬币游戏【kmp+高斯消元】

    有点神,按照1444的做法肯定会挂 注意到它的概率是相同的,所以可以简化状态 详见http://www.cnblogs.com/candy99/p/6701221.html https://www.c ...

  7. BZOJ 4820 [SDOI2017] 硬币游戏

    Description 周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数多谁胜利.大家纷纷觉得这个游戏非常符合同学们的特色,但只是扔硬币实在是太单调了.同学们觉得要加强趣味性,所以要找 ...

  8. 4820: [Sdoi2017]硬币游戏

    4820: [Sdoi2017]硬币游戏 链接 分析: 期望dp+高斯消元. 首先可以建出AC自动机,Xi表示经过节点i的期望次数,然后高斯消元,这样点的个数太多,复杂度太大.但是AC自动机上末尾节点 ...

  9. BZOJ 3105: [cqoi2013]新Nim游戏 [高斯消元XOR 线性基]

    以后我也要用传送门! 题意:一些数,选择一个权值最大的异或和不为0的集合 终于有点明白线性基是什么了...等会再整理 求一个权值最大的线性无关子集 线性无关子集满足拟阵的性质,贪心选择权值最大的,用高 ...

随机推荐

  1. JavaWeb_客户端相对/绝对路径和服务器端路径

    客户端的绝对路径和相对路径 相对路径:相对与某个基准目录的路径,在同一根目录下各子目录文件之间的相互引用, 绝对路径:指目录下的绝对位置,直接到的目标位置 @charset "UTF-8&q ...

  2. nginx负载均衡 之集群概念与负载均衡

    集群介绍 为什么要用集群

  3. python3笔记十:python数据类型-Tuple元组

    一:学习内容 元组概念 元组创建.访问.删除 元组操作 元组方法 二:元组概念 1.本质:是一种有序集合 2.特点:与列表非常相似.一旦初始化就不能修改.使用小括号 三:元组创建 1.创建空元组 tu ...

  4. go GTK msys2

    1 下载安装msys2 MSYS2 (Minimal SYStem 2) 是一个MSYS的独立改写版本,主要用于 shell 命令行开发环境.同时它也是一个在Cygwin (POSIX 兼容性层) 和 ...

  5. SELECT list is not in GROUP BY clause and contains nonaggregated

    安装了mysql5.7,用group by 查询时抛出如下异常 SQLSTATE[42000]: Syntax error or access violation: 1055 Expression # ...

  6. 异步上传&预览图片-不压缩图片

    本例使用ajaxFileUpload异步上传预览图片 <bean id="multipartResolver" class="org.springframework ...

  7. Dubbo Admin 控制台

    下载项目源码 https://github.com/apache/incubator-dubbo-admin 解压后配置 zookeeper 地址,路径为 dubbo-admin-server/src ...

  8. yum源相关

    yum软件仓库默认配置文件/etc/yum.conf,此文件定义了yum在线下载的rpm包存放位置及下载后是否保存. [root@localhost ~]# head /etc/yum.conf[ma ...

  9. ajaxform和ajaxgird中添加数据

    ajaxform添加数据 ajaxform.setRecord(response.getAjaxDataWrap("dataWrapBill").getData()); ajaxg ...

  10. leetcode 105从前序与中序遍历序列构造二叉树

    方法一:直接使用复制的数据递归:O(n)时间,O(n)空间,不计算递归栈空间: /** * Definition for a binary tree node. * struct TreeNode { ...