博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ 3456 城市规划
阅读量:6977 次
发布时间:2019-06-27

本文共 1895 字,大约阅读时间需要 6 分钟。

题目链接

题解

nn个点构成的连通图数量为

g(n)=2(n2)g(n)=2(n2)
枚举
11号点所在连通块的大小,得到
gn=i=1n(n1i1)fignign=∑i=1n(n−1i−1)fign−i
因此
2(n2)=i=1n(n1)!(i1)!(ni)!figni2(n2)=∑i=1n(n−1)!(i−1)!(n−i)!fign−i
整理得到
2(n2)(n1)!=i=1nfi(i1)!2(ni2)(ni)!2(n2)(n−1)!=∑i=1nfi(i−1)!2(n−i2)(n−i)!
F=fi(i1)!xi,G=2(i2)i!xi,T=2(i2)(i1)!xiF=fi(i−1)!xi,G=2(i2)i!xi,T=2(i2)(i−1)!xi
则有
T=FG(modxn+1)T=FG(modxn+1)
F=TG1(modxn+1)F=TG−1(modxn+1)
先求出
GG
TT,求
GG的逆元,卷积一遍就可以得到
FF,得到
FF之后很容易得到
fifi

代码

#include 
#include
#include
int read(){ int x=0,f=1; char ch=getchar(); while((ch<'0')||(ch>'9')) { if(ch=='-') { f=-f; } ch=getchar(); } while((ch>='0')&&(ch<='9')) { x=x*10+ch-'0'; ch=getchar(); } return x*f;}const int maxn=300000;const int mod=1004535809;const int G=3;int quickpow(int a,int b,int m){ int res=1; while(b) { if(b&1) { res=1ll*res*a%m; } a=1ll*a*a%m; b>>=1; } return res;}int plus(int a,int b,int m){ int res=a+b; if(res>=m) { res-=m; } return res;}int minus(int a,int b,int m){ int res=a-b; if(res<0) { res+=m; } return res;}int rev[maxn+10],tmp[maxn+10];int getrev(int n){ int m=1,len=0; while(m<=n) { m<<=1; ++len; } for(int i=0; i
>1]>>1)+((i&1)<<(len-1)); } return m;}int fft(int *s,int len){ for(int i=0; i
>1); ++k) { int x=s[j+k],y=1ll*g*s[j+k+(i>>1)]%mod; s[j+k]=plus(x,y,mod); s[j+k+(i>>1)]=minus(x,y,mod); g=1ll*g*gn%mod; } } } return 0;}int getinv(int *s,int *v,int deg){ if(deg==1) { v[0]=quickpow(s[0],mod-2,mod); return 0; } getinv(s,v,(deg+1)>>1); int p=getrev(deg<<1); for(int i=0; i

转载于:https://www.cnblogs.com/Canopus-wym/p/10376151.html

你可能感兴趣的文章
浅谈SQL Server中统计对于查询的影响
查看>>
WF4 Beta,RC版文章总结
查看>>
WPF 与Surface 2.0 SDK 亲密接触–LibraryContainer 篇
查看>>
C# 对应 Oracle 存储过程 的 SYS_REFCURSOR 应该 传入什么类型的参数?
查看>>
Unity3D移植到自己的Android程序
查看>>
【转】用示例说明索引数据块中出现热块的场景,并给出解决方案
查看>>
HDU 2034 人见人爱A-B
查看>>
【AngularJS】—— 12 独立作用域
查看>>
使用工作集(Working Set)整理项目
查看>>
MailMail、RegeX等程序的云端版
查看>>
[Erlang 0072] Erlang XML处理解决方案
查看>>
从C#到Objective-C,循序渐进学习苹果开发(7)--使用FMDB对Sqlite数据库进行操作
查看>>
mmap学习
查看>>
X3D中Profile如何翻译
查看>>
7.14. revision
查看>>
第 175 章 Open Source Requirements Management Tool
查看>>
CentOS7安装配置redis-3.0.0
查看>>
SQL server 专业词汇
查看>>
Selenium2+python自动化25-js处理日历控件(修改readonly属性)
查看>>
ArcGIS制图之Sub Points点抽稀
查看>>