L3-005 垃圾箱分布 (30分)（最短路）

L3-005 垃圾箱分布 (30分)
大家倒垃圾的时候，都希望垃圾箱距离自己比较近，但是谁都不愿意守着垃圾箱住。所以垃圾箱的位置必须选在到所有居民点的最短距离最长的地方，同时还要保证每个居民点都在距离它一个不太远的范围内。

现给定一个居民区的地图，以及若干垃圾箱的候选地点，请你推荐最合适的地点。如果解不唯一，则输出到所有居民点的平均距离最短的那个解。如果这样的解还是不唯一，则输出编号最小的地点。

输入格式：
输入第一行给出4个正整数：N（≤10​3​​ ）是居民点的个数；M（≤10）是垃圾箱候选地点的个数；K（≤10​4）是居民点和垃圾箱候选地点之间的道路的条数；D。​S是居民点与垃圾箱之间不能超过的最大距离。所有的居民点从1到N编号，所有的垃圾箱候选地点从G1到GM编号。

随后K行，每行按下列格式描述一条道路：

P1 P2 Dist

其中P1和P2是道路两端点的编号，端点可以是居民点，也可以是垃圾箱候选点。Dist是道路的长度，是一个正整数。

输出格式：
首先在第一行输出最佳候选地点的编号。然后在第二行输出该地点到所有居民点的最小距离和平均距离。数字间以空格分隔，保留小数点后1位。如果解不存在，则输出No Solution。

输入样例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

开始看题看到“最小最大”，以为要用二分，有点不敢想象pta也出二分了，看完题并不是，是最短路，看着有点麻烦，不过最近练最短路了。还没写就吃午饭了，边吃边想这个好麻烦啊。。。
反正这题不算难，就是麻烦，开始还定义二维数组了，也没必要。这题的关键点就是垃圾候选点也可以通过，开始不敢确定，手写模拟一遍才确定的，剩下的也没啥，输入的时候麻烦一点，跑n次最短路就好了，最后再比较一下。段错误一次开大内存就对了。

#include<stdio.h>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<stack>
using namespace std;
const int N=50005;
const int inf=0x3f3f3f3f;
int n,head[N],num=0,dis[N];
struct node{
	int to,next,wight;
}s[N];
void add(int u,int v,int w){
	s[++num].next=head[u];
	s[num].to=v;
	s[num].wight=w;
	head[u]=num;
}
void spfa(int x){
	int v[N]={0};
	memset(dis,inf,sizeof(dis));
	dis[x]=0;
	queue<int>q;
	q.push(x);
	while(!q.empty()){
		int t=q.front();
		q.pop();
		for(int i=head[t];i;i=s[i].next){
			int to=s[i].to;
			if(dis[to]>dis[t]+s[i].wight){
				dis[to]=dis[t]+s[i].wight;
				if(!v[to]){
					q.push(to);
				}
			}
		}
	}
}
int main() {
	int j,i,m,x,y,z,d,ans=0,k;
	char a[5],b[5];
	scanf("%d %d %d %d",&n,&m,&k,&d);
	while(k--){
		scanf("%s %s %d",&a,&b,&z);
		if(z<=d){
			x=y=0;
			if(a[0]!='G'){
				for(i=0;a[i];i++){
					x=10*x+(a[i]-'0');
				}
			} else {
				if(strlen(a)==3) x=n+10;
				else x=n+(a[1]-'0');
			}
			
			if(b[0]!='G'){
				for(i=0;b[i];i++){
					y=10*y+(b[i]-'0');
				}
			} else {
				if(strlen(b)==3) y=n+10;
				else y=n+(b[1]-'0');
			}
			
			add(x,y,z);
			add(y,x,z);
		}
	}
	int ans1=-inf,ans2=inf;
	for(i=1;i<=m;i++){
		spfa(n+i);
		int sum=0,mi=inf;
		for(j=1;j<=n;j++){
			if(dis[j]>d) {
				sum=inf;break;
			}
			sum+=dis[j];
			mi=min(mi,dis[j]);
		}
		if(sum!=inf){
			if((mi>ans1)||(mi==ans1&&sum<ans2)){
				ans=i;
				ans2=sum;
				ans1=mi;
			}
		}
	}
	if(ans2==inf) puts("No Solution");
	else {
		printf("G%d\n",ans);
		printf("%.1f %.1f",(double)ans1,(double)ans2/n);
	}
	return 0;
}