题目链接:POJ2112
题意:农场主John 将他的K(1≤K≤30)个挤奶器运到牧场,在那里有C(1≤C≤200)头奶牛,在奶牛和挤奶器之间有一组不同长度的路。K个挤奶器的位置用1~K的编号标明,奶牛的位置用K+1~K+C 的编号标明。每台挤奶器每天最多能为M(1≤M≤15)头奶牛挤奶。编写程序,寻找一个方案,安排每头奶牛到某个挤奶器挤奶,并使得C 头奶牛需要走的所有路程中的最大路程最小。每个测试数据中至少有一个安排方案。每条奶牛到挤奶器有多条路。
题解:看了好久才懂“C 头奶牛需要走的所有路程中的最大路程最小”这句话的意思。每头奶牛到挤奶器之间都有一个最短距离的路程,现在要求 所有奶牛的最短路程中的最大值最小。 先用Floyd求出各个实体之间的最短距离,然后二分枚举最短距离,用最大流求出最大能完成的奶牛个数,小于总的奶牛个数,说明距离太小,大于说明距离可以继续缩小。
最大流建图:源点和奶牛相连,容量为1,机器和汇点相连,容量为M,奶牛和机器之间,连接 距离小于当前枚举距离的边,容量为1。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int MAX=200+30+10;
const int INF=0x3f3f3f3f;
int dis[MAX][MAX];
int g[MAX][MAX];
int level[MAX];
int N;
int S,E;
int k,c,m;
int BFS()//构建层次网络
{
memset(level,0xff,sizeof(level));
queue<int> que;
que.push(S);
level[S]=0;
while(!que.empty())
{
int u=que.front();
que.pop();
for(int v=0;v<N;v++)
{
if(level[v]<0&&g[u][v]>0)//存在路径而且没有被标记层次
{
level[v]=level[u]+1;
que.push(v);
}
}
}
if(level[E]>0) return 1;
return 0;
}
int DFS(int x,int low)
{
if(x==E) return low;
for(int i=0;i<N;i++)
{
int a;
if(g[x][i]>0&&level[i]==level[x]+1&&(a=DFS(i,min(low,g[x][i]))))
{
g[x][i]-=a;
g[i][x]+=a;
return a;
}
}
return 0;
}
int Dinic()
{
int ans=0;
while(BFS())
{
int tans;
while(tans=DFS(S,INF)) ans+=tans;
}
return ans;
}
void floyd()
{
for(int m=1;m<=k+c;m++)
{
for(int i=1;i<=k+c;i++)
{
for(int j=1;j<=k+c;j++)
{
dis[i][j]=min(dis[i][j],dis[i][m]+dis[m][j]);
}
}
}
}
void build(int maxdis)
{
memset(g,0,sizeof(g));
for(int i=1;i<=k;i++)
{
g[i][E]=m;
}
for(int i=k+1;i<=k+c;i++)
{
g[S][i]=1;
}
for(int i=k+1;i<=k+c;i++)
{
for(int j=1;j<=k;j++)
{
if(dis[i][j]<=maxdis) g[i][j]=1;
}
}
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
while(scanf("%d%d%d",&k,&c,&m)!=EOF)
{
for(int i=1;i<=k+c;i++)
{
for(int j=1;j<=k+c;j++)
{
scanf("%d",&dis[i][j]);
if((dis[i][j]==0&&i!=j)) dis[i][j]=INF;
}
}
floyd();
N=k+c+2;S=0;E=N-1;
int l=0,r=200*MAX;
while(l<r)
{
int mid=(l+r)/2;
build(mid);
/* for(int i=0;i<N;i++)
{
for(int j=0;j<N;j++)
{
cout<<g[i][j]<<" ";
}
cout<<endl;
}*/
int ans=Dinic();
//cout<<ans<<endl;
if(ans>=c)r=mid;
else l=mid+1;
// cout<<l<<" "<<r<<endl;
}
cout<<r<<endl;
}
return 0;
}
学了二分图,又用二分图多重匹配写了一遍。依旧要用floyd跑一遍最短距离,然后二分最大距离,然后在牛和机器之间建立边,匹配的时候如果加一个小于距离的限制条件。
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAX=200+30+10;
const int MAXC=200+5;
const int MAXK=30+5;
const int INF=0x3f3f3f3f;
int map[MAX][MAX];
int g[MAXC][MAXK];
int k,c,m;
int mid;
void floyd(int n)
{
for(int k=0;k<n;k++)
{
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
map[i][j]=min(map[i][j],map[i][k]+map[k][j]);
}
}
}
}
int match[MAXK][MAXC];
int vis[MAXK];
bool dfs(int u)
{
for(int v=0;v<k;v++)
{
if(!vis[v]&&g[u][v]<=mid)
{
vis[v]=1;
if(match[v][0]<m)
{
match[v][++match[v][0]]=u;
return true;
}
for(int i=1;i<=m;i++)
{
if(dfs(match[v][i]))
{
match[v][i]=u;
return true;
}
}
}
}
return false;
}
int hungray()
{
int res=0;
memset(match,0,sizeof(match));
for(int i=0;i<c;i++)
{
memset(vis,0,sizeof(vis));
if(dfs(i)) res++;
}
return res;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
while(scanf("%d%d%d",&k,&c,&m)!=EOF)
{
memset(map,0x3f,sizeof(map));
for(int i=0;i<k+c;i++)
{
for(int j=0;j<k+c;j++)
{
scanf("%d",&map[i][j]);
if(map[i][j]==0&&i!=j) map[i][j]=INF;
}
}
floyd(k+c);
/* for(int i=0;i<k+c;i++)
{
for(int j=0;j<k+c;j++)
cout<<map[i][j]<<" ";
cout<<endl;
}*/
memset(g,0x3f,sizeof(g));
for(int i=0;i<c;i++)
{
for(int j=0;j<k;j++)
{
g[i][j]=map[k+i][j];
// cout<<g[k+i][j]<<" ";
}
// cout<<endl;
}
int l=0,r=50000;
while(l<r)
{
mid=(l+r)/2;
if(hungray()>=c) r=mid;
else l=mid+1;
}
printf("%d\n",r);
}
return 0;
}