网络流24题 太空飞行计划SPJ（最大流）

【题目描述】

W 教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合 E={E1，E2，…，Em}，和进行这些实验需要使用的全部仪器的集合 I={I1， I2，…In}。 实验 Ej 需要用到的仪器是 I 的子集 Rj∈I。配置仪器 Ik 的费用为 Ck 美元。实验 Ej 的赞助商已同意为该实验结果支付 Pj 美元。W 教授的任务是找出一个有效算法， 确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大。这里净收益是指进行实验所获得的全部收入与配置仪器的全部费用的差额。 
            对于给定的实验和仪器配置情况，编程找出净收益最大的试验计划，输出最大收益值。

【输入格式】

第 1 行有 2 个正整数 m 和 n 。m 是实验数，n是仪器数。接下来的 m行，每行是一个实验的有关数据。第一个数赞助商同意支付该实验的费用；接着是该实验需要用到的若干仪器的编号。最后一行的 n个数是配置每个仪器的费用。

【输出格式】

输出一个整数，即最大的净收益值。

【输入样例】

2 3

10 1 2

25 2 3

5 6 7

【输出样例】

17

【题目分析】

本题难点不在于建图，而在于读入。。。对于s向每个实验连一条边，边权为可赞助最大费用，每个实验向所需仪器连边，边权设为正无穷，每台仪器向t连边，边权为仪器费用，跑一遍最大流即为最小花费，用总收益-最小花费即为最大收益。

【代码】

#include<bits/stdc++.h>
using namespace std;
const int MAXN=3010;
const int INF=0x3f3f3f3f;

int cnt,n,m,s,t;
int head[MAXN],depth[MAXN],cur[MAXN];
int nxt[MAXN],to[MAXN],w[MAXN];
queue<int> q;
int vis[MAXN];
int ans[MAXN];

bool bfs()
{
	while(!q.empty())
	  q.pop();
	memset(depth,0,sizeof(depth));
	memset(vis,0,sizeof(vis));
	depth[s]=1;
	vis[s]=1;
	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]==0&&w[i]!=0)
			{
				depth[v]=depth[u]+1;
				vis[v]=1;
				q.push(v);
			}
		}
	}
	if(depth[t]==0)
	  return false;
	return true;
}

int dfs(int u,int dist)
{
	if(u==t)
	  return dist;
	for(int &i=cur[u];i!=-1;i=nxt[i])
	{
		int v=to[i];
		if(depth[v]==depth[u]+1&&w[i]>0)
		{
			int di=dfs(v,min(dist,w[i]));
			if(di>0)
			{
				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;
}

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 main()
{
	memset(head,-1,sizeof(head));
	memset(nxt,-1,sizeof(nxt));
	scanf("%d%d",&n,&m);
	s=0,t=n+m+1;
	int sum=0;
	for(int i=1;i<=n;i++)
    {
        int x;
        scanf("%d",&x);
        sum+=x;
        add(0,i,x);
        char ch=getchar();
        while((ch=getchar())!='\n')
        {
            x=ch-'0';
            while((ch=getchar())&&ch>='0'&&ch<='9')
                x=x*10+ch-'0';
            add(i,x+n,INF);
            if(ch=='\n')break;
        }
    }
    for(int i=1;i<=m;i++)
    {
        int x;
        scanf("%d",&x);
        add(i+n,t,x);
    }
    int tmp=dinic();
    printf("%d\n",sum-tmp);
	return 0;
}

 

转载于:https://www.cnblogs.com/Ishtar/p/10010848.html