CCF 201712-4 行车路线（100分）

声明：CCF提交后得100分，但程序存在bug，过段时间修改

试题编号：	201712-4
试题名称：	行车路线
时间限制：	1.0s
内存限制：	256.0MB
问题描述：	问题描述 　　小明和小芳出去乡村玩，小明负责开车，小芳来导航。 　　小芳将可能的道路分为大道和小道。大道比较好走，每走1公里小明会增加1的疲劳度。小道不好走，如果连续走小道，小明的疲劳值会快速增加，连续走s公里小明会增加s^2的疲劳度。 　　例如：有5个路口，1号路口到2号路口为小道，2号路口到3号路口为小道，3号路口到4号路口为大道，4号路口到5号路口为小道，相邻路口之间的距离都是2公里。如果小明从1号路口到5号路口，则总疲劳值为(2+2)^2+2+2^2=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。

问题连接：CCF 201712-4 行车路线

解题思路：使用dijkstra算法，但是每次更新结点的时候要考虑小路，这道题有重边，使用邻接矩阵的话有问题，开始的程序使用邻接矩阵没办法处理重边的问题，因此只得了50分，选择用vector实现邻接表这样就能解决这个问题，程序中sum[i]记录的是最短路径中以i结尾的连续小路的和 ，用这种方法写出的程序也只得了80分，后来发现出现的问题是程序各种变量都用的int类型，可能是溢出了（尤其是小路距离的乘积操作），后来换成long long 类型后得了100分

提交后得100分的C++程序：

#include<iostream>
#include<queue>
#include<cstring>

using namespace std;
typedef long long ll;
ll INF=1e18;
const int N=505;

bool vis[N];//记录结点i是否被访问 
ll dist[N];//记录最短路径
ll sum[N];//sum[i]记录最短路径中以i结尾的连续小路的长度和 

struct Edge{//边 
	int v,t;//到结点v的距离为w ,t是道路的类型 
	ll w;
	Edge(int v,int w,int t):v(v),w(w),t(t){}
}; 
struct Node{
	int u;//结点的编号，源点到结点的最短距离 
	ll w;
	Node(){}
	Node(int u,ll w):u(u),w(w){}
	bool operator<(const Node &a)const//使用优先队列，因此把小于重载为大于含义 
	{
		return w>a.w;
	} 
};
vector<Edge>g[N];//图的邻接表形式 

void dijkstra(int s,int n)
{
	priority_queue<Node>q;
	memset(vis,false,sizeof(vis));
	memset(sum,0,sizeof(sum));
	for(int i=1;i<=n;i++)
	  dist[i]=INF;
	dist[0]=dist[s]=0;
	q.push(Node(s,0));
	while(!q.empty()){
		Node e=q.top();
		q.pop();
		int u=e.u;
		if(!vis[u]){//如果此结点没有被访问 
			vis[u]=true;
			int num=g[u].size();//与结点u相连的结点个数 
			for(int i=0;i<num;i++){//如果结点已经访问了就跳过 
			    int v=g[u][i].v;
				if(vis[v])
				  continue;
				ll cost=0;
				if(g[u][i].t){//如果u到i是小道 
				  ll temp=sum[u]+g[u][i].w;
				  cost=dist[u]-sum[u]*sum[u]+temp*temp;
			    }
				else
				  cost=dist[u]+g[u][i].w;
				if(cost<dist[v]){//如果距离比原来小就更新 
					dist[v]=cost;
					//更新最短路中以i结尾的连续小路长度 
					if(g[u][i].t)
					  sum[v]=sum[u]+g[u][i].w;
					else
					  sum[v]=0;
					q.push(Node(v,cost));
				}
			}
		}
	}
}

int main()
{
	int n,m,t,a,b,c;
	scanf("%d%d",&n,&m);
	while(m--){
		scanf("%d%d%d%d",&t,&a,&b,&c);
		g[a].push_back(Edge(b,c,t));
		g[b].push_back(Edge(a,c,t));
	}
	dijkstra(1,n);
	printf("%lld\n",dist[n]);
	return 0;
}