POJ - 2349 Arctic Network （最小生成树分块）

题目：http://poj.org/problem?id=2349

题意：

s个卫星频道，p个哨兵，频道之间的哨兵通话无距离限制，没有频道的哨兵通话距离为Di，自由分配卫星频道给相应哨兵，使得通话限制的最大距离Di最小

分析：

因为所有哨兵都要连通，且距离要求最小，所以先求出哨兵距离为权值的最小生成树

因为最小生成树中，每少一条边，则整体多一棵树，所以将最小生成树的边最大到小减去s-1个（用卫星频道通话），则剩下的最大距离为所求D

代码：

#include <stdio.h>
#include <iostream>
#include <string.h>
#include <string>
#include <math.h>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <vector>

using namespace std;

#define MAX 2020
#define MIN -1e+10
#define INF 0x7f7f7f7f

int n, m, i, j;
double lowcost[MAX];
double Edge[MAX][MAX];

void Prim()
{
int ans = 0;
double a[MAX];
for(i = 1; i<=n; i++)
{
lowcost[i] = Edge[1][i];
}
lowcost[1] = -1;
for(i = 1; i<n; i++)
{
int v = -1, mint = INF;
for(j = 1; j<=n; j++)
{
if(lowcost[j]!=-1 && lowcost[j] < mint)
{
v = j; mint = lowcost[j];
}
}
if(j != -1)
{
//ans += lowcost[v];
a[i] = lowcost[v];
lowcost[v] = -1;
for(j = 1; j<=n; j++)
{
if(lowcost[j] !=-1 && lowcost[j] > Edge[v][j])
{
lowcost[j] = Edge[v][j];
}
}
}
else break;
}
sort(a + 1, a + i);
for(j = 0; j<m-1; j++) i--;
printf("%.2lf\n", a[i-1]);
}

double x[MAX], y[MAX];

int main()
{
freopen("a.txt", "r", stdin);

int t;
scanf("%d", &t);
while(t--)
{
scanf("%d%d", &m, &n);
for(i = 1; i<=n; i++)
{
scanf("%lf%lf", &x[i], &y[i]);
}
for(i = 1; i<n; i++)
{
for(j = i+1; j<=n; j++)
{
Edge[i][j] = Edge[j][i] = sqrt((x[i]-x[j])*(x[i]-x[j]) + (y[i]-y[j])*(y[i]-y[j]));
}
}
Prim();
}
return 0;
}