luogu P4103 [HEOI2014]大工程 虚树 + 树形 DP
Description
Input
第一行 n 表示点数。
Output
输出 q 行,每行三个数分别表示代价和,最小代价,最大代价。
题解:建出来虚树后就不是很难了
#include<bits/stdc++.h>
#define setIO(s) freopen(s".in","r",stdin), freopen(s".out","w",stdout)
#define maxn 2000001
#define inf 1000000000
#define ll long long
using namespace std;
vector<int>G[maxn];
int edges,tim,root,top;
int hd[maxn], to[maxn<<1], val[maxn<<1], nex[maxn<<1];
int dep[maxn],Top[maxn],hson[maxn],siz[maxn],dfn[maxn],fa[maxn],arr[maxn],S[maxn],mk[maxn];
inline void add(int u,int v)
{
nex[++edges]=hd[u],hd[u]=edges,to[edges]=v,val[edges]=1;
}
void dfs1(int u,int ff)
{
fa[u]=ff,siz[u]=1,dep[u]=dep[ff]+1,dfn[u]=++tim;
for(int i=hd[u];i;i=nex[i])
{
int v=to[i];
if(v==ff) continue;
dfs1(v,u);
siz[u]+=siz[v];
if(siz[v]>siz[hson[u]]) hson[u]=v;
}
}
void dfs2(int u,int tp)
{
Top[u]=tp;
if(hson[u]) dfs2(hson[u],tp);
for(int i=hd[u];i;i=nex[i])
{
int v=to[i];
if(v==fa[u]||v==hson[u]) continue;
dfs2(v,v);
}
}
inline int LCA(int x,int y)
{
while(Top[x]!=Top[y])
{
dep[Top[x]] > dep[Top[y]] ? x = fa[Top[x]] : y = fa[Top[y]];
}
return dep[x] < dep[y] ? x : y;
}
inline int getdis(int x,int y)
{
return dep[x] + dep[y] - (dep[LCA(x,y)] << 1);
}
inline void addvir(int u,int v)
{
G[u].push_back(v);
}
void insert(int x)
{
if(top<=1) { S[++top]=x; return; }
int lca=LCA(x, S[top]);
if(lca == S[top]) { S[++top] = x; return; }
while(top > 1 && dep[S[top - 1]] >= dep[lca]) addvir(S[top - 1], S[top]), --top;
if(lca != S[top]) addvir(lca, S[top]), S[top] = lca;
S[++top] = x;
}
bool cmp(int i,int j)
{
return dfn[i] < dfn[j];
}
ll ans=0,a1,a2;
int size[maxn],d1[maxn],d2[maxn],dmin1[maxn],k,dmin2[maxn];
void DP(int x)
{
size[x]=mk[x];
d1[x]=d2[x]=0;
if(!mk[x]) d1[x]=d2[x]=-inf;
dmin1[x]=dmin2[x]=inf;
if(mk[x]) dmin1[x]=0;
for(int i=0;i<G[x].size();++i)
{
int v = G[x][i],w = dep[G[x][i]] - dep[x];
DP(v);
if(mk[v])
{
if(w <= dmin1[x]) dmin2[x]=dmin1[x], dmin1[x]=w;
else if(w < dmin2[x]) dmin2[x]=w;
}
else
{
if(w + dmin1[v] <= dmin1[x]) dmin2[x]=dmin1[x], dmin1[x]=w + dmin1[v];
else if(w + dmin1[v] < dmin2[x]) dmin2[x] = w + dmin1[v];
}
int curd=w+d1[v];
if(curd >= d1[x])
{
d2[x]=d1[x], d1[x]=curd;
}
else if(curd > d2[x])
{
d2[x] = curd;
}
ans+=1ll*size[v]*w*(k-size[v]),size[x]+=size[v];
}
a1=max(a1, 1ll*(d1[x] + d2[x]));
a2=min(a2, 1ll*(dmin1[x] + dmin2[x]));
}
void init(int x)
{
d1[x]=d2[x]=0;
dmin1[x]=dmin2[x]=inf;
for(int i=0;i<G[x].size();++i) init(G[x][i]);
G[x].clear();
}
inline void work()
{
scanf("%d",&k);
for(int i=1;i<=k;++i) scanf("%d",&arr[i]);
for(int i=1;i<=k;++i) mk[arr[i]] = 1;
sort(arr+1,arr+1+k,cmp);
top=S[0]=root=ans=0;
if(arr[1]!=1) S[top=1]=1;
for(int i=1;i<=k;++i) insert(arr[i]);
while(top > 1) addvir(S[top-1], S[top]),--top;
a1=-inf, a2=inf, DP(1);
printf("%lld %lld %lld\n",ans,a2,a1);
init(1);
for(int i=1;i<=k;++i) mk[arr[i]]=0;
}
int main()
{
// setIO("input");
int n;
scanf("%d",&n);
for(int i=1;i<n;++i)
{
int a,b;
scanf("%d%d",&a,&b);
add(a,b), add(b,a);
}
dfs1(1,0),dfs2(1,1);
int Q;
scanf("%d",&Q);
for(int i=1;i<=Q;++i) work();
return 0;
}
luogu P4103 [HEOI2014]大工程 虚树 + 树形 DP的更多相关文章
- 洛谷P4103 [HEOI2014]大工程(虚树 树形dp)
题意 链接 Sol 虚树. 首先建出虚树,然后直接树形dp就行了. 最大最小值直接维护子树内到该节点的最大值,然后合并两棵子树的时候更新一下答案. 任意两点的路径和可以考虑每条边两边的贡献,\(d[x ...
- BZOJ.3611.[HEOI2014]大工程(虚树 树形DP)
题目链接 要求的和.最大值.最小值好像都可以通过O(n)的树形DP做,总询问点数<=2n. 于是建虚树就可以了.具体DP见DP()函数,维护三个值sum[],mx[],mn[]. sum[]要开 ...
- bzoj 3611(洛谷 4103) [Heoi2014]大工程——虚树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3611 https://www.luogu.org/problemnew/show/P4103 ...
- bzoj 3611: [Heoi2014]大工程 虚树
题目: 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶点上. 在 2 个国家 a,b 之间建一条新通道需要的代价为树上 ...
- BZOJ 3611 [Heoi2014]大工程 ——虚树
虚树第二题.... 同BZOJ2286 #include <map> #include <cmath> #include <queue> #include < ...
- bzoj 3611[Heoi2014]大工程 虚树+dp
题意: 给一棵树 每次选 k 个关键点,然后在它们两两之间 新建 C(k,2)条 新通道. 求: 1.这些新通道的代价和 2.这些新通道中代价最小的是多少 3.这些新通道中代价最大的是多少 分析:较常 ...
- P4103 [HEOI2014]大工程
题目 P4103 [HEOI2014]大工程 化简题目:在树上选定\(k\)个点,求两两路径和,最大的一组路径,最小的一组路径 做法 关键点不多,建个虚树跑一边就好了 \(sum_i\)为\(i\)子 ...
- 【BZOJ-3572】世界树 虚树 + 树形DP
3572: [Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 1084 Solved: 611[Submit][Status ...
- 【BZOJ-2286】消耗战 虚树 + 树形DP
2286: [Sdoi2011消耗战 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 2120 Solved: 752[Submit][Status] ...
随机推荐
- 疑难杂症——bash: /dev/null: Permission denied
描述 在使用 Devstack 的时候需要时常切换用户su stack,此时会触发错误: root@mickeyfan-dev:~# su stack bash: /dev/null: Permiss ...
- 阶段1 语言基础+高级_1-3-Java语言高级_1-常用API_1_第4节 ArrayList集合_16-ArrayList练习一_存储随机数
循环6次就是6.fori 循环子在外部+1就是得到的1到33的数字 list.fori遍历集合 自动生for循环的代码
- MySQL 修改密码,
ALTER USER 'root'@'localhost' IDENTIFIED BY 'admin';
- 【Linux】limits.conf 不重启就生效或者不生效的原因
前阵子,我要用到使LInux的文件打开数为65534个,而且需要永久生效,于是将配置写到了: vim /etc/security/limits.conf * soft nofile 65534* ha ...
- Web前端开发 --》 如何实现页面同时在移动端和pc端的兼容问题
很简单,只需要在html文件中对你引入的css进行一个类似于媒体查询的操作 <!DOCTYPE html> <html lang="en"> <hea ...
- SpringBoot(九) -- SpringBoot与数据访问
一.简介 对于数据访问层,无论是SQL还是NOSQL,Spring Boot默认采用整合Spring Data的方式进行统一处理,添加大量自动配置,屏蔽了很多设置.引入各种xxxTemplate,xx ...
- [Luogu P3825] [NOI2017] 游戏 (2-SAT)
[Luogu P3825] [NOI2017] 游戏 (2-SAT) 题面 题面较长,略 分析 看到这些约束,应该想到这是类似2-SAT的问题.但是x地图很麻烦,因为k-SAT问题在k>2的时候 ...
- 生日蛋糕 (poj1190) (dfs剪枝)
[题目描述] 7月17日是Mr.W的生日,ACM-THU为此要制作一个体积为Nπ的M层生日蛋糕,每层都是一个圆柱体. 设从下往上数第i(1 <= i <= M)层蛋糕是半径为Ri, 高度为 ...
- 最长回文子序列/最长回文子串(DP,马拉车)
字符子串和字符子序列的区别 字符字串指的是字符串中连续的n个字符:如palindrome中,pa,alind,drome等都属于它的字串 而字符子序列指的是字符串中不一定连续但先后顺序一致的n个字符: ...
- 06-File-文件
文件 长久保存信息的一种数据信息集合 常用操作 打开关闭(文件一旦打开,需要关闭操作) 读写内容 查找 open函数 open函数负责打开文件,带有很多参数 第一个参数: 必须有,文件的路径和名称 m ...