题目:

洛谷2624

分析:

本文中所有的 “树” 都是带标号的。

介绍一种把树变成一个序列的工具:Prufer 序列。

对于一棵 \(n\) 个结点的树,每次选出一个叶子(度数为 \(1\) 的结点),将唯一的那个与它相连的点标号加入 Prufer 序列末尾,然后删去这个叶子及其所连的边,直到最后剩下两个点和一条边。由于每次删且仅删一个点和一条边,所以 Prufer 序列长度为 \(n-2\) 。点 \(a\) 在序列中每次出现都意味着一条与它相连的边被删去了,一直删到 \(a\) 度数为 \(1\) (叶子)为止。所以任意一点在 Prufer 序列中的出现次数就是它的度数减 \(1\) 。

再扯远点。把一个 Prufer 序列还原成树的方式是:定义集合 \(V=\{i|1\leq i \leq n\}\) 每次选出 \(V\) 中编号最小的没有在当前序列出现过的点 \(u\),与 Prufer 序列第一个元素 \(v\) 相连,然后从 \(V\) 中删除 \(u\),并删除 Prufer 序列的第一个元素。最后 \(V\) 中剩余两个元素,将它们相连。

按照上述方法,任意一个由整数 \(1\) 到 \(n\) 组成的序列都可以生成一棵 \(n\) 个结点的树,且任意一棵 \(n\) 个结点的树也可以生成一个这样的序列,所以序列和树一一对应。由此可得推论:\(n\) 个结点的树有 \(n^{n-2}\) 种。

从 Prufer 序列的角度考虑可以解决很多树计数的问题。回到本题。设 \(sum=\sum_{d_i\neq -1} di\) 和 \(num=\sum_{d_i\neq -1} 1\) (即:有 \(num\) 个结点被指定度数,它们的度数之和为 \(sum\) ),则问题被转化为:求一个长为 \(n-2\) 个序列,使其中对于任意 \(i(d_i\neq -1)\) 满足 \(i\) 的出现次数是 \(d_i-1\) 。首先从序列种选出 \(sum-num\) 个位置来放置这些有限制的数,根据可重集合排列公式 \(\frac{(\sum a_i)!}{\prod (a_i!)}\) ,方案数是:

\[C_{n-2}^{sum-num}\frac{(sum-num)!}{\prod_{d_i\neq -1} ((d_i-1)!)}
\]

然后剩下 \(n-2-(sum-num)\) 个位置可以从 \(n-num\) 个数中任取,所以最终答案是:

\[C_{n-2}^{sum-num}\frac{(sum-num)!}{\prod_{d_i\neq -1} ((d_i-1)!)}\cdot (n-num)^{(n-2-(sum-num))}
\]

