POJ1459Power Network（电网）

POJ1459Power Network（电网）

题目描述：
一个电网包含一些结点（电站、消费者、调度站），这些结点通过电线连接。每个结点u 可能被供给s(u)的电能，s(u)≥0；同时也可能产生p(u)的电能，0≤p(u)≤pmax(u)；站点u 还有可能消费c(u)电能，0≤c(u)≤min( s(u), cmax(u) )；可能传输d(u)的电能，d(u) = s(u) + p(u) - c(u)。以上这些量存在以下限制关系：对每个电站，c(u) = 0；对每个消费者，p(u) = 0；对每个调度站，p(u) = c(u) = 0。在电网中两个结点u 和v 之间最多有一条电线连接。从结点u 到结点v 传输L(u,v)的电能，0≤L(u,v)≤Lmax(u,v)。定义Con为c(u)的总和，表示电网中消费电能的总和。本题的目的是求Con的最大值。电网的一个例子如图6.24 所示。在图(a)中，电站结点u 的标记”x/y”代表p(u) = x、pmax(u) =y。消费者结点u 的标记”x/y”代表c(u) = x、cmax(u)= y。每条电线所对应的边(u,v)，其标记”x/y”代表L(u,v) = x、Lmax(u,v) = y。在图(b)中，消费的最大电能Con = 6，图(a)列出了在此状态下各个站点的s(u)、p(u)、c(u)和d(u)。注意，如图(b)所示的电网中，电能的流动还存在其他状态，但消费的电能总和不超过6。

输入描述：
输入文件中包含多个测试数据。每个测试数据描述了一个电网。每个测试数据的第1 行为4个整数：n np nc m，其中，0≤n≤100，代表结点数目；0≤np≤n，代表电站数目；0≤nc≤n，代表消费者数目；0≤m≤n2，代表传输电线的数目。接下来有m 个三元组，(u,v)z，其中u 和v为结点序号（结点序号从0 开始计起），0≤z≤1000，代表Lmax(u,v)的值。接下来有np 个二元组，(u)z，其中u 为电站结点的序号，0≤z≤10000，代表pmax(u)的值；每个测试数据的最后是nc 个二元组，(u)z，其中u 为消费者结点的序号，0≤z≤10000，代表cmax(u)的值。所有数据都是整数。除三元组(u,v)z 和二元组(u)z 中不含空格外，输入文件中其他位置允许出现空格。测试数据一直到文件尾。

 

 

网络流问题

/*摘抄自：https://comzyh.com/blog/archives/568/*/ 


#include <cstdio>
#include<cstring>
#include<iostream>
#include<vector>
#include<queue>
using namespace std;
int N,NP,NC,M;
struct Edge
{
    int u,v,cap;
    Edge(){}
    Edge(int u,int v,int cap):u(u),v(v),cap(cap) {}
} es[150*150];
int R,S,T;
vector<int>tab[109]; // 边集
int dis[109];
int current[109];
void addedge(int u,int v,int cap)
{
    tab[u].push_back(R);
    es[R++]=Edge(u, v, cap); // 正向边
    tab[v].push_back(R);
    es[R++]=Edge(v,u,0); // 反向边容量为0
    // 正向边下标通过异或就得到反向边下标, 2 ^ 1 == 3 ; 3 ^ 1 == 2
}
int BFS()
{
    queue<int>q;
    q.push(S);
    memset(dis,0x3f,sizeof(dis));
    dis[S]=0;
    while (!q.empty())
    {
        int h=q.front();
        q.pop();
        for (int i=0;i<tab[h].size();i++)
        {
            Edge&e=es[tab[h][i]];
            if (e.cap>0&&dis[e.v]==0x3f3f3f3f)
            {
                dis[e.v]=dis[h]+1;
                q.push(e.v);
            }
        }
    }
    return dis[T]<0x3f3f3f3f; // 返回是否能够到达汇点
}
int dinic(int x, int maxflow)
{
    if (x==T)
        return maxflow;
    // i = current[x] 当前弧优化
    for (int i=current[x];i<tab[x].size();i++)
    {
        current[x]=i;
        Edge&e=es[tab[x][i]];
        if (dis[e.v]==dis[x]+1&&e.cap>0)
        {
            int flow=dinic(e.v,min(maxflow,e.cap));
            if (flow)
            {
                e.cap-=flow; // 正向边流量降低
                es[tab[x][i]^1].cap+=flow; // 反向边流量增加
                return flow;
            }
        }
    }
    return 0; // 找不到增广路 退出
}
int DINIC()
{
    int ans=0;

    while (BFS()) // 建立分层图
    {
        int flow;
        memset(current,0,sizeof(current)); // BFS后应当清空当前弧数组
        while (flow=dinic(S,0x3f3f3f3f)) // 一次BFS可以进行多次增广
            ans+=flow;
    }
    return ans;
}
int main()
{
    while(scanf("%d%d%d%d",&N,&NP,&NC,&M)!=EOF)
    {
        R=0;
        S=N;
        T=N+1;
        for(int i=0;i<=T;i++)
            tab[i].clear();
        for(int i=0;i<M;i++)
        {
            int u,v,cap;
            scanf("(%d,%d)%d",&u,&v,&cap);
            addedge(u,v,cap);
        }
        for(int i=0;i<NP;i++)
        {
            int u,p;
            scanf("(%d)%d",&u,&p);
            addedge(S,u,p);
        } 
        for(int i=0;i<NC;i++)
        {
            int u, c;
            scanf("(%d)%d",&u,&c);
            addedge(u,T,c);
        }

        printf("%d\n",DINIC());
    }
    return 0;
}