【BZOJ】3572: [Hnoi2014]世界树
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3572
算是细节题了吧。。
构出虚树,考虑z正反DP两次求出虚树中每个点最近的议事处是哪一个点,然后枚举虚树中的每一条边${X->Y}$,对于两点间在原树中的路径,显然存在一个分界点使得自分界点之上的所有点归最靠近$X$的议事处管辖,之下的点归最靠经$Y$的议事处管辖,还有一些没有考虑过的点,另外统计一下即可。
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
#define maxn 310000
#define llg long long
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
llg n,m,f[maxn][],deep[maxn],dfn[maxn],cnt,size[maxn],s[maxn],top;
llg a[maxn],b[maxn],bel[maxn],c[maxn],rem[maxn],ans[maxn]; vector<llg>g[maxn]; bool cmp(llg x,llg y) {return dfn[x]<dfn[y];} void link(llg x,llg y){g[x].push_back(y);} void dfs(llg x,llg fa)
{
dfn[x]=++cnt;
llg w=g[x].size(),v;
size[x]=;
for (llg i=;i<w;i++)
{
v=g[x][i];
if (v==fa) continue;
deep[v]=deep[x]+;
f[v][]=x;
dfs(v,x);
size[x]+=size[v];
}
} void make_f()
{
for (llg j=;j<=;j++)
for (llg i=;i<=n;i++)
f[i][j]=f[f[i][j-]][j-];
} llg lca(llg x,llg y)
{
if (deep[x]<deep[y]) swap(x,y);
for (llg i=;i>=;i--)
if (deep[f[x][i]]>=deep[y])
x=f[x][i];
if (x==y) return x;
for (llg i=;i>=;i--)
if (f[x][i]!=f[y][i])
{
x=f[x][i];
y=f[y][i];
}
return f[x][];
} llg dis(llg x,llg y) {return deep[x]+deep[y]-*deep[lca(x,y)];} void init()
{
cin>>n;
llg x,y;
for (llg i=;i<n;i++)
{
scanf("%lld%lld",&x,&y);
g[x].push_back(y),g[y].push_back(x);
}
deep[]=;
dfs(,-);
make_f();
} void build_tree()
{
cnt=top=;
sort(a+,a+m+,cmp);
if (bel[]!=) s[++top]=;
for (llg i=;i<=m;i++)
{
llg x=a[i],fa=;
while (top)
{
fa=lca(s[top],x);
if (top> && deep[fa]<deep[s[top-]])
{
link(s[top-],s[top]);
top--;
}
else
if (deep[fa]<deep[s[top]]){ link(fa,s[top--]); break;}
else
break; }
if (s[top]!=fa) s[++top]=fa;
s[++top]=x;
}
while (--top) link(s[top],s[top+]);
} void dfs1(llg x,llg fa)
{
c[++cnt]=x;
rem[x]=size[x];
llg w=g[x].size(),v;
for (llg i=;i<w;i++)
{
v=g[x][i];
dfs1(v,x);
if (v==fa) continue;
if (!bel[v]) continue;//why
llg d1=dis(bel[v],x),d2=dis(bel[x],x);
if ((d1==d2 && bel[v]<bel[x]) || (d1<d2) || !bel[x]) bel[x]=bel[v];
}
} void dfs2(llg x,llg fa)
{
llg w=g[x].size(),v;
for (llg i=;i<w;i++)
{
v=g[x][i];
if (v==fa) continue;
llg d1=dis(bel[x],v),d2=dis(bel[v],v);
if ((d1==d2 && bel[v]>bel[x]) || (d1<d2) || !bel[v])
bel[v]=bel[x];
dfs2(v,x);
}
} void work(llg x,llg y)
{
llg son=y,mid=y;
for (llg i=;i>=;i--)
if (deep[f[son][i]]>deep[x]) son=f[son][i];
rem[x]-=size[son];
if (bel[x]==bel[y]) {ans[bel[x]]+=size[son]-size[y]; return;}
for (llg i=;i>=;i--)
{
llg ne=f[mid][i];
if (deep[ne]<=deep[x]) continue;
llg d1=dis(bel[x],ne),d2=dis(bel[y],ne);
if (d1>d2 || (d1==d2 && bel[y]<bel[x])) mid=ne;
}
ans[bel[x]]+=size[son]-size[mid];
ans[bel[y]]+=size[mid]-size[y];
} int main()
{
yyj("tree");
init();
llg T; cin>>T;
for (llg i=;i<=n;i++)
{
while (!g[i].empty()) g[i].pop_back();
//std::vector(c).swap(vec);
}
while (T--)
{
scanf("%lld",&m);
for (llg i=;i<=m;i++)
{
scanf("%lld",&a[i]);
b[i]=a[i]; bel[a[i]]=a[i];
}
build_tree();
dfs1(,-); dfs2(,-);
for (llg i=;i<=cnt;i++)
{
llg w=g[c[i]].size();
llg la=-;
for (llg j=;j<w;j++)
{
if (g[c[i]][j]==la) continue;
la=g[c[i]][j];
work(c[i],g[c[i]][j]);
}
}
for (llg i=;i<=cnt;i++) ans[bel[c[i]]]+=rem[c[i]];
for (llg i=;i<=m;i++) printf("%lld ",ans[b[i]]);
printf("\n");
for (llg i=;i<=cnt;i++)
{
ans[c[i]]=bel[c[i]]=rem[c[i]]=;
while (!g[c[i]].empty()) g[c[i]].pop_back();
}
}
return ;
}
【BZOJ】3572: [Hnoi2014]世界树的更多相关文章
- BZOJ 3572: [Hnoi2014]世界树
BZOJ 3572: [Hnoi2014]世界树 标签(空格分隔): OI-BZOJ OI-虚数 OI-树形dp OI-倍增 Time Limit: 20 Sec Memory Limit: 512 ...
- bzoj 3572: [Hnoi2014]世界树 虚树 && AC500
3572: [Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 520 Solved: 300[Submit][Status] ...
- bzoj 3572 [Hnoi2014]世界树(虚树+DP)
3572: [Hnoi2014]世界树 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 645 Solved: 362[Submit][Status] ...
- bzoj 3572 [Hnoi2014]世界树——虚树
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3572 关于虚树:https://www.cnblogs.com/zzqsblog/p/556 ...
- BZOJ 3572: [Hnoi2014]世界树 虚树 树形dp
https://www.lydsy.com/JudgeOnline/problem.php?id=3572 http://hzwer.com/6804.html 写的时候参考了hzwer的代码,不会写 ...
- bzoj 3572: [Hnoi2014]世界树 虚树
题目: Description 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森,在他们的信条里,公平是使世界树能够生生 ...
- BZOJ 3572 [HNOI2014]世界树 (虚树+DP)
题面:BZOJ传送门 洛谷传送门 题目大意:略 细节贼多的虚树$DP$ 先考虑只有一次询问的情况 一个节点$x$可能被它子树内的一个到x距离最小的特殊点管辖,还可能被管辖fa[x]的特殊点管辖 跑两次 ...
- BZOJ 3572: [Hnoi2014]世界树 [虚树 DP 倍增]
传送门 题意: 一棵树,多次询问,给出$m$个点,求有几个点到给定点最近 写了一晚上... 当然要建虚树了,但是怎么$DP$啊 大爷题解传送门 我们先求出到虚树上某个点最近的关键点 然后枚举所有的边$ ...
- 【BZOJ】3572: [Hnoi2014]世界树 虚树+倍增
[题意]给定n个点的树,m次询问,每次给定ki个特殊点,一个点会被最近的特殊点控制,询问每个特殊点控制多少点.n,m,Σki<=300000. [算法]虚树+倍增 [题解]★参考:thy_asd ...
随机推荐
- Java程序员秋招面经大合集(BAT美团网易小米华为中兴等)
Cvte提前批 阿里内推 便利蜂内推 小米内推 金山wps内推 多益网络 拼多多学霸批 搜狗校招 涂鸦移动 中国电信it研发中心 中兴 华为 苏宁内推 美团内推 百度 腾讯 招商银行信用卡 招银网络科 ...
- django框架基础
所有的Web应用本质上就是一个socket服务端,而用户的浏览器就是一个socket客户端. 这样我们就可以自己实现Web框架了. 最简单的web框架 import socket sk = socke ...
- jq如何获取服务器的时间
$.ajax({type:"OPTIONS",url:"/",complete:function(x){alert(x.getResponseHeader(&q ...
- Caused by: java.lang.ClassNotFoundException: Illegal access: this web application instance has been stopped already. Could not load [org.jboss.netty.util.internal.ByteBufferUtil]. The following stack
Caused by: java.lang.ClassNotFoundException: Illegal access: this web application instance has been ...
- oracle闪回、闪回数据归档Flashback Data Archive (Oracle Total Recall)的真正强大之处、11gR2增强以及合理使用
oracle的闪回很早就出来了,准确的说一直以来应该都较少被真正用户广为使用,除了dba和极少部分开发人员偶尔用于逻辑出错.误删恢复之外,较少被用于产生更有价值的用途. 各种闪回表flashback ...
- 【题解】Luogu CF915E Physical Education Lessons
原题传送门:CF915E Physical Education Lessons 前置芝士:珂朵莉树 窝博客里对珂朵莉树的介绍 没什么好说的自己看看吧 这道题很简单啊 每个操作就是区间赋值,顺带把总和修 ...
- docker 实践
https://doc.yonyoucloud.com/doc/docker_practice/etcd/etcdctl.html 启动http restful API docker批量映射端口 怎么 ...
- oracle基础——内存管理、优化
内存图解: 自动管理:11g:AMM 10g:ASMM SGA(system global area):由所有服务进程和后台进程共享 PGA(program global area): 由每个服务 ...
- topcoder srm 686 div1
problem1 link 左括号和右括号较少的一种不会大于20.假设左括号少.设$f[i][mask][k]$表示处理了前$i$个字符,其中留下的字符以$k$开头($k=0$表示'(',$k=1$表 ...
- FireMonkey 源码学习(4)
(4)DoDrawLayout DoDrawLayout函数的源代码分析如下: procedure TTextLayoutNG.DoDrawLayout(const ACanvas: TCanvas) ...