【BZOJ4820】【SDOI2017】硬币游戏
Description
Solution
设当前走出了一个不匹配任何字符串的串\(S\)。
若在\(S\)后随机增添\(m\)个字符,单看这些字符而言,这\(m\)个字符匹配到每个玩家的字符串的概率是相同的,记为\(P\)。
问题在于,对于每个字符串来说,并不是所有情况下一定要通过新增添\(m\)个字符才能匹配到自己,有可能加到中途时,就已经与\(S\)的某个后缀组成了自己,又或者是与\(S\)的某个后缀组成了别的字符串,早该停止了。
但是,对于每个串,不管每种情况中途该不该停下,我们计算出每种情况继续增加满\(m\)个字符并匹配到自己的概率,它们的概率之和还是\(P\)。
记每个人\(i\)成功被匹配到的概率是\(f_i\)(答案的定义)。
现在枚举对于一个人\(i\),在新加\(m\)个字符尝试匹配自己时,所有中途应该停下的情况。
枚举另一个人\(j\),如果\(j\)的\(len\)后缀与\(i\)的\(len\)前缀相同,则配合\(S\)的随机性,可能出现了这种情况:
这时候早就该停了,但为了凑齐\(P\),要计算在这种情况下继续匹配完全时所需的概率。继续匹配完\(i\)的子串,则还需要\((\frac 1 2)^{m-len}\)的概率。因此,这种情况对总和的贡献有\((\frac 1 2)^{m-len}f_j\)。当然,\(i\)和\(j\)之间不止有一个\(len\)满足条件,应找出所有符合描述的\(len\),累加\(f_j\)的贡献系数,记最终\(f_j\)贡献系数为\(a_j=\sum(\frac 1 2)^{m-len}\)。
其实\(j\)可以等于\(i\),这代表着提前匹配到自己的情况。
我们可以列出等式:
\]
总的来说,每一个人获胜的概率之和应该是1,因此有等式
\]
算上\(P\),我们拿到了一个\(n+1\)个未知数的\(n+1\)条方程,高斯消元解决即可,尽管我们并不需要知道\(P\)的具体取值。
Code
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
const int N=305;
int n,m;
char str[N][N];
double a[N][N],x[N];
int nex[N*2];
double mi2[N];
double kmp(int x,int y){
static char b[N*2];
for(int i=1;i<=m;i++) b[i]=str[x][i],b[m+i]=str[y][i];
nex[1]=0;
for(int i=2,j;i<=m*2;i++){
j=nex[i-1];
while(j&&b[j+1]!=b[i]) j=nex[j];
if(b[j+1]==b[i]) nex[i]=j+1;
else nex[i]=0;
}
double res=0;
for(int i=m*2;i;i=nex[i])
if(i<=m) res+=mi2[m-i];
return res;
}
void fill_matrix(){
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++)
a[i][j]=kmp(i,j);
a[i][n+1]=-1;
}
for(int i=1;i<=n;i++) a[n+1][i]=1;
a[n+1][n+2]=1;
}
void gaussian(int n){
int best;
for(int i=1;i<=n;i++){
best=i;
for(int j=i+1;j<=n;j++)
if(fabs(a[j][i])>fabs(a[best][i])) best=j;
if(best!=i)
for(int j=i;j<=n+1;j++) swap(a[i][j],a[best][j]);
for(int j=i+1;j<=n;j++){
double t=a[j][i]/a[i][i];
for(int k=i;k<=n+1;k++)
a[j][k]-=a[i][k]*t;
}
}
for(int i=n;i>=1;i--){
for(int j=i+1;j<=n;j++) a[i][n+1]-=a[i][j]*x[j];
x[i]=a[i][n+1]/a[i][i];
}
}
int main(){
scanf("%d%d",&n,&m);
mi2[0]=1;
for(int i=1;i<=m;i++) mi2[i]=mi2[i-1]*0.50000000000;
for(int i=1;i<=n;i++) scanf("%s",str[i]+1);
fill_matrix();
gaussian(n+1);
for(int i=1;i<=n;i++) printf("%.10lf\n",x[i]);
return 0;
}
【BZOJ4820】【SDOI2017】硬币游戏的更多相关文章
- BZOJ4820 Sdoi2017 硬币游戏 【概率期望】【高斯消元】【KMP】*
BZOJ4820 Sdoi2017 硬币游戏 Description 周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数多谁胜利.大家纷纷觉得这个游戏非常符合同学们的特色,但只是扔硬币实 ...
- [bzoj4820][Sdoi2017]硬币游戏
来自FallDream的博客,未经允许,请勿转载,谢谢. 周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数多谁胜利.大家纷纷觉得这个游戏非常符合同学们的特色,但只是扔硬币实在是太单调了 ...
- [BZOJ4820][SDOI2017]硬币游戏(高斯消元+KMP)
比较神的一道题,正解比较难以理解. 首先不难得出一个(nm)^3的算法,对所有串建AC自动机,将在每个点停止的概率作为未知数做高斯消元即可. 可以证明,AC自动机上所有不是模式串终止节点的点可以看成一 ...
- BZOJ4820 SDOI2017硬币游戏(概率期望+高斯消元+kmp)
容易想到的做法是建出AC自动机,高斯消元.然而自动机上节点数量是nm的. 注意到我们要求的变量只有n个,考虑将其他不用求的节点合并为一个变量.这个变量即表示随机生成一个串,其不包含任何一个模板串的概率 ...
- 【BZOJ4820】[SDOI2017]硬币游戏(高斯消元)
[BZOJ4820][SDOI2017]硬币游戏(高斯消元) 题面 BZOJ 洛谷 题解 第一眼的感觉就是构\(AC\)自动机之后直接高斯消元算概率,这样子似乎就是\(BZOJ1444\)了.然而点数 ...
- 【BZOJ4820】[Sdoi2017]硬币游戏 AC自动机+概率DP+高斯消元
[BZOJ4820][Sdoi2017]硬币游戏 Description 周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数多谁胜利.大家纷纷觉得这个游戏非常符合同学们的特色,但只是扔硬 ...
- BZOJ:4820: [Sdoi2017]硬币游戏&&BZOJ:1444: [Jsoi2009]有趣的游戏(高斯消元求概率)
1444: [Jsoi2009]有趣的游戏 4820: [Sdoi2017]硬币游戏 这两道题都是关于不断随机生成字符后求出现给定字符串的概率的问题. 第一题数据范围较小,将串建成AC自动机以后,以A ...
- [Sdoi2017]硬币游戏 [高斯消元 KMP]
[Sdoi2017]硬币游戏 题意:硬币序列,H T等概率出现,\(n \le 300\)个人猜了一个长为$ m \le 300$的字符串,出现即获胜游戏结束.求每个人获胜概率 考场用了[1444: ...
- 4820: [Sdoi2017]硬币游戏
4820: [Sdoi2017]硬币游戏 链接 分析: 期望dp+高斯消元. 首先可以建出AC自动机,Xi表示经过节点i的期望次数,然后高斯消元,这样点的个数太多,复杂度太大.但是AC自动机上末尾节点 ...
- [BZOJ 4820] [SDOI2017] 硬币游戏(高斯消元+概率论+字符串hash)
[BZOJ 4820] [SDOI2017] 硬币游戏(高斯消元+概率论+字符串hash) 题面 扔很多次硬币后,用H表示正面朝上,用T表示反面朝上,会得到一个硬币序列.比如HTT表示第一次正面朝上, ...
随机推荐
- 使用maven&&make-distribution.sh编译打包spark源码
1>基础环境准备: jdk1.8.0_101 maven 3.3.9scala2.11.8 安装好上述软件,配置好环境变量,并检查是否生效. 2>配置maven:intellij idea ...
- 2019第十届蓝桥杯 E题 迷宫
/*输入 30 50 01010101001011001001010110010110100100001000101010 00001000100000101010010000100000001001 ...
- Hyperledger Fabric(v1.2.0)代码分析1——channel创建
Hyperledger Fabric(v1.2.0)代码分析1--channel创建 0. e2e_cli Hyperledger Fabric提供了一个e2e的例子,该例中创建了一个基础的区块链网络 ...
- learning of a previous team
作为一个软件工程团队,离不开下面三个要素:支持,即分享.责任和合作. 分享是出色技术团队的另一个关键要素,它是团队的基石之一.只有通过分享,团队才有可能实现1+1 > 2这种效应,分享也是让团 ...
- 奔跑吧DKY——团队Scrum冲刺阶段-Day 3
今日完成任务 各个成员今日完成的任务(如果完成的任务为开发或测试任务,需给出对应的Github代码签入记录截图:如果完成的任务为调研任务,需给出对应的调研总结博客链接:如果完成的任务为学习技术任务,需 ...
- json反序列化对象
这个是同事研究的wcf中中根据type类型反序列化json的示例 /// <summary> /// json转对象 /// </summary> /// <param ...
- 索引超出了数组界限。 在 System.Collections.Generic.Dictionary`2.Resize
博问:Dictionary 超出了数组界限 异常: Exception type: IndexOutOfRangeException Exception message: 索引超出了数组界限. 在 S ...
- Week2-作业1
第一章:引用:如果一架民用飞机上有一个功能,用户使用它的概率是百万分之一,你还要做这个功能么? 选择之后,这个功能是什么呢?谜底是飞机的安全功能. 个人认为,飞机的安全功能这个 ...
- angularJS1笔记-(14)-自定义指令(scope)
index.html: <!DOCTYPE html> <html lang="en"> <head> <meta charset=&qu ...
- .NET 类库研究必备参考 扣丁格鲁
.NET 类库的强大让我们很轻松的解决常见问题,作为一个好专研的程序员,为了更上一层楼,研究CLR的基础类库实现是快速稳定的捷径. 一般场景下,采用 Reflector可以反射出.NET 的部分实现出 ...