[ 网络流 树上启发式合并 主席树 ] BZOJ3681
来源:互联网 发布:同花顺软件开发怎么样 编辑:程序博客网 时间:2024/05/17 06:30
题解
#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 1e9vector<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;}
阅读全文
0 0
- [ 网络流 树上启发式合并 主席树 ] BZOJ3681
- 【BZOJ3681】Arietta,主席树优化网络流
- BZOJ 3123: [Sdoi2013]森林 启发式合并 树上主席树
- [bzoj3123][洛谷P3302] [SDOI2013]森林(树上主席树+倍增lca+启发式合并)
- 树上启发式合并
- bzoj3123 森林 主席树&启发式合并
- 【bzoj3123】【sdoi2013】【森林】【启发式合并+主席树】
- BZOJ 3123 主席树 启发式合并
- BZOJ 3123 森林 主席树启发式合并
- bzoj 3123(主席树+启发式合并)
- BZOJ 3123 LCA + 主席树 + 启发式合并
- 【算法】树上启发式合并算法
- [ 树上启发式合并 线段树 单调栈 ] BZOJ5040
- BZOJ 3123 [Sdoi2013]森林 主席树+启发式合并
- BZOJ 3123: [Sdoi2013]森林|主席树|启发式合并
- 3123: [Sdoi2013]森林 主席树+启发式合并+LCA
- [主席树 启发式合并] BZOJ 3123 [Sdoi2013]森林
- [BZOJ3123][Sdoi2013]森林(主席树启发式合并)
- Java JVM:垃圾回收(GC 在什么时候,对什么东西,做了什么事情)
- 运维生涯-开始篇(2)
- Mybatis
- JavaScript 事件类型
- 请使用socket相关函数(非curl)实现如下功能:构造一个post请求,发送到指定httpserver的指定端口的指定请求路径(如http://www.example.com:8080/test)
- [ 网络流 树上启发式合并 主席树 ] BZOJ3681
- 【1701H1】【穆晨】【171011】连续第一天总结
- 贝叶斯分类——朴素贝叶斯算法
- JVM堆与栈
- 洛谷 1273 有线电视网
- 数据库
- 蓝桥杯——剪邮票
- maven私有库的安装、配置和使用
- 作业(10/1-10/3)