HDU 4862 Jump （费用流）

题目链接：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;
}