题目传送门

Description

周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数多谁胜利。

大家纷纷觉得这个游戏非常符合同学们的特色,但只是扔硬币实在是太单调了。

同学们觉得要加强趣味性,所以要找一个同学扔很多很多次硬币,其他同学记录下正反面情况。

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

但扔到什么时候停止呢?大家提议,选出 $ n $ 个同学, 每个同学猜一个长度为 $ m $ 的序列,当某一个同学猜的序列在硬币序列中出现时,就不再扔硬币了,并且这个同学胜利。为了保证只有一个同学胜利,同学们猜的 $ n $ 个序列两两不同。

很快,$ n $个同学猜好序列,然后进入了紧张而又刺激的扔硬币环节。你想知道,如果硬币正反面朝上的概率相同,每个同学胜利的概率是多少。

对于 $ 10% $ 的数据,$ 1 \leq n, m \leq 3 $;

对于 $ 40% $ 的数据,$ 1 \leq n, m \leq 18 $;

对于另外 $ 20%$ 的数据,$ n = 2 $;

对于 $ 100% $ 的数据,$ 1 \leq n, m \leq 300 $。

Solution

考虑使用概率生成函数。

我们设 \([x^n]G(x)\) 表示在 \(n\) 次操作后仍未结束的概率,\([x^n]F_i(x)\) 表示在 \(n\) 次操作后第 \(i\) 个人获胜的概率,\(a_{i,j,k}\) 表示第 \(i\) 个人的 \([1,k]\) 与第 \(j\) 个人的 \([m-k+1,m]\) 是否相同。

那我们可以得到:

\[G(x)+\sum_{i=1}^{n} F_i(x)=1+x·G(x)
\]
\[G(x)(\frac{1}{2})^mx^m=\sum_{j=1}^{n}\sum_{k=1}^{m}a_{i,j,k}(\frac{1}{2})^{m-k}x^{m-k}F_j(x)
\]

因为你最后只需要求出 \(F_i(1)\) ,所以可以把 \(x\) 都当成 \(1\) ,然后第一个式子就变为了:

\[\sum_{i=1}^{n} F_i(1)=1
\]

(似乎很显然的不需要第一个式子)

然后你就有 \(n+1\) 个方程,可以直接高斯消元了。

Code

#include <bits/stdc++.h>

using namespace std;

#define Int register int

#define MAXN 305

template <typename T> inline void read (T &t){t = 0;char c = getchar();int f = 1;while (c < '0' || c > '9'){if (c == '-') f = -f;c = getchar();}while (c >= '0' && c <= '9'){t = (t << 3) + (t << 1) + c - '0';c = getchar();} t *= f;}

template <typename T,typename ... Args> inline void read (T &t,Args&... args){read (t);read (args...);}

template <typename T> inline void write (T x){if (x < 0){x = -x;putchar ('-');}if (x > 9) write (x / 10);putchar (x % 10 + '0');}

template <typename T> inline void chkmax (T &a,T b){a = max (a,b);}

template <typename T> inline void chkmin (T &a,T b){a = min (a,b);}

int n,m;

char S[MAXN][MAXN];

bool a[MAXN][MAXN][MAXN];

double mat[MAXN][MAXN],pw2[MAXN],ans[MAXN];

#define seed 1000000007

#define ull unsigned long long

ull pw[MAXN],pre[MAXN][MAXN],suf[MAXN][MAXN];

void gauss (){

	int up = n + 1;

	for (Int i = 1;i <= up;++ i){

		int ind = i;

		for (Int j = i;j <= up;++ j) if (abs (mat[j][i]) >= 1e-7){ind = j;break;}

		if (ind ^ i) swap (mat[i],mat[ind]);

		for (Int j = i + 1;j <= up;++ j){

			double det = mat[j][i] / mat[i][i];

			for (Int k = i;k <= up + 1;++ k) mat[j][k] -= mat[i][k] * det;

		}

	}

	for (Int i = up;i >= 0;-- i){

		for (Int j = i + 1;j <= up;++ j) mat[i][up + 1] -= mat[i][j] * ans[j];

		ans[i] = mat[i][up + 1] / mat[i][i];

	}

}

signed main (){

	read (n,m),pw[0] = pw2[0] = 1;

	for (Int i = 1;i <= m;++ i) pw[i] = pw[i - 1] * seed,pw2[i] = pw2[i - 1] * 2;

	for (Int i = 1;i <= n;++ i) scanf ("%s",S[i] + 1);

	for (Int i = 1;i <= n;++ i){

		for (Int j = 1;j <= m;++ j) pre[i][j] = pre[i][j - 1] * seed + (S[i][j] == 'T');

		for (Int j = m;j >= 1;-- j) suf[i][j] = suf[i][j + 1] + (S[i][j] == 'T') * pw[m - j];

	}

	for (Int i = 1;i <= n;++ i) 

		for (Int j = 1;j <= n;++ j)

			for (Int k = 1;k <= m;++ k)

				a[i][j][k] = (pre[i][k] == suf[j][m - k + 1]);

	for (Int i = 1;i <= n;++ i){

		for (Int j = 1;j <= n;++ j)

			for (Int k = 1;k <= m;++ k)

				if (a[i][j][k]) mat[i][j] += pw2[k];

		mat[i][n + 1] = -1;

	}

	for (Int i = 1;i <= n;++ i) mat[n + 1][i] = 1;mat[n + 1][n + 2] = 1;

	gauss ();for (Int i = 1;i <= n;++ i) printf ("%.6f\n",ans[i]);

	return 0;

}

