BZOJ 3438: 小M的作物 最大权闭合子图

3438: 小M的作物

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 1403  Solved: 607
[Submit][Status][Discuss]

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。

最大权闭合子图

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;

typedef long long ll;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=3010,inf=0X3f3f3f3f;

int n,D,F,S=0,T=N-1;

int ecnt=1,last[N];
struct EDGE{int to,nt,val;}e[N*N];
inline void readd(int u,int v,int val)
{e[++ecnt]=(EDGE){v,last[u],val};last[u]=ecnt;}
inline void add(int u,int v,int val)
{readd(u,v,val);readd(v,u,0);}

int d[N],q[N];

bool bfs()
{
	memset(d,0,sizeof(d));
	register int head=0,tail=1,u,i;
	q[head]=0;d[0]=1;
	while(head<tail)
	{
		u=q[head++];
		for(i=last[u];i;i=e[i].nt)if(e[i].val&&!d[e[i].to])
		{d[e[i].to]=d[u]+1;q[tail++]=e[i].to;}
	}
	return d[T];
}

int dfs(int u,int lim)
{
	if(u==T|!lim)return lim;
	register int i,res=0,tmp;
	for(i=last[u];i;i=e[i].nt)
	if(e[i].val&&d[u]+1==d[e[i].to])
	{
		tmp=dfs(e[i].to,min(e[i].val,lim));
		lim-=tmp;res+=tmp;e[i].val-=tmp;e[i^1].val+=tmp;
		if(!tmp)d[e[i].to]=-1;if(!lim)break;
	}
	return res;
}

ll flow;

void dinic()
{while(bfs())flow+=dfs(S,inf);}

int main()
{
	n=read();
	register int u,val,c1,c2,k,i;ll ans=0;
	for(i=1;i<=n;++i){val=read();add(S,i,val);ans+=val;}
	for(i=1;i<=n;++i){val=read();add(i,T,val);ans+=val;}
	int m=read();
	for(i=1;i<=m;++i)
	{
		k=read();c1=read();c2=read();ans+=c1+c2;
		add(i+m+n,T,c2);add(S,i+n,c1);
		while(k--){u=read();add(u,i+m+n,inf);add(i+n,u,inf);}
	}
	dinic();
	print(int(ans-flow));puts("");
	return 0;
}
/*
3
4 2 1
2 3 2
1
2 3 2 1 2

11
*/