BZOJ3631松鼠的新家

3631: [JLOI2014]松鼠的新家
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 1152 Solved: 586
Description
松鼠的新家是一棵树,前几天刚刚装修了新家,新家有n个房间,并且有n-1根树枝连接,每个房间都可以相互到达,且俩个房间之间的路线都是唯一的。天哪,他居然真的住在“树”上。松鼠想邀请小熊维尼前来参观,并且还指定一份参观指南,他希望维尼能够按照他的指南顺序,先去a1,再去a2,……,最后到an,去参观新家。
可是这样会导致维尼重复走很多房间,懒惰的维尼不听地推辞。可是松鼠告诉他,每走到一个房间,他就可以从房间拿一块糖果吃。维尼是个馋家伙,立马就答应了。
现在松鼠希望知道为了保证维尼有糖果吃,他需要在每一个房间各放至少多少个糖果。因为松鼠参观指南上的最后一个房间an是餐厅,餐厅里他准备了丰盛的大餐,所以当维尼在参观的最后到达餐厅时就不需要再拿糖果吃了。
Input
第一行一个整数n,表示房间个数
第二行n个整数,依次描述a1-an
接下来n-1行,每行两个整数x,y,表示标号x和y的两个房间之间有树枝相连。
Output
一共n行,第i行输出标号为i的房间至少需要放多少个糖果,才能让维尼有糖果吃。
Sample Input
5
1 4 5 3 2
1 2
2 4
2 3
4 5
Sample Output
1
2
1
2
1
HINT
2<= n <=300000
第一次遇到链剖但是没有询问的题。。突然就一下子没思路了。。
其实只要打个标记就可以了。。
链剖裸题。。
附上本蒟蒻的代码:

#include
#include
using namespace std;
#define MAXN 300001
int n,cnt,sz,a[MAXN],h[MAXN],father[MAXN][19],pos[MAXN],belong[MAXN],deep[MAXN],size[MAXN],delta[MAXN<<2];
bool vis[MAXN];
struct data
{
    int to,next;
}edge[MAXN<<1];
struct kx
{
    int value;
}node[MAXN<<2];

int read()
{
    int w=0,c=1; char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-') c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      w=w*10+ch-'0',ch=getchar();
    return w*c;
}

void add(int u,int v)
{
    cnt++,edge[cnt].next=h[u],h[u]=cnt,edge[cnt].to=v;
}

void dfs1(int x)
{
    int i;
    vis[x]=true,size[x]=1;
    for (i=1;i<=18;i++)
      {
        if (deep[x]<(1<break;
        father[x][i]=father[father[x][i-1]][i-1];
      }
    for (i=h[x];i;i=edge[i].next)
      {
        if (vis[edge[i].to]) continue;
        deep[edge[i].to]=deep[x]+1,father[edge[i].to][0]=x;
        dfs1(edge[i].to),size[x]+=size[edge[i].to];
      }
}

void dfs2(int x,int chain)
{
    int k=0,i;
    pos[x]=++sz,belong[x]=chain;
    for (i=h[x];i;i=edge[i].next)
      if (deep[edge[i].to]>deep[x] && size[edge[i].to]>size[k])
        k=edge[i].to;
    if (!k) return;
    dfs2(k,chain);
    for (i=h[x];i;i=edge[i].next)
      if (deep[edge[i].to]>deep[x] && k!=edge[i].to)
        dfs2(edge[i].to,edge[i].to);
}

void build(int s,int l,int r)
{
    if (l==r) return;
    build(s*2,l,(l+r)/2),build(s*2+1,(l+r)/2+1,r);
}

void paint(int s,int z)
{
    node[s].value+=z,delta[s]+=z;
}

void pushdown(int s)
{
    paint(s*2,delta[s]),paint(s*2+1,delta[s]),delta[s]=0;
}

int query(int s,int l,int r,int x)
{
    int mid=(l+r)/2;
    if (l==r) return node[s].value;
    pushdown(s);
    if (x<=mid) return query(s*2,l,mid,x);
    return query(s*2+1,mid+1,r,x);
}

void insert(int s,int l,int r,int x,int y,int z)
{
    int mid=(l+r)/2;
    if (l==x && y==r)
      {
        paint(s,z);
        return;
      }
    pushdown(s);
    if (y<=mid) insert(s*2,l,mid,x,y,z);
    else if (x>mid) insert(s*2+1,mid+1,r,x,y,z);
    else insert(s*2,l,mid,x,mid,z),insert(s*2+1,mid+1,r,mid+1,y,z);
}

void solveinsert(int x,int y,int z)
{
    for (;belong[x]!=belong[y];x=father[belong[x]][0])
      {
        if (deep[belong[x]]y]]) swap(x,y);
        insert(1,1,n,pos[belong[x]],pos[x],z);
      }
    if (deep[x]y]) swap(x,y);
    insert(1,1,n,pos[y],pos[x],z);
}

int main()
{
    int i,x,y;
    n=read();
    for (i=1;i<=n;i++) a[i]=read();
    for (i=1;ix=read(),y=read(),add(x,y),add(y,x);
    dfs1(1),dfs2(1,1),build(1,1,n);
    for (i=2;i<=n;i++) solveinsert(a[i-1],a[i],1);
    for (i=1;i<=n;i++)
      if (i!=a[1]) printf("%d\n",query(1,1,n,pos[i])-1);
      else printf("%d\n",query(1,1,n,pos[i]));
    return 0;
}

你可能感兴趣的:(线段树,树链剖分)