很不幸发现这题没模数,也看不懂什么分解质因数的做法,乖乖写高精度吧……(其实写这篇博客的主要目的是记录高精度板子qwq

代码:

注意一定要把\(\prod_{d_i\neq -1}((d_i-1)!)\)先全部乘起来再除!!!

算 \(n\) 遍除法会 TLE !!!

(我就不幸因为这个 T 了一发……

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cctype>
using namespace std; namespace zyt
{
typedef long long ll;
const int N = 1e3 + 10;
class long_long_long
{
private:
typedef long_long_long Tril;
const static int DIGIT = 4, BASE = 10000, LEN = (3e3 + 100) / DIGIT;
int data[LEN], len;
public:
long_long_long(const int x = 0)
{
memset(data, 0, sizeof(data));
len = 0;
data[0] = x;
while (data[len])
{
data[len + 1] = data[len] / BASE;
data[len++] %= BASE;
}
}
bool operator == (const Tril &b) const
{
if (len != b.len)
return false;
for (int i = 0; i < len; i++)
if (data[i] != b.data[i])
return false;
return true;
}
bool operator < (const Tril &b) const
{
if (len != b.len)
return len < b.len;
for (int i = len - 1; i >= 0; i--)
if (data[i] != b.data[i])
return data[i] < b.data[i];
return false;
}
Tril operator + (const Tril &b) const
{
Tril ans;
ans.len = max(len, b.len);
for (int i = 0; i < ans.len; i++)
{
ans.data[i] += data[i] + b.data[i];
ans.data[i + 1] += ans.data[i] / BASE;
ans.data[i] %= BASE;
}
if (ans.data[ans.len])
++ans.len;
return ans;
}
Tril operator - (const Tril &b) const
{
Tril ans;
ans.len = len;
for (int i = 0; i < len; i++)
{
ans.data[i] += data[i] - b.data[i];
if (ans.data[i] < 0)
{
--ans.data[i + 1];
ans.data[i] += BASE;
}
}
while (ans.len && !ans.data[ans.len - 1])
--ans.len;
return ans;
}
Tril operator * (const Tril &b) const
{
static ll tmp[LEN];
Tril ans;
if (*this == 0 || b == 0)
return ans;
ans.len = len + b.len - 1;
memset(tmp, 0, sizeof(ll[len + b.len + 2]));
for (int i = 0; i < len; i++)
for (int j = 0; j < b.len; j++)
tmp[i + j] += data[i] * b.data[j];
for (int i = 0; i < ans.len; i++)
{
tmp[i + 1] += tmp[i] / BASE;
ans.data[i] = tmp[i] % BASE;
}
while (tmp[ans.len])
{
tmp[ans.len + 1] += tmp[ans.len] / BASE;
ans.data[ans.len] = tmp[ans.len] % BASE;
++ans.len;
}
return ans;
}
Tril operator / (const Tril &b) const
{
Tril ans, rest;
for (int i = len - 1; i >= 0; i--)
{
rest = rest * BASE + data[i];
if (rest < b)
continue;
int l = 1, r = BASE - 1, anss;
while (l <= r)
{
int mid = (l + r) >> 1;
if (rest < b * mid)
r = mid - 1;
else
l = mid + 1, anss = mid;
}
ans.data[i] = anss;
rest = rest - b * anss;
if (!ans.len)
ans.len = i + 1;
}
return ans;
}
void print()
{
if (!len)
printf("0");
else
{
printf("%d", data[len - 1]);
for (int i = len - 2; i >= 0; i--)
printf("%0*d", DIGIT, data[i]);
}
}
};
typedef long_long_long Tril;
Tril ans, fac[N];
int n;
Tril C(const int n, const int m)
{
return fac[n] / fac[m] / fac[n - m];
}
Tril power(Tril a, int b)
{
Tril ans = 1;
while (b)
{
if (b & 1)
ans = ans * a;
a = a * a;
b >>= 1;
}
return ans;
}
int d[N];
int work()
{
int sum = 0, num = 0;
scanf("%d", &n);
fac[0] = 1;
for (int i = 1; i <= n; i++)
{
fac[i] = fac[i - 1] * i;
scanf("%d", &d[i]);
if (d[i] == -1)
continue;
sum += d[i], ++num;
}
ans = fac[sum - num];
for (int i = 1; i <= n; i++)
if (d[i] > 0)
ans = ans / fac[d[i] - 1];
ans = ans * C(n - 2, sum - num);
ans = ans * power(n - num, n - 2 - (sum - num));
ans.print();
return 0;
}
}
int main()
{
return zyt::work();
}

【洛谷2624_BZOJ1005】[HNOI2008] 明明的烦恼(Prufer序列_高精度_组合数学)的更多相关文章

  1. 【bzoj1005】[HNOI2008]明明的烦恼 Prufer序列+高精度

    题目描述 给出标号为1到N的点,以及某些点最终的度数,允许在任意两点间连线,可产生多少棵度数满足要求的树? 输入 第一行为N(0 < N < = 1000),接下来N行,第i+1行给出第i ...

  2. bzoj1005: [HNOI2008]明明的烦恼 prufer序列

    https://www.lydsy.com/JudgeOnline/problem.php?id=1005 给出标号为1到N的点,以及某些点最终的度数,允许在任意两点间连线,可产生多少棵度数满足要求的 ...

  3. bzoj 1005: [HNOI2008]明明的烦恼 prufer编号&&生成树计数

    1005: [HNOI2008]明明的烦恼 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 2248  Solved: 898[Submit][Statu ...

  4. BZOJ 1005 [HNOI2008]明明的烦恼 (Prufer编码 + 组合数学 + 高精度)

    1005: [HNOI2008]明明的烦恼 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 5786  Solved: 2263[Submit][Stat ...

  5. bzoj 1005 [HNOI2008] 明明的烦恼 (prufer编码)

    [HNOI2008]明明的烦恼 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 5907  Solved: 2305[Submit][Status][Di ...

  6. BZOJ 1005: [HNOI2008]明明的烦恼 Purfer序列 大数

    1005: [HNOI2008]明明的烦恼 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://www.lydsy.com/JudgeOnline/ ...

  7. BZOJ 1005 [HNOI2008]明明的烦恼 purfer序列,排列组合

    1005: [HNOI2008]明明的烦恼 Description 自从明明学了树的结构,就对奇怪的树产生了兴趣......给出标号为1到N的点,以及某些点最终的度数,允许在任意两点间连线,可产生多少 ...

  8. BZOJ.1005.[HNOI2008]明明的烦恼(Prufer 高精 排列组合)

    题目链接 若点数确定那么ans = (n-2)!/[(d1-1)!(d2-1)!...(dn-1)!] 现在把那些不确定的点一起考虑(假设有m个),它们在Prufer序列中总出现数就是left=n-2 ...

  9. BZOJ 1005 [HNOI2008]明明的烦恼 ★(Prufer数列)

    题意 N个点,有些点有度数限制,问这些点可以构成几棵不同的树. 思路 [Prufer数列] Prufer数列是无根树的一种数列.在组合数学中,Prufer数列是由一个对于顶点标过号的树转化来的数列,点 ...

  10. Luogu P2624 [HNOI2008]明明的烦恼 Prufer+组合+高精

    好的我把标准版过了... 设$ r_i$为$i$的度数 首先,我们设 $ sum = \Sigma r_i-1$,$ tot $ 为所有能够确定度数的点 所以我们有 $ C ^ {sum} _{n-2 ...

随机推荐

  1. python 全栈之路

    目录 Python 全栈之路 一. Python 1. Python基础知识部分 2. Python -函数 3. Python - 模块 4. Python - 面对对象 5. Python - 文 ...

  2. 51nod1485 字母排序

    [题解] 开26棵线段数,记录区间内每种字母的出现次数,修改的时候就用区间设置为一个数操作即可.同时也有平衡树做 #include<cstdio> #include<algorith ...

  3. POJ - 2007 极角排序(Java 实现)

    POJ 2007 将所有的点按逆时针输出 import java.io.*; import java.util.*; public class Main { static class Point im ...

  4. 2.8 补充:shell脚本执行方法

    bash shell 脚本的方法有多种,现在作个小结.假设我们编写好的shell脚本的文件名为hello.sh,文件位置在/data/shell目录中并已有执行权限.   方法一:切换到shell脚本 ...

  5. tomcat上传图片失败

    今天,突然遇到个奇怪的问题,tomcat上传图片时好时坏,经查线上四台服务有一台服务器硬盘满了. 解决一下硬盘空间的问题,有好使了,那么图片上传通过流写到远程对象存储服务中,并没有落盘和硬盘满有哪些关 ...

  6. 为什么Java使用System.getenv()获取刚刚设置的环境变量时为空

    场景: 在Ubuntu下已经启动了Eclipse,然后通过终端设置了环境变量(export $ENV=123),然后通过System.getenv()获取时显示为空. 解释: 环境变量仅在进程树下方, ...

  7. 源码分析-react2-根节点渲染

    //FiberNode{ alternate : '通过该属性和后面的切片进行比较', child : '改切片的子切片', firstEffect : '当前要加入的切片', stateNode : ...

  8. [C++]_[获取Utf8字符串的字符个数和子字符串]

    场景: 1.有时候须要统计utf8字符串的个数,单纯统计字节个数是不行的. 2.有时候也须要获取从某个位置開始的n个连续字符用于显示或计算. static int GetUtf8LetterNumbe ...

  9. OpenCV2马拉松第25圈——直线拟合与RANSAC算法

    计算机视觉讨论群162501053 转载请注明:http://blog.csdn.net/abcd1992719g/article/details/28118095 收入囊中 最小二乘法(least ...

  10. PHPthinking官方论坛招募版主

    时间飞逝.就在昨天,我们PHPthinking的官方论坛刚刚上线了我们自己的论坛! 欢迎大家注冊账号,活跃在论坛的大家庭中,我们会及时关注论坛公布的全部内容.在开发学习的过程中,遇到的不论什么问题,有 ...