CCF认证 201712-4 行车路线
思路：好难啊，去学习大佬的解法了，结果吭吭哧哧的还出现很多不bug。首先用的是是spfa的算法。但是由于处理小路的时候用来floyd,所以这个时间复杂度也不小，当我把floy的j=i+1改成1开始就超时了。g0[i][j]保存着小路i-j的总长 然后分大路小路两种情况更新最短路。（参考博客）
问题描述
小明和小芳出去乡村玩，小明负责开车，小芳来导航。
小芳将可能的道路分为大道和小道。大道比较好走，每走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。
#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn=500+10;
typedef long long ll;
const ll inf=1e18;
bool vis[maxn];
int que[maxn];
ll dis[maxn],dis0[maxn];
ll g[maxn][maxn],g0[maxn][maxn];
//81
void spfa(int start,int n){
int front=0,rear=0;//手动队列
for(int i=1;i<=n;i++)
{
if(i==start)
{
que[rear++]=i;
vis[i]=true;
dis[i]=dis0[i]=0;
}
else
{
vis[i]=false;dis[i]=dis0[i]=inf;
}
}
while(front!=rear)
{
int u=que[front++];//cout<<u<<endl;
vis[u]=false;//cout<<u<<" "<<dis[u]<<endl;
if(front>=maxn) front=0;//front%=maxn;//超过限度 front=front%maxn
for(int i=1;i<=n;i++)
{
ll v=g[u][i];//大陆的花费
if(dis[i]>dis[u]+v)
{
dis[i]=dis[u]+v;
if(!vis[i])
{
vis[i]=true;
que[rear++]=i;
if(rear>=maxn) rear=0;//rear%=maxn;//rear%=maxn
}
}
if(dis[i]>dis0[u]+v)
{
dis[i]=dis0[u]+v;
if(!vis[i])
{
vis[i]=true;
que[rear++]=i;
if(rear>=maxn) rear=0;//rear%=maxn;
}
}
if(g0[u][i]!=inf)//可由小路到达
{
v=g0[u][i]*g0[u][i];
if(dis0[i]>dis[u]+v)//走小路，前驱为小路时不能走小路
{
dis0[i]=dis[u]+v;
if(!vis[i])
{
vis[i]=true;
que[rear++]=i;
if(rear>=maxn) rear=0;
}
}
}
}
}
}
//zuiduan
int main()
{
ll a,b,c,d;
int n,m;
scanf("%d%d",&n,&m);
for(int i=0;i<=n;i++)
for(int j=0;j<=n;j++) g[i][j]=g0[i][j]=inf;
for(int i=0;i<m;i++)
{
scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
if(a&&g0[b][c]>d) g0[b][c]=g0[c][b]=d;//小路
if(!a&&g[b][c]>d) g[b][c]=g[c][b]=d;//大路
//cout<<g[b][c]<<" "<<g0[b][c]<<endl;
}
for(int k=1;k<=n;k++)
{
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
{
if(k==i||j==k) continue;
if(g0[i][j]>g0[i][k]+g0[k][j]) g0[i][j]=g0[i][k]+g0[k][j];
}
}
spfa(1,n);
printf("%lld\n",min(dis[n],dis0[n]));
return 0;
}