CF #77 Div 1 C ,Div 2 D

        题目意思是对于一个给定的图，求从给定起点坐的士到给定终点的最小花费，并且的士的送达距离有限。

        该题目可以用最短路径来做。首先对原图求各个顶点到其余所有点的最短距离，这是因为人可以在任意地方下的士并且换乘其他的士。然后保存在距离上符合要求的点，这里用一个二维数组来保存（见代码），求最短路径是为了找出一条路径，保证对其中的点求“花费”最短路径时不再有其他条件的限制。接下来就是对保存有路径的二维数组求一次“花费”最短路径即可。

        之前第一次做的时候是用了构图的方法，先用dfs+回溯构图，即距离差在根节点的距离要求以内时，这两个节点即可连边，边的权值是的士的花费。但可能dfs处理得很笨，产生了不少的冗余边，结果只过了前14的tests。接下来的这份代码是看了牛人的，最短路的算法是spfa。

 

       

#include<cstdio>
#include<cstring>
#include<iostream>
#include<string>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstdlib>
using namespace std;
const int M=20010;
const int N=1100;
const __int64 inf=(__int64)1<<58;

struct Node
{
	int v;
	int next;
	__int64 w;

}edge[M];
int head[N],num;
vector<int> que;
__int64 dis[N],a[N],b[N];
int map[N][N];
bool vis[N];
int n,m,s,t;

void init()
{
	for(int i=0;i<=n+2;i++)
		head[i]=-1;
	num=0;
}

void addege(int u,int v,__int64 w)
{
	edge[num].v=v;
	edge[num].w=w;
	edge[num].next=head[u];
	head[u]=num++;
}

void solve(int s)
{
	memset(dis,60,sizeof(dis));
	dis[s]=0;
	memset(vis,0,sizeof(vis));
	vis[s]=1;
	que.clear();
	que.push_back(s);
	for(int i=0;i<que.size();i++)
	{
		int u=que[i];
		for(int j=head[u];j!=-1;j=edge[j].next)
		{
			int v=edge[j].v;
			if(dis[v]<=dis[u]+edge[j].w) continue;
			dis[v]=dis[u]+edge[j].w;
			if(vis[v]) continue;
			vis[v]=1;
			que.push_back(v);
		}
		vis[u]=0;
	}
	map[s][0]=0;
	for(int i=1;i<=n;i++)
		if(i!=s && dis[i]<=a[s]) map[s][++map[s][0]]=i;
}

int main()
{
	while(scanf("%d%d",&n,&m)==2)
	{
		scanf("%d%d",&s,&t);
		init();
		int u,v;
		__int64 c;
		for(int i=1;i<=m;i++)
		{
			scanf("%d%d%I64d",&u,&v,&c);
			addege(u,v,c);
			addege(v,u,c);
		}
		for(int i=1;i<=n;i++)
			scanf("%I64d%I64d",&a[i],&b[i]);
		for(int i=1;i<=n;i++)
			solve(i);
		memset(dis,60,sizeof(dis));
		dis[s]=0;
		memset(vis,0,sizeof(vis));
		vis[s]=1;
		que.clear();
		que.push_back(s);
		for(int i=0;i<que.size();i++)
		{
			int u=que[i];
			for(int i=1;i<=map[u][0];i++)
			{
				int v=map[u][i];
				if(dis[v]<=dis[u]+b[u]) continue;
				dis[v]=dis[u]+b[u];
				if(vis[v]) continue;
				vis[v]=1;
				que.push_back(v);
			}
			vis[u]=0;
		}
		if(dis[t]>(1ll<<60))
			dis[t]=-1;
		printf("%I64d\n",dis[t]);
	}
	return 0;
}