ZOJ 3910 Market

Market

---

Time Limit: 2 Seconds      Memory Limit: 65536 KB

---

There's a fruit market in Byteland. The salesmen there only sell apples.

There are n salesmen in the fruit market and the i-th salesman will sell at most w_i apples. Every salesman has an immediate manager p_i except one salesman who is the boss of the market. A salesmanA is said to be the superior of another salesman B if at least one of the followings is true:

• Salesman A is the immediate manager of salesman B.
• Salesman B has an immediate manager salesman C such that salesman A is the superior of salesman C.

The market will not have a managerial cycle. That is, there will not exist a salesman who is the superior of his/her own immediate manager.

We will call salesman x a subordinate of another salesman y, if either y is an immediate manager of x, or the immediate manager of x is a subordinate to salesman y. In particular, subordinates of the boss are all other salesmen of the market. Let the degree of the boss be 0. Then if the degree of i-th salesman is k, the immediate subordinates of i-th salesman will have degree k + 1.

Today, m buyers come to market for apples. The i-th buyer will buy at most c_i apples only from the x_i-th salesman and his subordinates whose degree is no larger than x_i-th salesman's degree plus d_i.

The boss wants to know how many apples can be sold in salesmen's best effort (i.e. the maximum number).

Input

There are multiple test cases. The first line of input contains an integer T indicating the number of test cases. For each test case:

The first line contains two integers n and m (1 ≤ n, m ≤ 10000) — the number of salesmen and the number of buyers.

The second line contains n integers w₁, w₂, ..., w_n (1 ≤ w_i ≤ 10⁵). Every w_i denotes the number of apples that i-th salesman can sell.

The next line contains n integers p_i (1 ≤ p_i ≤ n or p_i = -1). Every p_i denotes the immediate manager for the i-th salesman. If p_i is -1, that means that the i-th salesman does not have an immediate manager.

Each of the next m lines contains three integers c_i, x_i and d_i (1 ≤ c_i ≤ 10⁵, 1 ≤ x_i ≤ n, 0 ≤ d_i ≤ n) — the information of i-th buyer.

It is guaranteed that the total number of salesmen in the input doesn't exceed 10⁵, and the total number of buyers also doesn't exceed 10⁵. The number of test cases in the input doesn't exceed 500.

Output

For each test case, output a single integer denoting the maximum number of apples can be sold.

Sample Input

1
4 2
1 2 3 4
-1 1 2 3
3 2 1
5 1 1

Sample Output

6

---

Author: LIN, Xi
Source: ZOJ Monthly, October 2015

解题：网络流，关键是，由于点多，导致边更多，直接建图，爆内存。。。

