CCF-CSP 201712-4 行车路线(100分)

问题描述
　　小明和小芳出去乡村玩，小明负责开车，小芳来导航。
　　小芳将可能的道路分为大道和小道。大道比较好走，每走1公里小明会增加1的疲劳度。小道不好走，如果连续走小道，小明的疲劳值会快速增加，连续走s公里小明会增加s2的疲劳度。
　　例如：有5个路口，1号路口到2号路口为小道，2号路口到3号路口为小道，3号路口到4号路口为大道，4号路口到5号路口为小道，相邻路口之间的距离都是2公里。如果小明从1号路口到5号路口，则总疲劳值为(2+2)2+2+22=16+2+4=22。
　　现在小芳拿到了地图，请帮助她规划一个开车的路线，使得按这个路线开车小明的疲劳度最小。
输入格式:
　　输入的第一行包含两个整数n, m，分别表示路口的数量和道路的数量。路口由1至n编号，小明需要开车从1号路口到n号路口。
　　接下来m行描述道路，每行包含四个整数t, a, b, c，表示一条类型为t，连接a与b两个路口，长度为c公里的双向道路。其中t为0表示大道，t为1表示小道。保证1号路口和n号路口是连通的。
输出格式:
　　输出一个整数，表示最优路线下小明的疲劳度。
样例输入
6 7
1 1 2 3
1 2 3 2
0 1 3 30
0 3 4 20
0 4 5 30
1 3 5 6
1 5 6 1
样例输出
76
样例说明
　　从1走小道到2，再走小道到3，疲劳度为52=25；然后从3走大道经过4到达5，疲劳度为20+30=50；最后从5走小道到6，疲劳度为1。总共为76。
数据规模和约定
　　对于30%的评测用例，1 ≤ n ≤ 8，1 ≤ m ≤ 10；
　　对于另外20%的评测用例，不存在小道；
　　对于另外20%的评测用例，所有的小道不相交；
　　对于所有评测用例，1 ≤ n ≤ 500，1 ≤ m ≤ 105，1 ≤ a, b ≤ n，t是0或1，c ≤ 105。保证答案不超过106。

思路：Dijtstra算法，与最短路径不同，我们要将路径转化为疲劳值，方法与最短路径类似，每次选取疲劳值最小的节点，然后更新其与其他相连接的节点的疲劳值，如果相连接道路为大道则直接相加判断即可，如果是小道，则需向前迭代查找小道(这里需要开一个pre[]数组用来存储该节点最短路径上的前驱节点），将总小道长的平方算出来+大道长进行判断，这里选择用优先队列，可以加快查找疲劳值最小的节点，注意小道长的平方可能会爆int，需要开longlong，由于点的个数比较少，所以可以采用邻接矩阵的存储方式，要注意会有重边，在输入时进行处理。

#include<bits/stdc++.h>
#define N 510
#define INF 0x3f3f3f3f
using namespace std;
int n,m,vis[N],pre[N],g[N][N],type[N][N],dist[N];
typedef pair<int,int>PII;
void Dij(int u){
	memset(dist,0x3f,sizeof(dist));
	dist[u]=0;//u到自身的距离为0
	priority_queue<PII,vector<PII>,greater<PII> >heap;//创建小根堆 
	heap.push({0,u});
	while(!heap.empty())
	{
		PII t=heap.top();
		heap.pop();
		if(vis[t.second]==1) continue;//为了避免重复使用
		vis[t.second]=1;
		//cout<<"ok"<<endl;
		for(int j=1;j<=n;j++){
			if(g[t.second][j]!=INF){ 
			    if(type[t.second][j]==0) //大道
				{
					if(t.first+g[t.second][j]<dist[j]){
						dist[j]=t.first+g[t.second][j];
						pre[j]=t.second;
						heap.push({dist[j],j});
						//cout<<"0"<<" j="<<j<<" "<<"dist[j]="<<dist[j]<<endl;
					}
				}
				else //小道
				{ 
				    long long s=g[t.second][j];
				    int p=pre[t.second];
				    int q=t.second;
				    while(type[p][q]==1&&p){//迭代向前寻找小道,直到为大道为止 
				    	s+=g[p][q];
				    	q=p;
				    	p=pre[p];
					}
					if(s*s+dist[q]<dist[j]){
						dist[j]=s*s+dist[q];
						pre[j]=t.second;
						heap.push({dist[j],j});
						//cout<<"1 "<<"j="<<j<<" "<<"dist[j]="<<dist[j]<<endl;
					}
				}
			}
		}
	}
}
int main(){
	scanf("%d%d",&n,&m);
	memset(g,0x3f,sizeof(g));
	for(int i=1;i<=m;i++){
		int t,u,v,w;
		scanf("%d%d%d%d",&t,&u,&v,&w);
		//注意重边的情况,要通过比较选出最优道路 
		if(g[u][v]!=INF){
			if(type[u][v]==0&&t==0){
				g[u][v]=min(w,g[u][v]);
				g[v][u]=min(w,g[v][u]);
			}
			else if(type[u][v]==1&&t==0){
				g[u][v]=g[v][u]=w;
				type[u][v]=type[v][u]=t;
			}
			else if(type[u][v]==1&&t==1){
				if(w<g[u][v]){
					g[u][v]=g[v][u]=w;
				}
			}
		}
		else{
			g[u][v]=g[v][u]=w;
		    type[u][v]=type[v][u]=t;
		}

	}
	Dij(1);
	printf("%d",dist[n]);
	return 0;
}