CF855G. Harry Vs Voldemort(边双,并查集,dp)

CF855G. Harry Vs Voldemort

Solution

考虑每增加一条边都会把路径上的边双都连成一个大边双,考虑合并xxxy=faxy = fa_xy=fax 这两个边双的贡献,分类讨论:

  1. 选取三个同边双内的点。
  2. 选取在同一个边双内选两个点,剩下一个在其他边双内。
  3. 选取三个来自不同边双的点。

第一个的贡献即为Aszx3A_{sz_x}^3Aszx3
第二个的贡献即为2Aszx2(n−szx)2A_{sz_x}^2(n - sz_x)2Aszx2(nszx)
第三个比较麻烦:

  • 首先我们去掉在一二两种中重复的部分,这部分可能是(x,y,z)(x,y,z)(x,y,z)(y,x,z)(y,x,z)(y,x,z)(z,x,y)(z,x,y)(z,x,y)(z,y,x)(z,y,x)(z,y,x)
  • 然后我们发现合并(x,y)(x,y)(x,y)之后,就可以从本来xxx子树中的点ppp开始走到yyy再回走到xxx的另一个子树中的点qqq,即(p,y,q)(p,y,q)(p,y,q)(q,y,p)(q,y,p)(q,y,p)都可以选择,这一部分是新多出来的,这部分相当于xxx的不同子树内的点对两两可达,于是我们再维护一个hx=∑yszy2h_x = \sum_{y}sz_y^2hx=yszy2即可快速统计贡献。(对于yyy这一部分的贡献同理)

具体统计方法见CodeCodeCode

并查集维护边双联通分量即可,时间复杂度O(nlgn)O(nlgn)O(nlgn)

Code

#include <bits/stdc++.h>using namespace std;template<typename T> inline bool upmin(T &x, T y) { return y < x ? x = y, 1 : 0; }
template<typename T> inline bool upmax(T &x, T y) { return x < y ? x = y, 1 : 0; }#define MP(A,B) make_pair(A,B)
#define PB(A) push_back(A)
#define SIZE(A) ((int)A.size())
#define LEN(A) ((int)A.length())
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define fi first
#define se second
#define int lltypedef long long ll;
typedef unsigned long long ull;
typedef long double lod;
typedef pair<int, int> PR;
typedef vector<int> VI; const lod eps = 1e-11;
const lod pi = acos(-1);
const int mods = 998244353;
const int oo = 1 << 30;
const ll loo = 1ll << 62;
const int MAXN = 600005;
const int INF = 0x3f3f3f3f;//1061109567
/*--------------------------------------------------------------------*/
inline int read() {int f = 1, x = 0; char c = getchar();while (c < '0' || c > '9') { if (c == '-') f = -1; c = getchar(); }while (c >= '0' && c <= '9') { x = (x << 3) + (x << 1) + (c ^ 48); c = getchar(); }return x * f;
}
vector<int> e[MAXN];
int fa[MAXN], f[MAXN], dep[MAXN];
ll sz[MAXN], num[MAXN], g[MAXN], h[MAXN], ans = 0, n;
int find(int x) { return f[x] == x ? f[x] : f[x] = find(f[x]); }
void dfs(int x, int father) {sz[x] = 0, fa[x] = father, dep[x] = dep[father] + 1;for (auto v : e[x]) if (v != father) dfs(v, x);for (auto v : e[x]) {if (v == father) continue;ans += sz[x] * sz[v] * 2;g[x] += g[v] + sz[v];h[x] += sz[v] * sz[v];sz[x] += sz[v]; }++ sz[x];h[x] += (n - sz[x]) * (n - sz[x]);ans += (n - sz[x]) * (sz[x] - 1) * 2;
}
void merge(int x, int y) {ans -= num[x] * (num[x] - 1) * (num[x] - 2); //part 1 xans -= num[y] * (num[y] - 1) * (num[y] - 2); //part 1 yans -= num[x] * (num[x] - 1) * (n - num[x]) * 2; //part 2 xans -= num[y] * (num[y] - 1) * (n - num[y]) * 2; //part 2 yans -= (sz[x] - num[x]) * num[x] * num[y] * 2 + (n - sz[x] - num[y]) * num[x] * num[y] * 2; //part 3.1ans += num[y] * ((sz[x] - num[x]) * (sz[x] - num[x]) - (h[x] - (n - sz[x]) * (n - sz[x]))); //part 3.2 xans += num[x] * ((n - sz[x] - num[y]) * (n - sz[x] - num[y]) - (h[y] - sz[x] * sz[x])); //part 3.2 yf[x] = y, num[y] += num[x], h[y] += h[x] - sz[x] * sz[x] - (n - sz[x]) * (n - sz[x]);ans += num[y] * (num[y] - 1) * (num[y] - 2); //part 1 new ans += num[y] * (num[y] - 1) * (n - num[y]) * 2; //part 2 new 
}
signed main() {n = read();for (int i = 1, u, v; i < n ; ++ i) u = read(), v = read(), e[u].PB(v), e[v].PB(u);for (int i = 1; i <= n ; ++ i) f[i] = i, num[i] = 1;dfs(1, 0);printf("%lld\n", ans);int Case = read();while (Case --) {int u = read(), v = read(), U = find(u), V = find(v);while (U != V) {if (dep[U] < dep[V]) swap(U, V);merge(U, find(fa[U]));U = find(U);}	printf("%lld\n", ans);}return 0;
}

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

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

