P2762 太空飞行计划问题

题目

题目

思路

显然网络瘤。
建超级源，超级汇。
超级源连线实验，流量为获利。
实验连器材，流量inf。
器材连超级汇，流量为费用。
考虑最小割。
则总收益-最大流为答案。
code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<cmath>
#include<queue>
#include<deque>
using namespace std;
//When I wrote this code,God and I unterstood what was I doing 
long long n1,n2,n3,m,s,t,x,y,tot=2,head[41001],dep[41001],l,r,u[41001],r2;
long long w,ans;
long long mn(long long x,long long y)
{
	return (x>y?y:x);
}
struct f{
	long long to,net;
	long long w;
} a[1000011];
void add(long long x,long long y,long long w)
{
	a[tot].to=y,a[tot].w=w,a[tot].net=head[x],head[x]=tot++;
	a[tot].to=x,a[tot].w=0,a[tot].net=head[y],head[y]=tot++;
	return;
}
bool bfs()
{
	memset(dep,0,sizeof(dep));
	l=0,r=0;
	u[r++]=s;
	dep[s]=1;
	while (l<r)
	{
		r2=r;
		for (long long i=l;i<r;i++)
		{
			for (long long j=head[u[i]];j;j=a[j].net)
			{
				if (a[j].w!=0&&dep[a[j].to]==0)
				{
					u[r2++]=a[j].to;
					dep[a[j].to]=dep[u[i]]+1;
				}
			}
		}
		l=r,r=r2;
	}
	return dep[t];
}
long long dfs(long long d,long long in)
{
	if (d==t) return in;
	long long out=0;
	long long uw=dep[d]+1;
	for (long long j=head[d];j&&in;j=a[j].net)
	{
		if (a[j].w==0||dep[a[j].to]!=uw) continue;
		long long s=dfs(a[j].to,mn(in,a[j].w));
		out+=s,in-=s;
		a[j].w-=s,a[j^1].w+=s;
	}
	if (out==0)
	{
		dep[d]=0;
	}
	return out;
}
char tools[10000];
int main()
{
	int sum=0;
	scanf("%lld%lld",&n1,&n2);
	s=n1+n2+1,t=n1+n2+2;
	for (int i=1;i<=n1;i++)
	{
		scanf("%lld",&x);
		add(s,i,x);
		sum+=x;
		memset(tools,0,sizeof(tools));
		cin.getline(tools,10000);
		int ulen=0,tool;
		while (sscanf(tools+ulen,"%d",&tool)==1)
		{
			add(i,n1+tool,1e8);
		    if (tool==0) ulen++;
	    	else
			{
    		    while (tool)
				{
        		    tool/=10;
        	    	ulen++;
	    	    }
    		}
    		ulen++;
		}
	}
	for (int i=1;i<=n2;i++)
	{
		scanf("%d",&x);
		add(n1+i,t,x);
	}
	while (bfs())
	{
		ans+=dfs(s,1e8);
	}
	for (int i=1;i<=n1;i++)
	{
		if (dep[i]) cout<<i<<' ';
	}
	cout<<endl;
	for (int i=1;i<=n2;i++)
	{
		if (dep[i+n1]) cout<<i<<' ';
	}
	cout<<endl;
	cout<<sum-ans;
    return 0;
}
//Now,only God know
//但行好事，莫问前程