第一次学prufer数列相关,一开始觉得可能可以树形DP但发现复杂度大的我无法想象
然后发现是prufer数列
具体相关的知识可以自行去百度,下面是我搜到的一些资料
将树转化成Prufer数列的方法
一种生成Prufer序列的方法是迭代删点,直到原图仅剩两个点。对于一棵顶点已经经过编号的树T,顶点的编号为{1,2,...,n},在第i步时,移去所有叶子节点(度为1的顶点)中标号最小的顶点和相连的边,并把与它相邻的点的编号加入Prufer序列中,重复以上步骤直到原图仅剩2个顶点。
例子
以右边的树为例子,首先在所有叶子节点中编号最小的点是2,和它相邻的点的编号是3,将3加入序列并删除编号为2的点。接下来删除的点是4,5被加入序列,然后删除5,1被加入序列,1被删除,3被加入序列,此时原图仅剩两个点(即3和6),Prufer序列构建完成,为{3,5,1,3}
将Prufer数列转化成树的方法
设{a1,a2,..an-2}为一棵有n个节点的树的Prufer序列,另建一个集合G含有元素{1..n},找出集合中最小的未在Prufer序列中出现过的数,将该点与Prufer序列中首项连一条边,并将该点和Prufer序列首项删除,重复操作n-2次,将集合中剩余的两个点之间连边即可。
例子
仍为上面的树,Prufer序列为{3,5,1,3},开始时G={1,2,3,4,5,6},未出现的编号最小的点是2,将2和3连边,并删去Prufer序列首项和G中的2。接下来连的边为{4,5},{1,5},{1,3},此时集合G中仅剩3和6,在3和6之间连边,原树恢复
然后我们发现一个prufer数列非常nice的性质:
就是一个节点在prufer数列中出现的次数是这个节点度数减一
很好证明,如果这个点是叶子节点那么不会出现的啊,即为0
如果不是叶子节点,那么每一次他的叶子节点被删时,他的度数减一,减到为一时,即为叶子节点了,所以出现次数为度数减一
然后排列组合一算就可以了
发现题目中n很小,直接质因数分解一下就可以了,不需要高精度
然后的话,通过prufer的定义,我们会发现,到最后prufer数列的中的数一共有n-2个
至于prufer数列的个数一共有n-2的证明我补上一下,其实非常显然,我每加一个点必须先删掉一个点
那么非常显然的是我一共要删n-2个点,那么prufer序列中的数一共有n-2个
那么所有的数的度数减一的和不为n-2就没有方案了
下面是代码
#include<iostream>
#include<algorithm>
using namespace std;
int n,tot,p[2][151],x;
long long ans(1);
void cal(int x,int k){
for(int i=2;;i++){
while(x%i==0){
p[k][i]++,x/=i;
}
if(x==1){
break;
}
}
}//直接把数字质因数分解记录进数组里面
long long qpow(int x,int y){
if(y==1){
return x;
}else{
long long fanhui=qpow(x,y/2);
if(y%2==1){
fanhui*=x;
}
return fanhui;
}
}//快速幂
int main(){
cin>>n;
for(int i=1;i<=n;i++){
cin>>x;
if(x>2){
for(int j=1;j<=x-1;j++){
cal(j,1);
}
}//当x>2时,拆开记录一个p数组中
if(x>1){
tot+=x-1;
}//记所有点的度数减一然后加起来
}
if(n==1){
cout<<!x<<endl;
return 0;
}//特判只有一个节点的情况
if(tot!=n-2){
cout<<0<<endl;
return 0;
}//所有点的度数减一之和不为n-2的话就无解
for(int i=2;i<=tot;i++){
cal(i,0);
}//将tot!拆开存到另一个p数组里
for(int i=1;i<=150;i++){
if(p[0][i]-p[1][i]>0){//若某一因子多了,那么算进答案中
ans*=(long long)qpow(i,p[0][i]-p[1][i]);
}
}//算答案
cout<<ans<<endl;
return 0;
}
/*
in:
4
2 1 2 1
out:
2
*/