题目大意:

n个城市构成一个树 m支军队 每只军队守卫 在xi到yi的最短路径上的城市

q个重要人物从vi出发 找到离根最近的点使从vi到这个点上所有路径都可以被至少ki个军队完全覆盖

输出每个答案的点的深度

思路:

对于每个军队 可以拆成两个链

在深度较大的节点的权值线段树上把深度较低的点+1

然后合并线段树

查询的时候查询子树线段树中第k小的

 #include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#define inf 2139062143
#define ll long long
#define MAXN 200100
using namespace std;
inline int read()
{
int x=,f=;char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') f=-;ch=getchar();}
while(isdigit(ch)) {x=x*+ch-'';ch=getchar();}
return x*f;
}
int n,m,f[MAXN][],dep[MAXN],hsh[MAXN],sz[MAXN],val[MAXN<<];
int nxt[MAXN<<],fst[MAXN],to[MAXN<<],cnt,rt[MAXN],ls[MAXN<<],rs[MAXN<<];
void add(int u,int v) {nxt[++cnt]=fst[u],fst[u]=cnt,to[cnt]=v;}
void dfs(int x)
{
hsh[x]=++cnt,sz[x]=;
for(int i=;(<<i)<=dep[x];i++) f[x][i]=f[f[x][i-]][i-];
for(int i=fst[x];i;i=nxt[i])
if(to[i]!=f[x][]){dep[to[i]]=dep[x]+,f[to[i]][]=x;dfs(to[i]);sz[x]+=sz[to[i]];}
}
int lca(int u,int v)
{
if(dep[u]<dep[v]) swap(u,v);
int t=dep[u]-dep[v];
for(int i=;i<=;i++)
if((<<i)&t) u=f[u][i];
if(u==v) return u;
for(int i=;i>=;i--)
if(f[u][i]!=f[v][i]) u=f[u][i],v=f[v][i];
return f[u][];
}
void inst(int &k,int l,int r,int x)
{
if(!k) k=++cnt;
val[k]++;
if(l==r) return ;
int mid=(l+r)>>;
if(x<=mid) inst(ls[k],l,mid,x);
else inst(rs[k],mid+,r,x);
}
void merge(int &a,int b)
{
if(!a) {a=b;return ;}
if(!b) return ;
val[a]+=val[b];
merge(ls[a],ls[b]);
merge(rs[a],rs[b]);
}
int query(int k1,int k2,int l,int r,int a,int b)
{
if(!k1) return ;
if(l==a&&r==b) return val[k1]-val[k2];
int mid=(l+r)>>;
if(b<=mid) return query(ls[k1],ls[k2],l,mid,a,b);
else if(a>mid) return query(rs[k1],rs[k2],mid+,r,a,b);
else return query(ls[k1],ls[k2],l,mid,a,mid)+query(rs[k1],rs[k2],mid+,r,mid+,b);
}
int ans(int k1,int k2,int l,int r,int x)
{
//cout<<l<<" "<<r<<" "<<val[k1]-val[k2]<<" "<<x<<endl;
if(l==r) return l;
int mid=(l+r)>>;
if(x<=val[ls[k1]]-val[ls[k2]]) return ans(ls[k1],ls[k2],l,mid,x);
else return ans(rs[k1],rs[k2],mid+,r,x-val[ls[k1]]+val[ls[k2]]);
}
int main()
{
n=read(),m=read();int a,b,c;
for(int i=;i<n;i++) {a=read(),b=read();add(a,b);add(b,a);}
cnt=;dfs();cnt=;
while(m--)
{
a=read(),b=read(),c=lca(a,b);
inst(rt[hsh[a]],,n,dep[c]+);
inst(rt[hsh[b]],,n,dep[c]+);
}
for(int i=;i<=n;i++) merge(rt[i],rt[i-]);
m=read();
while(m--)
{
a=read(),b=read();
//cout<<query(rt[hsh[a]+sz[a]-1],rt[hsh[a]-1],1,n,1,dep[a]+1)<<endl;
if(query(rt[hsh[a]+sz[a]-],rt[hsh[a]-],,n,,dep[a]+)<b) puts("");
else printf("%d\n",dep[a]+-ans(rt[hsh[a]+sz[a]-],rt[hsh[a]-],,n,b));
}
}

最新文章

  1. 基于jQuery的Validate表单验证
  2. Unity3D 之脚本架构,优雅地管理你的代码
  3. html上下结构(上部固定高度,下部平铺)
  4. texlive2015+texstudio
  5. Jquery中css()方法获取边框长度
  6. Python进阶04 函数的参数对应
  7. javascript背景淡入淡出
  8. html19-----视频,音乐的插入
  9. HTTP Authorization
  10. apache Alias使用问题
  11. iOS 推送全解析,你不可不知的所有 Tips!
  12. 京东商品及评论爬虫(selenium)
  13. Maven 学习总结 (四)之 测试
  14. 微信小程序上拉下拉刷新
  15. OpenCV 学习笔记 06 图像检索以及基于图像描述符的搜索
  16. node.js初识06
  17. HDU 2051 Bitset
  18. 基于asp.net + easyui框架,一步步学习easyui-datagrid——实现分页和搜索(二)
  19. 同样的输入,为什么Objects.hash()方法返回的hash值每次不一样?
  20. POJ 1389 Area of Simple Polygons 扫描线+线段树面积并

热门文章

  1. 你知道你常用的dos和linux命令吗?
  2. assert.notDeepEqual()
  3. mysql查询排名
  4. JS 去除字符串空格
  5. 3D标签云
  6. 九度oj 题目1490:字符串链接
  7. CodeForces788B 欧拉路
  8. hdu -1251 统计难题(字典树水题)
  9. cogs——1364. 聚会
  10. Elasticsearch自定义客户端(TransportClient)资源池