码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj 3653 谈笑风生 —— 主席树

时间:2018-10-29 10:25:37      阅读:178      评论:0      收藏:0      [点我收藏+]

标签:ref   targe   putc   dfs   iostream   题目   update   href   www   

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3653

对于一个 (a,b,c),分成 b 是 a 的祖先和 b 在 a 子树里两部分;

第一部分 b 可以选 min(dep[a]-1,k) 个,c 可以选 siz[a]-1 个,乘起来即可;

第二部分,要求 ∑siz[u]-1,其中 u 是 a 子树内到 a 距离不超过 k 的点;

考虑到子树在 dfs 序上就是一段区间,距离范围又可以看作 dep 范围,所以在 dfs 序上建主席树,就可以查询对应 dep 的和;

别忘了 long long。

代码如下:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define mid ((l+r)>>1)
using namespace std;
typedef long long ll;
int const xn=3e5+5,xm=xn*20;
int n,hd[xn],ct,to[xn<<1],nxt[xn<<1],tim,dfn[xn],siz[xn],id[xn],dep[xn];
int cnt,ls[xm],rs[xm],rt[xn];
ll sum[xm];
int rd()
{
  int ret=0,f=1; char ch=getchar();
  while(ch<0||ch>9){if(ch==-)f=0; ch=getchar();}
  while(ch>=0&&ch<=9)ret=(ret<<3)+(ret<<1)+ch-0,ch=getchar();
  return f?ret:-ret;
}
int gt[20];
void wr(ll x)
{
  if(!x){putchar(0); return;}
  if(x<0){putchar(-); x=-x;}
  int t=0;
  while(x)gt[++t]=x%10,x/=10;
  for(int i=t;i;i--)putchar(gt[i]+0);
}
void add(int x,int y){to[++ct]=y; nxt[ct]=hd[x]; hd[x]=ct;}
void dfs(int x,int fa)
{
  dfn[x]=++tim; id[tim]=x; siz[x]=1; dep[x]=dep[fa]+1;
  for(int i=hd[x],u;i;i=nxt[i])
    {
      if((u=to[i])==fa)continue;
      dfs(u,x); siz[x]+=siz[u];
    }
}
void update(int &x,int y,int l,int r,int pos,int v)
{
  x=++cnt; ls[x]=ls[y]; rs[x]=rs[y]; sum[x]=sum[y]+v;
  if(l==r)return;
  if(pos<=mid)update(ls[x],ls[y],l,mid,pos,v);
  else update(rs[x],rs[y],mid+1,r,pos,v);
  sum[x]=sum[ls[x]]+sum[rs[x]];
}
ll query(int x,int y,int l,int r,int L,int R)
{
  if(l>=L&&r<=R)return sum[y]-sum[x];
  ll ret=0;
  if(mid>=L)ret+=query(ls[x],ls[y],l,mid,L,R);
  if(mid<R)ret+=query(rs[x],rs[y],mid+1,r,L,R);
  return ret;
}
int main()
{
  n=rd(); int q=rd();
  for(int i=1,x,y;i<n;i++)
    {
      x=rd(); y=rd();
      add(x,y); add(y,x);
    }
  dfs(1,0);
  for(int i=1;i<=n;i++)
    update(rt[i],rt[i-1],1,n,dep[id[i]],siz[id[i]]-1);
  for(int i=1,p,k;i<=q;i++)
    {
      p=rd(); k=rd();
      ll ans=(ll)min(dep[p]-1,k)*(siz[p]-1);
      ans+=query(rt[dfn[p]],rt[dfn[p]+siz[p]-1],1,n,dep[p]+1,dep[p]+k);
      wr(ans); puts("");
    }
  return 0;
}

 

bzoj 3653 谈笑风生 —— 主席树

标签:ref   targe   putc   dfs   iostream   题目   update   href   www   

原文地址:https://www.cnblogs.com/Zinn/p/9868577.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!