BZOJ 3438 小M的作物

Description
小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
只有一行，包括一个整数，表示最大收益
Sample Input
3

4 2 1

2 3 2

1

2 3 2 1 2
Sample Output
11

样例解释A耕地种1，2，B耕地种3，收益4+2+3+2=11。

1<=k< n<= 1000,0 < m < = 1000 保证所有数据及结果不超过2*10^9。
HINT

Source
Kpmcup#0 By Greens

吼题哇！
窝只想到了最小割，对于这个集合究竟如何建边一脸懵逼。
首先从s向i(1<=i<=n)连容量为 ai 的边，再从i向t连容量为 bi 的边。
之后想到最大权闭合子图，还有二分图的思想。
将集合拆点，将s向I连容量为 c1i 的边，向属于集合的点连容量为inf的边。对于另一半同理。累加总收益，减去最小割就好了。

#include
#include
#include
#include
using namespace std;
const int N=1005;
const int inf=1e9+7;
int n,m,s,t,cnt=1,ans,hd[3*N],a[N],b[N],d[3*N];
queue<int>q;
struct edge
{
    int to,nxt,f;
}v[2*N*N];
void addedge(int x,int y,int z)
{
    v[++cnt].to=y,v[cnt].f=z;
    v[cnt].nxt=hd[x],hd[x]=cnt;
}
void addedges(int x,int y,int z)
{
    addedge(x,y,z),addedge(y,x,0);
}
bool bfs()
{
    memset(d,0,sizeof(d));
    q.push(s);
    d[s]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=hd[u];i;i=v[i].nxt)
            if(v[i].f&&!d[v[i].to])
            {
                d[v[i].to]=d[u]+1;
                q.push(v[i].to);
            }
    }
    return d[t];
}
int dfs(int u,int lft)
{
    if(u==t||lft==0)
        return lft;
    int r=lft;
    for(int i=hd[u];i;i=v[i].nxt)
        if(r&&v[i].f&&d[v[i].to]==d[u]+1)
        {
            int w=dfs(v[i].to,min(r,v[i].f));
            v[i].f-=w,v[i^1].f+=w,r-=w;
            if(!r)
                return lft;
        }
    if(r==lft)
        d[u]=0;
    return lft-r;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)
        scanf("%d",&b[i]);
    scanf("%d",&m);
    s=0,t=n+2*m+1;
    for(int i=1;i<=n;i++)
    {
        ans+=a[i]+b[i];
        addedges(s,i,a[i]),addedges(i,t,b[i]);
    }
    for(int i=1;i<=m;i++)
    {
        int x,y,z,w;
        scanf("%d%d%d",&x,&y,&z);
        addedges(s,n+i,y),addedges(n+m+i,t,z);
        ans+=y+z;
        while(x--)
        {
            scanf("%d",&w);
            addedges(n+i,w,inf),addedges(w,n+m+i,inf);
        }
    }
    while(bfs())
        ans-=dfs(s,inf);
    printf("%d\n",ans);
    return 0;
}