PAT甲级真题1072 加油站

加油站的建造位置必须使加油站与距离它最近的房屋的距离尽可能远。

与此同时，它还必须保证所有房屋都在其服务范围内。

现在，给出了城市地图和加油站的几个候选位置，请你提供最佳建议。

如果有多个解决方案，请输出选取位置与所有房屋的平均距离最小的解决方案。

如果这样的解决方案仍然不是唯一的，请输出选取位置编号最小的解决方案。

输入格式
第一行包含四个整数 NN，房屋总数，MM，加油站的候选位置总数，KK，连接房屋或加油站的道路总数，DsDs 加油站的最大服务范围。

所有房屋的编号从 11 到 NN，所有加油站侯选位置编号从 G1 到 GM。

接下来 KK 行，每行格式如下：

P1 P2 Dist
其中，P1 和 P2 表示一条 无向 道路连接的两个房屋或加油站侯选位置的编号，Dist 是道路长度，这是一个整数。

输出格式
第一行输出所选位置的编号。

第二行输出加油站与距离其最近的房屋之间的距离以及与所有房屋之间的平均距离，精确到小数后一位。

如果解决方案不存在，则输出 No Solution。

数据范围
1≤N≤1031≤N≤103,
1≤M≤101≤M≤10,
1≤K≤1041≤K≤104,
1≤Ds≤2001≤Ds≤200,
1≤Dist≤50001≤Dist≤5000
输入样例1：
4 3 11 5
1 2 2
1 4 2
1 G1 4
1 G2 3
2 3 2
2 G2 1
3 4 2
3 G3 2
4 G1 3
G2 G1 1
G3 G2 2
输出样例1：
G1
2.0 3.3
输入样例2：
2 1 2 10
1 G1 9
2 G1 20
输出样例2：
No Solution

 #include<bits/stdc++.h>
using namespace std;
const int N=1020,INF = 0x3f3f3f3f;;
int n,m,k,d;
int dist[N],st[N],g[N][N];
int getid(string a){
	int id=0,i=0,len=a.length();
	while(i<len){
		if(a[i]!='G'){
		id=id*10+(a[i]-'0');
	}
		i++;
	}
	if(a[0]=='G')return n+id;
	else return id;
}
void dijiskstra(int s,int&mind,int&sumd){
	memset(dist,0x3f,sizeof dist);
	memset(st,0,sizeof st);
	dist[s]=0;
	for(int i=1;i<=n+m;i++){
		int t=-1;
	  for(int j=1;j<=n+m;j++){
	  	if(!st[j]&&(t==-1||dist[j]<dist[t]))
	  	    t=j;
	  }
	if(t==-1)break;
	st[t]=1;
	for(int v=1;v<=n+m;v++){
		dist[v]=min(dist[v],dist[t]+g[t][v]);
	}
  }
    for(int i=1;i<=n;i++) { //枚举居民房 
       if(dist[i]>d){
       	mind=-INF;
       	return ;//如果有一个居民房距离超过d，程序结束 
    }
  } 
    mind = INF, sumd = 0;
    for (int i = 1; i <= n; i ++ )
    {
        mind = min(mind, dist[i]);
        sumd += dist[i];
    } 
}
int main(){
	cin>>n>>m>>k>>d;
	memset(g, 0x3f, sizeof g);//别忘了初始化图 
	while(k--){
		string a,b;
		int c;
		cin>>a>>b>>c;
		int x=getid(a);
		int y=getid(b);
		g[x][y]=g[y][x]=min(g[x][y], c);	
	}
	int res=-1,mind=0,sumd=INF;
	for(int i=n+1;i<=n+m;i++){
		int d1,d2;
		dijiskstra(i,d1,d2);
		
		if(d1>mind)mind=d1,sumd=d2,res=i; 
		else if(d1==mind&&d2<sumd)sumd=d2,res=i;//距离相等，总和较小 	
	} 	
	if (res == -1) puts("No Solution");
    else printf("G%d\n%.1f %.1f\n", res - n, (double)mind, (double)sumd / n + 1e-8);
	
  return 0;
}