题解 「SDOI2017」硬币游戏的更多相关文章

  1. @loj - 2004@ 「SDOI2017」硬币游戏

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 周末同学们非常无聊,有人提议,咱们扔硬币玩吧,谁扔的硬币正面次数 ...

  2. 【LOJ 2004】「SDOI2017」硬币游戏

    LOJ 2004 100pts 首先我们肯定要建AC自动机的.. 那么这题就肯定是个AC自动机上\(dp\). 所以想想状态. 首先如果我们把状态设成这样行不行: \(dp(i)\)表示匹配到了i节点 ...

  3. 「SDOI2017」硬币游戏

    题目链接 问题分析 首先一个显然的做法就是建出AC自动机,然后高斯消元.但是这样的复杂度是\(O(n^3m^3)\)的. 我们发现其实只需要求AC自动机上\(n\)个状态的概率,而其余的概率是没有用的 ...

  4. 【LOJ】#2067. 「SDOI2016」硬币游戏

    题解 c一样的就是一个独立的游戏 我们对于2和3的指数 sg[i][j] 表示\(c \cdot 2^i \cdot 3^j\)的棋子,只有这个硬币是反面,翻转的硬币是正面的sg值 枚举sg函数所有可 ...

  5. 题解 「HDU6403」卡片游戏

    link Description 桌面上摊开着一些卡牌,这是她平时很爱玩的一个游戏.如今卡牌还在,她却不在我身边.不知不觉,我翻开了卡牌,回忆起了当时一起玩卡牌的那段时间. 每张卡牌的正面与反面都各有 ...

  6. 「SDOI2017」树点涂色 解题报告

    「SDOI2017」树点涂色 我sb的不行了 其实一开始有一个类似动态dp的想法 每个点维护到lct树上到最浅点的颜色段数,然后维护一个\(mx_{0,1}\)也就是是否用虚儿子的最大颜色 用个set ...

  7. LibreOJ 2003. 「SDOI2017」新生舞会 基础01分数规划 最大权匹配

    #2003. 「SDOI2017」新生舞会 内存限制:256 MiB时间限制:1500 ms标准输入输出 题目类型:传统评测方式:文本比较 上传者: 匿名 提交提交记录统计讨论测试数据   题目描述 ...

  8. AC日记——「SDOI2017」序列计数 LibreOJ 2002

    「SDOI2017」序列计数 思路: 矩阵快速幂: 代码: #include <bits/stdc++.h> using namespace std; #define mod 201704 ...

  9. 【LOJ】#2269. 「SDOI2017」切树游戏

    题解 把所有的数组一开始就FWT好然后再IFWT回去可以减小常数 从13s跑到0.7s-- 可以参照immortalCO的论文,感受一下毒瘤的动态动态DP 就是用数据结构维护线性递推的矩阵的乘积 由于 ...

随机推荐

  1. redis《三》连接池配置参数

    参数 值 setTestWhileIdle() 在空闲时检查有效性 true setMinEvictableIdleTimeMillis() 连接最小空闲时间 1800000L setTimeBetw ...

  2. Layui引起的对前端的一次记录

    前言 首先会做这次记录,也是因为自己也是第一次去接触这个框架,以前总是听说,并没有去用过.这次出于实习的原因,去学习了一下Layui这个"面向后端开发者的框架".其次,此篇记录仅供 ...

  3. 【SpringMVC】文件上传与下载、拦截器、异常处理器

    文件下载 使用ResponseEntity实现下载文件的功能 index.html <!DOCTYPE html> <html lang="en" xmlns:t ...

  4. AOP联盟通知类型和Spring编写代理半自动

    一.cglib功能更强大 二.Spring核心jar包 三.AOP联盟通知 三.代码实现Spring半自动代理 1.环绕通知的切面 2.bean.xml配置 3.创建bean容器,获取bean,即已经 ...

  5. PreparedStatement预编译的sql执行对象

    一.预编译,防sql注入 其中,设置参数值占位符索引从1开始:在由sql 连接对象创建 sql执行对象时候传入参数sql语句,在执行对象在执行方法时候就不用再传入sql语句: 数据库索引一般是从1开始 ...

  6. XML解析——Jsoup解析器

    一.Jsoup解析器快速入门案例 Docement对象,文本对象,包含着各个Dom树结构 1.引入Jsoup解析器的jar包放在lib文件夹下后,写java代码 其中, 二.Jsoup对象 1.Jso ...

  7. 数据导出生成Excel附件使用POI的HSSFWorkbook对象

    比较常用的实现Java导入.导出Excel的技术有两种Jakarta POI和Java Excel.Jakarta POI 是一套用于访问微软格式文档的Java API.Jakarta POI有很多组 ...

  8. MySQL实战45讲(01--05)-笔记

    目录 MySQL复习 01 | 基础架构:一条SQL查询语句是如何执行的? 连接器 查询缓存 分析器 优化器 执行器 02 | 日志系统:一条SQL更新语句是如何执行的? 重要的日志模块:redo l ...

  9. 【转】shell脚本中echo显示内容带颜色的实现方法

    shell脚本中echo显示内容带颜色的实现方法 shell脚本里使用echo输出颜色  

  10. Php实现简易购物商城系统

    实现功能: 1.系统功能模块包括: 1)登陆注册模块 包括验证码.找回密码.注册模块中要使用Ajax判断用户名是否已经存在,使用正则表达式判断电子邮件.手机号和用户密码的格式是否合法. 2)用户管理模 ...