JZOJ 5033. 【NOI2017模拟3.28】A
2024-10-21 10:01:24
A
题面
思路
非常抽象地让你构造树,很容易想到 \(prufer\) 序列(如果你会的话)
说明一下:\(prufer\) 序列可以唯一确定一颗树的形态
若树的节点个数为 \(n\),那么 \(prufer\) 序列长度为 \(n-2\) ,且一个节点出现的个数为它的度数减一(不要问我为什么,因为 \(prufer\) 序列就是这样的)
那么我们就考虑 \(dp\) 了
设 \(f_{i,j,k}\) 表示考虑前 \(i\) 个数,选出 \(j\) 个数,当前 \(prufer\) 序列长度为 \(k\)。
为何要设 \(k\) ?因为一个节点在 \(prufer\) 序列中出现可能不止一次
考虑转移: \(f_{i,j,k} = \sum_{l=1}^{\min(a_i-1,k)}\binom{k}{l}f_{i-1,j-1,k-l}+f_{i-1,j,k}\)
\(f_{i-1,j,k}\) 意思是第 \(i\) 位不选
选的话,\(l\) 枚举选多少个,\(\binom{k}{l}\) 表示选了之后放到序列中的方案数
那么答案如何计算?
\(ans_x=\sum_{j=1}^x\binom{n-j}{x-j}f_{n,j,x-2}\)
意思是考虑 \(prufer\) 序列中数的种数,用 \(j\) 个数凑出长为 \(x-2\) 的序列。
因为叶子节点不会出现在序列中,所以我们再从剩下 \(n-j\) 个数中选出还差的 \(x-j\) 个数
\(Code\)
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
typedef long long LL;
const int N = 55;
const LL P = 1e9 + 7;
LL f[N][N][N] , fac[N];
int a[N] , n , T;
inline LL fpow(LL x , LL y)
{
LL res = 1;
while (y)
{
if (y & 1) res = res * x % P;
y >>= 1 , x = x * x % P;
}
return res;
}
inline LL C(int n , int m){return fac[n] * fpow(fac[m] * fac[n - m] % P , P - 2) % P;}
int main()
{
freopen("a.in" , "r" , stdin);
freopen("a.out" , "w" , stdout);
fac[0] = 1;
for(register int i = 1; i <= 52; i++) fac[i] = (i * 1LL * fac[i - 1]) % P;
scanf("%d" , &T);
while (T--)
{
scanf("%d" , &n);
for(register int i = 1; i <= n; i++) scanf("%d" , &a[i]);
memset(f , 0 , sizeof f);
f[0][0][0] = 1;
for(register int i = 1; i <= n; i++)
for(register int j = 0; j <= i; j++)
for(register int k = j; k <= n - 2; k++)
{
f[i][j][k] = f[i - 1][j][k];
if (j != 0) for(register int l = 1; l <= min(a[i] - 1 , k); l++)
f[i][j][k] = (f[i][j][k] + f[i - 1][j - 1][k - l] * C(k , l)) % P;
}
printf("%lld " , (LL)n);
LL ans;
for(register int x = 2; x <= n; x++)
{
ans = 0;
for(register int j = 0; j <= x; j++) ans = (ans + f[n][j][x - 2] * C(n - j , x - j) % P) % P;
printf("%lld " , ans);
}
printf("\n");
}
}
最新文章
- PLSQL查询字段为科学计数法,修正显示
- Wireshark插件编写
- C# params关键字
- SQL SERVER 创建作业
- linux mysql开启远程链接
- sqlsevrer中output的用法
- qsort函数、sort函数 (精心整理篇)
- 未在本地计算机上注册“Microsoft.Jet.OLEDB.4.0”提供程序。
- Smtp协议与Pop3协议的简单实现
- 基于visual Studio2013解决C语言竞赛题之0614递归大元素
- 教你用SVG画出一条龙
- [poj3687]Labeling Balls_拓扑排序
- 算法面试题(python)——如何找出数组中出现一次的数
- A - Points and Segments CodeForces - 429E
- linux vue项目+npm run build + nginx
- Android 单元测试覆盖率计算
- IntelliJ IDEA安装主题详细步骤
- node学习笔记3——文件操作fs
- Hibernate实体映射转换列值
- 【转载】屏幕坐标向3维坐标的转化-DXUT的CD3DArcBall类
热门文章
- C#调用WPS转换文档到PDF的的实现代码。
- IOS移动端 -webkit-overflow-scrollin属性造成的问题
- Kubernetes(k8s)存储管理之数据卷volumes(一):volumes的引入和emptyDir数据卷
- Docker容器入门到精通
- Service层和Dao层的一些自我理解(╥╯^╰╥)(╥╯^╰╥)(学了这么久,这玩意儿似懂非懂的)
- Redis 如何批量设置过期时间?PIPLINE的使用
- 低代码开发平台YonBuilder移动开发,开发阅读APP教程
- 一文告诉你AVM中设置字体的方法
- [OpenCV实战]3 透明斗篷
- react 高效高质量搭建后台系统 系列 —— antd和样式