hdu1853 Cyclic Tour (KM算法求最小费用流）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1853

题解：对于环上的每个点，入度=出度，将每个点拆成入点和出点，每个入点和一个出点匹配即可。用km算法求最小权值和。

#include <stdio.h>    
#include <string.h>   
#include <math.h>  
#define INF 0x3f3f3f3f
#define MAXN 102

int w[MAXN][MAXN],match[MAXN],n;  
int lx[MAXN],ly[MAXN],slack[MAXN];  
int visitx[MAXN],visity[MAXN];

int Scan()      
{      
	char ch;      
	int ret=0;      
	while((ch=getchar())<'0'||ch>'9');      
	while(ch>='0'&&ch<='9')      
	{      
		ret=ret*10+(ch-'0');      
		ch=getchar();      
	}      
	return ret;      
}  

int find(int x)
{
	int i,temp;
	visitx[x]=1;
	for(i=1;i<=n;++i)
	{
		if(visity[i])
			continue;
		temp=lx[x]+ly[i]-w[x][i];
		if(temp==0)
		{
			visity[i]=1;
			if(match[i]==-1||find(match[i]))
			{
				match[i]=x;//找到增广轨
				return 1;
			}
		}
		else if(slack[i]>temp)
		{//不在相等子图中slack 取最小的
			slack[i]=temp;
		}
	}
	return 0;
}

void KM()
{
	int i,j,d;
	memset(ly,0,sizeof(ly));
	memset(match,-1,sizeof(match));
	for(i=1;i<=n;++i)
	{//lx初始化为与它关联边中最大的 
		lx[i]=w[i][1];
		for(j=2;j<=n;++j)
			if(w[i][j]>lx[i])
				lx[i]=w[i][j];
	}
	for(i=1;i<=n;++i)
	{
		for(j=1;j<=n;++j)
			slack[j]=INF;
		while(1)
		{
			memset(visitx,0,sizeof(visitx));
			memset(visity,0,sizeof(visity));
			if(find(i))
				break;
			d=INF;
			for(j=1;j<=n;++j)
			{
				if(!visity[j]&&d>slack[j])
					d=slack[j];
			}
			for(j=1;j<=n;++j)
			{
				if(visitx[j])
					lx[j]-=d;
				if(visity[j])
					ly[j]+=d;
				else
					slack[j]-=d;
			}
		}
	}
}

int main()
{
	int i,j,m,x,y,val,ans;
	while(scanf("%d%d",&n,&m)!=EOF)
	{
		for(i=1;i<=n;++i)
		{
			for(j=1;j<=n;++j)
				w[i][j]=~INF+1;
		}
		while(m--)
		{
			x=Scan();
			y=Scan();
			val=Scan();
			val=~val+1;//取反
			if(val>w[x][y])//有重边
				w[x][y]=val;
		}
		KM();
		ans=0;
		val=1;
		for(i=1;i<=n;++i)
		{
			if(match[i]==-1||w[match[i]][i]==(~INF+1))
			{
				val=0;
				break;
			}
			ans+=w[match[i]][i];
		}
		if(val)
			printf("%d\n",(~ans+1));
		else
			printf("-1\n");
	}
	return 0;
}