题意:一个有n个结点的树,设它的结点分别为v1, v2, …, vn,
已知第i个结点vi的度数为di,问满足这样的条件的不同的树有多少棵。
其中1<=n<=150,输入数据保证满足条件的树不超过10^17个。
思路:
Martix67:
一个有趣的推广是,n个节点的度依次为D1, D2, …, Dn的无根树共有(n-2)! / [ (D1-1)!(D2-1)!..(Dn-1)! ]个,因为此时Prüfer编码中的数字i恰好出现Di-1次。
POPOQQQ:
注意此题无解需要输出0
当n!=1&&d[i]==0时 输出0
当Σ(d[i]-1)!=n-2时输出0
1 var a,b,prime,cnt:array[0..1000]of longint; 2 n,m,s,i:longint; 3 ans:int64; 4 flag:boolean; 5 6 function mult(x,y:longint):int64; 7 var tmp:int64; 8 begin 9 mult:=1; tmp:=x; 10 while y>0 do 11 begin 12 if y and 1=1 then mult:=mult*tmp; 13 tmp:=tmp*tmp; 14 y:=y>>1; 15 end; 16 end; 17 18 function ins(x,y:longint):longint; 19 var i,k:longint; 20 begin 21 for i:=1 to m do 22 begin 23 k:=x; 24 while k>0 do 25 begin 26 cnt[i]:=cnt[i]+k div prime[i]*y; 27 k:=k div prime[i]; 28 end; 29 end; 30 end; 31 32 function su(x:longint):boolean; 33 var i:longint; 34 begin 35 if x=2 then exit(true); 36 for i:=2 to trunc(sqrt(x)) do 37 if x mod i=0 then exit(false); 38 exit(true); 39 end; 40 41 begin 42 assign(input,'bzoj1211.in'); reset(input); 43 // assign(output,'bzoj1211.out'); rewrite(output); 44 readln(n); 45 for i:=2 to 150 do 46 if su(i) then begin inc(m); prime[m]:=i; end; 47 for i:=1 to n do 48 begin 49 read(a[i]); 50 inc(b[a[i]]); 51 s:=s+a[i]-1; 52 end; 53 if n-2<>s then 54 begin 55 writeln(0); 56 exit; 57 end; 58 flag:=true; 59 for i:=1 to n do 60 if a[i]=0 then begin flag:=false; break; end; 61 if (n>1)and( not flag) then 62 begin 63 writeln(0); 64 exit; 65 end; 66 67 ins(n-2,1); 68 for i:=1 to n do 69 if b[i]-1>0 then ins(b[i]-1,-1); 70 ans:=1; 71 for i:=1 to m do ans:=ans*mult(prime[i],cnt[i]); 72 writeln(ans); 73 74 close(input); 75 //close(output); 76 end.