题面

题解

删一条边、加一条边,相当于把一个子树折下来,然后嫁接在一个点上,

那么最优的情况肯定是接在根上,对吧,很好理解吧

那么这个拆下来的子树大小就不能超过n/2。

我们用son[]来表示每个点为根的子树大小,

如果一个点x可以改造后变成重心,那么要么它本来就是重心,要么它最多只有一个儿子y的son[y]大于n/2,并且y的子树大小可以通过改造变得<=n/2。

要改造一个儿子的子树,最优的方法就是减去里面最大的小于等于n/2的子子树,我们用dp1[]来表示一个点的子树里这样一个子子树的大小,然后暴力判就是了。

由于每个点的父亲也要考虑,所以用换根DP,一个dp2[]表示对于除了点的子树以外的部分可以剪掉的最大的部分。

转移方程:

  1. dp1[x] = max{son[y]*2 <= n ? son[y] : dp1[y]} //y 为 x 的儿子
  2. dp2[x] = (n - son[x])*2 <= n ? (n - son[x]) : max(dp2[father[x]],max{son[y]*2 <= n ? son[y] : dp1[y]}) //y 为 x 的父亲除了 x 以外的其他儿子
  3. /*
  4. 这里的dp2[x]后面不能暴力枚举y,只能处理出前缀最大值和后缀最大值来算
  5. */

注意,只能改造一次。

CODE

zxy A了! orz or2 orz

  1. #include<cstdio>
  2. #include<cstring>
  3. #include<iostream>
  4. //-----------F1
  5. using namespace std;
  6. #include<algorithm>
  7. #include<cmath>
  8. //-----------F2
  9. #include<vector>
  10. #include<stack>
  11. #include<queue>
  12. #include<map>
  13. #define MAXN 400005
  14. #define LL long long
  15. #define lowbit(x) (-(x) & (x))
  16. #define ENDL putchar('\n')
  17. //#pragma GCC optimize(2)
  18. //#pragma G++ optimize(3)
  19. //#define int LL
  20. char char_read_before = 1;
  21. inline int read() {
  22. int f = 1,x = 0;char s = char_read_before;
  23. while(s < '0' || s > '9') {if(s == '-') f = -1;s = getchar();}
  24. while(s >= '0' && s <= '9') {x = x * 10 - '0' + s;s = getchar();}
  25. char_read_before = s;return x * f;
  26. }
  27. LL zxy = 1000000007ll; // 用来膜的
  28. inline LL qkpow(LL a,LL b) {
  29. LL res = 1;
  30. while(b>0) {
  31. if(b & 1) res = res * a % zxy;
  32. a = a * a % zxy;
  33. b >>= 1;
  34. }
  35. return res;
  36. }
  37. int n,m,i,j,s,o,k;
  38. struct it{
  39. int v,w;
  40. it(){v=w=0;}
  41. it(int V,int W){v = V;w = W;}
  42. };
  43. vector<int> g[MAXN];
  44. vector<int> pre[MAXN],suf[MAXN];
  45. int son[MAXN];
  46. int dp1[MAXN],dp2[MAXN];
  47. int ans[MAXN];
  48. inline void dfs1(int x,int fa) {
  49. dp1[x] = 0;
  50. int sum = 0;
  51. pre[x].push_back(0);
  52. son[x] = 1;
  53. stack<int> sf;
  54. for(int i = 0;i < g[x].size();i ++) {
  55. if(g[x][i] != fa) {
  56. dfs1(g[x][i],x);
  57. son[x] += son[g[x][i]];
  58. dp1[x] = max(dp1[x],son[g[x][i]]*2 <= n ? son[g[x][i]] : dp1[g[x][i]]);
  59. sum = max(sum,son[g[x][i]]*2 <= n ? son[g[x][i]] : dp1[g[x][i]]);
  60. }
  61. pre[x].push_back(sum);
  62. }
  63. sum = 0;
  64. sf.push(0);
  65. for(int i = (int)g[x].size() - 1;i >= 0;i --) {
  66. if(g[x][i] != fa) {
  67. sum = max(sum,son[g[x][i]]*2 <= n ? son[g[x][i]] : dp1[g[x][i]]);
  68. }
  69. sf.push(sum);
  70. }
  71. while(!sf.empty()) suf[x].push_back(sf.top()),sf.pop();
  72. return ;
  73. }
  74. inline void dfs2(int x,int fa,int ll,int rr) {
  75. ans[x] = 1;
  76. int ct = 1;
  77. if(fa) {
  78. dp2[x] = (n - son[x])*2 <= n ? (n - son[x]) : max(dp2[fa],max(pre[fa][ll],suf[fa][rr]));
  79. if((n - son[x])*2 > n) {
  80. if((n - son[x] - dp2[x])*2 > n) ans[x] = 0;
  81. ct = 0;
  82. }
  83. }
  84. else dp2[x] = 0;
  85. for(int i = 0;i < g[x].size();i ++) {
  86. if(g[x][i] != fa) {
  87. dfs2(g[x][i],x,i,i+1);
  88. if(son[g[x][i]]*2 > n) {
  89. if(ct == 0 || (son[g[x][i]] - dp1[g[x][i]])*2 > n) ans[x] = 0;
  90. ct = 0;
  91. }
  92. }
  93. }
  94. return ;
  95. }
  96. signed main() {
  97. n = read();
  98. for(int i = 2;i <= n;i ++) {
  99. s = read();o = read();
  100. g[s].push_back(o);
  101. g[o].push_back(s);
  102. }
  103. dfs1(1,0);
  104. dfs2(1,0,0,0);
  105. for(int i = 1;i <= n;i ++) {
  106. printf("%d ",ans[i]);
  107. }
  108. return 0;
  109. }

