BZOJ 3456: 城市规划 [多项式求逆元 DP]
题意:
求出n个点的简单(无重边无自环)无向连通图数目.
方案数mod 1004535809(479 * 2 ^ 21 + 1)即可.
n<=130000
DP求方案
g(n) n个点所有图的方案数 显然2C(n,2)=2n(n-1)
f(n) n个点连通图的方案数
然后枚举第一个点所在连通块的点数
g(n)=∑i=1..n-1{C(n-1,i-1)*f(i)*g(n-i)}
代入g(n) 两边同除(n-1)!消掉那个组合数上面那块,就变成了卷积的形式
我不写了直接看Miskcoo的公式啦 http://blog.miskcoo.com/2015/05/bzoj-3456
然后C(x)=A(x)*B(x)
A(x)=C(x)*B(x)-1
放在mod (x>n) 意义下求逆元就行了 因为需要的是a[n]
多项式求逆元
去看Miskcoo的教程吧 http://blog.miskcoo.com/2015/05/polynomial-inverse
简单的思路就是知道A(x) mod (x[n/2]) 下的逆元求mod (xn) 下的逆元
方法就是两个同余的式子写出来一减,两边平方再同乘A(x) 再移项
说一点关于意义的理解吧:
A(x)=Q(x)B(x)+R(x) degR<degB
A(x)Ξ0 (mod xn) 就是说A(x)的0..n-1项系数都是0
A(x)B(x)Ξ1 (mod xn) 它们每一项都有xn,否则不可能余数只有1;所以也有xn/2;
注意:
1.最后要乘(n-1)! 不要乘(n-1)
2.多项式求逆元每次长度都不确定,不能先预处理二进制反转
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=3e5+;
inline int read(){
char c=getchar();int x=,f=;
while(c<''||c>''){if(c=='-')f=-;c=getchar();}
while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
return x*f;
}
int P=,MOD=P;
ll Pow(ll a,ll b,ll MOD){
ll ans=;
for(;b;b>>=,a=a*a%MOD)
if(b&) ans=ans*a%MOD;
return ans;
}
struct NTT{
int n,rev[N];
ll g;
void ini(int m){
n=;
while(n<m) n<<=;
/*
int k=0;
while((1<<k)<n) k++;
for(int i=0;i<n;i++){
int t=0;
for(int j=0;j<k;j++) if(i&(1<<j)) t|=(1<<(k-j-1));
rev[i]=t;
}
*/
g=;
}
void transform(int *a,int flag,int n){
int k=;
while((<<k)<n) k++;
for(int i=;i<n;i++){
int t=;
for(int j=;j<k;j++) if(i&(<<j)) t|=(<<(k-j-));
if(t<i) swap(a[i],a[t]);
} for(int l=;l<=n;l<<=){
int m=l>>;
ll wn=Pow(g,flag==?(P-)/l:P--(P-)/l,P);
for(int *p=a;p!=a+n;p+=l){
ll w=;
for(int k=;k<m;k++){
ll t=w*p[k+m]%P;
p[k+m]=(p[k]-t+P)%P;
p[k]=(p[k]+t)%P;
w=w*wn%P;
}
}
}
if(flag==-){
ll inv=Pow(n,P-,P);
for(int i=;i<n;i++) a[i]=a[i]*inv%P;
}
}
int c[N];
void test(int *a,int n){for(int i=;i<n;i++) printf("%d ",a[i]);puts("");}
void polyInv(int deg,int *a,int *b){
if(deg==) b[]=Pow(a[],P-,P);
else{
polyInv((deg+)>>,a,b);
int n=;
while(n< deg<<) n<<=;
copy(a,a+deg,c);
fill(c+deg,c+n,);
transform(c,,n);
transform(b,,n);
for(int i=;i<n;i++)
b[i]=(ll)b[i]*(-(ll)b[i]*c[i]%P+P)%P;
transform(b,-,n);
fill(b+deg,b+n,);
}
}
}fft;
int n,inv[N],invFac[N],poc[N],A[N],B[N],C[N];
void getInv(int n){
inv[]=invFac[]=;
for(int i=;i<=n;i++){
if(i!=) inv[i]=-(ll)P/i*inv[P%i]%P;
if(inv[i]<) inv[i]+=P;
invFac[i]=(ll)invFac[i-]*inv[i]%P;
}
}
int main(){
//freopen("in","r",stdin);
n=read();
fft.ini(n);
getInv(n); poc[]=poc[]=;
for(int i=;i<=n;i++) poc[i]=Pow(,(ll)i*(i-)>>%(P-),P);
for(int i=;i<=n;i++) B[i]=(ll)poc[i]*invFac[i]%P;
for(int i=;i<=n;i++) C[i]=(ll)poc[i]*invFac[i-]%P;//printf("CC %d\n",C[i]);
fft.polyInv(fft.n,B,A); fft.n<<=;
fft.transform(A,,fft.n);
fft.transform(C,,fft.n);
for(int i=;i<fft.n;i++) A[i]=(ll)A[i]*C[i]%P;//,printf("ABC %d %d %d\n",i,A[i],C[i]); fft.transform(A,-,fft.n);
printf("%lld",(ll)A[n]*Pow(invFac[n-],P-,P)%P);
}
BZOJ 3456: 城市规划 [多项式求逆元 DP]的更多相关文章
- BZOJ 3456: 城市规划 [多项式求逆元 组合数学 | 生成函数 多项式求ln]
3456: 城市规划 题意:n个点组成的无向连通图个数 以前做过,今天复习一下 令\(f[n]\)为n个点的无向连通图个数 n个点的完全图个数为\(2^{\binom{n}{2}}\) 和Bell数的 ...
- BZOJ 3456: 城市规划 多项式求逆
Description 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了. 刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接 ...
- bzoj 3456 城市规划 多项式求逆+分治FFT
城市规划 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 1091 Solved: 629[Submit][Status][Discuss] Desc ...
- BZOJ 3456 城市规划 (组合计数、DP、FFT)
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=3456 著名的多项式练习题,做法也很多,终于切掉了纪念 首先求一波递推式: 令\(F(n ...
- bzoj 3456: 城市规划【NTT+多项式求逆】
参考:http://blog.miskcoo.com/2015/05/bzoj-3456 首先推出递推式(上面的blog讲的挺清楚的),大概过程是正难则反,设g为n个点的简单(无重边无自环)无向图数目 ...
- [BZOJ 3456]城市规划(cdq分治+FFT)
[BZOJ 3456]城市规划(cdq分治+FFT) 题面 求有标号n个点无向连通图数目. 分析 设\(f(i)\)表示\(i\)个点组成的无向连通图数量,\(g(i)\)表示\(i\)个点的图的数量 ...
- 51NOD 1258 序列求和 V4 [任意模数fft 多项式求逆元 伯努利数]
1258 序列求和 V4 题意:求\(S_m(n) = \sum_{i=1}^n i^m \mod 10^9+7\),多组数据,\(T \le 500, n \le 10^{18}, k \le 50 ...
- 多项式求逆元详解+模板 【洛谷P4238】多项式求逆
概述 多项式求逆元是一个非常重要的知识点,许多多项式操作都需要用到该算法,包括多项式取模,除法,开跟,求ln,求exp,快速幂.用快速傅里叶变换和倍增法可以在$O(n log n)$的时间复杂度下求出 ...
- BZOJ 3456: 城市规划 与 多项式求逆算法介绍(多项式求逆, dp)
题面 求有 \(n\) 个点的无向有标号连通图个数 . \((1 \le n \le 1.3 * 10^5)\) 题解 首先考虑 dp ... 直接算可行的方案数 , 容易算重复 . 我们用总方案数减 ...
随机推荐
- 去除FineReport报表点击单元格时出现的黑框
选择模版->模版web属性 添加加载结束事件,具体代码如下 $("td").bind("click",function(){return false;}) ...
- Coins、Tokens、山寨币:区别在哪里
[译] Coins, Tokens & Altcoins: What's the Difference? coin:指"正宗的"数字货币,比如比特币 tokens:比如用 ...
- Java数据结构和算法(十二)——2-3-4树
通过前面的介绍,我们知道在二叉树中,每个节点只有一个数据项,最多有两个子节点.如果允许每个节点可以有更多的数据项和更多的子节点,就是多叉树.本篇博客我们将介绍的——2-3-4树,它是一种多叉树,它的每 ...
- OpenStack Horizon创建虚拟机时增加虚拟机OS用户
背景 通过OpenStack的Horizon使用镜像创建虚拟机(以Ubuntu为例),如果不知道镜像的用户名和密码,在创建好虚拟机之后,无法登录虚拟机的OS.因此,我们需要一种方法,创建虚拟机时增加用 ...
- html日历(3)
<!DOCTYPE html> <html> <head> <meta http-equiv="Content-Type" content ...
- Navicat连接阿里云Mysql遇到的的坑
连上去那一刻,心态真是起飞了
- 数据库01创建表和DML语言
楼主用的数据库时mysql,用的时navacat for mysql. 数据库层面: 1.显示所有数据库 show databases; 2.创建数据库,名字叫lyh,编码为utf-8 create ...
- SDP(2):ScalikeJDBC-Connection Pool Configuration
scalikeJDBC可以通过配置文件来设置连接池及全局系统参数.对配置文件的解析是通过TypesafeConfig工具库实现的.默认加载classpath下的application.conf,app ...
- 通过Xshell如何从Linux服务器下载文件(亲测可行)
到网上下载lrzsz安装包,这里以lrzsz-0.12.20.tar.gz为例 2 打开终端 cd 到安装包所在目录 tar zxvf lrzsz-0.12.20.tar.gz 解压安装包 3 进入解 ...
- addslashes() 函数返回在预定义字符之前添加反斜杠的字符串
. 预定义字符是: 单引号(') 双引号(") 反斜杠(\) NULL 提示:该函数可用于为存储在数据库中的字符串以及数据库查询语句准备字符串. 注释:默认地,PHP 对所有的 GET.PO ...