bzoj 3681: Arietta

题意：

所有的 n 个音符形成一棵由音符 C ( 1 号节点) 构成的有根树,每一个音符有一个音高 Hi 。
Arietta 有 m 个力度,第 i 个力度能弹出 Di 节点的子树中,音高在 [Li,Ri] 中的任意一个音符。
为了乐曲的和谐,Arietta 最多会弹奏第 i 个力度 Ti 次。
Arietta 想知道她最多能弹出多少个音符。

题解：

显然一种暴力的做法是每个力度向它可以弹出的点连边，然后跑网络流，这样时空都会爆炸。
一个点连向一段区间的建图可以考虑线段树优化建图，但这是在一棵树上，所以要建主席树。
那么就dsu+主席树优化建图跑网络流即可。
code：

#include<queue>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;
const int inf=1<<28;
struct Node{
    int y,c,next,other;
}a[2300010];int len=0,last[630010];
int H[630010],st,ed;
struct lr{
    int l,r,t;
    lr() {}
    lr(int _l,int _r,int _t) {l=_l;r=_r;t=_t;}
};
vector<lr> v[10010];
struct trnode{
    int lc,rc;
}tr[630010];int tot=0,root[10010];
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void insert(int x,int y)
{
    a[++len].y=y;
    a[len].next=last[x];last[x]=len;
}
void ins(int x,int y,int c)
{
    a[++len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;
    a[++len].y=x;a[len].c=0;
    a[len].next=last[y];last[y]=len;
    a[len].other=len-1;a[len-1].other=len;
}
queue<int> q;
bool bt_h()
{
    memset(H,0,sizeof(H));
    while(!q.empty()) q.pop();
    q.push(st);H[st]=1;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=last[x];i;i=a[i].next)
        {
            int y=a[i].y;
            if(H[y]==0&&a[i].c>0) H[y]=H[x]+1,q.push(y);
        }
    }
    return H[ed]!=0;
}
int findflow(int x,int f)
{
    if(x==ed) return f;
    int t,ans=0;
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        if(H[y]==H[x]+1&&a[i].c>0&&ans<f)
        {
            ans+=(t=findflow(y,min(a[i].c,f-ans)));
            a[i].c-=t;a[a[i].other].c+=t;
        }
    }
    if(ans==0) H[x]=0;
    return ans;
}
int n,m,h[10010];
int size[10010],son[10010],Son;
void pre_node(int x)
{
    size[x]=1;
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        pre_node(y);size[x]+=size[y];
        if(size[son[x]]<size[y]) son[x]=y;
    }
}
void update(int &x,int fr,int l,int r,int k,int c)
{
    x=++tot;tr[x]=tr[fr];
    if(fr!=0) ins(x,fr,inf);
    if(l==r) {ins(x,c+n,inf);return;}
    int mid=(l+r)/2;
    if(k<=mid) update(tr[x].lc,tr[fr].lc,l,mid,k,c),ins(x,tr[x].lc,inf);
    else update(tr[x].rc,tr[fr].rc,mid+1,r,k,c),ins(x,tr[x].rc,inf);
}
void solve(int x,int &rt)
{
    update(rt,rt,1,n,h[x],x);
    for(int i=last[x];i;i=a[i].next)
        if(a[i].y!=Son) solve(a[i].y,rt);
}
void make(int x,int l,int r,int fl,int fr)
{
    if(l==fl&&r==fr) {ins(tot,x,inf);return;}
    int mid=(l+r)/2;
    if(fr<=mid) make(tr[x].lc,l,mid,fl,fr);
    else if(fl>mid) make(tr[x].rc,mid+1,r,fl,fr);
    else make(tr[x].lc,l,mid,fl,mid),make(tr[x].rc,mid+1,r,mid+1,fr);
}
void build(int x,int rt)
{
    for(int i=0;i!=v[x].size();i++)
    {
        int l=v[x][i].l,r=v[x][i].r,t=v[x][i].t;
        ins(st,++tot,t);
        make(rt,1,n,l,r);
    }
}
void dfs(int x)
{
    for(int i=last[x];i;i=a[i].next)
        if(a[i].y!=son[x]) dfs(a[i].y);
    if(son[x]) dfs(son[x]),Son=son[x];
    root[x]=root[son[x]];
    solve(x,root[x]);Son=0;
    build(x,root[x]);
}
int main()
{
    n=read();m=read();
    for(int i=2;i<=n;i++)
    {
        int fa;fa=read();
        insert(fa,i);
    }
    for(int i=1;i<=n;i++) h[i]=read();
    for(int i=1;i<=m;i++)
    {
        int l,r,x,t;l=read();r=read();x=read();t=read();
        v[x].push_back(lr(l,r,t));
    }
    tot=n+n+2;st=n+n+1;ed=n+n+2;
    for(int i=1;i<=n;i++) ins(n+i,ed,1);
    pre_node(1);dfs(1);
    int ans=0;
    while(bt_h()) ans+=findflow(st,inf);
    printf("%d",ans);
}