感谢Claris的帮助，在花费了一两天的时间，终于搞定了

 #include <bits/stdc++.h>
 using namespace std;
 const int N = ,M = ,INF = 0x3f3f3f3f;
 int n,m;
 struct arc{
     int to,flow,next;
     arc(int x = ,int y = ,int z = -){
         to = x;
         flow = y;
         next = z;
     }
 }e[];
 int cur[M],h[M],gap[M],head[M],S,T,tot;
 void addedge(int u,int v,int flow){
     if(!u || !v) return;
     e[tot] = arc(v,flow,head[u]);
     head[u] = tot++;
     e[tot] = arc(u,,head[v]);
     head[v] = tot++;
 }
 int sap(int u,int low){
     if(u == T) return low;
     int ret = ;
     for(int &i = cur[u]; ~i; i = e[i].next){
         if(e[i].flow && h[u] == h[e[i].to] + ){
             int a = sap(e[i].to,min(low,e[i].flow));
             e[i].flow -= a;
             e[i^].flow += a;
             low -= a;
             ret += a;
             if(!low) return ret;
         }
     }
     if(!(--gap[h[u]])) h[S] = T;
     gap[++h[u]]++;
     cur[u] = head[u];
     return ret;
 }
 namespace REDUCE {
 int cnt,tot,leaf[N],dep[N],root[N],l[M],r[M],head[N];
 int hd[N],num;
 struct arc {
     int to,next;
     arc(int x = ,int y = -) {
         to = x;
         next = y;
     }
 } e[N];
 struct QU {
     int to,L,R,next;
     QU(int to = ,int L = ,int R = ,int nxt = -) {
         this->to = to;
         this->L = L;
         this->R = R;
         this->next = nxt;
     }
 } Q[N];
 void add(int u,int v) {
     e[tot] = arc(v,head[u]);
     head[u] = tot++;
 }
 void addask(int u,int v,int L,int R) {
     Q[num] = QU(v,L,R,hd[u]);
     hd[u] = num++;
 }
 int merge(int x,int y,int L,int R) {
     if(!x) return y;
     if(!y) return x;
     int z = ++cnt;
     if(L == R) {
         addedge(z,x,INF);
         addedge(z,y,INF);
         return z;
     }
     int mid = (L + R)>>;
     addedge(z,l[z] = merge(l[x],l[y],L,mid),INF);
     addedge(z,r[z] = merge(r[x],r[y],mid + ,R),INF);
     return z;
 }
 int build(int L,int R,int pos) {
     int x = ++cnt;
     if(L == R) return x;
     int mid = (L + R)>>;
     if(pos <= mid) addedge(x,l[x] = build(L,mid,pos),INF);
     if(pos > mid) addedge(x,r[x] = build(mid+,R,pos),INF);
     return x;
 }
 void ask(int root,int L,int R,int lt,int rt,int node) {
     if(!root) return;
     if(lt <= L && rt >= R) {
         addedge(node,root,INF);
         return;
     }
     int mid = (L + R)>>;
     if(lt <= mid && l[root]) ask(l[root],L,mid,lt,rt,node);
     if(rt > mid && r[root]) ask(r[root],mid + ,R,lt,rt,node);
 }
 void dfs(int u,int depth) {
     dep[u] = depth;
     root[u] = build(,n,dep[u]);
     leaf[u] = cnt;
     for(int i = head[u]; ~i; i = e[i].next)
         dfs(e[i].to,depth + );
 }
 void dfs(int u) {
     for(int i = head[u]; ~i; i = e[i].next) {
         dfs(e[i].to);
         root[u] = merge(root[u],root[e[i].to],,n);
     }
     for(int i = hd[u]; ~i; i = Q[i].next)
         ask(root[u],,n,Q[i].L,min(n,Q[i].R),Q[i].to);
 }
 void init() {
     memset(head,-,sizeof head);
     num = tot = cnt = ;
     memset(hd,-,sizeof hd);
     memset(l,,sizeof l);
     memset(r,,sizeof r);
 }
 }
 int have[N],need[N],hs[N];
 int main() {
     int kase,u,v,w,rt;
     scanf("%d",&kase);
     memset(head,-,sizeof head);
     while(kase--) {
         scanf("%d%d",&n,&m);
 
         tot = ;
         REDUCE::init();
         for(int i = ; i <= n; ++i)
             scanf("%d",have + i);
         for(int i = ; i <= n; ++i) {
             scanf("%d",&u);
             if(~u) REDUCE::add(u,i);
             else rt = i;
         }
         REDUCE::dfs(rt,);
         for(int i = ; i <= m; ++i) {
             scanf("%d%d%d",need + i,&u,&v);
             REDUCE::addask(u,hs[i] = ++REDUCE::cnt,REDUCE::dep[u],v + REDUCE::dep[u]);
         }
         REDUCE::dfs(rt);
         S = ++REDUCE::cnt;
         T = ++REDUCE::cnt;
         for(int i = ; i <= m; ++i)
             addedge(S,hs[i],need[i]);
         for(int i = ; i <= n; ++i)
             addedge(REDUCE::leaf[i],T,have[i]);
         int maxflow = ;
         gap[] = T;
         while(h[S] < T) maxflow += sap(S,INF);
         printf("%d\n",maxflow);
         for(int i = ; i <= T; i++) {
             h[i] = gap[i]=;
             head[i] = -;
         }
     }
     return ;
 }