P4014 分配问题

题目描述:
有n件工作要分配给n个人做,第i个人做第j件工作产生的效益为Cij,试设计一个将n件工作分配给n个人做的分配方案,使产生的总效益最大/最小

输出格式:
两行分别输出最小总效益和最大总效益


题解:
比较裸的费用流,建一个超级源点连接每一个人,边权为1,费用为0,同样的方式将每一个工作都连向超级汇点,边权为1,费用为0,而人与工作直接连接,边权为1,费用为Cij,保证每个点都只能被使用一次

因为EK算法总的来说就是每次通过spfa找到花费最小的可行流,然后进行增广,所以我们跑最大费用最大流的时候就只需要将spfa改成最长路就OK啦


AC代码:

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define LL long long
const int MAXN = 2000+50;
const int MAXM = 2e6+50;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
int n,m,s,t,tot=1,head[MAXN],nxt[MAXM],w[MAXM],to[MAXM],co[MAXM];
int dis[MAXN],vis[MAXN],pre[MAXN],flow[MAXN];
inline void ade(int u,int v,int ww,int cost){
    to[++tot]=v; w[tot]=ww; co[tot]=cost; nxt[tot]=head[u]; head[u]=tot;
}
inline void add(int u,int v,int w,int cost){ ade(u,v,w,cost); ade(v,u,0,-cost); }
inline int spfa_min(){
    for(int i=0;i<=t;i++) dis[i]=INF;
    queue<int> que; que.push(s); dis[s]=0; vis[s]=1; flow[s]=INF;
    while(!que.empty()){
        int u=que.front(); que.pop(); vis[u]=0;
        for(int i=head[u];i;i=nxt[i]){
            int v=to[i];
            if(w[i] && dis[u]+co[i]<dis[v]){
                dis[v]=dis[u]+co[i]; pre[v]=i;
                flow[v]=min(flow[u],w[i]);
                if(!vis[v]) vis[v]=1,que.push(v);
            }
        }
    }
    return dis[t]!=INF;
}
inline int spfa_max(){
    for(int i=0;i<=t;i++) dis[i]=-INF;
    queue<int> que; que.push(s); dis[s]=0; vis[s]=1; flow[s]=INF;
    while(!que.empty()){
        int u=que.front(); que.pop(); vis[u]=0;
        for(int i=head[u];i;i=nxt[i]){
            int v=to[i];
            if(w[i] && dis[u]+co[i]>dis[v]){
                dis[v]=dis[u]+co[i]; pre[v]=i;
                flow[v]=min(flow[u],w[i]);
                if(!vis[v]) vis[v]=1,que.push(v);
            }
        }
    }
    return dis[t]!=-INF;
}
signed main(){
#ifndef ONLINE_JUDGE
    freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
#endif // ONLINE_JUDGE
    scanf("%d",&n); s=0,t=2*n+1;
    for(int i=1;i<=n;i++) add(s,i,1,0);
    for(int i=1;i<=n;i++) add(i+n,t,1,0);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            int x; scanf("%d",&x);
            add(i,j+n,1,x);
        }
    }
    int res_min=0,res_max=0;
    while(spfa_min()){
        res_min += dis[t]*flow[t];
        int x=t;
        while(x!=s){
            int i=pre[x];
            w[i]-=flow[t]; w[i^1]+=flow[t]; x=to[i^1];
        }
    }
    printf("%d\n",res_min);
    for(int i=2;i<=tot;i+=2) w[i]+=w[i^1],w[i^1]=0;
    while(spfa_max()){
        res_max += dis[t]*flow[t];
        int x=t;
        while(x!=s){
            int i=pre[x];
            w[i]-=flow[t]; w[i^1]+=flow[t]; x=to[i^1];
        }
    }
    printf("%d\n",res_max);
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值