P4233-射命丸文的笔记【NTT,多项式求逆】
时间:2021-08-08
本文章向大家介绍P4233-射命丸文的笔记【NTT,多项式求逆】,主要包括P4233-射命丸文的笔记【NTT,多项式求逆】使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
正题
题目链接:https://www.luogu.com.cn/problem/P4233
题目大意
随机选择一条有哈密顿回路的\(n\)个点的竞赛图,求选出图的哈密顿回路的期望个数。
对于每个\(n\in[1,N]\)求答案。
\(1\leq N\leq 10^5\)
解题思路
竟然自己推出来了泪目( Ĭ ^ Ĭ )
如果是统计所以的哈密顿回路个数是一个很简单的题目,我们可以求出\(n\)的一个圆排列表示一条回路,然后剩下的边随便排即可。也就是\((n-1)!\times 2^{\frac{n(n-1)}{2}-n}\)条哈密顿路,但是因为求的是期望所以我们还得求出有哈密顿回路的竞赛图个数。
这个是问题所在,我们可以考虑用城市规划的推法,设\(f_i\)表示\(i\)个点有哈密顿回路的竞赛图个数。
那么有
\[2^{\frac{n(n-1)}2}=2\sum_{i=0}^{n-1}2^{\frac{i(i-1)}{2}}f_{n-i}\binom{n}{i}
\]
但是注意\(n=0\)的时候要特别处理算出来为\(1\)。
化一下式子有
\[2^{\frac{n(n-1)}2}=2\sum_{i=0}^{n-1}2^{\frac{i(i-1)}{2}}f_{n-i}\frac{n!}{i!(n-i)!}
\]
\[\frac{2^{\frac{n(n-1)}2}}{n!}=\sum_{i=0}^{n-1}\frac{2^{\frac{i(i-1)}{2}}}{i!}\frac{2f_{n-i}}{(n-i)!}
\]
设\(F=\sum_{i=0}^{\infty}\frac{2f_i}{i!},G=\sum_{i=0}^{\infty}\frac{2^{\frac{i(i-1)}{2}}}{i!}\),那么有
\[G=FG+1\Rightarrow F=\frac{G-1}{G}
\]
上多项式求逆就可以求出\(f\)了。
时间复杂度\(O(n\log n)\)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=131072,M=N<<1,P=998244353;
ll n,fac[M],G[M],H[M],r[M],tmp[M];
ll power(ll x,ll b){
ll ans=1;
while(b){
if(b&1)ans=ans*x%P;
x=x*x%P;b>>=1;
}
return ans;
}
void NTT(ll *f,ll n,ll op){
for(ll i=0;i<n;i++)
if(i<r[i])swap(f[i],f[r[i]]);
for(ll p=2;p<=n;p<<=1){
ll len=(p>>1),tmp=power(3,(P-1)/p);
if(op==-1)tmp=power(tmp,P-2);
for(ll k=0;k<n;k+=p){
ll buf=1;
for(ll i=k;i<k+len;i++){
ll tt=buf*f[i+len]%P;
f[i+len]=(f[i]-tt+P)%P;
f[i]=(f[i]+tt)%P;
buf=buf*tmp%P;
}
}
}
if(op==-1){
ll invn=power(n,P-2);
for(ll i=0;i<n;i++)
f[i]=f[i]*invn%P;
}
return;
}
void GetInv(ll n,ll *f,ll *g){
if(!n)
{g[0]=power(f[0],P-2);return;}
GetInv(n>>1,f,g);ll m=n<<1;
for(ll i=0;i<n;i++)tmp[i]=f[i];
for(ll i=0;i<m;i++)r[i]=(r[i>>1]>>1)|((i&1)?(m>>1):0);
NTT(tmp,m,1);NTT(g,m,1);
for(ll i=0;i<m;i++)
g[i]=(2*g[i]-tmp[i]*g[i]%P*g[i]%P+P)%P;
NTT(g,m,-1);
for(ll i=n;i<m;i++)g[i]=0;
return;
}
signed main()
{
scanf("%lld",&n);fac[0]=1;
for(ll i=1;i<N;i++)fac[i]=fac[i-1]*i%P;
for(ll i=0;i<N;i++)G[i]=power(2,i*(i-1)/2ll)*power(fac[i],P-2)%P;
GetInv(N,G,H);G[0]--;
NTT(G,M,1);NTT(H,M,1);
for(ll i=0;i<M;i++)G[i]=G[i]*H[i]%P;
NTT(G,M,-1);
for(ll i=1;i<=n;i++){
if(i==1){puts("1");continue;}
G[i]=G[i]*fac[i]%P;
if(!G[i]){puts("-1");continue;}
ll ans=fac[i-1]*power(2,i*(i-1)/2ll-i)%P;
printf("%d\n",ans*power(G[i],P-2)%P);
}
return 0;
}
原文地址:https://www.cnblogs.com/QuantAsk/p/15116484.html
- Playrix Codescapes Cup (Codeforces Round #413, rated, Div. 1 + Div. 2)(A.暴力,B.优先队列,C.dp乱搞)
- [libvirt][nginx]libvirt文档访问速度提高的小技巧
- 图论----同构图(详解)
- 基数排序与桶排序,计数排序【详解】
- SG函数和SG定理【详解】
- 密码学经典之生日悖论与生日攻击【详解】
- POJ 1659 Frogs' Neighborhood(可图性判定—Havel-Hakimi定理)【超详解】
- BZOJ 1192: [HNOI2006]鬼谷子的钱袋(新生必做的水题)
- 快速傅里叶变换(FFT)算法【详解】
- Codeforces Round #416 (Div. 2)(A,思维题,暴力,B,思维题,暴力)
- 作为程序员的你在外行人眼里是一个怎样的群体?
- 高斯消元法(Gauss Elimination)【超详解&模板】
- [快学Python3]读写Excel - openpyxl库
- HDU 2147 kiki's game(规律,博弈)
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法