Description
求n个点的带标号连通简单图的个数。答案对1004535809取模。
n<=130000
Solution
根据套路,首先我们需要求出n个点的带标号简单图的个数,设为g[n]
显然我们有g[n]=2^(n*(n-1)/2)
考虑设答案为f[n],要怎么求出?
枚举1号点所在的联通块的大小,我们可以知道
g[n]=∑i=1nf[i]Ci−1n−1g[n−i]
因为f是保证联通的所以不会算重
把g带入,C拆开
2n∗(n+1)/2=∑i=1nf[i](n−1)!(i−1)!(n−i)!2(n−i)∗(n−i−1)/2
两边同除(n-1)!
2n∗(n+1)/2(n−1)!=∑i=1nf[i](i−1)!2(n−i)∗(n−i−1)/2(n−i)!
可以发现这是一个卷积的形式
那么我们就直接上FFT就好了
需要多项式求逆
第一次写写的很丑QwQ
Code
#include <cstdio>
#include <cstring>
#include <algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
typedef double db;
ll read() {
char ch;
for(ch=getchar();ch<'0'||ch>'9';ch=getchar());
int x=ch-'0';
for(ch=getchar();ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0';
return x;
}
void write(int x) {
if (!x) {puts("0");return;}
char ch[20];int tot=0;
for(;x;x/=10) ch[++tot]=x%10+'0';
fd(i,tot,1) putchar(ch[i]);
puts("");
}
const int Mo=1004535809,D=3,N=5*1e5+5;
int pwr(int x,ll y) {
int z=1;
for(;y;y/=2,x=(ll)x*x%Mo)
if (y&1) z=(ll)z*x%Mo;
return z;
}
int n,len,lg,inv,t[N],a[N],b[N],c[N],W[N];
int p[N],h[N],g[N],f[N];
void pre(int len) {
W[0]=1;W[1]=pwr(D,(Mo-1)/len);
for(int i=2;i<=len;i++) W[i]=(ll)W[i-1]*W[1]%Mo;
inv=pwr(len,Mo-2);
}
void DFT(int *a,int len,int lg,int flag) {
for(int i=0;i<len;i++) {
int p=0;
for(int j=i,k=0;k<lg;k++,j/=2) p=(p<<1)+(j&1);
t[p]=a[i];
}
for(int m=2;m<=len;m<<=1) {
int half=m>>1,times=len/m;
for(int i=0;i<half;i++) {
int w=(flag==1)?W[i*times]:W[len-i*times];
for(int j=i;j<len;j+=m) {
int u=t[j],v=(ll)t[j+half]*w%Mo;
t[j]=(u+v)%Mo;
t[j+half]=(u-v+Mo)%Mo;
}
}
}
for(int i=0;i<len;i++) a[i]=t[i];
if (flag==-1) for(int i=0;i<len;i++) a[i]=(ll)a[i]*inv%Mo;
}
void NTT(int *A,int *B,int *C,int len,int lg) {
for(int i=0;i<len;i++) a[i]=A[i],b[i]=B[i];
DFT(a,len,lg,1);DFT(b,len,lg,1);
for(int i=0;i<len;i++) c[i]=(ll)a[i]*b[i]%Mo;
DFT(c,len,lg,-1);
for(int i=0;i<len;i++) C[i]=c[i];
}
void solve(int *A,int *B) {
for(len=1,lg=0;len<=n*2;len*=2,lg++) {
if (len==1) {B[0]=pwr(A[0],Mo-2);continue;}
for(int i=0;i<len;i++) p[i]=A[i];
pre(len*2);
NTT(p,B,h,len*2,lg+1);
for(int i=0;i<len*2;i++) h[i]=Mo-h[i];
(h[0]+=2)%=Mo;
NTT(h,B,B,len*2,lg+1);
for(int i=len;i<len*2;i++) B[i]=0;
}
}
int fact[N],Inv[N],F[N],G[N],Inv_G[N],H[N];
int main() {
scanf("%d",&n);
fact[0]=1;fo(i,1,n) fact[i]=(ll)fact[i-1]*i%Mo;
Inv[n]=pwr(fact[n],Mo-2);fd(i,n-1,0) Inv[i]=(ll)Inv[i+1]*(i+1)%Mo;
fo(i,1,n) F[i]=(ll)pwr(2,(ll)i*(i-1)/2)*Inv[i-1]%Mo;
fo(i,0,n) G[i]=(ll)pwr(2,(ll)i*(i-1)/2)*Inv[i]%Mo;
solve(G,Inv_G);
pre(len);NTT(F,Inv_G,H,len,lg);
printf("%d\n",(ll)H[n]*fact[n-1]%Mo);
}
UPD:听说一次倍增只需要3次DFT QwQ
结果我常数直接*2压倒性劣势