dijkstra模板+ 缩点-bupt新oj-187- Find the Fastest Server

先放一下dijkstra的模板。

不要问我为什么最短路这么基础的东西都要用模板！！！因为我图论弱啊！

typedef pair<int, int> P;

int V;
int d[SIZE_D];

void dijkstra(int s)
{
V = N;
priority_queue<P,vector<P>,greater<P> > que;
//fill(d,d + V, INF);
//为何楼上的fill没有用啊
memset(d,-1,sizeof(d));//第一次修改
d[s] = 0;
que.push(P(0,s));
while (!que.empty()){
P p = que.top();que.pop();
int ver = p.second;
if (d[ver] != -1 && d[ver] < p.first) continue;//第二次修改
for (int i = head[ver]; i != -1; i = pra[i].next){
int u = pra[i].to, cost = pra[i].w;
if (ser[u] != 0) continue;
if (d[u] == -1 || d[u] > d[ver] + cost){
d[u] = d[ver] + cost;
que.push(P(d[u],u));//第五次修改
}
}
}
}
//无法到达的点d是-1

题目链接：http://code.bupt.edu.cn/problem/p/187/

讲一下这道题目的题意

题意：有10^4个点，10^5条边，每条边有一个延时，告诉你10^3个服务站。给你10^4个查询点，问每个查询点到最快能到的服务站的延时是多少。注意：一共有40组样例！！！

解答：

1、我一开始用每输入一个服务站就bfs更新所有点的最小距离d【i】的方法，我算了复杂度是10^3 × 10^ 4 大约10的七次方，时限一秒可以过，但是有40组样例啊！10的8次方啊！

2、然后是用二分算一个距离，bfs看看在这个距离里面有没有服务站，有就跳出来，复杂度我不懂算啊！因为bfs随时可以跳出啊！但我觉得比方法一快啊！可是结果比方法一慢了一倍啊！

3、最后就是正解啦！

因为这道题里面，server没有区别，所以可以缩成一个点。这是一个好思路啊！

无差别的点可以缩成一个点啊！！！！！

然后跑一遍dijkstra。把这个超级点作为起点，记录d【i】，然后每次查询就是O（1）的了。

还可以用spfa一试，待会试试看，再来贴代码。

4、还有一件事情啊！那个样例输出是错的啊！！！

不过刚刚让田田把样例改对了。我应该早点跟他说的，刚刚交代码每次都交两遍，wa的次数要除以2了。不过现在的样例是对的了真是造福后人啊！

5、不知道为什么我用fill函数不能把数组变成INF.我查了一下http://www.cplusplus.com/reference/algorithm/fill/ 明明就是在头文件algorithm里面的啊！

粗心的傻逼错：

1、无向图啊！双向边啊！

2、图可能不连通啊！要输出-1啊！！！

3、变量居然打错了orz简直傻逼。。。

放一下dijkstra的代码

#include<bits/stdc++.h>
#define SIZE_D 20005
#define SIZE_B 500005
#define INF 0x3f3f3f3f
using namespace std;
int N,M;
int e,head[SIZE_D];
struct pp
{
int to,w,next;
}pra[SIZE_B];
void init()
{
e = 0;
memset(head,-1,sizeof(head));
}
int addedge(int x,int y,int z)
{
pra[e].to = y;
pra[e].w = z;
pra[e].next = head[x];
head[x] = e++;
}

bool ser[SIZE_D];

typedef pair<int, int> P;

int V;
int d[SIZE_D];

void dijkstra(int s)
{
V = N;
priority_queue<P,vector<P>,greater<P> > que;
//fill(d,d + V, INF);
//为何楼上的fill没有用啊
memset(d,-1,sizeof(d));//第一次修改
d[s] = 0;
que.push(P(0,s));
while (!que.empty()){
P p = que.top();que.pop();
int ver = p.second;
if (d[ver] != -1 && d[ver] < p.first) continue;//第二次修改
for (int i = head[ver]; i != -1; i = pra[i].next){
int u = pra[i].to, cost = pra[i].w;
if (ser[u] != 0) continue;
if (d[u] == -1 || d[u] > d[ver] + cost){
d[u] = d[ver] + cost;
que.push(P(d[u],u));//第五次修改
}
}
}
}
//无法到达的点d是-1

void copyc(int oldd,int newd)
{
for (int i = head[oldd]; i != -1; i = pra[i].next){
addedge(newd,pra[i].to,pra[i].w);
addedge(pra[i].to,newd,pra[i].w);//第三次修改
}
}
int main()
{
//freopen("input.txt","r",stdin);
int T;
scanf("%d",&T);
while (T--){
init();
scanf("%d %d",&N,&M);
for (int i = 0; i < M; i++){
int temp1,temp2,temp3;
scanf("%d %d %d",&temp1,&temp2,&temp3);
addedge(temp1,temp2,temp3);
addedge(temp2,temp1,temp3);
}
int k;
scanf("%d",&k);
memset(ser,0, sizeof(ser));
int super;
for (int i = 1; i <= k; i++){
int tempv;
scanf("%d",&tempv);
ser[tempv] = 1;
if (i == 1)
super = tempv;
else
copyc(tempv,super);
}
dijkstra(super);
int Q;
scanf("%d",&Q);
for (int i = 0; i < Q; i++){
int tempq;
scanf("%d",&tempq);
if (ser[tempq] == 1 ){
puts("0");
}else{
printf("%d\n",d[tempq]);
}
}
puts("");
}
return 0;

}