CCF认证 2017-12 行车路线

　　小明和小芳出去乡村玩，小明负责开车，小芳来导航。
　　小芳将可能的道路分为大道和小道。大道比较好走，每走1公里小明会增加1的疲劳度。小道不好走，如果连续走小道，小明的疲劳值会快速增加，连续走 s公里小明会增加 s ²的疲劳度。
　　例如：有5个路口，1号路口到2号路口为小道，2号路口到3号路口为小道，3号路口到4号路口为大道，4号路口到5号路口为小道，相邻路口之间的距离都是2公里。如果小明从1号路口到5号路口，则总疲劳值为(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，疲劳度为5 ²=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 ≤ 10⁵，1 ≤ a, b ≤ n，t是0或1，c ≤ 10⁵。保证答案不超过10⁶。

Dij算法变形 

在模板算权值的时候 加入小道判断 如果是小道 则diss = dist[u] - small*small + (small + len)*(small + len)

small为 通过small小道长度走到点u，len为 u到v的距离

#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#include<cmath>
using namespace std;
const int maxn = 20005;
const int maxm = 200000;
const int inf = 0x3f3f3f3f;
int tot, head[maxn], n;
struct edge
{
	int v, len, next,size;
}e[maxm * 2];

void init()
{
	tot = 0;
	memset(head, -1, sizeof(head));
}

void addedge(int u, int v, int len,int size)
{
	e[tot] = { v,len,head[u],size };
	head[u] = tot++;
}

struct node
{
	int p, num;
}a[20005];

struct dnode
{
	long long p, len, small;
	bool operator<(const dnode &T) const
	{
		return len>T.len;
	}
};
bool vis[maxn];
long long dist[maxn];
long long dij(int st, int en)
{
	memset(vis, 0, sizeof(vis));
	for (int i = 0; i <= n; i++)
		dist[i] = inf;
	priority_queue<dnode> que;
	while (!que.empty())
		que.pop();
	que.push({ st,0,0 });
	dist[st] = 0;
	int ans = -1;
	while (!que.empty())
	{
		long long u = que.top().p,small = que.top().small;
		if (u == en) {
			ans = en;
			break;
		}
		que.pop();
		if (vis[u])
			continue;
		vis[u] = 1;
		for (int i = head[u]; i != -1; i = e[i].next)
		{
			long long v = e[i].v;
			long long len = e[i].len;
			if (e[i].size == 0) {
				if (!vis[v] && dist[v] > dist[u] + len)
				{
					dist[v] = dist[u] + len;
					que.push({ v,dist[v],0 });
				}
			}
			//和模板不同在于多了 else中的判断小道还是大路
			else {
				long long diss = dist[u] - small*small + (small + len)*(small + len);
				if (!vis[v] && dist[v] > diss) {
					dist[v] = diss;
					que.push({ v,dist[v],small+len });
				}
			}
		}
	}
	if (ans == en)return dist[ans];
	else return -1;
}
int main()
{
	int  m, u, v, d, a;
	while (cin >> n >> m)
	{
		init();
		for (int i = 1; i <= m; i++)
		{
			scanf("%d%d%d%d", &a, &u, &v, &d);
			addedge(u, v, d, a);
			addedge(v, u, d, a);
		}
		cout << dij(1, n) << endl;
	}
	return 0;
}