http://acm.nefu.edu.cn/JudgeOnline/problemshow.php?problem_id=476
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 |
#include<cstdio> #include<iostream> using namespace std; const int oo=1e9; /**oo 表示无穷大*/ const int mm=111111; /**mm 表示边的最大数量,记住要是原图的两倍,在加边的时候都是双向的*/ const int mn=999; /**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]; /**head 节点的链表头,work 用于算法中的临时链表头,dis 计算距离*/ /**初始化链表及图的信息*/ 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 Dinic_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; } return ret; } int a[1005][1005],c[10005],b[10005]; int main() { int n,m; while(~scanf("%d%d",&n,&m)) { int sum=0; for(int i=1;i<=n;i++) { scanf("%d%d",&a[i][0],&c[i]); sum+=a[i][0]; for(int j=1;j<=c[i];j++) scanf("%d",&a[i][j]); } prepare(n+m+2,0,n+m+1); for(int j=1;j<=m;j++) scanf("%d",&b[j]); for(int i=n+1;i<=n+m;i++) addedge(i,dest,b[i-n]); for(int i=1;i<=n;i++) { addedge(src,i,a[i][0]); for(int j=1;j<=c[i];j++) addedge(i,a[i][j]+n,b[a[i][j]]); } int ans=Dinic_flow(); printf("%d\n",sum-ans); } return 0; }
有关最大权闭合图: