【CF923E】Perpetual Subtraction【生成函数】【积分推式子】【NTT卷积】

题意:有一个整数x∈[0,n]x\in[0,n]x[0,n],取iii的概率为pip_ipi。执行mmm次操作,每次把xxx等概率变成[0,x][0,x][0,x]中的一个整数,求操作完后等于每个数的概率。模998244353998244353998244353

n≤105,m≤1018n\leq 10^5,m\leq10^{18}n105,m1018

显然有一个dp式子

fi,j=∑k=jnfi−1,kk+1f_{i,j}=\sum_{k=j}^n\frac{f_{i-1,k}}{k+1}fi,j=k=jnk+1fi1,k

考虑生成函数

fi(x)=∑j=0n∑k=jn[xk]fi−1(x)k+1xjf_i(x)=\sum_{j=0}^n\sum_{k=j}^n\frac{[x^k]f_{i-1}(x)}{k+1}x^jfi(x)=j=0nk=jnk+1[xk]fi1(x)xj

=∑k=0n[xk]fi−1(x)k+1∑j=0kxj=\sum_{k=0}^n\frac{[x^k]f_{i-1}(x)}{k+1}\sum_{j=0}^kx^j=k=0nk+1[xk]fi1(x)j=0kxj

等比数列求和

=∑k=0n[xk]fi−1(x)k+1xk+1−1x−1=\sum_{k=0}^n\frac{[x^k]f_{i-1}(x)}{k+1}\frac{x^{k+1}-1}{x-1}=k=0nk+1[xk]fi1(x)x1xk+11

似乎推不动了

但是有奥妙重重的一步:

=∑k=0n[xk]fi−1(x)x−1xk+1−1k+1=\sum_{k=0}^n\frac{[x^k]f_{i-1}(x)}{x-1}\frac{x^{k+1}-1}{k+1}=k=0nx1[xk]fi1(x)k+1xk+11

看起来没啥用,但这样把右边写成了积分的形式

=∑k=0n[xk]fi−1(x)x−1∫1xtkdt=\sum_{k=0}^n\frac{[x^k]f_{i-1}(x)}{x-1}\int_1^xt^kdt=k=0nx1[xk]fi1(x)1xtkdt

然后盯着这个式子:

∑k=0n[xk]fi−1(x)∫1xtkdt\sum_{k=0}^n[x^k]f_{i-1}(x)\int_1^xt^kdtk=0n[xk]fi1(x)1xtkdt

ttt可能看不习惯,换个更通用的写法:

∑k=0n[xk]f(x)∫LRxkdx\sum_{k=0}^n[x^k]f(x) \int_L^Rx^kdxk=0n[xk]f(x)LRxkdx

相当于把多项式拆成单项式,每个单项式对区间[L,R][L,R][L,R]求定积分,最后还贴心地帮你乘上了系数

那不就是对多项式求积分吗

∫LRf(x)dx\int_L^Rf(x)dxLRf(x)dx

代回原式:

fi(x)=1x−1∫1xfi−1(t)dtf_i(x)=\frac{1}{x-1}\int_1^xf_{i-1}(t)dtfi(x)=x111xfi1(t)dt

然而这个积分下界是111,不好处理

然后又是奥妙重重的一步:

gi(x)=fi(x+1)g_i(x)=f_i(x+1)gi(x)=fi(x+1)

gi(x)=1x∫1x+1fi−1(t)dtg_i(x)=\frac{1}{x}\int_1^{x+1}f_{i-1}(t)dtgi(x)=x11x+1fi1(t)dt

=1x∫0xfi−1(t+1)dt=\frac{1}{x}\int_0^{x}f_{i-1}(t+1)dt=x10xfi1(t+1)dt

=1x∫0xgi−1(t)dt=\frac{1}{x}\int_0^{x}g_{i-1}(t)dt=x10xgi1(t)dt

那就gi−1g_{i-1}gi1的不定积分在xxx处的取值,容易得出

