[網絡流24題] 太空飛行計劃
★★☆ 輸入文件：shuttle.in 輸出文件：shuttle.out 簡單對比
時間限制：1 s 內存限制：128 MB
【問題描述】
W 教授正在為國家航天中心計劃一系列的太空飛行。每次太空飛行可進行一系列商業性實驗而獲取利潤。現已確定了一個可供選擇的實驗集合E={E1,E2,…,Em}，和進行這些實驗需要使用的全部儀器的集合I={ I1, I2,…,In }。實驗Ej 需要用到的儀器是I的子集Rj∈I。配置儀器Ik 的費用為ck 美元。實驗Ej 的贊助商已同意為該實驗結果支付pj 美元。W教授的任務是找出一個有效算法，確定在一次太空飛行中要進行哪些實驗并因此而配置哪些儀器才能使太空飛行的凈收益最大。這里凈收益是指進行實驗所獲得的全部收入與配置儀器的全部費用的差額。
【編程任務】
對于給定的實驗和儀器配置情況，編程找出凈收益最大的試驗計劃。
【數據輸入】
第1行有2個正整數m和n(m,n <= 100)。m是實驗數，n是儀器數。接下來的m行，每行是一個實驗的有關數據。第一個數贊助商同意支付該實驗的費用；接著是該實驗需要用到的若干儀器的編號。最后一行的n個數是配置每個儀器的費用。
【結果輸出】
第1行是實驗編號；第2行是儀器編號；最后一行是凈收益。
【輸入文件示例】shuttle.in
2 3
10 1 2
25 2 3
5 6 7
【輸出文件示例】shuttle.out
1 2
1 2 3
17

/* 最大權閉合子圖的一般姿勢. 這題要輸出方案. 但是自己yy的輸出姿勢不對沒想明白... 覺得可能是殘余網絡不理想,有反向弧啥的. 但是我們通過反向弧總能找到方案. */ #include<iostream> #include<cstdio> #include<queue> #define MAXN 210 #define INF 1e9 using namespace std; int n,m,S,T,ans,cut=1,dis[MAXN],head[MAXN],c[MAXN],p[MAXN]; queue<int>q; bool b[MAXN]; struct data{int u,v,next,c;}e[MAXN*MAXN*6]; int read() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();return x*f; } void add(int u,int v,int c) {e[++cut].u=u;e[cut].v=v;e[cut].c=c;e[cut].next=head[u];head[u]=cut;e[++cut].u=v;e[cut].v=u;e[cut].c=0;e[cut].next=head[v];head[v]=cut; } bool bfs() {for(int i=S;i<=T;i++) dis[i]=-1;dis[S]=0;q.push(S);while(!q.empty()){int u=q.front();q.pop();b[u]=0;for(int i=head[u];i;i=e[i].next){int v=e[i].v;if(dis[v]==-1&&e[i].c){dis[v]=dis[u]+1;q.push(v); }}}return dis[T]!=-1; } int dfs(int u,int y) {if(u==T) return y;int rest=0;for(int i=head[u];i&&rest<y;i=e[i].next){int v=e[i].v;if(dis[v]==dis[u]+1&&e[i].c){int x=dfs(v,min(y-rest,e[i].c));e[i].c-=x;e[i^1].c+=x;rest+=x;}}if(!rest) dis[u]=-1;return rest; } void slove(int u) {for(int i=head[u];i;i=e[i].next)if(e[i].c&&!b[e[i].v]) b[e[i].v]=true,slove(e[i].v);return ; } void dinic() {while(bfs()) ans-=dfs(S,INF);slove(S);/*for(int i=head[S];i;i=e[i].next)if(e[i].c) b[e[i].v]=true;*/for(int i=1;i<=n;i++)if(b[i]) printf("%d ",i),b[i]=false;printf("\n");/*for(int i=head[T];i;i=e[i].next)if(e[i].c) b[e[i].v-n]=true;*/for(int i=n+1;i<=n+m;i++)if(b[i]) printf("%d ",i-n);printf("\n"); } int main() {freopen("shuttle.in","r",stdin);freopen("shuttle.out","w",stdout);int x,t=1;char ch;n=read(),m=read();S=0,T=n+m+1;for(int i=1;i<=n;i++){scanf("%d",&x),add(S,i,x),ans+=x;x=0;ch=getchar();while(ch!='\n'){if(ch>='0'&&ch<='9') x=x*10+ch-48;else if(x) add(i,x+n,INF),x=0;ch=getchar();}if(x) add(i,x+n,INF);}for(int i=1;i<=m;i++) p[i]=read();for(int i=1;i<=m;i++) add(i+n,T,p[i]);dinic();printf("%d",ans);return 0; }

轉載于:https://www.cnblogs.com/nancheng58/p/10068047.html

總結

以上是生活随笔為你收集整理的Cogs 727. [网络流24题] 太空飞行计划(最大权闭合子图)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。