题目:

洛谷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. Java 十二周总结

  2. Maven学习总结(29)——Maven项目的pom.xml中log4j2配置

        <dependency>         <groupId>org.apache.logging.log4j</groupId>         <a ...

  3. 泛型转换https://www.cnblogs.com/eason-chan/p/3633210.html

    import java.lang.reflect.ParameterizedType;import java.lang.reflect.Type;//总结1.st.getClass==Student. ...

  4. 初次使用Let's encrypt

    wget --no-check-certificate -O shadowsocks.sh https://raw.githubusercontent.com/teddysun/shadowsocks ...

  5. HDU 1274 递归拼接字符串

    题目大意: 根据所给的数字,表示其相连的字符的输出个数,或是下一个括号中的所有字符的输出个数 每一个相互对应的 '(' 和 ')' 中的所有字母均作为一组数据处理 在每一次dfs过程中都处理好这样一个 ...

  6. [luoguP1069] 细胞分裂(数论)

    传送门 分解质因数,不说了 这题坑了我2个多小时 教训 不熟悉位运算的优先级一定要加括号!!!! #include <cstdio> #include <iostream> # ...

  7. codevs2486 太鼓达人

    题目描述 Description 七夕祭上,Vani牵着cl的手,在明亮的灯光和欢乐的气氛中愉快地穿行.这时,在前面忽然出现了一台太鼓达人机台,而在机台前坐着的是刚刚被精英队伍成员XLk.Poet_s ...

  8. Java配置方式读取外部的资源配置文件

    通过@PropertySource可以指定读取的配置文件,通过@Value注解获取值,具体用法: package cn.qlq; import org.springframework.context. ...

  9. nyoj_127_星际之门(一)_201403282033

    星际之门(一) 时间限制:3000 ms  |  内存限制:65535 KB 难度:3   描述 公元3000年,子虚帝国统领着N个星系,原先它们是靠近光束飞船来进行旅行的,近来,X博士发明了星际之门 ...

  10. Ubuntu 16.04安装VLC播放器,替代系统默认播放器

    VLC播放器应该说是开源项目中最好的视频播放器,但功能不止于视频播放,还有视频直播等等.可以通过安装字幕插件搜索字母等. 安装步骤: 1.安装: sudo add-apt-repository ppa ...