网络流24题c语言解法,太空飞行计划问题(网络流24题，七)



太空飞行计划问题

题目分析：

是一道中问题，就不深入的分析了。就是叫你求给你M个实验，N个仪器。每个实验要用到N个仪器中的一些，而每个仪使用每个仪器都要花费一些钱。而完成每个实验会有一定的报酬，现在叫你求出最大的利润。

利润 = 总收益 - 总花费

算法分析：

最大权闭合子图。

建模分析：

我们可以根据最大闭合子图的算法思想，建立一个图。

1、我们建立两个超级点S,T。

2、对每个实验跟S链接一条容量为收入的边。

3、对每个一起跟T链接一条容量为花费的边。

4、对每个实验要用到的一起链接一条容量为无穷大的边。

C语言链表实现：

#include

#include

#include

#include

#include

using namespace std;

const int MAXN = 101,MAXM = MAXN*MAXN*2,INF = ~0U >> 1;

struct Edge{

Edge *next,*op;

int t,c;

}*V[MAXN],*P[MAXN],ES[MAXM],*Stae[MAXN];

int N,M,S,T,EC,Ans,Maxflow;

int Lv[MAXN],Stap[MAXN];

void Clear()

{

for(int i = 0;i < MAXN;++i)

V[i] = NULL,ES[i].next = NULL,ES[i].op = NULL;

EC = 0; Ans = 0; Maxflow = 0;

}

inline void addedge(int a,int b,int c)

{

ES[++EC].next = V[a];V[a] = ES+EC;V[a]->t = b;V[a]->c = c;

ES[++EC].next = V[b];V[b] = ES+EC;V[b]->t = a;V[b]->c = 0;

V[a]->op = V[b]; V[b]->op = V[a];

}

bool Dinic_Lable() //构造层次图

{

int head,tail,i,j;

Stap[head=tail=0] = S;

memset(Lv,-1,sizeof(Lv));

Lv[S] = 0;

while(head<=tail){

i = Stap[head++];

for(Edge *e = V[i];e;e = e->next){

j = e->t;

if(e->c&&Lv[j]==-1){

Lv[j] = Lv[i] + 1;

if(j==T)

return true;

Stap[++tail] = j;

}

}

}

return false;

}

void Dinic_Augment() //增广

{

int i,j,delta,Stop;

for(i = S;i <= T;++i)

P[i] = V[i];

Stap[Stop=1] = S;

while(Stop){

i = Stap[Stop];

if(i != T){

for(;P[i];P[i] = P[i]->next)

if(P[i]->c&&Lv[i]+1==Lv[j=P[i]->t])

break;

if(P[i]){

Stap[++Stop] = j;

Stae[Stop] = P[i];

}

else

Stop--,Lv[i] = -1;

}

else

{

delta = INF;

for(i = Stop;i >= 2;--i)

if(Stae[i]->c < delta)

delta = Stae[i]->c;

Maxflow += delta;

for(i = Stop;i >= 2;--i){

Stae[i]->c -= delta;

Stae[i]->op->c += delta;

if(Stae[i]->c == 0)

Stop = i-1;

}

}

}

}

void Dinic()

{

while(Dinic_Lable())

Dinic_Augment();

}

void init()

{

int i,a,c;

S = 0; T = M+N+1;

for(i = 1;i <= M;++i){

scanf("%d",&c);

addedge(S,i,c);

Ans += c;

// for(;;){

// int sum = 0;

// c = getchar();

// while(c == ' ')c = getchar();

// if(c=='\n') break;

// while(c!=' '&&c!='\n'){

// sum *= 10;

// sum += c - '0';

// c = getchar();

// } printf("sum == %d\n",sum);

// addedge(i,sum+M,INF);

// if(c=='\n') break;

// }

for (;;){

while((c=getchar())==' '); ungetc(c,stdin);

if (c==10 || c==13) break;

scanf("%d",&a);

addedge(i,a+M,INF);

}

}

for(i = 1;i <= N;++i){

scanf("%d",&c);

addedge(i+M,T,c);

}

}

void print()

{

for(int i = 1;i <= M;++i)

if(Lv[i]!=-1)

printf("%d ",i);

putchar('\n');

for(int i = M+1;i <= M+N;++i)

if(Lv[i]!=-1)

printf("%d ",i-M);

Ans -= Maxflow;

printf("\n%d\n",Ans);

}

int main()

{

// freopen("Input.txt","r",stdin);

while(~scanf("%d%d",&M,&N)){

Clear();

init();

Dinic();

print();

}

return 0;

}

C++容器实现：

#include

#include

#include

#include

#include

#include

#include

using namespace std;

const int MAXN = 1e3 + 5,INF = ~0U >> 2;

struct Edge{

int from,to,cap,flow;

Edge(int f,int t,int c,int _f)

:from(f),to(t),cap(c),flow(_f){};

};

vector edges;

vector G[MAXN];

int N,M,S,T;

int d[MAXN],cur[MAXN];

bool vst[MAXN];

void Init()

{

for(int i = 0;i < MAXN;++i)

G[i].clear();

edges.clear();

}

inline void AddEdge(int from,int to,int cap)

{

edges.push_back(Edge(from,to,cap,0));

edges.push_back(Edge(to,from,0,0));

int sz = edges.size();

G[from].push_back(sz-2);

G[to].push_back(sz-1);

}

bool BFS()

{

memset(vst,false,sizeof(vst));

queue Q;

Q.push(S);

d[S] = 0;

vst[S] = true;

while(!Q.empty()){

int x = Q.front();

Q.pop();

for(int i = 0;i < (int)G[x].size();++i){

Edge& e = edges[G[x][i]];

if(!vst[e.to]&&e.cap>e.flow){

vst[e.to] = true;

d[e.to] = d[x] + 1;

Q.push(e.to);

}

}

}

return vst[T];

}

int DFS(int u,int a)

{

if(u==T||a==0)

return a;

int f,flow = 0;

for(int& i = cur[u];i < (int)G[u].size();++i){

Edge& e = edges[G[u][i]];

if(d[u]+1==d[e.to]&&(f=DFS(e.to,min(a,e.cap-e.flow)))>0){

e.flow += f;

edges[G[u][i]^1].flow -= f;

flow += f;

a -= f;

if(a==0)break;

}

}

return flow;

}

int Maxflow()

{

int flow = 0;

while(BFS()){

memset(cur,0,sizeof(cur));

flow += DFS(S,INF);

}

return flow;

}

void Solve()

{

int i,c,ans = 0;

S = 0; T = M+N+1;

for(i = 1;i <= M;++i){

scanf("%d",&c);

AddEdge(S,i,c);

ans += c;

string str;

getline(cin,str);

istringstream ssin(str);

while(ssin>>c)

AddEdge(i,c+M,INF);

}

for(i = 1;i <= N;++i){

scanf("%d",&c);

AddEdge(i+M,T,c);

}

int maxflow = Maxflow();

for(i = 1;i <= M;++i)

if(vst[i])

printf("%d ",i);

putchar('\n');

for(i = M+1;i <= M+N;++i)

if(vst[i])

printf("%d ",i-M);

printf("\n%d\n",ans-maxflow);

}

int main()

{

// freopen("Input.txt","r",stdin);

while(~scanf("%d%d",&M,&N)){

Init();

Solve();

}

return 0;

}