[ 网络流 树上启发式合并 主席树 ] BZOJ3681

题解

#include<bits/stdc++.h>
using namespace std;
inline char nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void Read(int& x){
    char c=nc();
    for(;c<'0'||c>'9';c=nc());
    for(x=0;c>='0'&&c<='9';x=(x<<3)+(x<<1)+c-48,c=nc());
}
#define N 10010
#define M 1300010
#define INF 1e9
vector<int>g[N];
struct Edge{
    int t,nx,c;
}e[2600010];
bool b[M];
int s,t,h[M],num=1,cur[M],d[M];
int Rt[N],cl[M],cr[M],cnt;
int i,j,k,n,m,p,a[N];
int x,y,z,sz[N];
inline void Add(int x,int y,int w){
    e[++num].t=y;e[num].c=w;e[num].nx=h[x];h[x]=num;
    e[++num].t=x;e[num].c=0;e[num].nx=h[y];h[y]=num;
}
inline void Update(int& x,int y,int f,int l,int r,int z){
    x=++cnt;cl[x]=cl[y];cr[x]=cr[y];
    if(f)Add(f,x,INF);
    if(y)Add(x,y,INF);
    if(l==r){
        Add(x,z+m,INF);return;
    }
    int Mid=l+r>>1;
    if(a[z]<=Mid)Update(cl[x],cl[y],x,l,Mid,z);else Update(cr[x],cr[y],x,Mid+1,r,z);
}
inline void Calc(int x,int& y){
    Update(y,y,0,1,n,x);
    for(int i=0;i<g[x].size();i++)Calc(g[x][i],y);
}
inline void Dfs(int x){
    sz[x]=1;int S=0;
    for(int i=0;i<g[x].size();i++){
        int v=g[x][i];
        Dfs(v);
        sz[x]+=sz[v];
        if(sz[v]>sz[S])S=v;
    }
    Update(Rt[x],Rt[S],0,1,n,x);
    for(int i=0;i<g[x].size();i++){
        int v=g[x][i];
        if(v==S)continue;
        Calc(v,Rt[x]);
    }
}
inline void Query(int x,int l,int r,int L,int R,int y){
    if(!x||l>R||r<L)return;
    if(l>=L&&r<=R){
        Add(y,x,INF);
        return;
    }
    int Mid=l+r>>1;
    Query(cl[x],l,Mid,L,R,y);Query(cr[x],Mid+1,r,L,R,y);
}
int q[M],l,r;
inline bool Bfs(){
    for(int i=1;i<=cnt;i++)b[i]=0;
    l=0;q[r=1]=s;
    while(++l<=r){
        int x=q[l];
        for(int i=h[x];i;i=e[i].nx)
        if(e[i].c&&!b[e[i].t]){
            b[e[i].t]=1;d[e[i].t]=d[x]+1;q[++r]=e[i].t;
        }
    }
    return b[t];
}
inline int Dfs(int x,int a){
    if(!a||x==t)return a;
    int fl=0;
    for(int& i=cur[x];i;i=e[i].nx){
        int v=e[i].t,f;
        if(d[v]==d[x]+1&&(f=Dfs(v,min(a,e[i].c)))){
            fl+=f;
            e[i].c-=f;e[i^1].c+=f;
            if(a==fl)break;
        }
    }
    if(!fl)d[x]=-1;
    return fl;
}
inline int Dinic(){
    int Ans=0;b[s]=1;
    while(Bfs()){
        for(int i=0;i<=cnt;i++)cur[i]=h[i];
        Ans+=Dfs(s,INF);
    }
    return Ans;
}
int main(){
    Read(n);Read(m);
    for(i=2;i<=n;i++)Read(x),g[x].push_back(i);
    for(i=1;i<=n;i++)Read(a[i]);
    cnt=n+m;
    Dfs(1);
    for(i=1;i<=m;i++)Read(x),Read(y),Read(z),Read(k),Add(0,i,k),Query(Rt[z],1,n,x,y,i);
    t=++cnt;
    for(i=1;i<=n;i++)Add(i+m,t,1);
    printf("%d\n",Dinic());
    return 0;
}