題目

Source

http://www.tsinsen.com/A1493

Description

剛剛解決完電力網(wǎng)絡(luò)的問(wèn)題, 阿貍又被領(lǐng)導(dǎo)的任務(wù)給難住了.
剛才說(shuō)過(guò), 阿貍的國(guó)家有n個(gè)城市, 現(xiàn)在國(guó)家需要在某些城市對(duì)之間建立一些貿(mào)易路線, 使得整個(gè)國(guó)家的任意兩個(gè)城市都直接或間接的連通.
為了省錢, 每?jī)蓚€(gè)城市之間最多只能有一條直接的貿(mào)易路徑. 對(duì)于兩個(gè)建立路線的方案, 如果存在一個(gè)城市對(duì), 在兩個(gè)方案中是否建立路線不一樣, 那么這兩個(gè)方案就是不同的, 否則就是相同的. 現(xiàn)在你需要求出一共有多少不同的方案.
好了, 這就是困擾阿貍的問(wèn)題. 換句話說(shuō), 你需要求出n個(gè)點(diǎn)的簡(jiǎn)單(無(wú)重邊無(wú)自環(huán))無(wú)向連通圖數(shù)目.
由于這個(gè)數(shù)字可能非常大, 你只需要輸出方案數(shù)mod 1004535809(479 * 2 ^ 21 + 1)即可.

Input

僅一行一個(gè)整數(shù)n(<=130000)

Output

僅一行一個(gè)整數(shù), 為方案數(shù) mod 1004535809.

Sample Input

3
4
100000

Sample Output

4
38
829847355

?

分析

樓教主男人八題也又一題是和這題一樣，n個(gè)有標(biāo)號(hào)點(diǎn)能構(gòu)成的簡(jiǎn)單無(wú)向連通圖數(shù)，不過(guò)這題題目數(shù)據(jù)大很多。做法如下：

• 補(bǔ)集轉(zhuǎn)化，設(shè)$f[n]$表示n個(gè)帶標(biāo)號(hào)的點(diǎn)能構(gòu)成的簡(jiǎn)單圖數(shù)，而只要利用$f[n]$減去不連通的數(shù)目就是要求的答案了。其中$f[n]=2^{n(n-1)/2}$，因?yàn)橥耆珗D有$n(n-1)/2$條邊，對(duì)于每條邊選或不選來(lái)構(gòu)成一張簡(jiǎn)單圖。

• $dp[n]$表示n個(gè)帶標(biāo)號(hào)的點(diǎn)能構(gòu)成的簡(jiǎn)單無(wú)向連通圖數(shù)
• 考慮通過(guò)固定一個(gè)點(diǎn)，枚舉這個(gè)點(diǎn)所在連通塊的大小來(lái)轉(zhuǎn)移：$dp[n]\ =\ f[n]-\sum_{i=1}^{n-1}C_{n-1}^{i-1}dp[i]f[n-i]$
• 整理成卷積形式：$d[n]\ =\ f[n]-(n-1)!\sum_{i=1}^{n-1}((i-1)!)^{-1}d[i] \times ((n-i)!)^{-1}f[n-i]$
• 最后就是用分治FFT來(lái)做了。

?

代碼

#include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define MAXN 262144//const long long P=50000000001507329LL; // 190734863287 * 2 ^ 18 + 1 const int P=1004535809; // 479 * 2 ^ 21 + 1 //const int P=998244353; // 119 * 2 ^ 23 + 1 const int G=3;long long mul(long long x,long long y){return x*y%P; } long long qpow(long long x,long long k){long long ret=1;while(k){if(k&1) ret=mul(ret,x);k>>=1;x=mul(x,x);}return ret; }long long wn[25]; void getwn(){for(int i=1; i<=21; ++i){int t=1<<i;wn[i]=qpow(G,(P-1)/t);} }int len; void NTT(long long y[],int op){for(int i=1,j=len>>1,k; i<len-1; ++i){if(i<j) swap(y[i],y[j]);k=len>>1;while(j>=k){j-=k;k>>=1;}if(j<k) j+=k;}int id=0;for(int h=2; h<=len; h<<=1) {++id;for(int i=0; i<len; i+=h){long long w=1;for(int j=i; j<i+(h>>1); ++j){long long u=y[j],t=mul(y[j+h/2],w);y[j]=u+t;if(y[j]>=P) y[j]-=P;y[j+h/2]=u-t+P;if(y[j+h/2]>=P) y[j+h/2]-=P;w=mul(w,wn[id]);}}}if(op==-1){for(int i=1; i<len/2; ++i) swap(y[i],y[len-i]);long long inv=qpow(len,P-2);for(int i=0; i<len; ++i) y[i]=mul(y[i],inv);} } void Convolution(long long A[],long long B[],int n){for(len=1; len<(n<<1); len<<=1);for(int i=n; i<len; ++i){A[i]=B[i]=0;}NTT(A,1); NTT(B,1);for(int i=0; i<len; ++i){A[i]=mul(A[i],B[i]);}NTT(A,-1); }long long A[MAXN],B[MAXN];long long d[MAXN],f[MAXN]={1},fact[MAXN]={1},fact_inv[MAXN]={1}; /*d[n] = f[n]-fact[n-1]*Σ(fact_inv[i-1]*d[i]*f[n-i]*fact_inv[n-i]) */ void cdq(int l,int r){if(l==r) return;int mid=l+r>>1;cdq(l,mid);for(int i=l; i<=mid; ++i){A[i-l]=mul(fact_inv[i-1],d[i]);}for(int i=mid-l+1; i<=r-l; ++i) A[i]=0;for(int i=0; i<=r-l; ++i){B[i]=mul(fact_inv[i],f[i]);}Convolution(A,B,r-l+1);for(int i=mid+1; i<=r; ++i){d[i]-=mul(fact[i-1],A[i-l]);if(d[i]<0) d[i]+=P;}cdq(mid+1,r); }int main(){getwn();for(int i=1; i<=130000; ++i){fact[i]=mul(fact[i-1],i);fact_inv[i]=qpow(fact[i],P-2);d[i]=f[i]=qpow(2,(i-1LL)*i/2);}cdq(1,130000);int n;while(~scanf("%d",&n)){printf("%lld\n",d[n]);}return 0; }

?

轉(zhuǎn)載于:https://www.cnblogs.com/WABoss/p/5929692.html

總結(jié)

以上是生活随笔為你收集整理的Tsinsen A1493 城市规划（DP + CDQ分治 + NTT）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。