相关文章

xBIM(eXtensible Building Information Modelling)可扩展的建筑信息模型

一、xBIM 简介BIM&#xff08;Building Information Modelling&#xff09;建筑信息模型&#xff0c;xBIM&#xff08;eXtensible Building Information Modelling&#xff09;可扩展的建筑信息模型。它是一个.NET 开源软件开发BIM工具包&#xff0c;支持BuildingSmart数据模型…

Codeforces Global Round 12 C1 C2. Errich-Tac-Toe 思维构造 好题

传送门 题意&#xff1a; 给了如下规则&#xff0c;上面三个只要出现一个情况就是非平局&#xff0c;现在给你个字符矩阵&#xff0c;让后其中XXX字符有KKK个(hardhardhard版本XXX和OOO一共KKK个)&#xff0c;每次操作可以将XXX变成OOO&#xff0c;OOO变成XXX&#xff0c;用不…

AGC019D - Shift and Flip(枚举)

AGC019D - Shift and Flip 很久之前WAWAWA的题&#xff0c;终于补掉了。。。这题细节是真的烦。 Solution 这题数据范围较小&#xff0c;于是我们枚举最终AAA与BBB的哪一个字符开始匹配&#xff0c;设这个位置为SSS。 然后考虑分顺时针/逆时针转到SSS两种情况讨论。 以逆时…

译 | .NET Core 3.0 Preview 6 已发布

原文&#xff1a;Richard Lander翻译&#xff1a;Edi Wang今天&#xff0c;我们宣布 .NET Core 3.0 Preview 6。它的更新包括编译程序集以改进启动、使用链接器和 EventPipe 改进优化应用程序的大小。我们还在 ARM64 上发布了针对 Alpine 的新 Docker 镜像。立即在 Windows、ma…

Codeforces Global Round 12 D. Rating Compression 思维 + 贪心

传送门 题意&#xff1a; 给一个长度为nnn的数组aaa&#xff0c;定义一个数组bbb&#xff0c;且bjminj<i<jk−1aib_jmin_{j<i<jk-1}a_ibj​minj<i<jk−1​ai​&#xff0c;比如j3j3j3的时候&#xff0c;a[1,3,4,5,2]a[1,3,4,5,2]a[1,3,4,5,2]&#xff0c;b3…

CF528C. Data Center Drama(欧拉回路,构造)

CF528C. Data Center Drama Solution 容易发现&#xff0c;加边后的图必然满足所有点的度为偶数&#xff0c;并且总的边数是偶数&#xff0c;这启发我们使用欧拉回路。 设欧拉回路为vk1et1vk2et2vk3...vk∣E∣1v_{k_1}e_{t_1}v_{k_2}e_{t_2}v_{k_3}...v_{k_{|E|1}}vk1​​et…

庆祝.NET Core悄然崛起:免费送50本优秀技术书籍,请笑纳!!

为了庆祝 .NET社区蓬勃发展&#xff0c;今天特地联合几位.NET大佬给大家带来50本高质量技术书籍。.NET Core自开源以来&#xff0c;依托微软强大的科技实力日趋完善&#xff0c;再也不是10年前那种封闭状态。强大的跨平台&#xff0c;强大的IDE&#xff0c;强大的语言特性&…

P2371 [国家集训队]墨墨的等式 同余最短路

传送门 题意&#xff1a; 思路&#xff1a; 一个同于最短路的板子题&#xff0c;初始的时候值为0&#xff0c;所以设dis[0]0dis[0]0dis[0]0&#xff0c;让后选择一个最小的a[i]a[i]a[i]作为basebasebase&#xff0c;跑一遍同余最短路就好啦。跑完dis[i]dis[i]dis[i]表示在模b…

CF997E. Good Subsegments(线段树,单调栈)

CF997E. Good Subsegments Description 给定一个序列&#xff0c;多次询问一个区间内的连续段个数。 n,q≤2105n,q \leq 2\times 10^{5}n,q≤2105 Solution 算得上是经典题了吧。 Part one 先考虑求全部的连续段个数。 相当于求区间内mx−mn1−len0mx - mn 1 - len 0m…