Centroids (换根DP)的更多相关文章

  1. [BZOJ4379][POI2015]Modernizacja autostrady[树的直径+换根dp]

    题意 给定一棵 \(n\) 个节点的树,可以断掉一条边再连接任意两个点,询问新构成的树的直径的最小和最大值. \(n\leq 5\times 10^5\) . 分析 记断掉一条边之后两棵树的直径为 \ ...

  2. 2018.10.15 NOIP训练 水流成河(换根dp)

    传送门 换根dp入门题. 貌似李煜东的书上讲过? 不记得了. 先推出以1为根时的答案. 然后考虑向儿子转移. 我们记f[p]f[p]f[p]表示原树中以ppp为根的子树的答案. g[p]g[p]g[p ...

  3. 换根DP+树的直径【洛谷P3761】 [TJOI2017]城市

    P3761 [TJOI2017]城市 题目描述 从加里敦大学城市规划专业毕业的小明来到了一个地区城市规划局工作.这个地区一共有ri座城市,<-1条高速公路,保证了任意两运城市之间都可以通过高速公 ...

  4. 小奇的仓库:换根dp

    一道很好的换根dp题.考场上现场yy十分愉快 给定树,求每个点的到其它所有点的距离异或上m之后的值,n=100000,m<=16 只能线性复杂度求解,m又小得奇怪.或者带一个log像kx一样打一 ...

  5. 国家集训队 Crash 的文明世界(第二类斯特林数+换根dp)

    题意 ​ 题目链接:https://www.luogu.org/problem/P4827 ​ 给定一棵 \(n\) 个节点的树和一个常数 \(k\) ,对于树上的每一个节点 \(i\) ,求出 \( ...

  6. Acesrc and Travel(2019年杭电多校第八场06+HDU6662+换根dp)

    题目链接 传送门 题意 两个绝顶聪明的人在树上玩博弈,规则是轮流选择下一个要到达的点,每达到一个点时,先手和后手分别获得\(a_i,b_i\)(到达这个点时两个人都会获得)的权值,已经经过的点无法再次 ...

  7. bzoj 3566: [SHOI2014]概率充电器 数学期望+换根dp

    题意:给定一颗树,树上每个点通电概率为 $q[i]$%,每条边通电的概率为 $p[i]$%,求期望充入电的点的个数. 期望在任何时候都具有线性性,所以可以分别求每个点通电的概率(这种情况下期望=概率 ...

  8. codeforces1156D 0-1-Tree 换根dp

    题目传送门 题意: 给定一棵n个点的边权为0或1的树,一条合法的路径(x,y)(x≠y)满足,从x走到y,一旦经过边权为1的边,就不能再经过边权为0的边,求有多少边满足条件? 思路: 首先,这道题也可 ...

  9. [Bzoj3743][Coci2015] Kamp【换根Dp】

    Online Judge:Bzoj3743 Label:换根Dp,维护最长/次长链 题目描述 一颗树n个点,n-1条边,经过每条边都要花费一定的时间,任意两个点都是联通的. 有K个人(分布在K个不同的 ...

随机推荐

  1. c++ 辗转相除(动图)

    #include<iostream> #include<cstdio> #include<iomanip> #include<cstring> usin ...

  2. 2022年最强八股文《码出八股文-斩出offer线》

    宝剑锋从磨砺出,梅花香自苦寒来,大家好,我是小码哥 整理好的八股文终于完成了,希望看完对大家面试有所收获! 目录: 基础篇 javaOOP面试题 java集合/泛型面试题 java异常面试题 java ...

  3. SAP APO - 简介

    SAP高级计划和优化(SAP APO)是SAP SCM中的关键模块之一,在供应链流程中控制供应网络计划,备件计划,TP / VS和需求计划. 它可以帮助组织在动态环境中管理其供应链流程. SAP AP ...

  4. 【JS】两数之和

    给定一个整数数组 nums 和一个整数目标值 target,请你在该数组中找出 和为目标值 target 的那 两个 整数,并返回它们的数组下标. 你可以假设每种输入只会对应一个答案.但是,数组中同一 ...

  5. rhel安装程序

    Linux下软件分类     rpm软件包,包管理器 yum     deb软件包,包管理器 apt     源代码软件包            一般为".tar.gz".&quo ...

  6. mt19937 用法

    老是忘记怎么用,自己写一个用作备忘录吧. 首先需要的头文件: #include <random> 或者是 #include <bits/stdc++.h> //万能头 yyds ...

  7. 10道不得不会的JavaEE面试题

    10道不得不会的 JavaEE 面试题 我是 JavaPub,专注于面试.副业,技术人的成长记录. 以下是 JavaEE 面试题,相信大家都会有种及眼熟又陌生的感觉.看过可能在短暂的面试后又马上忘记了 ...

  8. Tapdata 在“疫”线:携手张家港市卫健委争分夺秒实时抗疫

      "抗疫两年以来最困难的时期,是漫长冬夜还是倒春寒?"--国家传染病医学中心主任张文宏 于3月14日凌晨   "等到疫情结束了,我一定要--",常怀这样的期许 ...

  9. 【一知半解】AQS

    什么是AbstractQueuedSynchronizer(AQS) 字面意思是抽象队列同步器,使用一个voliate修饰的int类型的同步状态,通过一个FIFO队列完成资源获取的排队工作,把每个参与 ...

  10. NOI-Online一日游(小炒肉心态爆炸记)

    提高组 早上起来赶来机房,逃掉数竟,耶!(光凯:?? T1 It takes me a long time想T1 努力了一个小时的成果 成功码出暴力代码,时间复杂度\(O(n^2)\) 㵘到15分 刑 ...