1. 程式人生 > >BZOJ2467 [中山市選2010]生成樹

BZOJ2467 [中山市選2010]生成樹

不同 clu ron inpu 表示 現在 輸出 void ++

Description

有一種圖形叫做五角形圈。一個五角形圈的中心有1個由n個頂點和n條邊組成的圈。在中心的這個n邊圈的每一條邊同時也是某一個五角形的一條邊,一共有n個不同的五角形。這些五角形只在五角形圈的中心的圈上有公共的頂點。如圖0所示是一個4-五角形圈。 現在給定一個n五角形圈,你的任務就是求出n五角形圈的不同生成樹的數目。還記得什麽是圖的生成樹嗎?一個圖的生成樹是保留原圖的所有頂點以及頂點的數目減去一這麽多條邊,從而生成的一棵樹。 註意:在給定的n五角形圈中所有頂點均視為不同的頂點。

Input

輸入包含多組測試數據。第一行包含一個正整數T,表示測試數據數目。每組測試數據包含一個整數n

( 2<=N<=100),代表你需要求解的五角形圈中心的邊數。

Output

對每一組測試數據,輸出一行包含一個整數x,表示n五角形圈的生成樹數目模2007之後的結果。

Sample Input

1
2

Sample Output

40 矩陣樹定理裸題 不過數據再大一點也可以做,因為 $ans_n=n*4*5^{n-1}$
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5
#include<cmath> 6 using namespace std; 7 int n,cnt,a[501][501],ans,Mod=2007; 8 void guass() 9 {int i,j,k; 10 n=cnt-1; 11 ans=1; 12 for (i=1;i<=n;i++) 13 { 14 for (j=1;j<=n;j++) 15 { 16 a[i][j]=(a[i][j]+Mod)%Mod; 17 } 18 } 19 for (i=1;i<=n;i++)
20 { 21 for (j=i+1;j<=n;j++) 22 { 23 while (a[j][i]) 24 { 25 int t=a[i][i]/a[j][i]; 26 for (k=i;k<=n;k++) 27 { 28 a[i][k]=(a[i][k]-1ll*t*a[j][k]%Mod+Mod)%Mod; 29 swap(a[i][k],a[j][k]); 30 } 31 ans*=-1; 32 } 33 } 34 ans=1ll*ans*a[i][i]%Mod; 35 } 36 ans=(ans+Mod)%Mod; 37 } 38 int main() 39 {int i,j,T; 40 cin>>T; 41 while (T--) 42 { 43 cin>>n;cnt=0; 44 memset(a,0,sizeof(a)); 45 for (i=1;i<=n-1;i++) 46 { 47 for (j=1;j<=4;j++) 48 { 49 ++cnt; 50 a[cnt][cnt+1]--;a[cnt][cnt]++; 51 a[cnt+1][cnt]--;a[cnt+1][cnt+1]++; 52 } 53 ++cnt; 54 a[cnt][cnt-4]--;a[cnt][cnt]++; 55 a[cnt-4][cnt]--;a[cnt-4][cnt-4]++; 56 cnt--; 57 } 58 for (j=1;j<=3;j++) 59 { 60 ++cnt; 61 a[cnt][cnt+1]--;a[cnt][cnt]++; 62 a[cnt+1][cnt]--;a[cnt+1][cnt+1]++; 63 } 64 ++cnt; 65 a[1][cnt]--;a[1][1]++; 66 a[cnt][1]--;a[cnt][cnt]++; 67 a[1][cnt-3]--;a[1][1]++; 68 a[cnt-3][1]--;a[cnt-3][cnt-3]++; 69 guass(); 70 cout<<ans<<endl; 71 } 72 }

BZOJ2467 [中山市選2010]生成樹