ASP.NET Core WebAPI帮助页--Swagger简单使用1.0

1、什么是Swagger&#xff1f;Swagger是一个规范且完整的框架&#xff0c;提供描述、生产、消费和可视化RESTful API&#xff0c;它是为了解决Web API生成有用文档和帮助页的问题。2、为啥选用swagger?1&#xff09;它具有交互式文档、客户端SDK生成和API可发现性等优点。2&am…

Codeforces Global Round 12 E. Capitalism 差分约束

传送门 题意&#xff1a; 思路&#xff1a; 一开始被题意迷惑了&#xff0c;没看出来差分约束&#xff0c;老菜鸡啦。首先看到ajai1a_ja_i1aj​ai​1可以把aia_iai​分成奇偶&#xff0c;让后这个图就变成一个二分图了。再考虑如何连边&#xff1a; (1) 对于b1b1b1的情况&…

ARC114E - Paper Cutting 2(组合数学,概率与期望)

ARC114E - Paper Cutting 2 Solution 考场上时间不够&#xff0c;没刚出来QAQ。 做法和官解本质相同&#xff0c;只是官解运用期望的线性性直接导出答案&#xff0c;而这里是对于所有方案统计贡献在除以方案数&#xff0c;从期望的定义上计算答案。可能稍显复杂。 Part one…

8个月打磨,一份送给程序员的「分布式系统」合集

这里是Z哥的个人公众号每周五早8点 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「75」篇原创敬上整理好的文章目录在文末&#xff0c;可直接拉到最后是的&#xff0c;这份礼物最佳受众是程序员。但是&#xff0c;如果你不是程序员&#xff0c;相信这些能使你…

HDU - 6071 Lazy Running 同余最短路 + 分层

传送门 题意&#xff1a; 给定四个点构成一个环&#xff0c;给出四个点之间的距离&#xff0c;让后从222号点出发&#xff0c;最终回到222号点&#xff0c;求经过的距离>k>k>k的最小距离。 思路&#xff1a; 由于从222开始&#xff0c;最终在222结束&#xff0c;所以…

CF1398F

CF1398F Solution 我又来贡献暴力做法了。。。听说两只log不可能过1e6? 有一个显然的想法是&#xff1a; 我们先预处理一个aia_iai​表示第iii个位置的最长后缀长度&#xff0c;满足该后缀中不同时存在0和1。 假设我们要求xixixi的答案&#xff0c;那么一个位置jjj可以作…

【翻译】无需安装Python,就可以在.NET里调用Python库

原文地址&#xff1a;https://henon.wordpress.com/2019/06/05/using-python-libraries-in-net-without-a-python-installation/pythonnet这个屌爆的项目的出现&#xff0c;使得我们可以用一种新的方式&#xff0c;让C#可以和Python之间进行互操作。但是它的设置和部署可能有点…

Codeforces Round #706 (Div. 2) E. Garden of the Sun 思维构造

传送门 题意&#xff1a; 给你一个nmnmnm的矩阵&#xff0c;其中包含字符′.′.′.′和′X′X′X′&#xff0c;你可以将任何′.′.′.′改成′X′X′X′&#xff0c;现在问你能否通过修改一些′.′.′.′来使′X′X′X′联通且不存在环。保证原本的′X′X′X′没有任何两个相…

CF1067E Random Forest Rank(树形dp,概率与期望,线性代数)

CF1067E Random Forest Rank Solution 考虑树的邻接矩阵的秩的意义&#xff0c;不难发现相当于每个点找一个“代表”&#xff0c;对于一个点xxx&#xff0c;它的“代表”为一个与xxx相邻的结点&#xff0c;要求保证所有点的“代表”不重复&#xff0c;求能找到“代表”的点的…

Ocelot(六)- 架构图

简介Ocelot是一个用.NET Core实现并且开源的API网关&#xff0c;它功能强大&#xff0c;包括了&#xff1a;路由、请求聚合、服务发现、认证、鉴权、限流熔断、并内置了负载均衡器与Service Fabric、k8s 集成。这些功能只都只需要简单的配置即可完成。架构图Ocelot的目标是使用…

P5170 【模板】类欧几里得算法(类欧)

P5170 【模板】类欧几里得算法 Description 要求在O(lgn)O(lgn)O(lgn)的时间内求出&#xff1a; ∑i0n⌊aibc⌋\sum_{i 0}^n{\lfloor\frac{aib}{c}\rfloor}∑i0n​⌊caib​⌋ ∑i0ni⌊aibc⌋\sum_{i 0}^n{i\lfloor\frac{aib}{c}\rfloor }∑i0n​i⌊caib​⌋ ∑i0n⌊aibc⌋2…