1 /*BFS+压入重标法*/  2   3   4 //Memory  Time   5 //384K    860MS    6   7 #include<iostream>  8 using namespace std;  9  10 const int inf=10001; 11  12 int n;  //总节点数 13 int np; //电站数 14 int nc; //用户数 15 int line;  //线路数 16 int cap[102][102];  //弧(u,v)的容量 17 int flow[102][102];  //弧(u,v)的流量 18 bool vist[102];   //标记点v是否已标号 19 int s,t;  //超级源，超级汇 20  21 class info   //当前点v的标记信息 22 { 23 public: 24     int pre;  //当前点v的前驱u 25     int lv;  //l(v) 26     int nei[101];  //当前结点直接指向的邻居结点 27     int pn;  //邻居结点的指针 28 }node[102]; 29  30 int min(int a,int b) 31 { 32     return a<b?a:b; 33 } 34  35 void back(void) 36 { 37     int x=t; 38     while(x!=s) 39     { 40         flow[ node[x].pre ][x] += node[t].lv;  //改进增广路 41         x=node[x].pre; 42  43     } 44     return; 45 } 46  47 bool bfs(void) 48 { 49     memset(vist,false,sizeof(vist)); 50     node[s].pre=-1; 51     node[s].lv=inf; 52     vist[s]=true; 53  54     int queue[102]; 55     int head=0; 56     int tail=0; 57     queue[tail++]=s; 58  59     while(head<=tail-1)  //注意，这是再也找不到增广路的结束条件 60     { 61         int x=queue[head]; 62         int y; 63         for(int i=0;i<node[x].pn;i++) 64         { 65             y=node[x].nei[i]; 66             if(!vist[y] && flow[x][y]<cap[x][y])   //搜索的目标要求是 未标记 & 非饱和弧 67             { 68                 queue[tail++]=y; 69  70                 vist[y]=true; 71                 node[y].pre=x; 72                 node[y].lv=min( node[x].lv , cap[x][y]-flow[x][y] ); 73             } 74             if(vist[t])  //当超级汇被标记 75                 break; 76         } 77         if(!vist[t]) 78             head++; 79         else 80             return true;  //搜索到一条增广路 81     } 82     return false; 83 } 84  85 int main(int i,int j,int u,int v,int z,char temp) 86 { 87     while(cin>>n>>np>>nc>>line) 88     { 89         /*Initial*/ 90  91         s=n; 92         t=n+1; 93         for(i=0;i<n+1;i++) 94             node[i].pn=0; 95  96         /*Input & Structure Maps*/ 97  98         for(i=1;i<=line;i++) 99         {100             cin>>temp>>u>>temp>>v>>temp>>z;101             if(u==v)102                 continue;   //不需要环103             cap[u][v]=z;104             flow[u][v]=0;   //每条边的流量都初始化为0105             node[u].nei[ node[u].pn++ ]=v;106         }107         for(i=1;i<=np;i++)108         {109             cin>>temp>>v>>temp>>z;110             cap[s][v]=z;     //建立超级源，指向所有电站111             flow[s][v]=0;112             node[s].nei[ node[s].pn++ ]=v;113         }114         for(i=1;i<=nc;i++)115         {116             cin>>temp>>u>>temp>>z;117             cap[u][t]=z;     //建立超级汇，被所有用户指向118             flow[u][t]=0;119             node[u].nei[ node[u].pn++ ]=t;120         }121 122         /*标号法找增广轨*/123 124         while(true)125         {126             if(bfs())  //如果能搜到到增广路127                 back();  //从超级汇开始回溯,改进增广路128             else129             {130                 int max=0;        //输出最大流131                 for(i=0;i<node[s].pn;i++)132                     max+=flow[s][ node[s].nei[i] ];133                 cout<<max<<endl;134                 break;135             }136         }137     }138     return 0;139 }