博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
Bzoj2152/洛谷P2634 聪聪可可(点分治)
阅读量:4946 次
发布时间:2019-06-11

本文共 2076 字,大约阅读时间需要 6 分钟。

题面

题解

点分治套路走一波,考虑\(calc\)函数怎么写,存一下每条路径在\(\%3\)意义下的路径总数,假设为\(tot[i]\)\(\equiv i(mod\ 3)\),这时当前的贡献就是\(tot[0]^2+2\times tot[1]\times tot[2]\)

#include 
#include
#include
#include
using std::min; using std::max;using std::swap; using std::sort;using std::__gcd;typedef long long ll;template
void read(T &x) { int flag = 1; x = 0; char ch = getchar(); while(ch < '0' || ch > '9') { if(ch == '-') flag = -flag; ch = getchar(); } while(ch >= '0' && ch <= '9') x = x * 10 + ch - '0', ch = getchar(); x *= flag;}const int N = 2e4 + 10, Inf = 1e9 + 7;int n, m, k, d[N], Size, ans, tot[4];int cnt, from[N], to[N << 1], nxt[N << 1], dis[N << 1];int p, siz[N], tmp; bool vis[N];inline void addEdge(int u, int v, int w) { to[++cnt] = v, nxt[cnt] = from[u], dis[cnt] = w, from[u] = cnt;}void getrt(int u, int f) { siz[u] = 1; int max_part = 0; for(int i = from[u]; i; i = nxt[i]) { int v = to[i]; if(v == f || vis[v]) continue; getrt(v, u), siz[u] += siz[v]; max_part = max(max_part, siz[v]); } max_part = max(max_part, Size - siz[u]); if(max_part < tmp) tmp = max_part, p = u;}void getpoi(int x, int y, int f) { ++tot[y]; for(int i = from[x]; i; i = nxt[i]) { int v = to[i]; if(v == f || vis[v]) continue; getpoi(v, (y + dis[i]) % 3, x); }}void calc (int x, int y, int dd) { memset(tot, 0, sizeof tot); getpoi(x, y % 3, 0); ans += dd * (tot[0] * tot[0] + 2 * tot[1] * tot[2]);}void doit(int x) { tmp = Inf, getrt(x, 0), vis[p] = true; calc(p, 0, 1); for(int i = from[p]; i; i = nxt[i]) { int v = to[i]; if(vis[v]) continue; calc(v, dis[i], -1); Size = siz[v], doit(v); }}int main () { read(n), Size = n; for(int i = 1, u, v, w; i < n; ++i) { read(u), read(v), read(w); addEdge(u, v, w), addEdge(v, u, w); } doit(1); int m = n * n, gg = __gcd(m, ans); printf("%d/%d\n", ans / gg, m / gg); return 0;}

转载于:https://www.cnblogs.com/water-mi/p/10304382.html

你可能感兴趣的文章
foreach遍历数组、数组的转置与方阵的迹
查看>>
Still unable to dial persistent://blog.csdn.net:80 after 3 attempts
查看>>
HTML超文本标记语言(九)——表单输入类型
查看>>
基于busybox制作mini2440根文件系统及使用nfs挂载
查看>>
信道容量及信道编码原理学习
查看>>
浅谈独立特征(independent features)、潜在特征(underlying features)提取、以及它们在网络安全中的应用...
查看>>
从随机过程的熵率和马尔科夫稳态过程引出的一些思考 - 人生逃不过一场马尔科夫稳态...
查看>>
《A First Course in Abstract Algebra with Applications》-chaper1-数论-关于素数
查看>>
算法笔记_145:拓扑排序的应用(Java)
查看>>
JS获取农历日期
查看>>
PHP中的HTTP协议
查看>>
Java WebService入门实例
查看>>
css样式之补充
查看>>
结构与联合
查看>>
BUPT复试专题—众数(2014)
查看>>
20145316 《信息安全系统设计基础》第十四周学习总结
查看>>
Liferay7 BPM门户开发之18: 理解ServiceContext
查看>>
Intel Galileo development documentation
查看>>
EV: Workaround to Allow Only One Instance or Window of outlook
查看>>
数据校验,
查看>>