题目链接:http://acm.split.hdu.edu.cn/showproblem.php?pid=4862
题解:K次路径覆盖的问题,可以通过费用流来解决。
首先通过二分图来实现最小路径覆盖,具体思路为:
因为一条路径上,只有一个点的入度为零,一个点的出度为零,那么我们可以通过将一个店拆分为出度和入度,然后左半部分向右半部分连边,然后源点向左半部分连边,右半部分向汇点连边,然后通过最大流可以求出整个图的入度(出度),然后通过总点数减去入度,可以得到路径的数量。
如何实现K个路径呢?
在左半部分点集中加一个点Q,与源点见边(k,0),然后Q向右半部分点集建边(1,0),然后跑费用流,如果流量等于n*m说明,能通过k条路径覆盖所有点集。
花费可以通过费用流来求出。
#include<bits/stdc++.h>
using namespace std;
//最小费用最大流,求最大费用只需要取相反数,结果取相反数即可。
//点的总数为N,点的编号0~N-1
const int MAXN=505;
const int MAXM=100000;
const int INF=0x3f3f3f3f;
namespace MCMF
{
struct Edge
{
int to,next,cap,flow,cost;
}edge[MAXM];
int head[MAXN],tol;
int pre[MAXN],dis[MAXN];
bool vis[MAXN];
int N;//节点总个数,节点编号从0~N-1
void init(int n)
{
N=n;
tol=0;
memset(head,-1,sizeof(head));
}
void addedge(int u,int v,int cap,int cost)
{
edge[tol].to=v;edge[tol].cap=cap;edge[tol].cost=cost;edge[tol].flow=0;edge[tol].next=head[u];
head[u]=tol++;
edge[tol].to=u;edge[tol].cap=0;edge[tol].cost=-cost;edge[tol].flow=0;edge[tol].next=head[v];
head[v]=tol++;
}
bool spfa(int s,int t)
{
queue<int>q;
for(int i=0;i<N;i++)
{
dis[i]=INF;
vis[i]=false;
pre[i]=-1;
}
dis[s]=0;
vis[s]=true;
q.push(s);
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;
if(edge[i].cap>edge[i].flow&&dis[v]>dis[u]+edge[i].cost)
{
dis[v]=dis[u]+edge[i].cost;
pre[v]=i;
if(!vis[v])
{
vis[v]=true;
q.push(v);
}
}
}
}
if(pre[t]==-1)
return false;
else
return true;
}
//返回的是最大流,cost存的是最小费用
int minCostMaxflow(int s,int t,int &cost)
{
int flow=0;
cost=0;
while(spfa(s,t))
{
int Min=INF;
for(int i=pre[t];i!=-1;i=pre[edge[i^1].to])
{
if(Min>edge[i].cap-edge[i].flow)
Min=edge[i].cap-edge[i].flow;
}
for(int i=pre[t];i!=-1;i=pre[edge[i^1].to])
{
edge[i].flow+=Min;
edge[i^1].flow-=Min;
cost+=edge[i].cost*Min;
}
flow+=Min;
}
return flow;
}
}
using namespace MCMF;
char mp[15][15];
int sv[15][15];
int get(int i,int j,int ti,int tj)
{
int co=abs(ti-i)+abs(tj-j)-1;
if(sv[i][j]==sv[ti][tj])
co-=sv[i][j];
return co;
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int T;
scanf("%d",&T);
for(int _=1;_<=T;_++)
{
int n,m,k;
scanf("%d%d%d",&n,&m,&k);
init(n*m*2+3);
for(int i=0;i<n;i++)
{
scanf("%s",mp[i]);
}
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
{
sv[i][j]=mp[i][j]-'0';
}
}
addedge(n*m*2,n*m*2+1,k,0);
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
{
addedge(n*m*2,i*m+j,1,0);
addedge(n*m+i*m+j,n*m*2+2,1,0);
addedge(n*m*2+1,n*m+i*m+j,1,0);
}
}
for(int i=0;i<n;i++)
{
for(int j=0;j<m;j++)
{
for(int i1=i;i1<n;i1++)
{
for(int j1=j;j1<m;j1++)
{
if((i==i1&&j1>j)||(i1>i&&j==j1))
addedge(i*m+j,n*m+i1*m+j1,1,get(i,j,i1,j1));
}
}
}
}
int co;
int mx=minCostMaxflow(n*m*2,n*m*2+2,co);
//printf("%d\n",mx);
printf("Case %d : ",_);
if(mx!=n*m)
{
printf("-1\n");
}
else
{
printf("%d\n",-co);
}
}
return 0;
}