prufer编码(无根树的唯一表示)

树的编码:prufer序列初始为空。每次从树上选出一个编号最小的叶子节点,然后将与该叶子节点相邻的那个节点的编号写入prufer序列的末尾,之后从树上删掉这个叶子节点。循环这个步骤n-2次,最后得到一个长度为n-2的prufer序列
编码到树:首先,将每个节点的度数设为1加上该节点在prufer序列中出现的次数。然后以下循环执行n-2次。第i次循环,选择此时度数为1的编号最小的节点u,将其与此时prufer序列的第i个元素v连边,然后将u和v的度数都减去1。这n-2次执行完之后,仅剩下两个节点他们的度数都是1,将这两个点连边,这样就得到一个有n-1条边的树
性质:一个点的入度为d,那么它最多在编码中出现d-1次。编码共n-2长度,每个相同的数字出现d-1次,则满足给出每个点的度数要求能生成的不同的树(编码)为(n-2)!/((d_1-1)!(d_2-1)!×⋯×(d_n-1)!)

bzoj1211
应用prufer编码+阶乘的质因子分解

#include<cmath>
#include<cstdio>
const int maxn=155;
int n,tot,cnt;
int pri[maxn],d[maxn],num[maxn];
long long ans=1;
long long s[25];
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9') {if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9') {x*=10;x+=ch-'0';ch=getchar();}
    return x*f;
}
bool jud(int x)
{
    for(int i=2;i<=sqrt(x);i++)
    if(x%i==0) return 0;
    return 1;
}
void getpri()
{
    for(int i=2;i<=150;i++)
    if(jud(i)) pri[++cnt]=i;
}
void solve(long long x,int f)  //按照指数分解质因数 
{
    for(int i=1;i<=cnt;i++)
    {
        if(x<=1) return;
        while(x%pri[i]==0) {num[i]+=f;x/=pri[i];}
    }
}
int main()
{
    s[1]=1;
    for(int i=2;i<=22;i++) s[i]=s[i-1]*i;
    getpri();
    n=read();
    if(n==1)
    {
        int x=read();
        if(x==0) printf("1");
        else printf("0");
        return 0;
    }
    for(int i=1;i<=n;i++)
    {
        d[i]=read();
        if(d[i]==0) {printf("0");return 0;}
        d[i]--;
        tot+=d[i];
    }
    if(tot!=n-2) {printf("0"); return 0;}
    solve(s[n-2],1);  //计算阶乘并将结果分解质因数 
    for(int i=1;i<=n;i++) solve(s[d[i]],-1); //同上 
    for(int i=1;i<=cnt;i++)
        while(num[i]--) ans*=pri[i];  //统计结果 
    printf("%lld",ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值