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

【点分治】bzoj1316 树上的询问

时间:2015-03-11 10:49:12      阅读:122      评论:0      收藏:0      [点我收藏+]

标签:

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define MAXN 10001
#define INF 2147483647
typedef pair<int,int> Point;
int n,m,Ks[101];
bool vis[101];
int v[MAXN<<1],w[MAXN<<1],first[MAXN],next[MAXN<<1],en;
int dis[MAXN],En,last;
void AddEdge(const int &U,const int &V,const int &W)
{
    v[++en]=V;
    w[en]=W;
    next[en]=first[U];
    first[U]=en;
}
bool centroid[MAXN];
int size[MAXN];
int calc_sizes(int U,int Fa)
{
    int res=1;
    for(int i=first[U];i;i=next[i])
      if(v[i]!=Fa&&(!centroid[v[i]]))
        res+=calc_sizes(v[i],U);
    return size[U]=res;
}
Point calc_centroid(int U,int Fa,int nn)
{
    Point res=make_pair(INF,-1);
    int sum=1,maxv=0;
    for(int i=first[U];i;i=next[i])
      if(v[i]!=Fa&&(!centroid[v[i]]))
        {
          res=min(res,calc_centroid(v[i],U,nn));
          maxv=max(maxv,size[v[i]]);
          sum+=size[v[i]];
        }
    maxv=max(maxv,nn-sum);
    res=min(res,make_pair(maxv,U));
    return res;
}
void calc_dis(int U,int Fa,int d)
{
    dis[En++]=d;
    for(int i=first[U];i;i=next[i])
      if(v[i]!=Fa&&(!centroid[v[i]]))
        calc_dis(v[i],U,d+w[i]);
}
void calc_pairs()
{
    sort(dis,dis+last);
    for(int i=last;i<En;++i) 
      for(int j=1;j<=m;++j)
        {
          int *p=lower_bound(dis,dis+last,Ks[j]-dis[i]);
          if(*p==Ks[j]-dis[i]&&p!=dis+last)
            vis[j]=1;
        }
}
void solve(int U)
{
    calc_sizes(U,-1);
    int s=calc_centroid(U,-1,size[U]).second;
    centroid[s]=1;
    for(int i=first[s];i;i=next[i])
      if(!centroid[v[i]])
        solve(v[i]);
    En=0; dis[En++]=0;
    for(int i=first[s];i;i=next[i])
      if(!centroid[v[i]])
        {
          last=En;
          calc_dis(v[i],s,w[i]);
          calc_pairs();
        }
    centroid[s]=0;
}
int main()
{
    scanf("%d%d",&n,&m);
    int a,b,c;
    for(int i=1;i<n;++i)
      {
        scanf("%d%d%d",&a,&b,&c);
        AddEdge(a,b,c); AddEdge(b,a,c);
      }
    for(int i=1;i<=m;++i) scanf("%d",&Ks[i]);
    solve(1);
    for(int i=1;i<=m;++i) puts((Ks[i]==0||vis[i])?"Yes":"No");
    return 0;
}

【点分治】bzoj1316 树上的询问

标签:

原文地址:http://www.cnblogs.com/autsky-jadek/p/4329126.html

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