csu 1808 地铁 迪杰斯特拉+拆点建图

csu 1808

思路：2016年省赛的最短路题，去年省赛也是考了最短路，那么问题来了，今年会不会继续考最短路呢....这个题和一般的最短路不太一样，每个点还多了个路线，同一个点不同路线还得多花时间 |ci - cj | ，那么将点和路线结合（u，ci）看成一个点，前面输入的边容易转化成符合条件的边，但是问题来了，假设一个点u有1e5条路线，那么就有1e5个这样的点，如何在他们之间建立边呢，全部建立的话，就有1e10条边，（一个下午就被这个错误思想带偏了...），其实将这1e5个点按路线排序，只要在路线相邻的点建边就可以了，无需全部建边，这样就只有1e5条边，问题迎刃而解。网上还有很多思路是将边看做点，也去学学看....

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=4e5+10;
typedef long long LL;
const LL inf=1e14;
struct HeapNode
{
	LL d;
	int u;
	bool operator<(const HeapNode& rhs)const
	{
		return d>rhs.d;
	}
};
struct Edge
{
	int from,to;
	LL dist;
};
struct Dijkstra 
{
	int n,m;
	vector<Edge>edges;
	vector<int>G[maxn];
	bool done[maxn];
	LL d[maxn];
	int p[maxn]; 
	
	void init(int n)
	{
		this->n=n;
		for(int i=1;i<=n;i++)G[i].clear();
		edges.clear();
	}
	void AddEdge(int from,int to,LL dist)
	{
		edges.push_back((Edge){from,to,dist});
		m=edges.size();
		G[from].push_back(m-1);
	}
	void dij(vector<int>list)
	{
		priority_queue<HeapNode>Q;
		for(int i=1;i<=n;i++)d[i]=inf;
		for(int i=0;i<list.size();i++)
		{
			int u=list[i];
			d[u]=0;
			Q.push((HeapNode){0,u});
		}
		memset(done,0,sizeof(done));	
		while(!Q.empty())
		{
			HeapNode x=Q.top();Q.pop();
			int u=x.u;
			if(done[u])continue;
			done[u]=true;
			for(int i=0;i<G[u].size();i++)
			{
				Edge& e=edges[G[u][i]];
				if(d[e.to]>d[u]+e.dist)
				{
					d[e.to]=d[u]+e.dist;
					p[e.to]=G[u][i];
					Q.push((HeapNode){d[e.to],e.to});
				}
			}
		}
	}
}solver;
map<pair<int,int>,int>mp;
vector<int>c[maxn];
int tot;
void init(int n)
{
	for(int i=1;i<=n;i++)
	c[i].clear();
	mp.clear();
	tot=0;
}
int main()
{
	int n,m;
	while(~scanf("%d%d",&n,&m))
	{
		init(n);
		solver.init(maxn-2);
		for(int i=1;i<=m;i++)
		{
			int u,v,ci,ti;
			scanf("%d%d%d%d",&u,&v,&ci,&ti);
			if(!mp[make_pair(u,ci)])
			{
				c[u].push_back(ci);
				mp[make_pair(u,ci)]=++tot;
			}
			if(!mp[make_pair(v,ci)])
			{
				c[v].push_back(ci);
				mp[make_pair(v,ci)]=++tot;
			}
			u=mp[make_pair(u,ci)];
			v=mp[make_pair(v,ci)];
			solver.AddEdge(u,v,ti);
			solver.AddEdge(v,u,ti);
		}
		for(int i=1;i<=n;i++)
		sort(c[i].begin(),c[i].end());
		for(int i=1;i<=n;i++)
		{
			for(int j=0;j<c[i].size()-1;j++)
			{
				int u=mp[make_pair(i,c[i][j])];
				int v=mp[make_pair(i,c[i][j+1])];
				int w=abs(c[i][j]-c[i][j+1]);
				solver.AddEdge(u,v,w);
				solver.AddEdge(v,u,w);
			}
		}
		vector<int>list;
		for(int i=0;i<c[1].size();i++)
		{
			int u=mp[make_pair(1,c[1][i])];
			list.push_back(u);
		}
		solver.dij(list);
		LL ans=inf;
		for(int i=0;i<c[n].size();i++)
		{
			int u=mp[make_pair(n,c[n][i])];
			ans=min(ans,solver.d[u]);
		}
		printf("%lld\n",ans);
	}
}