【牛客挑战赛31D】雷的打字机

题目

首先看到这个出现长度至少为\(2\)的回文子串

这就等价于不能出现两个连续且相同的字符

于是我们用概率生成函数来搞

设\(g_i\)表示\(i\)次操作后游戏没有结束的概率，\(f_{i,j}\)表示\(i\)次操作之后出现最后两个字符都是\(j\)的概率，这样的话游戏就结束了

再定义\(f_i=\sum_{j=1}^nf_{i,j}\)

我们要求的东西显然就是\(F'(1)\)，即期望长度

那个非常有用的\(g_i=g_{i+1}+f_{i+1}\)还是成立的

于是

\[xG(x)+1=F(x)+G(x)\]

代入\(x=1\)，并两边求导

\[G(1)=F'(1)\]

显然对于任意的\(F_i(x)\)都会存在

\[G(x)p_i^2x^2=F_i(x)p_ix+F_i(x)\]

就是往一个非结束状态后面强塞两个字符\(i\)，代入\(x=1\)就有

\[G(1)p_i^2=F_i(1)(1+p_i)\]

\[F_i(1)=\frac{p_i^2}{1+p_i}G(1)\]

又因为\(\sum_{i=1}^nF_i(1)=1\)，于是\(\sum_{i=1}^n\frac{p_i^2}{1+p_i}G(1)=1\)，于是答案就是

\[\frac{1}{\sum_{i=1}^n\frac{p_i^2}{1+p_i}}\]

之后发现\(n\leq 10^7\)，带\(log\)求逆元显然不是很可行了

于是离线求一下逆元就可以了，就像十二省联考里的D2T3那样就好了

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define re register
const int mod=998244353;
const int maxn=1e7+5;
int n,ans;unsigned int seed;
int fac[maxn],p[maxn];
inline int Rand() {seed=(seed<<2)^(seed>>5)^(seed<<7)^(seed>>11);return seed%mod;
}
inline int ksm(int a,int b) {int S=1;for(;b;b>>=1,a=1ll*a*a%mod) if(b&1) S=1ll*S*a%mod;return S;
}
inline void Generate() { scanf("%d%u",&n,&seed);int sum=0; for(re int i=1;i<n;++i) p[i]=Rand(),sum=(sum+p[i])%mod;  p[n]=(mod+1-sum)%mod;
}
int main() {Generate();fac[0]=1;for(re int i=1;i<=n;i++) fac[i]=1ll*fac[i-1]*(p[i]+1)%mod;int Inv=ksm(fac[n],mod-2),inv;for(re int i=n;i;i--) {inv=1ll*Inv*fac[i-1]%mod;ans=(ans+1ll*p[i]*p[i]%mod*inv%mod)%mod;Inv=1ll*Inv*(p[i]+1)%mod;}printf("%d\n",ksm(ans,mod-2));return 0;
}

转载于:https://www.cnblogs.com/asuldb/p/10950678.html

【牛客挑战赛31D】雷的打字机相关推荐

牛客挑战赛47 C 条件（Floyd bitset优化）
牛客挑战赛47 C 条件思路:首先我们要两个图,一个是一定能到达的,一个是可能到达的,如果我们使用floyd (n^3)就有可能会超时,因为只要求询问能否到达,所以权值只有0和1,那我们可以使用bi ...
牛客挑战赛47 A 一道GCD问题
牛客挑战赛47 A 一道GCD问题思路参考牛客上的题解: 根据多维的更相减损术得gcd(x,y,z)=gcd(x,y−x,z−y)得 gcd(a1+k,a2+k,a3+k-,an+k)=gcd(a1 ...
牛客挑战赛47 D Lots of Edges（最短路+递归枚举子集）
牛客挑战赛47 D Lots of Edges 思路:点的权值最多只有(1<<17)-1(131071) ,那我们可以枚举终点的值来算最短路,每个点能连边的值都是固定的,可以通过递归枚举子 ...
牛客挑战赛42 A.小睿睿的数列
牛客挑战赛42 A.小睿睿的数列题目链接题目描述小睿睿给了你一个长度为n的数列,他想问你该数列中满足条件(区间内存在某个数是区间内所有数的公因数)的最长区间有多少个输入描述: 第一行 111 ...
牛客挑战赛36 - 纸飞机
题目链接:牛客挑战赛36 - 纸飞机题目描述直线上有n座山峰,第i座的高度为hi.从某座山峰上放飞一架纸飞机,它可以从左往右依次经过一系列高度严格递减的山头. 假设五座山峰的高度依次是3,4,3, ...
牛客挑战赛34 A 能天使的愿望（dp 分组背包）
链接:https://ac.nowcoder.com/acm/contest/2271/A 来源:牛客网题目描述出题人寄给大家的一些闲话:参加了CSP-J/S 2019 的同学,考的都怎么样啊?不 ...
牛客挑战赛34 A.能天使的愿望（分组背包）
链接:https://ac.nowcoder.com/acm/contest/2271/A 来源:牛客网时间限制:C/C++ 2秒,其他语言4秒空间限制:C/C++ 65536K,其他语言1310 ...
牛客挑战赛48--A and B
A-天使果冻链接:https://ac.nowcoder.com/acm/contest/11161/A 来源:牛客网题意:一个数组包含n个数,求前i个数的次大值是多少(倘若最大值有两个,则次大值 ...
【牛客挑战赛46】F.柠檬树
Description 这是一棵有 n 个柠檬的柠檬树,由 n-1 条枝条连接而成. 秋天了,柠檬都成熟了, 牛牛和牛妹准备选一些柠檬送给他们的朋友们. 对于每一个朋友,牛妹会选择第 l-r 个柠檬送 ...

【牛客挑战赛31D】雷的打字机

【牛客挑战赛31D】雷的打字机相关推荐

最新文章

热门文章

【牛客挑战赛31D】 雷的打字机

【牛客挑战赛31D】 雷的打字机相关推荐

最新文章

热门文章

【牛客挑战赛31D】雷的打字机

【牛客挑战赛31D】雷的打字机相关推荐