虚树(Bzoj3611: [Heoi2014]大工程)
题面
虚树
把跟询问有关的点拿出来建树,为了方便树\(DP\)
在\(LCA\)处要合并答案,那么把这些点的\(LCA\)也拿出来
做法:把点按\(dfs\)序排列,然后求出相邻两个点的\(LCA\),把这些点建一个虚树,维护一个栈就好了
Sol
虚树+树\(DP\)
# include <bits/stdc++.h>
# define IL inline
# define RG register
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
IL int Input(){
RG int x = 0, z = 1; RG char c = getchar();
for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
return x * z;
}
const int maxn(1e6 + 5);
const int inf(1e9);
int n, first1[maxn], cnt, first2[maxn], id[maxn], p[maxn];
ll g[maxn], f[maxn], num[maxn], sum, mx, mn, tot;
int dfn[maxn], size[maxn], son[maxn], fa[maxn];
int s[maxn], top[maxn], idx, deep[maxn];
struct Edge{
int to, next;
} e2[maxn << 1], e1[maxn << 1];
IL void Add1(RG int u, RG int v){
e1[cnt] = (Edge){v, first1[u]}, first1[u] = cnt++;
}
IL void Add2(RG int u, RG int v){
e2[cnt] = (Edge){v, first2[u]}, first2[u] = cnt++;
}
IL void Dfs1(RG int u){
size[u] = 1;
for(RG int e = first1[u]; e != -1; e = e1[e].next){
RG int v = e1[e].to;
if(!size[v]){
deep[v] = deep[u] + 1;
Dfs1(v), size[u] += size[v], fa[v] = u;
if(size[v] > size[son[u]]) son[u] = v;
}
}
}
IL void Dfs2(RG int u, RG int tp){
top[u] = tp, dfn[u] = ++idx;
if(son[u]) Dfs2(son[u], tp);
for(RG int e = first1[u]; e != -1; e = e1[e].next)
if(!dfn[e1[e].to]) Dfs2(e1[e].to, e1[e].to);
}
IL int LCA(RG int u, RG int v){
while(top[u] ^ top[v])
deep[top[u]] > deep[top[v]] ? u = fa[top[u]] : v = fa[top[v]];
return deep[u] > deep[v] ? v : u;
}
IL int Dis(RG int u, RG int v){
RG int lca = LCA(u, v);
return deep[u] + deep[v] - 2 * deep[lca];
}
IL int Cmp(RG int u, RG int v){
return dfn[u] < dfn[v];
}
IL void DP(RG int u){
g[u] = inf, f[u] = -inf;
if(num[u]) g[u] = f[u] = 0;
for(RG int e = first2[u]; e != -1; e = e2[e].next){
RG int v = e2[e].to, w = Dis(u, v);
DP(v), num[u] += num[v];
sum += (tot - num[v]) * num[v] * w;
mn = min(mn, g[u] + w + g[v]);
mx = max(mx, f[u] + w + f[v]);
g[u] = min(g[u], g[v] + w);
f[u] = max(f[u], f[v] + w);
}
}
int main(){
n = Input();
for(RG int i = 1; i <= n; ++i) first1[i] = first2[i] = -1;
for(RG int i = 1; i < n; ++i){
RG int u = Input(), v = Input();
Add1(u, v), Add1(v, u);
}
Dfs1(1), Dfs2(1, 1);
for(RG int q = Input(); q; --q){
RG int k = Input(); cnt = 0, tot = k;
for(RG int i = 1; i <= k; ++i) p[i] = Input(), num[p[i]] = 1;
sort(p + 1, p + k + 1, Cmp);
for(RG int i = 1, t = k; i < t; ++i) p[++k] = LCA(p[i], p[i + 1]);
sort(p + 1, p + k + 1, Cmp), k = unique(p + 1, p + k + 1) - p - 1;
for(RG int i = 1; i <= k; ++i) first2[p[i]] = -1;
RG int t = 0;
for(RG int i = 1; i <= k; ++i){
while(t && dfn[p[i]] >= dfn[s[t]] + size[s[t]]) --t;
if(t) Add2(s[t], p[i]);
s[++t] = p[i];
}
mx = sum = 0, mn = inf;
DP(p[1]);
printf("%lld %lld %lld\n", sum, mn, mx);
for(RG int i = 1; i <= k; ++i) num[p[i]] = 0;
}
return 0;
}
虚树(Bzoj3611: [Heoi2014]大工程)的更多相关文章
- [BZOJ3611][Heoi2014]大工程
[BZOJ3611][Heoi2014]大工程 试题描述 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶点上. 在 ...
- [Bzoj3611][Heoi2014]大工程(虚树)
3611: [Heoi2014]大工程 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 2000 Solved: 837[Submit][Status ...
- [BZOJ3611][Heoi2014]大工程(虚树上DP)
3611: [Heoi2014]大工程 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 2464 Solved: 1104[Submit][Statu ...
- BZOJ2286 [Sdoi2011]消耗战 和 BZOJ3611 [Heoi2014]大工程
2286: [Sdoi2011]消耗战 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 6371 Solved: 2496[Submit][Statu ...
- BZOJ3611:[HEOI2014]大工程(树形DP,虚树)
Description 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶点上. 在 2 个国家 a,b 之间建一条新通 ...
- BZOJ3611 [Heoi2014]大工程 【虚树】
题目 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶点上. 在 2 个国家 a,b 之间建一条新通道需要的代价为树上 a ...
- [BZOJ3611] [Heoi2014]大工程(DP + 虚树)
传送门 $dp[i][0]$表示节点i到子树中的所有点的距离之和 $dp[i][1]$表示节点i到子树中最近距离的点的距离 $dp[i][2]$表示节点i到子树中最远距离的点的距离 建好虚树后dp即可 ...
- BZOJ3611 HEOI2014大工程
先建虚树,然后统计答案. 对于这个两点间最大值和最小值的操作我参考了hzwer的代码. 建虚树时注意判自环 By:大奕哥 #include<bits/stdc++.h> using nam ...
- bzoj 3611 [Heoi2014]大工程(虚树+DP)
3611: [Heoi2014]大工程 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 408 Solved: 190[Submit][Status] ...
随机推荐
- Common xaml controls(补交作业)
Common xaml controls 常见的xaml控件: 先上一段代码,把他们基本都实现出来: <Grid Name="MyGrid"> <Button N ...
- 浅谈Android选项卡(四)
前面几篇介绍的选项的用法,基本上使用TabActivity.ViewPager.已经基本上满足开发需求了.但是这里再介绍一种小技巧,在有的时候,感觉使用前面的ViewPager和Fragment时候, ...
- 代码 | 自适应大邻域搜索系列之(4) - Solution定义和管理的代码实现解析
前言 上一篇讲解了destroy和repair方法的具体实现代码,好多读者都在喊酸爽和得劲儿--今天这篇就讲点简单的,关于solution的定义和管理的代码实现,让大家回回神吧--哈哈. 01 总体概 ...
- POJ 2209
#include<iostream> #include<stdio.h> #include<algorithm> #include<math.h> #d ...
- spring属性配置执行过程,单列和原型区别
Spring配置中,采用属性注入时,当创建IOC容器时,也直接创建对象,并且执行相对应的setter方法 Student.java package com.scope; public class St ...
- FlowPortal-BPM——离线审批(邮箱审批)配置
一.将系统文件复制到安装目录下 二.以管理员身份运行bat安装程序 三.开启邮件离线审批服务 四.创建数据库表(JAVA数据库 或 .Net数据库) 五.配置config文件(发件箱.收件箱.错误问题 ...
- Java导包问题
eclipse中,导包只可以导到当前层,不可以再导入包里面的包中的类 包目录如下: 导包如下 package def; import abc.*; public class CCC { public ...
- Java代码度量分析工具:Designite简介
前言 在Java面向对象课程的学习过程中,我们需要使用度量工具来分析自己程序的代码结构.此类的度量工具有许多,或以插件形式存在于各个IDE中,或以.jar包的形式供用户使用.在这里,笔者向大家简单的介 ...
- Bluetooth Lowe Energy
BTL---------- // Wikipedia --------The first review paper to read when you counterred a new filed . ...
- Visual Studio 跨平台開發實戰(1) - Hello Xamarin! (转帖)
前言 應用程式發展的腳步, 從來沒有停過. 從早期的Windows 應用程式, 到網路時代的web 應用程式, 再到近幾年相當盛行的行動裝置應用程式(Mobile Application), 身為C# ...