hdu 5253

老 Jack 有一片农田，以往几年都是靠天吃饭的。但是今年老天格外的不开眼，大旱。所以老 Jack 决定用管道将他的所有相邻的农田全部都串联起来，这样他就可以从远处引水过来进行灌溉了。当老 Jack 买完所有铺设在每块农田内部的管道的时候，老 Jack 遇到了新的难题，因为每一块农田的地势高度都不同，所以要想将两块农田的管道链接，老 Jack 就需要额外再购进跟这两块农田高度差相等长度的管道。 

现在给出老 Jack农田的数据，你需要告诉老 Jack 在保证所有农田全部可连通灌溉的情况下，最少还需要再购进多长的管道。另外，每块农田都是方形等大的，一块农田只能跟它上下左右四块相邻的农田相连通。

Input

第一行输入一个数字 T(T≤10) T(T≤10)，代表输入的样例组数 

输入包含若干组测试数据，处理到文件结束。每组测试数据占若干行，第一行两个正整数  N,M(1≤N,M≤1000) N,M(1≤N,M≤1000)，代表老 Jack 有N行*M列个农田。接下来 N 行，每行 M 个数字，代表每块农田的高度，农田的高度不会超过100。数字之间用空格分隔。 

Output

对于每组测试数据输出两行： 

第一行输出："Case #i:"。i代表第i组测试数据。 

第二行输出 1 个正整数，代表老 Jack 额外最少购进管道的长度。

Sample Input

2
4  3
9 12 4
7 8 56
32 32 43
21 12 12
2  3
34 56 56
12 23 4

Sample Output

Case #1:
82
Case #2:
74

  题意概括： 有一块n*m的农田，农田是高低不平的，现在给一个n*m的数阵分别代表每块农田的高度，jack想用最少的水管把 所有的农田都连接起来进行灌溉，问最少需要多少水管？

解题思路： 根据题意只有相邻的农田能进行连接，首先我们需要算出所有农田的高度差当作权值，然后地的编号和权值对应 存好即可。

代码1：

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define N 1100
#define M 1100
int m,n,a[N][M],l,i,j;
int f[N*M*2];
struct edge
{
	int u,v,w;
}e[N*M*2];
int abs(int a)
{
	if(a<0)
		a=-a;
	return a;
}
bool cmp(edge a,edge b)
{
	return a.w<b.w;
}
void Ufset()
{
	for(i=1;i<=m*n*2;i++)
	{
		f[i]=i;
	}
}
int getf(int u)
{
	if(f[u]==u)
	{
		return u;
	}
	else
	{
		f[u]=getf(f[u]);
		return f[u];
	}
}
int merge(int u,int v)
{
	int t1=getf(u);
	int t2=getf(v);
	if(t1!=t2)
	{
		f[t2]=t1;
		return 1;
	}
	return 0;
}
void Kruskal()
{
	int num=0,sum=0;
	Ufset();
	
	for(i=0;i<l;i++)
	{
		if(merge(e[i].u,e[i].v)!=0)
		{
			num++;
			sum+=e[i].w;
		
		}
		if(num==n*m-1)
			break;
	}
	printf("%d\n",sum);
}
int main()
{
	int t,g;
	scanf("%d",&t);
	g=0;
	while(t--)
	{
		g++;
		l=0;
		memset(a,0,sizeof(a));
		scanf("%d%d",&n,&m);
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				scanf("%d",&a[i][j]);
			}
		}
		
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				if(a[i][j+1]!=0)
				{
					e[l].u=(i-1)*m+j;
					e[l].v=(i-1)*m+1+j;
					e[l].w=abs(a[i][j]-a[i][j+1]);
					l++;
				}
				if(a[i+1][j]!=0)
				{
					e[l].u=j+(i-1)*m;
					e[l].v=j+(i)*m;
					e[l].w=abs(a[i][j]-a[i+1][j]);
					l++;
				}
				
			}
		}
	
		printf("Case #%d:\n",g);
		sort(e,e+l,cmp);
		Kruskal();
		
	}
	return 0;
}

代码2：

这个并查集进行了优化，能比第一个省些时间。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define N 1100
#define M 1100
int m,n,a[N][M],l,i,j;
int p[N*M*2];
struct edge
{
	int u,v,w;
}e[N*M*2];
int abs(int a)
{
	if(a<0)
		a=-a;
	return a;
}
bool cmp(edge a,edge b)
{
	return a.w<b.w;
}
void Ufset()
{
	for(i=1;i<=m*n;i++)
	{
		p[i]=-1;
	}
}
int Find(int x)
{
	int s;
	for(s=x;p[s]>=0;s=p[s]);
	while(s!=x)
	{
		int tmp=p[x];
		p[x]=s;
		x=tmp;
	}
	return s;
}
void Union(int R1,int R2)
{
	int r1=Find(R1);
	int r2=Find(R2);
	int tmp=p[r1]+p[r2];
	if(p[r1]>p[r2])
	{
		p[r1]=r2;
		p[r2]=tmp;
	}
	else
	{
		p[r2]=r1;
		p[r1]=tmp;
	}
}
void Kruskal()
{
	int num=0,sum=0;
	Ufset();
	
	for(i=0;i<l;i++)
	{
		if(Find(e[i].u)!=Find(e[i].v))
		{
			num++;
			sum+=e[i].w;
			Union(e[i].u,e[i].v);
		}
		if(num==n*m-1)
			break;
	}
	printf("%d\n",sum);
}
int main()
{
	int t,g;
	scanf("%d",&t);
	g=0;
	while(t--)
	{
		g++;
		l=0;
		memset(a,0,sizeof(a));
		scanf("%d%d",&n,&m);
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				scanf("%d",&a[i][j]);
			}
		}
		
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				if(a[i][j+1]!=0)
				{
					e[l].u=(i-1)*m+j;
					e[l].v=(i-1)*m+1+j;
					e[l].w=abs(a[i][j]-a[i][j+1]);
					l++;
				}
				if(a[i+1][j]!=0)
				{
					e[l].u=j+(i-1)*m;
					e[l].v=j+(i)*m;
					e[l].w=abs(a[i][j]-a[i+1][j]);
					l++;
				}
				
			}
		}
	
		printf("Case #%d:\n",g);
		sort(e,e+l,cmp);
		Kruskal();
		
	}
	return 0;
}