1. 程式人生 > >1459 Power Network解題報告(網路流最大流 超級源點建圖)

1459 Power Network解題報告(網路流最大流 超級源點建圖)

題目大意:

有三種結點,第一種是發電站,它能不會消耗電能;一種是使用者,他不會產生電能;再有一種就是一個我不認識的單詞,它既不會產生電也不會消耗電。每個點都有兩個屬性,產生電量,消耗電量;現在問你最多這個網路能消耗多少的電量。

分析:

網路流建圖:

超級源點連到所有的發電站,容量為各個發電站的產生電的能力;所有的使用者連線到超級匯點,容量為各個使用者的耗電能力;然後原圖中的各個邊保持不變,容量為所給最大可流通電量。

程式碼:

#include<iostream>
#include<stdio.h>
#include<queue>
#include<vector>
#include<string.h> #define inf 0x3f3f3f3f #define maxn 500 using namespace std; struct Edge { int from,to,cap,flow; Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f){} }; int n;int m;int np,nc; vector<Edge>edge; vector<int>G[maxn]; int a[maxn]; int p[maxn]; void init(int
n) { for(int i=0;i<n;i++)G[i].clear(); edge.clear(); } void add_edge(int from,int to,int cap) { edge.push_back(Edge(from,to,cap,0)); edge.push_back(Edge(to,from,0,0)); int lm=edge.size(); G[from].push_back(lm-2); G[to].push_back(lm-1); } int maxflow(int s,int t) { int
flow=0; while(1) { memset(a,0,sizeof(a)); queue<int>Q; Q.push(s); a[s]=inf; while(!Q.empty()) { int x=Q.front();Q.pop(); for(int i=0;i<G[x].size();i++) { Edge& e=edge[G[x][i]]; if(!a[e.to]&&e.cap>e.flow) { p[e.to]=G[x][i]; a[e.to]=min(a[x],e.cap-e.flow); Q.push(e.to); } } if(a[t])break; } if(!a[t])break; for(int u=t;u!=s;u=edge[p[u]].from) { edge[p[u]].flow+=a[t]; edge[p[u]^1].flow-=a[t]; } flow+=a[t]; } return flow; } char c; //int map[maxn][maxn]={0}; int main() { int fm=0,fnp=0,fnc=0;int la,lb,lc; while(cin>>n>>np>>nc>>m) { fm=0;fnp=0;fnc=0; //memset(map,0,sizeof(map)); memset(p,0,sizeof(p)); init(maxn); while(m!=0&&scanf("%c",&c)) { if(fm==m)break; if(c=='(') { scanf("%d,%d)%d",&la,&lb,&lc); //if(map[la][lb]<lc)map[la][lb]=lc; add_edge(la,lb,lc); fm++; } } /*for(int i=0;i<n;i++) { for(int j=i+1;j<n;j++) { if(map[i][j]) } }*/ while(np!=0&&scanf("%c",&c)) { if(fnp==np)break; if(c=='(') { scanf("%d)%d",&la,&lb); add_edge(n,la,lb); fnp++; } } while(nc!=0&&scanf("%c",&c)) { if(fnc==nc)break; if(c=='(') { scanf("%d)%d",&la,&lb); add_edge(la,n+1,lb); fnc++; } } cout<<maxflow(n,n+1)<<endl; } }