gi,j=1j+1gi−1,jg_{i,j}=\frac{1}{j+1}g_{i-1,j}gi,j=j+11gi1,j

这样如果知道g0g_0g0就可以轻松算出gmg_mgm

fffggg直接二项式定理暴拆再卷积一下就可以互相转换,问题解决

复杂度O(nlog⁡n)O(n\log n)O(nlogn)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#define MAXN 262144
using namespace std;
const int MOD=998244353;
typedef long long ll;
inline int read()
{int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans;
}
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD;p>>=1;}return ans;
}
#define inv(x) qpow(x,MOD-2)
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
int rt[2][24];
int r[MAXN],l,lim;
inline void init(){lim=1<<l;for (int i=0;i<lim;i++) r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));}
inline void NTT(int* a,int type)
{for (int i=0;i<lim;i++) if (i<r[i]) swap(a[i],a[r[i]]);for (int L=0;L<l;L++){int mid=1<<L,len=mid<<1;ll Wn=rt[type][L+1];for (int s=0;s<lim;s+=len)for (int k=0,w=1;k<mid;k++,w=w*Wn%MOD){int x=a[s+k],y=(ll)w*a[s+mid+k]%MOD;a[s+k]=add(x,y),a[s+mid+k]=dec(x,y);}}if (type){int t=inv(lim);for (int i=0;i<lim;i++) a[i]=(ll)a[i]*t%MOD;}
}
int fac[MAXN],finv[MAXN];
int a[MAXN];
int f[MAXN],g[MAXN];
int main()
{rt[0][23]=qpow(3,119);rt[1][23]=inv(rt[0][23]);for (int i=22;i>=0;i--){rt[0][i]=(ll)rt[0][i+1]*rt[0][i+1]%MOD;rt[1][i]=(ll)rt[1][i+1]*rt[1][i+1]%MOD;}int n=read();ll m;scanf("%lld",&m);m%=MOD-1;for (int i=0;i<=n;i++) a[i]=read();fac[0]=1;for (int i=1;i<=n;i++) fac[i]=(ll)fac[i-1]*i%MOD;finv[n]=inv(fac[n]);for (int i=n-1;i>=0;i--) finv[i]=(ll)finv[i+1]*(i+1)%MOD;for (int i=0;i<=n;i++) f[i]=(ll)a[i]*fac[i]%MOD;for (int i=0;i<=n;i++) g[i]=finv[n-i];for (l=0;(1<<l)<=(n<<1);++l);init();NTT(f,0);NTT(g,0);for (int i=0;i<lim;i++) f[i]=(ll)f[i]*g[i]%MOD;NTT(f,1);for (int i=0;i<=n;i++) a[i]=(ll)f[i+n]*finv[i]%MOD;for (int i=0;i<=n;i++) a[i]=(ll)a[i]*inv(qpow(i+1,m))%MOD;for (int i=0;i<=n;i++) f[i]=(ll)a[i]*fac[i]%MOD;for (int i=0;i<=n;i++) g[i]=(((n-i)&1)? MOD-finv[n-i]:finv[n-i]);for (int i=n+1;i<lim;i++) f[i]=g[i]=0;NTT(f,0);NTT(g,0);for (int i=0;i<lim;i++) f[i]=(ll)f[i]*g[i]%MOD;NTT(f,1);for (int i=0;i<=n;i++) a[i]=(ll)f[i+n]*finv[i]%MOD;for (int i=0;i<=n;i++) printf("%d ",a[i]);return 0;
}

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

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

相关文章

使用Jenkins部署.Net Core遇到的几个坑

一、相关环境和版本搞过CI/CD的同学一定吃过不少苦头&#xff0c;或者说遇到不少坑&#xff0c;但是对自动化的执着住挡不了前进的步伐&#xff0c;如果你缺少了运维这一块知识&#xff0c;那么你的流水线总是不那么完美&#xff0c;本文记录的是自己躺过的坑&#xff0c;希望对…

