jzoj5365-[GDOI2018模拟9.14]通信【线段树合并】

正题


题目大意

nnn个节点的一棵树,随机选择一个区间,求这个区间的点所构成的虚树的期望权值和。


解题思路

考虑每一条边的贡献,定义一边的点为黑点,一边的为白点,显然包含黑白的区间都会产生贡献。考虑减去没有贡献的,也就是对于连续的一段长度为lll的颜色相同的区间会产生贡献n(n+1)2\frac{n(n+1)}{2}2n(n+1)

用线段树维护然后合并即可。

时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=1e5+10,M=N*20,XJQ=1e9+7;
ll f(ll n)
{return n*(n-1)/2;}
struct Seg_Tree{ll cnt,ls[M],rs[M],lb[M],rb[M],lw[M],rw[M],ans[M];void PushUp(ll x,ll L,ll R){ll mid=(L+R)>>1;if(!ls[x])ls[x]=++cnt,lb[cnt]=rb[cnt]=(mid-L+1),ans[cnt]=f(mid-L+1);if(!rs[x])rs[x]=++cnt,lb[cnt]=rb[cnt]=(R-mid),ans[cnt]=f(R-mid);ans[x]=ans[ls[x]]+ans[rs[x]];ans[x]=ans[x]-f(rb[ls[x]])-f(lb[rs[x]])+f(rb[ls[x]]+lb[rs[x]]);ans[x]=ans[x]-f(rw[ls[x]])-f(lw[rs[x]])+f(rw[ls[x]]+lw[rs[x]]);lb[x]=(lb[ls[x]]==mid-L+1)?(lb[ls[x]]+lb[rs[x]]):lb[ls[x]];rb[x]=(rb[rs[x]]==R-mid)?(rb[rs[x]]+rb[ls[x]]):rb[rs[x]];lw[x]=(lw[ls[x]]==mid-L+1)?(lw[ls[x]]+lw[rs[x]]):lw[ls[x]];rw[x]=(rw[rs[x]]==R-mid)?(rw[rs[x]]+rw[ls[x]]):rw[rs[x]];}void Change(ll &x,ll l,ll r,ll pos){if(!x)x=++cnt;if(l==r){lw[x]=rw[x]=1;return;}ll mid=(l+r)>>1;if(pos<=mid)Change(ls[x],l,mid,pos);else Change(rs[x],mid+1,r,pos);PushUp(x,l,r);}ll Merge(ll x,ll y,ll l,ll r){if(!x||!y)return x+y;if(lb[x]==r-l+1)return y;if(lb[y]==r-l+1)return x;ll mid=(l+r)>>1;ls[x]=Merge(ls[x],ls[y],l,mid);rs[x]=Merge(rs[x],rs[y],mid+1,r);PushUp(x,l,r);return x;}
}T;
struct node{ll to,next;
}a[N*2];
ll n,tot,ls[N],rt[N],ans;
void addl(ll x,ll y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;return;
}
ll power(ll x,ll b){ll ans=1;x%=XJQ;while(b){if(b&1)ans=ans*x%XJQ;x=x*x%XJQ;b>>=1;}return ans;
}
void dfs(ll x,ll fa){T.Change(rt[x],1,n,x);for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(y==fa)continue;dfs(y,x);T.Merge(rt[x],rt[y],1,n);}ans+=f(n)-T.ans[rt[x]];return;
}
int main()
{freopen("communicate.in","r",stdin);freopen("communicate.out","w",stdout);scanf("%lld",&n);for(ll i=1;i<n;i++){ll x,y;scanf("%lld%lld",&x,&y);addl(x,y);addl(y,x);}dfs(1,1);ans=ans*2%XJQ;printf("%lld",ans*power(n*(n+1)/2,XJQ-2)%XJQ);
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/320534.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

HDU5129 - Yong Zheng's Death

HDU5129 - Yong Zhengs Death 做法&#xff1a;看题之后&#xff0c;就想到要建好trie树&#xff0c;然后输出非根节点数的平方&#xff0c;这样显然错了。因为会出现重算的情况&#xff0c;即&#xff1a;你有 ab,cde ,abc,de, abcd, e,这几个串那么 [ab, cde] [abc,de] [abcd…

牛客网【每日一题】3月27日 数学考试

题号 NC15553 名称 数学考试 来源 2018年长沙理工大学第十三届程序设计竞赛 链接&#xff1a;https://ac.nowcoder.com/acm/problem/15553 来源&#xff1a;牛客网 试题传送门 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;…

【SPFA】重建道路(jzoj 1212)

重建道路 jzoj 1212 题目大意&#xff1a; 有一个图&#xff0c;其中的一些路坏了&#xff0c;为了从a走到b&#xff0c;问最少修复多长的路 输入样例 3 2 1 2 1 2 3 2 1 1 2 1 3输出样例 1数据范围 2⩽N⩽1002\leqslant N\leqslant 1002⩽N⩽100 N−1⩽D⩽M⩽N∗(N−1)/…

ASP.NET CORE 根据环境变量支持多个 appsettings.json

0.背景在开发项目的过程当中&#xff0c;生产环境与调试环境的配置肯定是不一样的。拿个最简单的例子来说&#xff0c;比如连接字符串这种东西&#xff0c;调试环境肯定是不能连接生产数据库的。在之前的话&#xff0c;这种情况只能说是你 COPY 两个同名的配置文件来进行处理。…

GYM101933I - Intergalactic Bidding

GYM101933I - Intergalactic Bidding 题解&#xff1a;不考虑首先显然是个背包&#xff0c;一开始直接用set模拟&#xff0c;然后map存方案&#xff0c;这样会mle。发现物品的体积有的特殊性 only one participant was allowed to make a bid at a time, each participant was …

【结论】【dfs】费解的开关(joyoi-tyvj 1266)

费解的开关 joyoi-tyvj 1266 题目大意&#xff1a; 有5*5的一个图&#xff0c;每个点的数值是1或0&#xff0c;如果将一个点的数值取反&#xff0c;那这个点上下左右的点都会取反&#xff0c;现在问你将所有点都变为1最少要多少步&#xff0c;如果步数大于6或无法全变成1的话…

P3586-[POI2015]LOG【线段树】

正题 题目链接:https://www.luogu.com.cn/problem/P3586 题目大意 nnn个数&#xff0c;有操作 修改一个数假如每次选出ccc个正数让它们减去111&#xff0c;求能否进行操作sss次。 解题思路 如果有cntcntcnt个数大于sss那么有解当且仅当满足cnt≥ccnt\geq ccnt≥c或者小于sss…

小G的项链(Manacher)

我看网上也没有写这个题的&#xff0c;顺便写一下&#xff08;可能是大佬都觉得太简单了 &#xff09; 链接&#xff1a;牛客网 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 32768K&#xff0c;其他语言65536K 64bit IO Format: %lld 题目描述…

基于Win10极简SonarQube C#代码质量分析

博客有些好些时间未更新了&#xff0c;这几个月的时间里&#xff0c;离开了实习的公司、大学毕了业、来了新公司、转了户口&#xff0c;有点忙&#xff0c;最近总算稍微闲下来了&#xff0c;打算重新拾起博客&#xff0c;坚持写下去。言归正转&#xff0c;什么是SonarQube ?So…

2018年牛客多校算法寒假训练营练习比赛(第一场)C. 六子冲

2018年牛客多校算法寒假训练营练习比赛&#xff08;第一场&#xff09;C. 六子冲 题目链接 做法&#xff1a;模拟即可 #include <bits/stdc.h> #define P pair<int,int> #define fir first #define sec second #define rep(i,a,b) for(int ia;i<b;i) using name…

YbtOJ#20078-[NOIP2020模拟赛B组Day7]路径之和【分治,Flody】

正题 题目链接:http://noip.ybtoj.com.cn/contest/108/problem/4 题目大意 nnn个点&#xff0c;每次只封闭一个点&#xff0c;求剩下点对两两之间的最短路和。 解题思路 FlodyFlodyFlody中如果枚举kkk时不管某一个点&#xff0c;那么就可以求不经过那个点的路径。考虑分治求解…

大吉大利【牛客网】(牛客练习赛60)

传送 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 262144K&#xff0c;其他语言524288K 64bit IO Format:%lld 题目描述 输入描述: 第一行一个整数n. 第二行n个整数ai. 输出描述: 一个整数表示上述求和式的答案. 示例1 输入 5 1 2 3 4 5输出…

【差分】Tallest Cow(poj 3263/luogu 2879)

Tallest Cow poj 3263 luogu 2879 题目大意&#xff1a; 现在有n头牛&#xff0c;两头牛如果要相互看到&#xff0c;那他们之间的牛必须比他们两低&#xff0c;现在给出n&#xff0c;最高牛的位置和高度&#xff0c;和m对关系&#xff0c;要你求每头牛最高是多少 输入样例…

【项目管理】git和码云的使用

缘起说了那么多关于git和码云相关的事&#xff0c;一直都没给大伙讲解这个码云究竟是个啥玩意儿。今天就给大伙说说如何通过git和码云搭建属于自己的代码库。码云码云(GitOSC)是开源中国社区团队推出的基于Git的快速的、免费的、稳定的在线代码托管平台,不限制私有库和公有库数…

牛客网暑期ACM多校训练营(第十场)F.Rikka with Line Graph

牛客网暑期ACM多校训练营&#xff08;第十场&#xff09;F.Rikka with Line Graph 做法&#xff1a;\(G\) 中的对应原图两条边(a,b) (c,d)的最短路为:\[ w[a][b] w[c][d] 2* min(dis[a][c], dis[a][d], dis[b][c], dis[b][d])\] 其中\(dis[i][j]\)表示原图G中i 到 j 的最短路…

AT2165-[AGC006D]MedianPyramidHard【二分,贪心】

正题 题目链接:https://www.luogu.com.cn/problem/AT2165 题目大意 nnn行的一个金字塔&#xff0c;第iii层有2i−12\times i-12i−1个格子&#xff0c;每个格子是它左右下的中位数。给出最下面的2n−12n-12n−1个数&#xff0c;求第一行的数。 解题思路 我们先二分一个答案&a…

三角形周长和【牛客网】牛客网练习赛60

题目传送 时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 262144K&#xff0c;其他语言524288K 64bit IO Format:%lld 题目描述 输入描述: 输入描述 第一行一个整数表示n. 接下来n行每行两个整数x,y表示一个点. 输出描述: 输出一个整数表示周长…

【二分】递增(luogu 3902)

递增 luogu 3902 题目大意&#xff1a; 给出一个序列&#xff0c;问最少修改多少个数字可以使他单调递增 输入样例 3 1 3 2输出样例 1 数据范围 对于50% 的数据&#xff0c;N≤103N \le 10^3N≤103 对于100% 的数据&#xff0c;1≤N≤105,1≤Ai≤1091 \le N \le 10^5 , …

AspNetCore 基于AOP实现Polly的使用

前言 说起AOP&#xff0c;其实我们在做MVC/API 的时候应该没少接触&#xff0c;比如说各种的Fitter 就是典型的AOP了。本来在使用Polly的时候我最初的打算是使用过滤器来实现的&#xff0c;后来发现实现起来相当的困难&#xff0c;利用NetCore的中间以及过滤器去实现一个AOP的…

HDU5794 - A Simple Chess

HDU5794 - A Simple Chess 做法&#xff1a;首先的想法就是用总方案数减去&#xff0c;经过过障碍的方案数A。第一个思路就是容斥&#xff0c;但是显然不符合数据规模。另一个思路就是将障碍物从左上到右下排序&#xff0c;dp[i] 表示不经过前i-1个障碍&#xff0c;到达第i个障…