P3706-[SDOI2017]硬币游戏【高斯消元,字符串hash】

正题

题目链接:https://www.luogu.com.cn/problem/P3706


题目大意

给出 nnn 个长度为 mmmH/TH/TH/T 串。

开始一个空序列,每次随机在后面加一个 H/TH/TH/T ,求每个串第一次出现的概率。

1≤n,m≤3001\leq n,m\leq 3001n,m300


解题思路

数据范围显然不能在AC自动机上高斯消元,所以得考虑别的方法。

考虑一个很妙的做法,设一个状态NNN表示目前还没有匹配完成的串,然后考虑串A=HHTA=HHTA=HHT和串B=THHB=THHB=THH,那么在NNN后面直接插入一个AAA的概率就是p(N+A)=p(N)×2−mp(N+A)=p(N)\times 2^{-m}p(N+A)=p(N)×2m

但是考虑到有可能NNN先拼出了BBB然后再拼出AAA,此时考虑BBB的后缀对应AAA前缀的有H,HHH,HHH,HH。那么就有
p(N)×2−m=p(N+A)=p(A)+p(B)×2−2+p(B)×2−1p(N)\times 2^{-m}=p(N+A)=p(A)+p(B)\times 2^{-2}+p(B)\times 2^{-1}p(N)×2m=p(N+A)=p(A)+p(B)×22+p(B)×21

这样不难发现对于别的串如果它的一些后缀是这个串的前缀那么就会产生一些概率,用字符串hashhashhash匹配即可。

然后会发现还是少了一个方程,最后一个就是所有串的概率和为111就好了。

这样就有n+1n+1n+1个方程了。

时间复杂度:O(n2m+n3)O(n^2m+n^3)O(n2m+n3)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ull unsigned long long
using namespace std;
const int N=310;
const ull g=131;
int n,m;
double a[N][N],b[N],pw[N];
ull h[N][N],p[N];char s[N];
ull geth(int x,int l,int r)
{return h[x][r]-h[x][l-1]*p[r-l+1];}
void Gauss(int n){for(int i=1;i<=n;i++){int z=i;for(int j=i+1;j<=n;j++)if(a[j][i]>a[z][i])z=i;swap(a[i],a[z]);double x=a[i][i];b[i]/=x;for(int j=i;j<=n;j++)a[i][j]/=x;for(int j=i+1;j<=n;j++){double rate=-a[j][i];for(int k=i;k<=n;k++)a[j][k]+=rate*a[i][k];b[j]+=rate*b[i];}}for(int i=n;i>=1;i--){for(int j=1;j<i;j++){b[j]-=a[j][i]*b[i];a[j][i]=0;}}return;
}
int main()
{scanf("%d%d",&n,&m);p[0]=1;pw[0]=1;for(int i=1;i<=m;i++)pw[i]=pw[i-1]*0.5,p[i]=p[i-1]*g;for(int i=1;i<=n;i++){scanf("%s",s+1);for(int j=1;j<=m;j++)h[i][j]=h[i][j-1]*g+s[j];}for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)for(int k=1;k<=m;k++)if(geth(i,1,k)==geth(j,m-k+1,m))a[i][j]+=pw[m-k];for(int i=1;i<=n;i++)a[i][n+1]=-pw[m],a[n+1][i]=1;b[n+1]=1;Gauss(n+1);for(int i=1;i<=n;i++)printf("%.12lf\n",b[i]);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318304.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

YBTOJ:字符串匹配(KMP)

文章目录题目描述解析代码题目描述 解析 看了题解。。。 这题的关键在于可以变换匹配的一个充要条件&#xff1a; 每个字符与前一个相同字符的距离相同 这个搞出来之后就可以以它为关键字进行KMP了 注意&#xff01; 当与前一个字符的距离超过匹配长度时&#xff0c;是没有意义…

Strange Birthday Party CodeForces - 1471C

