题目链接：点击查看

题目大意：给出n个工人和n个工作，每个人做每一个工作的效率都是不同的，问如何分配能让效率最低/最高

题目分析：最小费用最大流和最大费用最大流的模板题，直接套模板跑答案就行了，没有任何细节需要注意，就是两个模板组合在一起时会有很多函数的功能都是相同的，确切来说修改一下spfa的内部实现然后复制一份一模一样的就好了，一开始可以先把题目给出的效率矩阵储存下来，在跑费用流之前重新建边就好了，不然在跑完第一次费用流后的图已经变为残余网络了，需要重新建边再跑一遍费用流就好了

感觉这个题毕竟是二分图最大权匹配，用KM跑应该能更快一些吧

代码：

费用流：40ms

#include<iostream>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<climits>
#include<cmath>
#include<cctype>
#include<stack>
#include<queue>
#include<list>
#include<vector>
#include<set>
#include<map>
#include<sstream>
using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=210;//点const int M=N*N;//边struct Edge
{int to,w,cost,next;
}edge[M];int head[N],cnt,maze[N][N],n,st=N-1,ed=st-1;void addedge(int u,int v,int w,int cost)
{edge[cnt].to=v;edge[cnt].w=w;edge[cnt].cost=cost;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].to=u;edge[cnt].w=0;edge[cnt].cost=-cost;edge[cnt].next=head[v];head[v]=cnt++;
}int d[N],incf[N],pre[N];bool vis[N];bool spfa1(int s,int t)
{memset(d,inf,sizeof(d));memset(vis,false,sizeof(vis));memset(pre,-1,sizeof(pre));queue<int>q;q.push(s);vis[s]=true;incf[s]=inf;d[s]=0;while(!q.empty()){int u=q.front();q.pop();vis[u]=false;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;int cost=edge[i].cost;if(!w)continue;if(d[v]>d[u]+cost){d[v]=d[u]+cost;pre[v]=i;incf[v]=min(incf[u],w);if(!vis[v]){vis[v]=true;q.push(v);}}}}return pre[t]!=-1;
}bool spfa2(int s,int t)
{memset(d,0xcf,sizeof(d));memset(vis,false,sizeof(vis));memset(pre,-1,sizeof(pre));queue<int>q;q.push(s);vis[s]=true;incf[s]=inf;d[s]=0;while(!q.empty()){int u=q.front();q.pop();vis[u]=false;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;int cost=edge[i].cost;if(!w)continue;if(d[v]<d[u]+cost){d[v]=d[u]+cost;pre[v]=i;incf[v]=min(incf[u],w);if(!vis[v]){vis[v]=true;q.push(v);}}}}return pre[t]!=-1;
}int update(int s,int t)
{int x=t;while(x!=s){int i=pre[x];edge[i].w-=incf[t];edge[i^1].w+=incf[t];x=edge[i^1].to;}return d[t]*incf[t];
}void init()
{memset(head,-1,sizeof(head));cnt=0;
}int solve1(int st,int ed)
{int ans=0;while(spfa1(st,ed))ans+=update(st,ed);return ans;
}int solve2(int st,int ed)
{int ans=0;while(spfa2(st,ed))ans+=update(st,ed);return ans;
}void build()
{init();for(int i=1;i<=n;i++){addedge(st,i,1,0);addedge(i+n,ed,1,0);}for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)addedge(i,j+n,1,maze[i][j]);
}int main()
{
//  freopen("input.txt","r",stdin);
//  ios::sync_with_stdio(false);init();scanf("%d",&n);for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)scanf("%d",&maze[i][j]);build();printf("%d\n",solve1(st,ed));build();printf("%d\n",solve2(st,ed));return 0;
}

KM：27ms

