BZOJ3438 小M的作物

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=3438

Description

 背景

    小M还是个特么喜欢玩MC的孩纸。。。

 描述

    小M在MC里开辟了两块巨大的耕地A和B（你可以认为容量是无穷），现在，小P有n中作物的种子，每种作物的种子有1个（就是可以种一棵作物）（用1...n编号），现在，第i种作物种植在A中种植可以获得ai的收益，在B中种植可以获得bi的收益，而且，现在还有这么一种神奇的现象，就是某些作物共同种在一块耕地中可以获得额外的收益，小M找到了规则中共有m种作物组合，第i个组合中的作物共同种在A中可以获得c1i的额外收益，共同总在B中可以获得c2i的额外收益，所以，小M很快的算出了种植的最大收益，但是他想要考考你，你能回答他这个问题么？

Input

    第一行包括一个整数n

    第二行包括n个整数，表示ai

    第三行包括n个整数，表示bi

    第四行包括一个整数m

    接下来m行，对于接下来的第i行：第一个整数ki，表示第i个作物组合中共有ki种作物，接下来两个整数c1i，c2i，接下来ki个整数，表示该组合中的作物编号。

Output

   只有一行，包括一个整数，表示最大收益

Orz Greens

Orz KPM

放官方题解跑

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define rep(i,l,r) for(int i=l; i<=r; i++)
#define clr(x,y) memset(x,y,sizeof(x))
#define travel(x) for(Edge *p=last[x]; p; p=p->pre)
using namespace std;
const int INF = 0x7fffffff;
const int maxn = 3010;
inline int read(){
    int ans = 0, f = 1;
    char c = getchar();
    for(; !isdigit(c); c = getchar())
    if (c == '-') f = -1;
    for(; isdigit(c); c = getchar())
    ans = ans * 10 + c - '0';
    return ans * f;
}
struct Edge{
    Edge *pre,*rev; int to,cost;
}edge[4010000],*last[maxn],*cur[maxn],*pt = edge;
int n,m,x,y,z,k,c1,c2,cnt,ans=0,d[maxn];
queue <int> q;
inline void addedge(int x,int y,int z){
    pt->pre = last[x]; pt->to = y; pt->cost = z; last[x] = pt++;
}
inline void add(int x,int y,int z){
    addedge(x,y,z); addedge(y,x,0); last[x]->rev = last[y]; last[y]->rev = last[x];
}
bool bfs(){
    while (!q.empty()) q.pop();
    clr(d,-1); d[0] = 0; q.push(0);
    while (!q.empty()){
        int now = q.front(); q.pop();
        travel(now){
            if (d[p->to] == -1 && p->cost > 0){
                d[p->to] = d[now] + 1;
                q.push(p->to);
                if (p->to == n + 1) return 1;
            }
        }
    }
    return 0;
}
int dfs(int x,int flow){
    if (x == n + 1 || (!flow)) return flow; int w = 0;
    for (Edge *p = cur[x]; p && w < flow; p = p->pre){
        if (d[p->to] == d[x] + 1 && p->cost > 0){
            int delta = dfs(p->to,min(flow-w,p->cost));
            p->cost -= delta; p->rev->cost += delta; w += delta;
            if (p->cost) cur[x] = p;
        }
    }
    if (w < flow) d[x] = -1;
    return w;
}
int main(){
    n = read(); clr(last,0); cnt = n + 2;
    rep(i,1,n) x = read(), ans += x, add(0,i,x);
    rep(i,1,n) x = read(), ans += x, add(i,n+1,x);
    m = read();
    rep(i,1,m){
        k = read(); c1 = read(); c2 = read();
        ans += c1 + c2;
        int u = cnt++, v = cnt++;
        add(0,u,c1); add(v,n+1,c2);
        rep(j,1,k){
            x = read();
            add(u,x,INF); add(x,v,INF);
        }
    }
    int mxflow = 0;
    while (bfs()){
        rep(i,0,cnt-1) cur[i] = last[i];
        mxflow += dfs(0,INF);
    }
    printf("%d
",ans-mxflow);
    return 0;
}