luoguP4103 [HEOI2014]大工程
题意
建出虚树DP。
设\(f[i]\)表示i的子树的第一问答案,\(minn[i]\)表示\(i\)的子树中到\(i\)最近的关键点,\(maxx[i]\)表示\(i\)的子树中到i距离最远的关键点,\(size[i]\)表示\(i\)子树中的关键点个数。
\(f[x]=\sum\limits_{y\in\ son_x}f[y]+size[y]*(tot-size[y])*w\)
\(minn\)和\(maxx\)就类似树的直径。
code:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1000010;
const ll inf=1e18;
int n,m,cnt,t,top,tot,tim;
int head[maxn],size[maxn],a[maxn],dep[maxn],sta[maxn],dfn[maxn];
int f[maxn][25];
ll ans1,ans2;
ll g[maxn],maxx[maxn],minn[maxn];
bool check[maxn];
struct edge{int to,nxt,dis;}e[maxn<<1];
inline int read()
{
char c=getchar();int res=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9')res=res*10+c-'0',c=getchar();
return res*f;
}
inline bool cmp(int x,int y){return dfn[x]<dfn[y];}
inline void add(int u,int v,int w)
{
e[++cnt].nxt=head[u];
head[u]=cnt;
e[cnt].to=v;
e[cnt].dis=w;
}
void dfs_pre(int x,int fa)
{
dep[x]=dep[fa]+1;dfn[x]=++tim;
for(int i=1;i<=t;i++)f[x][i]=f[f[x][i-1]][i-1];
for(int i=head[x];i;i=e[i].nxt)
{
int y=e[i].to;
if(y==fa)continue;
f[y][0]=x;dfs_pre(y,x);
}
}
inline int lca(int x,int y)
{
if(dep[x]>dep[y])swap(x,y);
for(int i=t;~i;i--)if(dep[f[y][i]]>=dep[x])y=f[y][i];
if(x==y)return x;
for(int i=t;~i;i--)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
return f[x][0];
}
inline int dis(int x,int y){return dep[x]+dep[y]-2*dep[lca(x,y)];}
void dfs(int x)
{
//cerr<<"test::"<<x<<endl;
g[x]=0;
maxx[x]=-inf,minn[x]=inf;size[x]=0;
if(check[x])size[x]=1,minn[x]=0,maxx[x]=0;
for(int i=head[x];i;i=e[i].nxt)
{
int y=e[i].to;
dfs(y);size[x]+=size[y];
g[x]+=g[y]+1ll*e[i].dis*size[y]*(tot-size[y]);
ans1=min(ans1,minn[x]+minn[y]+e[i].dis);
minn[x]=min(minn[x],minn[y]+e[i].dis);
ans2=max(ans2,maxx[x]+maxx[y]+e[i].dis);
maxx[x]=max(maxx[x],maxx[y]+e[i].dis);
}
}
void solve()
{
cnt=0;ans1=inf;ans2=-inf;
tot=read();
for(int i=1;i<=tot;i++)a[i]=read(),check[a[i]]=1;
sort(a+1,a+tot+1,cmp);
sta[top=1]=1,head[1]=0;
for(int i=1;i<=tot;i++)
{
if(a[i]==1)continue;
int x=lca(a[i],sta[top]);
if(x!=sta[top])
{
while(top>1&&dfn[sta[top-1]]>dfn[x])add(sta[top-1],sta[top],dis(sta[top-1],sta[top])),top--;
if(x!=sta[top-1])head[x]=0,add(x,sta[top],dis(x,sta[top])),sta[top]=x;
else add(x,sta[top],dis(x,sta[top])),top--;
}
head[a[i]]=0;sta[++top]=a[i];
}
for(int i=1;i<top;i++)add(sta[i],sta[i+1],dis(sta[i],sta[i+1]));
dfs(1);
printf("%lld %lld %lld\n",g[1],ans1,ans2);
for(int i=1;i<=tot;i++)check[a[i]]=0;
}
int main()
{
//freopen("test.in","r",stdin);
//freopen("test.out","w",stdout);
n=read();t=(int)log2(n)+1;
for(int i=1;i<n;i++)
{
int u=read(),v=read();
add(u,v,1),add(v,u,1);
}
dfs_pre(1,0);
m=read();
for(int i=1;i<=m;i++)solve();
return 0;
}
luoguP4103 [HEOI2014]大工程的更多相关文章
- [BZOJ3611][Heoi2014]大工程
[BZOJ3611][Heoi2014]大工程 试题描述 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶点上. 在 ...
- bzoj 3611 [Heoi2014]大工程(虚树+DP)
3611: [Heoi2014]大工程 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 408 Solved: 190[Submit][Status] ...
- 【LG4103】[HEOI2014]大工程
[LG4103][HEOI2014]大工程 题面 洛谷 题解 先建虚树,下面所有讨论均是在虚树上的. 对于第一问:直接统计所有树边对答案的贡献即可. 对于第\(2,3\)问:记\(f[x]\)表示在\ ...
- 3611: [Heoi2014]大工程
3611: [Heoi2014]大工程 链接 分析: 树形dp+虚树. 首先建立虚树,在虚树上dp. dp:sum[i]为i的子树中所有询问点之间的和.siz[i]为i的子树中有多少询问点,mn[i] ...
- P4103 [HEOI2014]大工程
题目 P4103 [HEOI2014]大工程 化简题目:在树上选定\(k\)个点,求两两路径和,最大的一组路径,最小的一组路径 做法 关键点不多,建个虚树跑一边就好了 \(sum_i\)为\(i\)子 ...
- BZOJ2286 [Sdoi2011]消耗战 和 BZOJ3611 [Heoi2014]大工程
2286: [Sdoi2011]消耗战 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 6371 Solved: 2496[Submit][Statu ...
- 【BZOJ3611】[Heoi2014]大工程 欧拉序+ST表+单调栈
[BZOJ3611][Heoi2014]大工程 Description 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶 ...
- [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 ...
随机推荐
- 爬取沪深a股数据
首先从东方财富网获取股票代码 再从网易财经下载股票历史数据 import requests import random from bs4 import BeautifulSoup as bs impo ...
- Leetcode 216. 组合总和 III
地址 https://leetcode-cn.com/problems/combination-sum-iii/ 找出所有相加之和为 n 的 k 个数的组合.组合中只允许含有 1 - 9 的正整数,并 ...
- 浅谈vue中的计算属性和侦听属性
计算属性 计算属性用于处理复杂的业务逻辑 计算属性具有依赖性,计算属性依赖 data中的初始值,只有当初始值改变的时候,计算属性才会再次计算 计算属性一般书写为一个函数,返回了一个值,这个值具有依赖性 ...
- C语言中,嵌套的if语句的一些经验...
double f(double x){double result;if(x<0)result=-x+sin(x);else if(x==0)result=0; else if(10> ...
- Codeforces Round #598 (Div. 3) A. Payment Without Change 水题
A. Payment Without Change You have a coins of value n and b coins of value 1. You always pay in exac ...
- 《Spring Cloud微服务 入门 实战与进阶》
很少在周末发文,还是由于昨晚刚收到实体书,还是耐不住性子马上发文了. 一年前,耗时半年多的时间,写出了我的第一本书<Spring Cloud微服务-全栈技术与案例解析>. 时至今日,一年的 ...
- Filter 原理
二.Filter 原理 2.1 Filter 概述 Filter(过滤器)是 DirectShow 中最基本的概念.DirectShow 是通过 Filter Graph 来管理 Filter 的.F ...
- LeetCode | 2 的幂
LeetCode 题库的第 231 题 —— 2 的幂 这题也是比较容易的一题,前提是找到规律即可.如果从 10 进制的角度观察 2 的幂次方,可能并不容易发现规律,那么可以从 2 进制的角度进行观察 ...
- 多模块springboot项目启动访问不了jsp页面
被springboot项目maven打包.启动.访问折腾的头都大了,一步一个坑,在解决了所有的问题之后,那种欣喜若狂的心情只有自己能体会,决定要好好记录一下了,废话不多说了,直接进入正题. 问题和排查 ...
- 转载:点云上实时三维目标检测的欧拉区域方案 ----Complex-YOLO
感觉是机器翻译,好多地方不通顺,凑合看看 原文名称:Complex-YOLO: An Euler-Region-Proposal for Real-time 3D Object Detection ...