Strange Birthday Party CodeForces - 1471C 题意&#xff1a; 我有n个朋友&#xff0c;商店有m种商品&#xff0c;这m种商品按序号价格从小到大排列&#xff0c;对于每一个朋友我给出一个序号k&#xff0c;我可以直接给朋友序号k的商品价格的金钱或给朋友买一个序号小于k的商…

TransactionScope事务处理方法介绍及.NET Core中的注意事项

作者&#xff1a;依乐祝原文链接&#xff1a;https://www.cnblogs.com/yilezhu/p/10170712.html今天在写CzarCms的UnitOfWork的使用使用到了这个TransactionScope事务&#xff0c;因此对它进行了相关资料的查阅并记录如下&#xff0c;希望对大伙在.NET Core中使用有所帮助。写在…

P4548-[CTSC2006]歌唱王国【概率生成函数,KMP】

正题 题目链接:https://www.luogu.com.cn/problem/P4548 题目大意 ttt次询问&#xff0c;给出一个长度为mmm的串SSS和一个空串TTT&#xff0c;每次在TTT后面随机加入1∼n1\sim n1∼n的字符&#xff0c;直到TTT中出现SSS为止&#xff0c;求期望次数。 1≤n≤105,t≤50,1≤m≤10…

YBTOJ:字符匹配(KMP)

文章目录题目描述解析代码题目描述 解析 显然应该要尝试套kmp的板子 关键是如何套 也就是那个判断匹配的条件是什么的问题 本题的关键是当kmp匹配时&#xff0c;匹配位之前的所有位大小关系的顺序都是匹配的&#xff0c;所以我们只需要看当前位即可 考虑对b预处理出3个数组&a…

历史上的今天(history)+ 勇者斗恶龙(dragon)

朋友们我来了&#xff0c;好久都没有更新了&#xff0c;手实在痒痒&#xff0c;不撸两道&#xff0c;内心过意不去 文章目录A&#xff1a;历史上的今天&#xff08;history&#xff09;题目题解代码实现B&#xff1a; 勇者斗恶龙&#xff08;dragon&#xff09;题目描述题解代…

Codeforces Round #694 Div. 2

Codeforces Round #694 Div. 2 CodeForces 1471A Strange Partition CodeForces 1471B Strange List CodeForces 1471C Strange Birthday Party 签到题 稍有难度 CodeForces 1471D Strange Definition CodeForces 1471E Strange Shuffle CodeForces 1471F Strange Housing

领域驱动设计,让程序员心中有码(四)

#领域驱动设计&#xff0c;让程序员心中有码&#xff08;四&#xff09;----------------------追忆有关分层的古老往事我一直认为&#xff0c;程序员也是艺术家&#xff0c;他们撰写的每一行代码&#xff0c;是献给这大好世界的优美诗篇。不同的人&#xff0c;写的代码也许风格…

YBTOJ洛谷P4551:最长异或路径(trie树)

洛谷传送门 文章目录题目描述解析代码题目描述 解析 本题关键就在于一点&#xff1a; 若把每个点的深度dep[i]定义为从根到节点边权的异或和 那么i到j的路径异或和可以表示为&#xff1a; dep[i] ^ dep[j] 首先要是i、j在不同子树上显然成立 如果他们在同一子树上&#xff0c;…

P5074-Eat the Trees【插头dp】

正题 题目链接:https://www.luogu.com.cn/problem/P5074 题目大意 给出一个nmn\times mnm的网格&#xff0c;有的必须铺线有的不能&#xff0c;铺成若干条闭合回路&#xff0c;求方案数。 1≤n,m≤121\leq n,m\leq 121≤n,m≤12 解题思路 考虑插头dpdpdp&#xff0c;因为可以…

初赛—算法复杂度

对于 \(T(n) a\times T(\frac{n}{b})c\times n^k\) 这样的递归关系&#xff0c;有这样的结论&#xff1a; \(a>b^k\) 则有 &#xff1a; \(T(n) O(c\times n^{\log_{b} a})\)\(ab^k\) 则有 &#xff1a; \(T(n) O(c\times n^k\times log n)\)\(a<b^k\) 则有 &#xff…

