Description
W教授正在为国家航天中心计划一系列的太空飞行。每次太空飞行可进行一系列商业性实验而获取利润。现已确定了一个可供选择的实验集合E={E1,E2,…,Em},和进行这些实验需要使用的全部仪器的集合I={I1,I2,…In}。实验Ej需要用到的仪器是I的子集。配置仪器Ik的费用为ck美元。实验Ej的赞助商已同意为该实验结果支付pj美元。W教授的任务是找出一个有效算法,确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大。这里净收益是指进行实验所获得的全部收入与配置仪器的全部费用的差额。
Input
多组数据输入.
每组输入第1行有2 个正整数m和n。m是实验数,n是仪器数。接下来的m 行,每行是一个实验的有关数据。第一个数赞助商同意支付该实验的费用;接下来是该实验需要的仪器数;接着是对应仪器的编号。最后一行的n个数是配置每个仪器的费用。
Output
每组输出最佳实验方案的净收益
Sample Input
2 3
10 2 1 2
25 2 2 3
5 6 7
Sample Output
17
题目出自nefu476
思路:这道题目是最大权闭合图,利润是实验的经费-最大权闭合图的权值(最小割)
源点连实验,权值为经费,实验连对应的仪器权值为无穷大,仪器连汇点,权值为仪器的费用,这样就可以建图了
#include<stdio.h>
#include<iostream>
#include <string.h>
using namespace std;
const int oo=1e9;
/**oo 表示无穷大*/
const int mm=99999999;
/**mm 表示边的最大数量,记住要是原图的两倍,在加边的时候都是双向的*/
const int mn=999;
/**mn 表示点的最大数量*/
struct sa
{
int val;
int num;
}yiqi[mn];
struct sa2
{
int val;
int need[mn];
}shiyan[mn];
int node,src,dest,edge;
/**node 表示节点数,src 表示源点,dest 表示汇点,edge 统计边数*/
int ver[mm],flow[mm],next[mm];
/**ver 边指向的节点,flow 边的容量 ,next 链表的下一条边*/
int head[mn],work[mn],dis[mn],q[mn];
void prepare(int _node, int _src,int _dest)
{
node=_node,src=_src,dest=_dest;
for(int i=0; i<node; ++i)head[i]=-1;
edge=0;
}
/**增加一条 u 到 v 容量为 c 的边*/
void addedge( int u, int v, int c)
{
ver[edge]=v,flow[edge]=c,next[edge]=head[u],head[u]=edge++;
ver[edge]=u,flow[edge]=0,next[edge]=head[v],head[v]=edge++;
}
/**广搜计算出每个点与源点的最短距离,如果不能到达汇点说明算法结束*/
bool Dinic_bfs()
{
int i,u,v,l,r=0;
for(i=0; i<node; ++i)dis[i]=-1;
dis[q[r++]=src]=0;
for(l=0; l<r; ++l)
for(i=head[u=q[l]]; i>=0; i=next[i])
if(flow[i]&&dis[v=ver[i]]<0)
{
/**这条边必须有剩余容量*/
dis[q[r++]=v]=dis[u]+1;
if(v==dest) return 1;
}
return 0;
}
/**寻找可行流的增广路算法,按节点的距离来找,加快速度*/
int Dinic_dfs( int u, int exp)
{
if(u==dest) return exp;
/**work 是临时链表头,这里用 i 引用它,这样寻找过的边不再寻找*/
for( int &i=work[u],v,tmp; i>=0; i=next[i])
if(flow[i]&&dis[v=ver[i]]==dis[u]+1&&(tmp=Dinic_dfs(v,min(exp,flow[i])))>0)
{
flow[i]-=tmp;
flow[i^1]+=tmp;
/**正反向边容量改变*/
return tmp;
}
return 0;
}
int Dicnic_flow()
{
int i,ret=0,delta;
while(Dinic_bfs())
{
for(i=0; i<node; ++i)work[i]=head[i];
while(delta=Dinic_dfs(src,oo))
{
ret+=delta;
// cout<<ret<<endl;
}
}
return ret;
}
int main()
{
int m,n,x;
ios::sync_with_stdio(false);
while(cin>>m>>n)
{
memset(shiyan,0,sizeof(shiyan));
memset(yiqi,0,sizeof(yiqi));
int ans=0;
prepare(n+m+2,0,n+m+1);
for(int i=1;i<=m;i++)
{
cin>>shiyan[i].val;
ans+=shiyan[i].val;
cin>>yiqi[i].num;
for(int j=1;j<=yiqi[i].num;j++)
{
cin>>shiyan[i].need[j];
}
}
for(int i=1;i<=n;i++)
{
cin>>yiqi[i].val;
}
for(int i=1;i<=m;i++)
{
addedge(0,i,shiyan[i].val);
for(int j=1;j<=yiqi[i].num;j++)
{
addedge(i,shiyan[i].need[j]+m,oo);
}
}
int flag[mn];
memset(flag,0,sizeof(flag));
for(int i=1;i<=n;i++)
{
/* for(int j=1;j<=yiqi[i].num;j++)//之前一直错在这块了,写的太麻烦了,错了还找不到,不需遍历然后再看该仪器是否出现再进行加边,直接把仪器到汇点连起来就行了
{
if(!flag[shiyan[i].need[j]])
{
addedge(shiyan[i].need[j]+m,n+m+1,yiqi[shiyan[i].need[j]].val);
flag[shiyan[i].need[j]]=1;
}
}*/
addedge(i+m,n+m+1,yiqi[i].val);
}
int sum=ans-Dicnic_flow();
cout<<sum<<endl;
}
return 0;
}