洛谷1361 小M的作物(最小割)

传送门

【题目分析】

第一眼看有A,B两个以为是二分图,根据题目各种想办法建边,感觉还是弄不成。

推了好久才发现求的就是一个最小割。

首先如果令m=0,那么可以建这样一个图:

当然贪心就可以了,最后选择的方案也一定是取ai和bi中的较大值,算出来也正是这张图的最小割,也是\Sigma^{n}_{i=1} a_{i}+b_{i}-最大流。

那如果m不为0呢?

我们举一个有两个点的例子:

a表示全部选种A地,b表示全部选种B地。

将a,b到两个点的容量设为INF,所以做最小割的时候一定不会去割INF的边,所以这个图的最小割就是这两种作物不选用的最小值,用总值减去即可得到最大贡献。

所以以此类推即可。

所以就是一道最小割嘛。

【代码~】

#include<bits/stdc++.h>
using namespace std;
const int MAXN=2010;
const int MAXM=1e6+10;
const int INF=0x3f3f3f3f;

int n,m,cnt,s,t,sum;
int head[MAXN],depth[MAXN],cur[MAXN];
int nxt[MAXM],to[MAXM],w[MAXM];
int a[MAXN],b[MAXN];

int Read(){
	int i=0,f=1;
	char c;
	for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void Add(int x,int y,int z){
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
	w[cnt]=z;
	cnt++;
}

void add(int x,int y,int z){
	Add(x,y,z);
	Add(y,x,0);
}

int bfs(){
	memset(depth,0,sizeof(depth));
	depth[s]=1;
	queue<int> q;
	q.push(s);
	while(!q.empty()){
		int u=q.front();
		q.pop();
		for(int i=head[u];i!=-1;i=nxt[i]){
			int v=to[i];
			if(!depth[v]&&w[i]){
				depth[v]=depth[u]+1;
				q.push(v);
				if(v==t)
				  return 1;
			}
		}
	}
	return 0;
}

int dfs(int u,int flow){
	if(u==t)
	  return flow;
	for(int &i=cur[u];i!=-1;i=nxt[i]){
		int v=to[i];
		if(depth[v]==depth[u]+1&&w[i]){
			int di=dfs(v,min(w[i],flow));
			if(di){
				w[i]-=di;
				w[i^1]+=di;
				return di;
			}
		}
	}
	return 0;
}

int dinic(){
	int ans=0;
	while(bfs()){
		for(int i=s;i<=t;++i)
		  cur[i]=head[i];
		while(int d=dfs(s,INF))
		  ans+=d;
	}
	return ans;
}

int main(){
	memset(head,-1,sizeof(head));
	n=Read();
	for(int i=1;i<=n;++i){
		a[i]=Read();
		sum+=a[i];
	}
	for(int i=1;i<=n;++i){
		b[i]=Read();
		sum+=b[i];
	}
	m=Read();
	s=0,t=n+m+m+1;
	for(int i=1;i<=n;++i){
		add(s,i,a[i]);
		add(i,t,b[i]);
	}
	for(int i=1;i<=m;++i){
		int k=Read();
		int c1=Read(),c2=Read();
		sum+=c1+c2;
		add(s,i+n,c1);
		add(i+n+m,t,c2);
		for(int j=1;j<=k;++j){
			int x=Read();
			add(i+n,x,INF);
			add(x,i+n+m,INF);
		}
	}
	cout<<sum-dinic();
	return 0;
}

 

posted @ 2018-12-07 21:02  Ishtar~  阅读(133)  评论(0编辑  收藏  举报