1414. 牛异或

题意&#xff1a; 区间长度为n的数组a&#xff0c;取连续一段子序列&#xff0c;问哪段子序列的异或和最大&#xff0c;如果存在多个这样的序列&#xff0c;那么选择序列末端整数对应的编号更小的那个序列。 如果仍然存在多个可选的序列&#xff0c;那么选择长度最短的那个序列…

基于Asp.Net Core打造轻量级内部服务治理RPC(一)

继之前的《Asp.Net Core Docker 搭建》文章末尾说过的&#xff0c;将陆续编写基于asp.net core 打造一个内部服务治理的rpc框架。不过前端时间较忙&#xff0c;所以搁置了一段时间。闲话不多说&#xff0c;下面就来讲讲为什么需要去做一个该框架&#xff0c;以及想法的来源和设…

P7046-「MCOI-03」诗韵【SAM,倍增,树状数组】

正题 题目链接:https://www.luogu.com.cn/problem/P7046 题目大意 给出一个长度为 nnn 的字符串&#xff0c;然后 mmm 次把它的一个子串加入集合。如果一个字符串在这个集合中作为字符串的后缀出现次数大于 kkk 那么这个字符串就会被计入贡献。 每次求计入贡献的字符串数和最…

YBTOJ:前缀询问(trie树)

文章目录题目描述解析代码题目描述 解析 &#xff08;没有做出来&#xff0c;这个ans的处理方式其实也不难想…qwq&#xff09; 考虑把T都作为模板串加入trie树 加入每个模板串自然就是按照i顺序的 所以我们在插入t的时候沿途标记一下 新出现的未标记的i的间隔就是当前的i与上…

DP专练1( [NOIP 2003]加分二叉树 + 太空梯 )

我们先慢慢来加分二叉树题目题解简单讲解前序//中序//后序遍历代码实现太空梯题目题解代码实现加分二叉树 题目 题解 简单讲解前序//中序//后序遍历 其实说白了&#xff0c;这个*序就是根root的遍历顺序 先序就是root–>left–>right 中序就是left–>root–>rig…

P1896 [SCOI2005]互不侵犯

P1896 [SCOI2005]互不侵犯 题目&#xff1a; 在NN的棋盘里面放K个国王&#xff0c;使他们互不攻击&#xff0c;共有多少种摆放方案。国王能攻击到它上下左右&#xff0c;以及左上左下右上右下八个方向上附近的各一个格子&#xff0c;共8个格子。 题解&#xff1a; 需要知道…

AT2368-[AGC013B]Hamiltonish Path【构造】

正题 题目链接:https://www.luogu.com.cn/problem/AT2368 题目大意 给出 nnn 个点 mmm 条边的一张无向图&#xff0c;然后求一条路径满足 路径长度不小于二。路径无交。对于所有的 xxx 与路径的端点相连&#xff0c;那么 xxx 在路径上。 1≤n,m≤1051\leq n,m\leq 10^51≤n,…

YBTOJ:运动积分(trie树)

文章目录题目描述解析代码题目描述 解析 做了巨长时间… 进行了一次刺激的阅读理解竞赛… 感谢whh dalao&#xff01; 那么让我们分析一下这道题 首先我们考虑单个求x选手的q值 不难发现 i 在第j天和x的大小关系只与a[x]与a[i]二进制下不同的最高位k有关 j的第k位与x相同时&am…

年终总结,我为什么离开舒适区?

当圣诞夜过去&#xff0c;也意味着这一年即将结束&#xff0c;迎来崭新的开始&#xff0c;一年时间既短暂&#xff0c;又漫长&#xff0c;当离开人生的舒适区&#xff0c;将迎来一个又一个的挑战。如果说一个人的优秀&#xff0c;取决于他天赋&#xff0c;以及他为之付出的额外…