Deltix Round, Spring 2021 E. Crypto Lights 组合数学 + 推公式

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个灯&#xff0c;每次可以打开一个灯&#xff0c;当连续的kkk个灯有至少两个灯开着的时候停止&#xff0c;问最终期望能打开多少灯。 思路&#xff1a; 由于不想打latexlatexlatex&#xff0c;所以…

【洛谷P5385】须臾幻境/【BZOJ3514】Codechef MARCH14 GERALD07加强版【LCT】【主席树】

题意&#xff1a;有nnn个点mmm条边&#xff0c;qqq次询问连接区间[L,R][L,R][L,R]中的边后的连通块个数。强制在线。 n,m,q≤2105n,m,q\leq 2\times10^5n,m,q≤2105 显然连通块个数n−任意一个生成森林的边数连通块个数n-任意一个生成森林的边数连通块个数n−任意一个生成森林…

.NET Core 仿魔兽世界密保卡实现

《魔兽世界》的老玩家都知道&#xff0c;密保卡曾经被用于登录验证&#xff0c;以保证账号安全。今天我用.NET Core模拟了一把密保卡&#xff08;也叫矩阵卡&#xff09;的实现&#xff0c;分享给大家。密保卡的原理这是一张典型的魔兽世界密保卡。序列号用于绑定游戏账号&…

Deltix Round, Spring 2021 D. Love-Hate 随机化 + sos dp(高维前缀和)

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个朋友&#xff0c;一共有mmm种货币&#xff0c;一个朋友最多喜欢ppp种&#xff0c;用二进制给出111代表喜欢&#xff0c;让你选出最多的一个货币集合使得至少有⌈n2⌉\left \lceil \frac{n}{2} \rig…

【UOJ207】共价大爷游长沙【LCT】【异或】【随机化】

传送门 题意&#xff1a;维护一棵无权树和一个路径集合SSS&#xff0c;支持以下操作&#xff1a; 断边连边在SSS加入中加入一条路径删除SSS中的一条路径询问是否SSS中的所有路径都经过了边(x,y)(x,y)(x,y) n≤105,q≤3105n\leq10^5,q\leq3\times10^5n≤105,q≤3105 给每条加…

浅谈C#泛型

一.为什么要提出泛型的概念我们在声明对象或者方法中&#xff0c;对象中成员变量的定义或者函数参数都传递都要指定具体的对象类型&#xff0c;但是有的时候参数的类型是变化的&#xff0c;但是实现的功能却又差不多&#xff0c;这个时候我们就想&#xff0c;是否存在一种东西可…

Codeforces Round #635 (Div. 1) C. Kaavi and Magic Spell 区间dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个串s,ts,ts,t&#xff0c;每次都可以从sss的开头拿一个字符放到AAA串的开头或结尾&#xff0c;问最终有多少种方案使得ttt是AAA的前缀&#xff0c;注意sss不必全部拿完。 m,n≤3000m,n\le3000m,n≤30…

【IOI2018】会议【笛卡尔树】【dp】【线段树】

题意&#xff1a;长度为nnn的序列&#xff0c;qqq次询问&#xff0c;每次给定一个区间&#xff0c;钦定区间中的一个位置xxx&#xff0c;使得区间所有点 与xxx之间的最大值&#xff08;含端点&#xff09; 之和 最小&#xff0c;输出最小值。 n,q≤7.5105n,q\leq7.5\times10^5n…

【半小时大话.net依赖注入】(一)理论基础+实战控制台程序实现AutoFac注入

第一章|理论基础实战控制台程序实现AutoFac注入第二章|AutoFac的常见使用套路第三章|实战Asp.Net Framework Web程序实现AutoFac注入第四章|实战Asp.Net Core自带DI实现依赖注入第五章|实战Asp.Net Core引入AutoFac的两种方式简介该系列共5篇文章&#xff0c;旨在以实战模式&am…

P4198 楼房重建 线段树 + 区间合并

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 题面有点问题&#xff0c;按照人类正常的理解来就好啦。 思路&#xff1a; 可以想到维护每个位置的一个斜率&#xff0c;模拟的话就是从第一个位置开始向后选&#xff0c;当某个位置斜率大于当前位置的时…