#include<iostream>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<climits>
#include<cmath>
#include<cctype>
#include<stack>
#include<queue>
#include<list>
#include<vector>
#include<set>
#include<map>
#include<sstream>
using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=210;int n;int la[N],lb[N];//顶标bool visa[N],visb[N];//访问标记int maze[N][N];//边权int match[N];//右部点匹配了哪一个左部点int upd[N];bool dfs(int x)
{visa[x]=true;//访问标记：x在交错树中for(int i=1;i<=n;i++){if(!visb[i]){if(la[x]+lb[i]-maze[x][i]==0)//相等子图{visb[i]=true;//访问标记：y在交错树中if(!match[i]||dfs(match[i])){match[i]=x;return true;}}elseupd[i]=min(upd[i],la[x]+lb[i]-maze[x][i]);}}return false;
} int KM1()
{memset(match,0,sizeof(match));for(int i=1;i<=n;i++){la[i]=-inf;lb[i]=0;for(int j=1;j<=n;j++)la[i]=max(la[i],maze[i][j]);}for(int i=1;i<=n;i++){while(1)//直到左部点找到匹配{memset(visa,false,sizeof(visa));memset(visb,false,sizeof(visb));memset(upd,inf,sizeof(upd));if(dfs(i))break;int delta=inf;for(int j=1;j<=n;j++)if(!visb[j])delta=min(delta,upd[j]);for(int j=1;j<=n;j++)//修改顶标{if(visa[j])la[j]-=delta;if(visb[j])lb[j]+=delta;}}}int ans=0;for(int i=1;i<=n;i++)ans+=maze[match[i]][i];return ans;
}int KM2()
{memset(match,0,sizeof(match));for(int i=1;i<=n;i++){la[i]=-inf;lb[i]=0;for(int j=1;j<=n;j++)la[i]=max(la[i],maze[i][j]);}for(int i=1;i<=n;i++){while(1)//直到左部点找到匹配{memset(visa,false,sizeof(visa));memset(visb,false,sizeof(visb));memset(upd,inf,sizeof(upd));if(dfs(i))break;int delta=inf;for(int j=1;j<=n;j++)if(!visb[j])delta=min(delta,upd[j]);for(int j=1;j<=n;j++)//修改顶标{if(visa[j])la[j]-=delta;if(visb[j])lb[j]+=delta;}}}int ans=0;for(int i=1;i<=n;i++)ans+=maze[match[i]][i];return ans;
}int main()
{
//  freopen("input.txt","r",stdin);
//  ios::sync_with_stdio(false);scanf("%d",&n);for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){int val;scanf("%d",&val);maze[i][j]=-val;}printf("%d\n",-KM1());for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)maze[i][j]*=-1;printf("%d\n",KM2());return 0;
}

洛谷 - P4014 分配问题(费用流/KM)相关推荐

洛谷 - P4015 运输问题(费用流)
题目链接:点击查看题目大意:有n个卖家和m个买家,每个卖家会卖ai个物品,每个买家会买bi个物品,每个卖家向每个卖家卖东西会有一定的代价,问如何匹配才能让代价最小/最大题目分析:和上一道题大同小异 ...
洛谷P3376 网络最大流
不是讲网络流,就是存个板子另外我的Dinic跑得比EK慢一倍可还行( 附两份比较好的教程,均来自洛谷日报 EK\sf \color{blue}EKEK Dinic\sf \color{blue}Di ...
洛谷3376 网络最大流
题目描述如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...
洛谷 P3128 [USACO15DEC]最大流Max Flow
题意简述给定一颗树,每次操作可以使两个点最短路上的点+1,求最大的点题解思路树上差分若操作u, v,则++f[u], ++f[v], --f[lca(u, v)], --f[father(lc ...
洛谷P2472-网络最大流（点的拆分）
题目链接-P2472 还是点的拆分,把所有点拆为两个点,其连接容量为点的数值,然后建图求解即可 #include <iostream> #include <cstdio> # ...
（费用流）洛谷P2469 [SDOI2010]星际竞速
洛谷P2469 [SDOI2010]星际竞速思路: 可以看出这是最小路径覆盖问题,最小路径覆盖=n−=n-=n−最大匹配数.所以这题就是最小费用最大流. 建图: 1.源点SSS向每个入点iii建(S ...
洛谷P4043 费用流
这题的建图方式可以类比洛谷P1251 我是由那个题才想到这么建的,由于每条边至少经过一次,我们又不清楚需要跑多少次,把边看成点,点与汇点相连,可是我们又不知道最大流应该是多少,直接这么连会发生错误.利 ...
洛谷P3980 志愿者招募——费用流
洛谷P3980 [NOI2008]志愿者招募--费用流题目介绍题目描述输入格式输出格式测试样例题解代码题目介绍题目描述链接: 传送门. 申奥成功后,布布经过不懈努力,终于成为奥组委 ...
洛谷 P4012 深海机器人问题【费用流】
题目链接:https://www.luogu.org/problemnew/show/P4012 洛谷 P4012 深海机器人问题输入输出样例输入样例#1: 1 1 2 2 1 2 3 4 5 6 ...

洛谷 - P4014 分配问题(费用流/KM)

费用流：40ms

KM：27ms

洛谷 - P4014 分配问题(费用流/KM)相关推荐

最新文章

热门文章