P4103 [HEOI2014]大工程
题目
P4103 [HEOI2014]大工程
化简题目:在树上选定\(k\)个点,求两两路径和,最大的一组路径,最小的一组路径
做法
关键点不多,建个虚树跑一边就好了
\(sum_i\)为\(i\)子树各关键点到根节点的距离和,\(small_i\)为其最小值,\(big_i\)为其最大值
My complete code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<vector>
using namespace std;
typedef long long LL;
inline LL Read(){
LL x(0),f(1); char c=getchar();
while(c<'0'||c>'9'){
if(c=='-')f=-1; c=getchar();
}
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-'0',c=getchar();
return x*f;
}
const LL maxn=2*1e6+9,inf=0x3f3f3f3f;
LL n,ans1,ans2,ans3;
LL a[maxn],sta[maxn];
struct Tree{
LL tot,num,head[maxn],dfn[maxn],dep[maxn],size[maxn],key[maxn],inc[maxn][25],fa[maxn];
LL sum[maxn],big[maxn],small[maxn];
struct node{
LL to,next,d;
}dis[maxn];
inline void Add(LL u,LL v,LL d){
dis[++num]=(node){v,head[u],d},head[u]=num;
}
void Fdfs(LL u){
dfn[u]=++tot;
inc[u][0]=fa[u];
for(LL i=1;i<=21;++i)
inc[u][i]=inc[inc[u][i-1]][i-1];
for(LL i=head[u];i;i=dis[i].next){
LL v(dis[i].to);
if(v==fa[u])
continue;
fa[v]=u,dep[v]=dep[u]+1,
Fdfs(v);
}
}
inline LL Lca(LL x,LL y){
if(dep[x]<dep[y])
swap(x,y);
for(LL i=20;i>=0;--i)
if(dep[inc[x][i]]>=dep[y])
x=inc[x][i];
if(x==y)
return x;
for(LL i=20;i>=0;--i)
if(inc[x][i]!=inc[y][i])
x=inc[x][i],y=inc[y][i];
return inc[x][0];
}
void Dfs(LL u){
size[u]=key[u],
sum[u]=big[u]=0,
small[u]=(key[u]==0)?inf:0;
for(LL i=head[u];i;i=dis[i].next){
LL v(dis[i].to),d(dis[i].d);
Dfs(v);
if(size[u]){
ans1+=size[u]*size[v]*d+size[u]*sum[v]+size[v]*sum[u],
ans2=min(ans2,small[u]+d+small[v]),
ans3=max(ans3,big[u]+d+big[v]);
}
sum[u]+=sum[v]+size[v]*d;
small[u]=min(small[u],small[v]+d),
big[u]=max(big[u],big[v]+d),
size[u]+=size[v];
}
head[u]=key[u]=0;
}
}T,X;
inline bool cmp1(LL x,LL y){
return T.dfn[x]<T.dfn[y];
}
inline void Solve(){
LL kase=Read();
while(kase--){
LL num(Read());
for(LL i=1;i<=num;++i)
a[i]=Read();
sort(a+1,a+1+num,cmp1);
LL tp; sta[tp=1]=1;
for(LL i=1;i<=num;++i){
X.key[a[i]]=1;
LL lca=T.Lca(a[i],sta[tp]);
while(T.dep[sta[tp]]>T.dep[lca])
if(T.dep[sta[tp-1]]<T.dep[lca])
X.Add(lca,sta[tp],T.dep[sta[tp]]-T.dep[lca]),sta[tp]=lca;
else
X.Add(sta[tp-1],sta[tp],T.dep[sta[tp]]-T.dep[sta[tp-1]]),--tp;
if(sta[tp]!=a[i])
sta[++tp]=a[i];
}
while(tp>1)
X.Add(sta[tp-1],sta[tp],T.dep[sta[tp]]-T.dep[sta[tp-1]]),--tp;
ans1=0,ans2=inf,ans3=0;
X.Dfs(1);
printf("%lld %lld %lld\n",ans1,ans2,ans3);
}
}
int main(){
n=Read();
for(LL i=1;i<n;++i){
LL u(Read()),v(Read());
T.Add(u,v,1),T.Add(v,u,1);
}
T.Fdfs(1);
Solve();
return 0;
}/*
*/
P4103 [HEOI2014]大工程的更多相关文章
- luogu P4103 [HEOI2014]大工程 虚树 + 树形 DP
Description 国家有一个大工程,要给一个非常大的交通网络里建一些新的通道. 我们这个国家位置非常特殊,可以看成是一个单位边权的树,城市位于顶点上. 在 2 个国家 a,b 之间建一条新通 ...
- 洛谷P4103 [HEOI2014]大工程(虚树 树形dp)
题意 链接 Sol 虚树. 首先建出虚树,然后直接树形dp就行了. 最大最小值直接维护子树内到该节点的最大值,然后合并两棵子树的时候更新一下答案. 任意两点的路径和可以考虑每条边两边的贡献,\(d[x ...
- [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] ...
- 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 ...
随机推荐
- mongoDB Java SDK CRUD操作例子
Example: package com.sdfwerwer.test; import java.net.UnknownHostException; import com.mongodb.Basi ...
- SQLServer中存储过程StoredProcedure创建及C#调用(转)
此文作为入门了解用,转自http://www.2cto.com/database/201502/378260.html 存储过程就是已经编译好的.优化过的放在数据库服务器中的一些SQL语句:可供应用程 ...
- Effective JavaScript Item 10 避免使用with
本系列作为Effective JavaScript的读书笔记. Item 9:避免使用withkeyword 重点: 设计withkeyword本来是为了让代码变简洁,可是却起到了相反的效果.比方: ...
- springboot整合docker部署(两种构建Docker镜像方式)
项目结构 package hello; import org.springframework.boot.SpringApplication; import org.springframework.bo ...
- awk.sed.grep三剑客详解
事前准备1.主机node1:172.16.133.112.作为实验的文件/etc/passwd /etc/fstab qinqin cp /etc/passwd . cp /etc/fstab . 一 ...
- unity开发android游戏(一)搭建Unity安卓开发环境
unity开发android游戏(一)搭建Unity安卓开发环境 分类: Unity2014-03-23 16:14 5626人阅读 评论(2) 收藏 举报 unity开发androidunity安卓 ...
- 《Lucene in Action 第二版》第三章节的学习总结----IndexSearcher以及Term和QueryParser
本章节告诉我们怎么用搜索.通过这章节的学习,虽然搜索的内部原理不清楚,但是至少应该学会简单的编写搜索程序了本章节,需要掌握如下几个主要API1.IndexSearcher类:搜索索引的门户,发起者. ...
- NFS网络文件系统服务(配置实战)
NFS网络文件系统服务(实战) NFS(Network File System)即网络文件系统,它允许网络中的计算机之间通过TCP/IP网络共享资源.让不同的主机系统(NFS的客户端)可以透明地读写位 ...
- grep -rl tttt /data/ 命令在 /data 目录下面搜寻包含tttt字符的命令
grep --help -R, -r, --recursive equivalent to --directories=recurse -l, --files-with-matches print o ...
- saltstack之定时管理
1.设置定时任务 /srv/salt/cron/ntpdate.sls /usr/sbin/ntpdate 10.31.10.3; /sbin/hwclock -w: cron.present: - ...