upc--小赵的盆栽（dfs序+线段树）

题目描述

小赵和小静逛街回来之后在路边看到了一家花店。钢铁小赵10分钟后买了一个绿油油盆栽送给了小静。回到了学校，小赵自顾很细心的照顾盆栽。由于没给小静买花，小静有些郁闷，所以小静决定给小赵制造一点麻烦。

为方便理解，这棵盆栽呈现出树的结构，每个分叉为一个节点，树枝为边。而且每个分叉上都有一片叶子，叶子刚买回来长度为0。为了方便，小赵给每个节点都进行了标号，根节点标号为1。

小赵每天会对所有节点喷上AC药剂，喷上AC药剂叶子会长长1mm。但小静每天会在小赵喷过药剂之后，选择对盆栽的一树枝上所有节点喷上WA药剂，WA药剂会让小赵今天喷的AC药剂无效。

m天后，小赵发现了这件事情，儒雅随和的小赵并没有生气，但是他想知道现在每个节点上叶子都长了多长。

输入

本题有多组测试数据，处理到文件结束。

每组数据的第一行输入两个整数，n代表盆栽的节点数，m代表小赵过了多少天才发现小静的恶作剧。第二行是n-1个数字，代表2~n节点的父亲节点编号，第三行是m个数字，代表小静每天选择的树枝的根节点标号。（n<=1e5, m<=1e5，节点编号为1~n）

输出

每组数据输出一行n个数字用空格隔开，代表每个叶子的长度。

样例输入

4 3
1 1 2
1 4 2
 

样例输出

2 1 2 0

思路

先用dfs把树转换成区间，接下来就变成了用线段树解决区间修改和单点查询的问题
（yysy,这算是真正意义上自己解决的第一道dfs序+线段树的问题，感觉值得纪念一下（o(￣▽￣)ブ）（其实最后被小于号卡了半小时，还是菜了/(ㄒoㄒ)/~~））

代码

#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<algorithm>
#include<iomanip>
#include<cmath>
#include<fstream>
#define X first
#define Y second
#define best 131 
#define INF 0x3f3f3f3f
#define P pair<int,int>
#define ls p<<1
#define rs p<<1|1
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const double eps=1e-5;
const double pai=acos(-1.0);
const int N=1e5+10;
const int maxn=110;
const int mod=1000000007;
int head[N],cnt,tot,in[N],out[N];
struct Node
{
    int to;
    int next;
}tr[N<<2];
struct node
{
    int l;
    int r;
    int sum;
    int lazy;
}t[N<<2];
void add(int x,int y)
{
    tr[cnt].to=y;
    tr[cnt].next=head[x];
    head[x]=cnt++;
}
void dfs(int x,int y)
{
    in[x]=++tot;
    for(int i=head[x];i!=-1;i=tr[i].next)
    {
        int end=tr[i].to;
        if(end==y) continue;
        dfs(end,x);
    }
    out[x]=tot;
}
void pushup(int p)
{
    t[p].sum=t[ls].sum+t[rs].sum;
}
void pushdown(int p)
{
    if(t[p].lazy)
    {
        t[ls].sum+=t[p].lazy*(t[ls].r-t[ls].l+1);
        t[rs].sum+=t[p].lazy*(t[rs].r-t[rs].l+1);
        t[ls].lazy+=t[p].lazy;
        t[rs].lazy+=t[p].lazy;
        t[p].lazy=0;
    }
}
void build(int p,int l,int r)
{
    t[p].l=l,t[p].r=r,t[p].lazy=0;
    if(l==r)
    {
        t[p].sum=0;
        return ;
    }
    int mid=(l+r)>>1;
    build(ls,l,mid);
    build(rs,mid+1,r);
    pushup(p);
}
void change(int p,int l,int r,int d)
{
    if(l<=t[p].l&&t[p].r<=r)
    {
        t[p].sum+=d*(t[p].r-t[p].l+1);
        t[p].lazy+=d;
        return ;
    }
    pushdown(p);
    int mid=(t[p].r+t[p].l)>>1;
    if(l<=mid) change(ls,l,r,d);
    if(r>mid) change(rs,l,r,d);
    pushup(p);
}
int ask(int p,int l)
{
    if(t[p].l==t[p].r)
    {
        return t[p].sum;
    }
    pushdown(p);
    int mid=(t[p].r+t[p].l)>>1;
    int ans;
    if(l<=mid) ans=ask(ls,l);
    else ans=ask(rs,l);
    return ans;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        cnt=tot=0;
        memset(head,-1,sizeof(head));
        memset(in,0,sizeof(in));
        memset(out,0,sizeof(out));
        for(int i=2;i<=n;i++)
        {
            int x;
            scanf("%d",&x);
            add(x,i);
            add(i,x);
        }
        dfs(1,0);
        build(1,1,n);
        for(int i=1;i<=m;i++)
        {
            int x;
            scanf("%d",&x);
            change(1,in[x],out[x],1);
        }
        for(int i=1;i<=n;i++)
        {
            printf("%d ",(m-(ask(1,in[i]))<0?0:m-(ask(1,in[i]))));
        }
        puts("");
    }
    return 0;
}