洛谷P4014分配问题——网络流24题

题目：https://www.luogu.org/problemnew/show/P4014

最大/小费用最大流裸题。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
queue<int>q;
int const inf=1e9;
int n,t,head[],ct=,pre[],mx,mn,incf[],dis[];
bool vis[];
struct N{
    int hd,to,next,w,v;
    N(int h=,int t=,int n=,int w=,int v=):hd(h),to(t),next(n),w(w),v(v) {}
}edge[],tmp[];
bool spfa()
{
    memset(dis,-,sizeof dis);
    memset(incf,0x3f,sizeof incf);
    memset(vis,,sizeof vis);
    memset(pre,,sizeof pre);
    while(q.size())q.pop();
    dis[]=;vis[]=;q.push();
    while(q.size())
    {
        int x=q.front();q.pop();vis[x]=;
        for(int i=head[x];i;i=edge[i].next)
        {
            int u=edge[i].to;
            if(dis[u]<dis[x]+edge[i].w&&edge[i].v)
            {
                dis[u]=dis[x]+edge[i].w;
                pre[u]=i;
                incf[u]=min(incf[x],edge[i].v);
                if(!vis[u])vis[u]=,q.push(u);
            }
        }
    }
    return pre[t];
}
bool spfa2()
{
    memset(dis,0x3f,sizeof dis);
    memset(incf,0x3f,sizeof incf);
    memset(vis,,sizeof vis);
    memset(pre,,sizeof pre);
    while(q.size())q.pop();
    dis[]=;vis[]=;q.push();
    while(q.size())
    {
        int x=q.front();q.pop();vis[x]=;
        for(int i=head[x];i;i=tmp[i].next)
        {
            int u=edge[i].to;
            if(dis[u]>dis[x]+tmp[i].w&&tmp[i].v)
            {
                dis[u]=dis[x]+tmp[i].w;
                pre[u]=i;
                incf[u]=min(incf[x],tmp[i].v);
                if(!vis[u])vis[u]=,q.push(u);
            }
        }
    }
    return pre[t];
}
void mcf()
{
    memcpy(tmp,edge,sizeof edge);
    while(spfa())
    {
        mx+=incf[t]*dis[t];
        for(int i=pre[t];i;i=pre[edge[i].hd])
        {
            edge[i].v-=incf[t];
            edge[i^].v+=incf[t];
        }
    }
    while(spfa2())
    {
        mn+=incf[t]*dis[t];
        for(int i=pre[t];i;i=pre[edge[i].hd])
        {
            tmp[i].v-=incf[t];
            tmp[i^].v+=incf[t];
        }
    }
}
int main()
{
    scanf("%d",&n);
    t=*n+;
    for(int i=;i<=n;i++)
    {
        edge[++ct]=(N){,i,head[],,};head[]=ct;
        edge[++ct]=(N){i,,head[i],,};head[i]=ct;
    }
    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
        {
            int x;
            scanf("%d",&x);
            edge[++ct]=(N){i,j+n,head[i],x,};head[i]=ct;
            edge[++ct]=(N){j+n,i,head[j+n],-x,};head[j+n]=ct;
        }
    for(int j=;j<=n;j++)
    {
        edge[++ct]=(N){j+n,t,head[j+n],,};head[j+n]=ct;
        edge[++ct]=(N){t,j+n,head[t],,};head[t]=ct;
    }
    mcf();
    printf("%d\n%d",mn,mx);
    return ;
}