#115. 無源匯有上下界可行流

描述

這是一道模板題。

n n n 個點，m m m 條邊，每條邊 e e e 有一個流量下界 lower(e) \text{lower}(e) lower(e) 和流量上界 upper(e) \text{upper}(e) upper(e)，求一種可行方案使得在所有點滿足流量平衡條件的前提下，所有邊滿足流量限制。

輸入格式

第一行兩個正整數 n n n、m m m。

之後的 m m m 行，每行四個整數 s s s、t t t、lower \text{lower} lower、upper \text{upper} upper。

輸出格式

如果無解，輸出一行 NO。

否則第一行輸出 YES，之後 m m m 行每行一個整數，表示每條邊的流量。

樣例

樣例輸入 1

4 6
1 2 1 2
2 3 1 2
3 4 1 2
4 1 1 2
1 3 1 2
4 2 1 2

樣例輸出 1

NO

樣例輸入 2

4 6
1 2 1 3
2 3 1 3
3 4 1 3
4 1 1 3
1 3 1 3
4 2 1 3

樣例輸出 2

YES
1
2
3
2
1
1

數據範圍與提示

1≤n≤200,1≤m≤10200 1 \leq n \leq 200, 1 \leq m \leq 10200 1≤n≤200,1≤m≤10200

顯示分類標簽

板子題，就不細將了，有空整理一下。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int MAXN=2000001;
inline char nc()
{
    static char buf[MAXN],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,MAXN,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
    char c=nc();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=nc();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=nc();}
    return x*f;
}
int n,m,s,t;
struct node
{
    int u,v,flow,nxt;
}edge[MAXN];
int head[MAXN],cur[MAXN],A[MAXN];
int num=0;
void AddEdge(int x,int y,int z)
{
    edge[num].u=x;
    edge[num].v=y;
    edge[num].flow=z;
    edge[num].nxt=head[x];
    head[x]=num++;
}
void add_edge(int x,int y,int z)
{
    AddEdge(x,y,z);
    AddEdge(y,x,0);
}
int deep[MAXN],L[MAXN];
bool BFS()
{
    memset(deep,0,sizeof(deep));
    deep[s]=1;
    queue<int>q;
    q.push(s);
    while(q.size()!=0)
    {
        int p=q.front();
        q.pop();
        for(int i=head[p];i!=-1;i=edge[i].nxt)
            if(!deep[edge[i].v]&&edge[i].flow)
            {
                deep[edge[i].v]=deep[edge[i].u]+1;q.push(edge[i].v);
                if(edge[i].v==t) return 1;
            }
                
    }
    return deep[t];
    
}
int DFS(int now,int nowflow)
{
    if(now==t||nowflow<=0)
        return nowflow;
    int totflow=0;
    for(int &i=cur[now];i!=-1;i=edge[i].nxt)
    {
        if(deep[edge[i].v]==deep[edge[i].u]+1&&edge[i].flow)
        {
            int canflow=DFS(edge[i].v,min(nowflow,edge[i].flow));
            edge[i].flow-=canflow;
            edge[i^1].flow+=canflow;
            totflow+=canflow;
            nowflow-=canflow;
            if(nowflow<=0)
                break;
        }
    }
    return totflow;
}
int Dinic()
{
    int ans=0;
    while(BFS())
    {
        for(int i=0;i<=n;i++)
            cur[i]=head[i];
        ans+=DFS(s,1e8);
    }
    return ans;
}
int main()
{
    #ifdef WIN32
    freopen("a.in","r",stdin);
    #else
    #endif
    n=read();m=read();s=0;t=n+1;
    memset(head,-1,sizeof(head));
    for(int i=1;i<=m;i++)
    {
        int x=read(),y=read(),lower=read(),upper=read();L[i-1]=lower;
        add_edge(x,y,upper-lower);A[x]-=lower;A[y]+=lower;
    }
    int sum=0;
    for(int i=1;i<=n;i++)
    {
        if(A[i]>0) sum+=A[i],add_edge(s,i,A[i]);
        else add_edge(i,t,-A[i]);
    }
    if(Dinic()!=sum) printf("NO");
    else
    {
        printf("YES\n");
        for(int i=0;i<m;i++)
            printf("%d\n",edge[i*2|1].flow+L[i]);
    }
    return  0;
}