JZOJ5444. 【NOIP2017提高A组冲刺11.2】救赎
来源:互联网 发布:ssh毕业设计2017源码 编辑:程序博客网 时间:2024/06/05 16:46
Description
“是的。”我回答,“我不会忘记你。在森林里我会一点点记起往日的世界。要记起的大概很多很多:各种人、各种场所、各种光、各种歌曲……”
——村上春树《世界尽头与冷酷仙境》
在没有心存在的世界尽头,音乐能够使小镇居民消散的心重新聚拢成形。作为镇子里唯一一个还残留着些许音乐记忆的人,我逐渐记起了往昔点滴……
记忆中有一棵无根树,有n个节点。
对于一棵有根树的每一个非叶子节点,我们都等概率选中其一个儿子节点作为偏好儿子。对于一条从父亲指向儿子的树边(u,v),如果v是u的偏好儿子,则称这条边为重边,否则为轻边。
我们定义一棵有根树的权值为其每一个节点到根路径上的轻边条数的和的期望值。
请对无根树每一个节点输出其为根的有根树的权值。答案模998244353。
Input
文件第一行是一个正整数n。
接下来n-1行,每行两个正整数(x,y)表示一条树边。
Output
输出文件共n行,每一行一个整数表示答案。
Sample Input
5
1 2
1 3
3 4
3 5
Sample Output
3
1
665496238
499122178
499122178
Data Constraint
对于10%的数据,保证n<=10。
对于30%的数据,保证n<=2000。
对于100%的数据,保证n<=10^5。
题解
求期望,
就是直接上套路。
设
求一个点为根的期望的复杂度是O(n)的,
如果要求每个点,那么总的复杂度就是O(
如果每次只改变根节点与它儿子的边,就是将根节点的某个儿子旋转到根节点上,这样修改只与两个点有关,其他都是不变的,这样复杂度就是O(n)的。
code
#include <iostream>#include <cstdio>#include <cmath>#include <cstring>#include <queue>#define N 200003#define ll long long#define mo 998244353using namespace std;char ch;void read(int& n){ n=0; ch=getchar(); while(ch<'0'||ch>'9')ch=getchar(); while('0'<=ch && ch<='9')n=(n<<1)+(n<<3)+ch-'0',ch=getchar();}void write(ll x){ if(x>9)write(x/10); putchar(x%10+48);}int nxt[N*2],b[N],to[N*2],tot,n,m,x,y;ll f[N],size[N],son[N],ans[N];int z[N],l,r,fa[N],top;void ins(int x,int y){ nxt[++tot]=b[x]; to[tot]=y; b[x]=tot; son[y]++;}ll ksm(ll x,int y){ ll s=1; while(y) { if(y&1)s=s*x%mo; x=x*x%mo; y>>=1; } return s;}void dfs(int x,int fa){ ll ss=ksm(son[x]-1,mo-2); for(int i=b[x];i;i=nxt[i]) if(to[i]!=fa) { dfs(to[i],x); f[x]=(f[x]+ss*f[to[i]]%mo+(son[x]-2)*ss%mo*(f[to[i]]+size[to[i]])%mo)%mo; size[x]+=size[to[i]]; } size[x]++;}void exchange(int x/*root*/,int y/*new root*/){ son[x]--;son[y]++; size[x]=size[y]=1; f[x]=f[y]=0; //reset x ll ss=ksm(son[x]-1,mo-2); for(int i=b[x];i;i=nxt[i]) if(to[i]!=y) { f[x]=(f[x]+ss*f[to[i]]%mo+(son[x]-2)*ss%mo*(f[to[i]]+size[to[i]])%mo)%mo; size[x]+=size[to[i]]; } //reset y ss=ksm(son[y]-1,mo-2); for(int i=b[y];i;i=nxt[i]) { f[y]=(f[y]+ss*f[to[i]]%mo+(son[y]-2)*ss%mo*(f[to[i]]+size[to[i]])%mo)%mo; size[y]+=size[to[i]]; }}void work(int x,int fa){ ans[x]=f[x]; for(int i=b[x];i;i=nxt[i]) if(to[i]!=fa) { exchange(x,to[i]); work(to[i],x); exchange(to[i],x); }}int main(){ freopen("redemption.in","r",stdin); freopen("redemption.out","w",stdout); read(n); for(int i=1;i<n;i++) read(x),read(y),ins(x,y),ins(y,x); memset(size,0,sizeof(size)); memset(f,0,sizeof(f)); son[1]++; dfs(1,0); work(1,0); for(int i=1;i<=n;i++) { write(ans[i]); putchar('\n'); } return 0;}
- Jzoj5444【NOIP2017提高A组冲刺11.2】救赎
- JZOJ5444. 【NOIP2017提高A组冲刺11.2】救赎
- 【NOIP2017提高A组冲刺11.2】救赎(数学期望)
- JZOJ 5444. 【NOIP2017提高A组冲刺11.2】救赎
- JZOJ 5444. 【NOIP2017提高A组冲刺11.2】救赎
- 【NOIP2017提高A组冲刺11.2】总结
- 【NOIP2017提高A组冲刺11.2】失格
- 【JZOJ 5445】【NOIP2017提高A组冲刺11.2】失格
- JZOJ5445. 【NOIP2017提高A组冲刺11.2】失格
- Jzoj5443【NOIP2017提高A组冲刺11.2】字典序
- JZOJ5443. 【NOIP2017提高A组冲刺11.2】字典序
- JZOJ 5443. 【NOIP2017提高A组冲刺11.2】字典序
- jzoj【NOIP2017提高A组冲刺11.2】失格
- [JZOJ5445]【NOIP2017提高A组冲刺11.2】失格
- JZOJ 5445. 【NOIP2017提高A组冲刺11.2】失格
- 【NOIP2017提高A组冲刺11.2】失格(乱搞)
- Jzoj5445【NOIP2017提高A组冲刺11.2】失格
- JZOJ 5445. 【NOIP2017提高A组冲刺11.2】失格
- Unity调试Android安装包
- Java线程池快速学习教程
- 微信APP下载页面
- 【TensorFlow实战计算机视觉】 之 猫狗大战
- JMeter简单压力测试
- JZOJ5444. 【NOIP2017提高A组冲刺11.2】救赎
- git过滤
- 练习使用消息队列
- BFS(广度优先搜索)
- 解决Nginx无法站外访问
- spring定时器--时间设置规则
- 常用的数组API
- Spark性能优化:数据倾斜调优
- python:读取视频,处理后,实时计算帧数fps