【[国家集训队] Crash 的文明世界】
先写一个五十分的思路吧
首先这道题有一个弱化版
相当于\(k=1\),于是就是一个非常简单的树形\(dp\)的\(up\ \ and\ \ down\)思想
但是我们现在要求的是这个柿子了
\]
感觉这个东西很组合数学啊,感觉这个柿子像是天生为二项式定理准备的
我们还是考虑树形\(dp\)
在第一遍\(up\)的时候,我们设\(dp[i][k]\)表示
\]
\(j\in{i}\)表示\(j\)在\(i\)子树内部
于是我们考虑一下化这个柿子
到达\(i\)肯定要先达到\(i\)的一个儿子,于是就有
\]
我们用二项式定理来将这个柿子展开
\]
后面两个\(\sum\)换一下位置
\]
之后就会惊奇的发现\(\sum_{t\in{j}}dis(t,j)^r\)就是\(dp[j][r]\),于是现在就有了
\]
这就是\(up\)的转移方程式,\(down\)的方程式也很好推
\(down\)的时候\(dp[i][k]\)表示的不仅局限于\(i\)的子树内部了,而是整棵树了
到达\(i\)首先要到达\(fa[i]\),于是就有
\]
\]
我们的要求不就是\(j\)不能来自于\(i\)内部吗,于是我们大力容斥就好了
我们把来自于\(i\)子树内部的答案减掉,于是就有了一个可以更新的柿子了
现在就有了一个我们就可以转移了,复杂度大概是\(O(nk^2)\)
代码
#include<iostream>
#include<cstdio>
#include<cstring>
#define re register
#define maxn 50005
#define int long long
const int mod=10007;
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
struct E
{
int v,nxt;
}e[maxn<<1];
int head[maxn],deep[maxn];
int n,m,num;
long long dp[maxn][151];
long long now[151];
long long c[151][151];
void dfs(int x)
{
dp[x][0]=1;
for(re int i=head[x];i;i=e[i].nxt)
if(!deep[e[i].v])
{
deep[e[i].v]=deep[x]+1;
dfs(e[i].v);
dp[x][0]+=dp[e[i].v][0];
for(re int k=1;k<=m;k++)
for(re int r=0;r<=k;r++)
dp[x][k]=(dp[x][k]+c[k][r]*dp[e[i].v][r])%mod;
}
}
void redfs(int x)
{
for(re int i=head[x];i;i=e[i].nxt)
if(deep[e[i].v]>deep[x])
{
memset(now,0,sizeof(now));
for(re int k=0;k<=m;k++)
{
now[k]=dp[x][k];
for(re int r=0;r<=k;r++)
now[k]=(now[k]-dp[e[i].v][r]*c[k][r]+mod)%mod;
}//先容斥,不能来自于e[i].v子树内部
for(re int k=0;k<=m;k++)
for(re int r=0;r<=k;r++)
dp[e[i].v][k]=(dp[e[i].v][k]+c[k][r]*now[r])%mod;//用容斥之后的答案来更新
redfs(e[i].v);
}
}
inline void add_edge(int x,int y)
{
e[++num].v=y;
e[num].nxt=head[x];
head[x]=num;
}
signed main()
{
n=read(),m=read();
int x,y;
for(re int i=1;i<n;i++)
x=read(),y=read(),add_edge(x,y),add_edge(y,x);
c[0][0]=1;
for(re int i=1;i<=m;i++) c[i][0]=c[i][i]=1;
for(re int i=1;i<=m;i++)
for(re int j=1;j<i;j++) c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
deep[1]=1;
dfs(1);
redfs(1);
for(re int i=1;i<=n;i++)
printf("%lld\n",(dp[i][m]%mod+mod)%mod);
return 0;
}
这是一个要T的复杂度,但是已经有了\(50\)分的好成绩
点分治据说可以做到\(O(nklogk)\),但是不会
我们继续组合数做吧
我们发现求\(x^k\)可以理解为把\(k\)个物品放到\(x\)个互不相同的盒子里,允许有盒子空着不放的方案数
于是我们可以写成\(x^k=\sum_{i=1}^kS(k,i)*C_{x}^i*i!\)
其中\(S(k,i)\)是第二类斯特林数,表示的是将\(k\)个球分到\(i\)个盒子里,这\(i\)个盒子没有差别,而且没有盒子是空的的方案数
\(C_{x}^i*i!\)其实就是排列数了,就相当于我们给\(i\)个盒子强行制造了差别
于是这个柿子可以理解为\(i\)枚举的是当前有几个盒子是有球的,之后通过加法原理合并了答案
其实我一开始觉得这里的\(\sum\)的上标应该写\(x\),好像也只有写\(x\)才满足组合的意义,之后发现自己非常naive
当\(k>x\)的时候显然是没有什么问题的了,因为\(\binom{x}{i}\)在\(x>i\)的时候取0,于是没有什么影响
当\(k<x\)的时候,如果上标取到比\(k\)大的数了,那么也会导致\(S(k,i)\)变成\(0\),于是写成\(k\)就可以了,在数值上没有什么影响
之后我们继续化柿子
\]
\]
\]
那么我们现在只需要求出\(\sum_{j=1}^n\binom{dis(t,j)}{i}\)就好了
我们都知道组合数有一个非常好的转移的方式就是\(\binom{n}{m}=\binom{n-1}{m-1}+\binom{n-1}{m}\)
于是
\]
现在是不是又可以用树形dp来转移了,因为到达\(t\)这个点还是要先到达\(t\)的儿子或者是父亲
于是我们设\(dp[x][k]=\sum_{j\in{x}}\binom{dis(x,i)}{k}\)
于是在\(up\)里的方程式就是
\]
\(down\)里我们还是要先容斥一下,方程和\(up\)类似
代码
#include<iostream>
#include<cstring>
#include<cstdio>
#define re register
#define maxn 50005
#define LL long long
const LL mod=10007;
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
struct node
{
int v,nxt;
}e[maxn<<1];
int head[maxn],deep[maxn];
LL dp[maxn][151],fac[151],s[151][151];
LL now[151];
int n,m,num;
inline void add_edge(int x,int y)
{
e[++num].v=y;
e[num].nxt=head[x];
head[x]=num;
}
void dfs(int x)
{
dp[x][0]=1;
for(re int i=head[x];i;i=e[i].nxt)
if(!deep[e[i].v])
{
deep[e[i].v]=deep[x]+1;
dfs(e[i].v);
dp[x][0]+=dp[e[i].v][0];
for(re int j=1;j<=m;j++)
dp[x][j]=(dp[x][j]+dp[e[i].v][j-1]+dp[e[i].v][j])%mod;
}
}
void redfs(int x)
{
for(re int i=head[x];i;i=e[i].nxt)
if(deep[e[i].v]>deep[x])
{
memset(now,0,sizeof(now));
now[0]=dp[x][0]-dp[e[i].v][0];
for(re int j=1;j<=m;j++)
now[j]=(dp[x][j]-dp[e[i].v][j]-dp[e[i].v][j-1]+mod)%mod;
//依旧是先容斥一遍
dp[e[i].v][0]+=now[0];
for(re int j=1;j<=m;j++)
dp[e[i].v][j]=(dp[e[i].v][j]+now[j]+now[j-1])%mod;
redfs(e[i].v);
}
}
int main()
{
n=read(),m=read();
int x,y;
for(re int i=1;i<n;i++)
x=read(),y=read(),add_edge(x,y),add_edge(y,x);
s[0][0]=1;
for(re int i=1;i<=m;i++)
s[i][1]=s[i][i]=1;
for(re int i=1;i<=m;i++)
for(re int j=1;j<i;j++)
s[i][j]=(s[i-1][j-1]+s[i-1][j]*j)%mod;//预处理第二类斯特林数
fac[0]=1;
for(re int i=1;i<=m;i++)
fac[i]=(fac[i-1]*i)%mod;//预处理阶乘
deep[1]=1;
dfs(1);
redfs(1);
for(re int i=1;i<=n;i++)
for(re int j=1;j<=m;j++)
dp[i][j]=(dp[i][j]%mod+mod)%mod;
for(re int i=1;i<=n;i++)
{
LL ans=0;
for(re int j=1;j<=m;j++)
ans=(ans+s[m][j]*fac[j]%mod*dp[i][j])%mod;
printf("%lld\n",ans);//统计答案
}
return 0;
}
【[国家集训队] Crash 的文明世界】的更多相关文章
- [国家集训队] Crash 的文明世界(第二类斯特林数)
题目 [国家集训队] Crash 的文明世界 前置 斯特林数\(\Longrightarrow\)斯特林数及反演总结 做法 \[\begin{aligned} ans_x&=\sum\limi ...
- [国家集训队] Crash的文明世界
Description 给定一棵 \(n\) 个点的树,对于每个点 \(i\) 求 \(S(i)=\sum\limits_{j=1}^n \operatorname{dist(i,j)}^k\) .\ ...
- [国家集训队] Crash 的文明世界
不错的树形$ DP$的题 可为什么我自带大常数啊$ cry$ 链接:here 题意:给定一棵$ n$个节点的树,边权为$ 1$,对于每个点$ x$求$ \sum\limits_{i=1}^n dist ...
- 洛谷P4827 [国家集训队] Crash 的文明世界 [斯特林数,组合数,DP]
传送门 思路 又见到这个\(k\)次方啦!按照套路,我们将它搞成斯特林数: \[ ans_x=\sum_{i=0}^k i!S(k,i)\sum_y {dis(x,y) \choose i} \] 前 ...
- P4827 [国家集训队] Crash 的文明世界
传送门:洛谷 题目大意:设$$S(i)=\sum_{j=1}^ndis(i,j)^k$$,求$S(1),S(2),\ldots,S(n)$. 数据范围:$n\leq 50000,k\leq 150$ ...
- 解题:国家集训队 Crash 的文明世界
题面 这种套着高次幂的统计问题一般都要用到第二类斯特林数和自然数幂的关系:$a^k=\sum\limits_{i=0}^{k}S_k^iC_a^i*i!$ 那么对于每个点$x$有: $ans_x=\s ...
- P4827 [国家集训队] Crash 的文明世界(第二类斯特林数+树形dp)
传送门 对于点\(u\),所求为\[\sum_{i=1}^ndis(i,u)^k\] 把后面那堆东西化成第二类斯特林数,有\[\sum_{i=1}^n\sum_{j=0}^kS(k,j)\times ...
- 国家集训队 Crash 的文明世界(第二类斯特林数+换根dp)
题意 题目链接:https://www.luogu.org/problem/P4827 给定一棵 \(n\) 个节点的树和一个常数 \(k\) ,对于树上的每一个节点 \(i\) ,求出 \( ...
- 洛谷 P4827 [国家集训队] Crash 的文明世界
题目描述 给你一棵 n 个点的树,对于树上的每个节点 i,求 \(\sum_{j=1}^ndis(i,j)^k\).其中 \(dis(i,j)\) 为两点在树上的距离. 输入格式 第一行两个整 ...
随机推荐
- PL/SQL之存储过程和触发器实例
1.Oracle存储过程实例 /*不带任何参数存储过程(输出系统日期)*/ CREATE OR REPLACE PROCEDURE output_date IS BEGIN DBMS_OUTPUT.P ...
- oAuth2.0及jwt介绍
oAuth2.0流程示意如下: 关于jwt介绍: 说明: 关于jwt简单说明一下,jwt即为json web token,是用来和服务端建立加密通信所使用的的一种“约定”,主要组成见上图即可.服务端一 ...
- MySQL常见常用的SQL优化
应尽量避免在where中使用!=或<>操作符.否则会进行全表查询 对于查询,避免全盘扫描,考虑在where或order by涉及到的列上建立索引 避免在where中进行null值判断,否则 ...
- Win10系统安装vmware workstation 12后没有桥接网卡怎么办
原文 vmware workstation是一款虚拟机计算机软件,可以同时运行不同的操作系统,然而有win10系统用户在安装vmware workstation 12之后,却发现网络连接里面没有桥接网 ...
- MySql的InnoDB存储引擎--索引
索引分类: 1.聚集索引:索引顺序与物理顺序一致. MySql 的 InnoDB 中,主键索引就是聚集索引.好处是,进行搜索的时候,因为索引和物理顺序一致,所以找数据的时候更快. 2.非聚集索引:索引 ...
- C# 后台处理图片的几种方式
第一种: 将上传图片直接保存到本地 var supportedTypes = new[] { "jpg", "jpeg", "png", & ...
- 从JSON中读取数据追加到HTML中
本文来自https://www.jianshu.com/p/04127d74d88c,并非本人原创,只是作为自己学习使用的资料,如有浏览者请点击地址自行到原作者页面浏览 有时候我们需要将json数据直 ...
- Django,ajax实现表格增删查改,Django内置分页功能。
1.工程目录 2.urls.py """Django_ajax URL Configuration The `urlpatterns` list routes URLs ...
- FCC的javascript初级算法题解答
FCC上的javascript基础算法题 前一阵子做的基础算法题,感觉做完后收获还蛮大的,现在将自己的做法总结出来,供大家参考讨论.基本上做到尽量简短有效,但有些算法还可以继续简化,比如第七题若采用正 ...
- python模拟登陆豆瓣——简单方法
学爬虫有一段时间了,前面没有总结又重装了系统,导致之前的代码和思考都没了..所以还是要及时整理总结备份.下面记录我模拟登陆豆瓣的方法,方法一登上了豆瓣,方法二重定向到了豆瓣中“我的喜欢”列表,获取了第 ...