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

/* 最大權(quán)閉合子圖的一般姿勢(shì). 這題要輸出方案. 但是自己yy的輸出姿勢(shì)不對(duì)沒(méi)想明白... 覺(jué)得可能是殘余網(wǎng)絡(luò)不理想,有反向弧啥的. 但是我們通過(guò)反向弧總能找到方案. */ #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; }

轉(zhuǎn)載于:https://www.cnblogs.com/nancheng58/p/10068047.html

總結(jié)

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

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。