标签：pr cnt return int ++ 计数 prufer

树的计数

已知每个点的度数，求满足条件的树的个数。

思路

对于无标号图的计数，考虑 \(prufer\) 序列。

注意到 \(prufer\) 序列的性质：在 \(prufer\) 序列中，顶点出现的次数就是其度数 \(-1\) 。于是，可以直接利用多重集的排列数得出结论：

\(\displaystyle ans = \frac{(n-2)!}{(d_1-1)(d_2-1)\dots (d_n-1)}\)

由于 \(n\le 150\) ，考虑对上下分解质因数进行计算。

代码：

#include <bits/stdc++.h>
using LL = long long;
const int N = 200;
int n, pr[N], c[N], cnt;
void split(int x)
{
    for (int i = 2; i * i <= x; i++)
    {
        if (x % i == 0)
        {
            int k = 1;
            while (pr[k] < i) k++;
            while (x % i == 0) x /= i, c[k]--;
        }
    }
    if (x > 1)
    {
        int k = 1;
        while (pr[k] < x) k++;
        c[k]--;
    }
}
void sie(int x)
{
    for (int i = 2; i * i <= x; i++)
        if (x % i == 0)
        {
            int k = 1;
            while (pr[k] < i) k++;
            if (pr[k] != i)
            {
                pr[++cnt] = i;
                while (x % i == 0) x /= i, c[cnt]++;
            }
            else
                while (x % i == 0) x /= i, c[k]++;
        }
    if (x > 1)
    {
        int k = 1;
        while (pr[k] < x) k++;
        if (pr[k] != x) pr[++cnt] = x, c[cnt] = 1;
        else c[k]++;
    } 
}
int main()
{
    memset(pr, 0x3f, sizeof(pr));
    scanf("%d", &n);
    int d, tot = 0;
    if (n == 1)
    {
        scanf("%d", &d);
        if (d == 0) return puts("1"), 0;
        else return puts("0"), 0;
    }   
    for (int i = 2; i <= n - 2; i++)
        sie(i);
    for (int i = 1; i <= n; i++)
    {
        scanf("%d", &d);
        if (d < 1 || d > n) return puts("0"), 0;
        tot += d - 1;
        for (int i = 2; i <= d - 1; i++) split(d - 1);
    }
    if (tot != n - 2) return puts("0"), 0;
    LL ans = 1;
    for (int i = 1; i <= cnt; i++) while (c[i]--) ans *= pr[i];
    std::cout << ans << "\n";
    return 0;
}

标签：pr,cnt,return,int,++,计数,prufer
来源： https://www.cnblogs.com/mklzc/p/16525013.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。