BZOJ 2648 世界树
题目传送门
分析:
喜 闻 乐 见 的虚树
但是建好虚树后的DP也非常的恶心
我们先考虑每个关键点的归哪个点管
先DFS一次计算儿子节点归属父亲
再DFS一次计算父亲节点归属儿子
然后然后我们对于虚树上的每条边计算一下
首先先找到分割点mid
那么向上归属的是红色部分
向下的是绿色部分
对于每条边都算一下就好了
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#define maxn 600005
#define INF 1ll<<50
using namespace std;
inline int getint()
{
int num=0,flag=1;char c;
while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
return num*flag;
}
int n;
int fir[maxn],nxt[maxn],to[maxn],cnt;
int f[maxn][21],dpt[maxn],sz[maxn];
int h[maxn],a[maxn];
int bl[maxn],F[maxn],ans[maxn];
int In[maxn],Out[maxn],cur;
int stk[maxn],top;
int pos[maxn],tot;
inline bool cmp(int x,int y){return In[x]<In[y];}
inline void newnode(int u,int v)
{to[++cnt]=v,nxt[cnt]=fir[u],fir[u]=cnt;}
inline void dfs(int u,int fa)
{
sz[u]=1;In[u]=++cur;
for(int i=fir[u];i;i=nxt[i])if(to[i]!=fa)
{
f[to[i]][0]=u,dpt[to[i]]=dpt[u]+1;
dfs(to[i],u);
sz[u]+=sz[to[i]];
}
Out[u]=cur;
}
inline int LCA(int u,int v)
{
if(dpt[u]<dpt[v])swap(u,v);
int d=dpt[u]-dpt[v];
for(int i=20;~i;i--)if(d&(1<<i))u=f[u][i];
if(u==v)return u;
for(int i=20;~i;i--)if(f[u][i]!=f[v][i])u=f[u][i],v=f[v][i];
return f[u][0];
}
inline int getdis(int u,int v)
{return dpt[u]+dpt[v]-2*dpt[LCA(u,v)];}
inline void dfs1(int u,int fa)
{
pos[++tot]=u;F[u]=sz[u];
for(int i=fir[u];i;i=nxt[i])if(to[i]!=fa)
{
dfs1(to[i],u);
if(!bl[to[i]])continue;
if(!bl[u]){bl[u]=bl[to[i]];continue;}
int tmp1=getdis(bl[to[i]],u),tmp2=getdis(bl[u],u);
if(tmp1<tmp2||(tmp1==tmp2&&bl[to[i]]<bl[u]))bl[u]=bl[to[i]];
}
}
inline void dfs2(int u,int fa)
{
for(int i=fir[u];i;i=nxt[i])if(to[i]!=fa)
{
int tmp1=getdis(bl[u],to[i]),tmp2=getdis(bl[to[i]],to[i]);
if(tmp1<tmp2||(tmp1==tmp2&&bl[u]<bl[to[i]]))bl[to[i]]=bl[u];
dfs2(to[i],u);
}
}
inline void solve(int fa,int u)
{
int son=u,mid=u;
for(int i=20;~i;i--)if(dpt[fa]+(1<<i)<dpt[son])son=f[son][i];
F[fa]-=sz[son];
if(bl[fa]==bl[u]){ans[bl[fa]]+=sz[son]-sz[u];return;}
for(int i=20;~i;i--)
{
int tmp=f[mid][i];
if(dpt[tmp]<=dpt[fa])continue;
int tmp1=getdis(tmp,bl[fa]);
int tmp2=getdis(tmp,bl[u]);
if(tmp2<tmp1||(tmp1==tmp2&&bl[u]<bl[fa]))mid=tmp;
}
ans[bl[fa]]+=sz[son]-sz[mid];
ans[bl[u]]+=sz[mid]-sz[u];
}
inline void solve()
{
int K=getint();top=0;
int tt=K;
for(int i=1;i<=K;i++)a[i]=h[i]=getint(),bl[h[i]]=h[i];h[++K]=1;
sort(h+1,h+K+1,cmp);
for(int i=K-1;i;i--)h[++K]=LCA(h[i],h[i+1]);
sort(h+1,h+K+1,cmp);K=unique(h+1,h+K+1)-h-1;
stk[++top]=h[1];
for(int i=2;i<=K;i++)
{
while(top&&Out[stk[top]]<In[h[i]])top--;
newnode(stk[top],h[i]);
stk[++top]=h[i];
}
dfs1(h[1],h[1]),dfs2(h[1],h[1]);
for(int i=1;i<=K;i++)
for(int j=fir[pos[i]];j;j=nxt[j])
solve(pos[i],to[j]);
for(int i=1;i<=K;i++)ans[bl[pos[i]]]+=F[pos[i]];
for(int i=1;i<=tt;i++)printf("%d%c",ans[a[i]],(i==tt)?'\n':' ');
for(int i=1;i<=K;i++)F[h[i]]=ans[h[i]]=fir[h[i]]=bl[h[i]]=0;
cnt=tot=0;
}
int main()
{
n=getint();
for(int i=1;i<n;i++)
{
int u=getint(),v=getint();
newnode(u,v),newnode(v,u);
}
dfs(1,1);
for(int j=1;j<=20;j++)for(int i=1;i<=n;i++)f[i][j]=f[f[i][j-1]][j-1];
memset(fir,0,sizeof fir);cnt=0;
int m=getint();
while(m--)solve();
}
BZOJ 2648 世界树的更多相关文章
- BZOJ 3572 世界树
Description 世界树是一棵无比巨大的树,它伸出的枝干构成了整个世界.在这里,生存着各种各样的种族和生灵,他们共同信奉着绝对公正公平的女神艾莉森,在他们的信条里,公平是使世界树能够生生不息.持 ...
- BZOJ 2648: SJY摆棋子
2648: SJY摆棋子 Time Limit: 20 Sec Memory Limit: 128 MBSubmit: 2968 Solved: 1011[Submit][Status][Disc ...
- BZOJ 2648: SJY摆棋子 kdtree
2648: SJY摆棋子 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=2648 Description 这天,SJY显得无聊.在家自己玩 ...
- BZOJ 2648/2716(SJY把件-KD_Tree)[Template:KD_Tree]
2648: SJY把件 Time Limit: 20 Sec Memory Limit: 128 MB Submit: 1180 Solved: 391 [id=2648" style= ...
- kd-tree注解 & bzoj 2648 & 2716 & 3053 解决问题的方法
[KD-TREE简介]于SYC1999大神"迷住"下一个.我开始接触这样的算法. 首先.这个概念大概能去百度百科.详细的实施.我在看RZZ的代码长大的. 我们能够想象在平面上有N个 ...
- BZOJ 2648 SJY摆棋子(KD Tree)
http://www.lydsy.com/JudgeOnline/problem.php?id=2648 题意: 思路: KDtree模板题. 参考自http://www.cnblogs.com/ra ...
- bzoj 2648: SJY摆棋子&&2716: [Violet 3]天使玩偶 --kdtree
2648: SJY摆棋子&&2716: [Violet 3]天使玩偶 Time Limit: 20 Sec Memory Limit: 128 MB Description 这天,S ...
- BZOJ 2648 kd-tree模板
学习了一下kd-tree的基本写法 http://blog.csdn.net/jiangshibiao/article/details/34144829 配合 http://www.bilibili. ...
- bzoj 2648 SJY摆棋子——KDtree
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2648 第一道KDtree! 学习资料:https://blog.csdn.net/zhl30 ...
随机推荐
- Comet OJ - Contest #5
Comet OJ - Contest #5 总有一天,我会拿掉给\(dyj\)的小裙子的. A 显然 \(ans = min(cnt_1/3,cnt_4/2,cnt5)\) B 我们可以感性理解一下, ...
- CF1045G AI robots
CF1045G AI robots 题目大意就不说了 这道题可以用CDQ分治做 但是,如何选择CDQ分治的维度一直是CDQ分治的难点所在 这道题我们有三种选择 1.让智商高的数智商低的 2.让看的近的 ...
- dotnet core 输出调试信息到 DebugView 软件
本文告诉大家如何在 dotnet core 输出调试信息到 DebugView 软件 在之前告诉小伙伴,如何在 WPF 输出调试信息到 DebugView 软件,请看文章 WPF 调试 获得追踪输出 ...
- js window.open 页面之间的通讯(不同源)
一:a页面1:打开b页面let isB= window.open('b.html','b'); 2:a页面发送消息给b页面 isB.postMessage("dsddfsdf&quo ...
- 16python的map函数,filter函数,reduce函数
map num_l = [1,6,8,9] def map_test(func,array): ret = [] for i in array: res = func(i) ret.append(re ...
- 一款类似loadRunner的优秀国产压力测试工具——kylinTOP测试与监控平台
市面上流行的压力/负载/性能测试工具多是来自国外,近年来国内的性能测试工具也如雨后春笋般崛起,但大部分产品是基于Jmeter开源内核包装起来的性能测试工具,其中也不乏佼佼者,如:kylinTOP测试与 ...
- $[SHOI2007]$ 园丁的烦恼 二维数点/树状数组
\(Sol\) 设一个矩阵的左上角为\((x_1,y_1)\),右下角为\((x_2,y_2)\),\(s_{x,y}\)是到\((1,1)\)二维前缀和,那么这个矩阵的答案显然是\(s_{x_2,y ...
- 使用这些idea插件让开发效率提高5倍
idea 有很多非常好用的插件,用好了这些插件能够极大的提高开发效率 插件用的好,bug 就追不上了我
- HDU3394 Railway 题解(边双连通分量)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3394 题目大意: 给定一个无向图,如果从一个点出发经过一些点和边能回到该点本身,那么一路走过来的这些点 ...
- linux入门系列2--CentOs图形界面操作及目录结构
上一篇文章"linux入门系列1--环境准备及linux安装"直观演示了虚拟机软件VMware和Centos操作系统的安装,按照文章一步一步操作,一定都可以安装成功.装好系统之后, ...