【IOI2018】狼人【Kruscal重构树】【主席树】

题意&#xff1a;nnn个点mmm条边的无向图&#xff0c;qqq次询问&#xff0c;每次给定s,t,L,Rs,t,L,Rs,t,L,R&#xff0c;判断是否存在一条sss到ttt的路径&#xff0c;使得路径上可以找到一点kkk,满足此路径s∼ks\sim ks∼k的部分标号都≥L\geq L≥L且k∼tk\sim tk∼t标号都≤R\…

基于 Docker Compose 实践 .NET Core 的现代化架构 2:在容器中集成 Skywalking APM

本章节代码已经上传至 https://github.com/siegrainwong/.NET-Core-with-Docker/tree/master/Part2系列大纲还是先介绍一下目录&#xff0c;这次我们讲第二篇&#xff1a;用 docker-compose 启动 WebApi 和 SQL Server在容器中集成 Skywalking APM通过 nginx-proxy 对 ES、Skyw…

【雅礼集训2017】字符串【后缀自动机】【数据分治】

题意&#xff1a;给定一个字符串SSS和mmm个区间[li,ri][l_i,r_i][li​,ri​],qqq次询问&#xff0c;每次给定长度为kkk的字符串www和区间[a,b][a,b][a,b]&#xff0c;求对于所有i∈[a,b]i\in[a,b]i∈[a,b]&#xff0c;www在[li,ri][l_i,r_i][li​,ri​]内的子串在SSS中出现次数…

AtCoder Regular Contest 100 E - Or Plus Max Sos dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为2n2^n2n的数组&#xff0c;让你对于所有的1≤k≤2n−11\le k\le 2^n-11≤k≤2n−1求最大的aiaj,0≤i<j≤2n−1,iorj≤ka_ia_j,0\le i<j\le2^n-1,i\ \ or \ \ j\le kai​aj​,0≤i<j≤…

使用Minikube部署本地Kubernetes集群(二十九)

前言 使用Minikube部署本地k8s集群相对比较简单&#xff0c;非常推荐将其用于本地k8s开发环境&#xff0c;唯一麻烦点的仅仅是网络问题。在本篇教程中&#xff0c;我们使用了国内的镜像来完成本地k8s集群的搭建。如果搭建过程中出现错误&#xff0c;建议使用mini…

Manthan, Codefest 19 (open for everyone, rated, Div. 1 + Div. 2) F. Bits And Pieces sosdp预处理超集

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤1e6,ai≤2e6n\le1e6,a_i\le2e6n≤1e6,ai​≤2e6 思路&#xff1a; 由于(aj&ak)(a_j \And a_k)(aj​&ak​)打的括号&#xff0c;所以应该放在一起考虑&#xff0c;现在我们可以枚举aia_iai​&…

【PKUWC2018】随机算法【状压dp】【组合计数】

题意&#xff1a;一张nnn个点mmm条边的无向无权图&#xff0c;求以下算法计算最大独立集的正确率&#xff1a;随机一个排列&#xff0c;依次考虑排列中每一个点&#xff0c;如果不与任何一个独立集中的点相邻则将其加入独立集。模998244353998244353998244353。 n≤20n\leq20n≤…

架构杂谈《一》

架构杂谈《一》从传统单体架构到服务化架构的发展历程典型的单体架构分为三个层级&#xff0c;Web层、业务逻辑层和数据存储层&#xff0c;每个层的指责分别如下&#xff1a;Web 层&#xff1a;负责与用户交互或者对外提供接口业务逻辑层&#xff1a;为了实现业务逻辑而设计的流…

hdu 3308 LCIS 线段树 + 区间合并

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 日常水一篇题解。 带修改的求区间连续的递增序列&#xff0c;我们考虑用线段树维护。 直接维护mlenmlenmlen是区间最长的递增序列&#xff0c;lslsls是从左端点开始的最长递增序列&#xf…