问题描述
小明和小芳出去乡村玩，小明负责开车，小芳来导航。
小芳将可能的道路分为大道和小道。大道比较好走，每走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。
现在小芳拿到了地图，请帮助她规划一个开车的路线，使得按这个路线开车小明的疲劳度最小。
小芳将可能的道路分为大道和小道。大道比较好走，每走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号路口是连通的。
接下来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
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。
对于另外20%的评测用例，不存在小道；
对于另外20%的评测用例，所有的小道不相交；
对于所有评测用例，1 ≤ n ≤ 500，1 ≤ m ≤ 105，1 ≤ a, b ≤ n，t是0或1，c ≤ 105。保证答案不超过106。
做法：把大路和小路分两个图算，先对小路跑了N次dij然后把这些小路边加入到大路边里面，为了区分二者我对小路边取了个相反数。然后最短路分两个，d[i][0/1]表示上一条边走的是0(大路)还是1(小路).
这个图边有点多所以T了，只有80分，懒得改了。
改的话我想的是别用堆优化了直接用O(N^2)的dij试试.
1 #include<bits/stdc++.h> 2 using namespace std; 3 #define pii pair<int,int> 4 #define mp make_pair 5 #define inf 0x3f3f3f3f 6 #define piii pair<int,pii> 7 const int maxn=505; 8 9 vector<pii>e1[maxn],e2[maxn]; 10 int dis1[maxn][maxn],dis2[maxn][maxn]; 11 int N,M; 12 bool vis[maxn]; 13 int d[maxn]; 14 void dij(int s,vector<pii> (&e)[maxn]){//值传递会怎样 15 memset(vis,0,sizeof(vis)); 16 memset(d,inf,sizeof(d)); 17 priority_queue<pii,vector<pii>,greater<pii> >q; 18 q.push(mp(0,s)); 19 d[s]=0; 20 while(!q.empty()){ 21 int u=q.top().second;q.pop(); 22 if(vis[u])continue; 23 vis[u]=1; 24 for(pii p:e[u]){ 25 int v=p.first,w=p.second; 26 if(d[v]>d[u]+w){ 27 d[v]=d[u]+w; 28 q.push(mp(d[v],v)); 29 } 30 } 31 } 32 } 33 int d2[maxn][2]; 34 bool vis2[maxn][2]; 35 void dij2(){ 36 memset(vis2,0,sizeof(vis2)); 37 memset(d2,inf,sizeof(d2)); 38 d2[1][0]=0; 39 priority_queue<piii,vector<piii>,greater<piii> >q; 40 q.push(mp(0,mp(1,0))); 41 while(!q.empty()){ 42 piii o=q.top();q.pop(); 43 int u=o.second.first; 44 int type=o.second.second; 45 if(vis2[u][type])continue; 46 vis2[u][type]=1; 47 for(pii p:e1[u]){ 48 int v=p.first,w=p.second; 49 if(w>=0){ 50 if(d2[v][0]>d2[u][type]+w){ 51 d2[v][0]=d2[u][type]+w; 52 q.push(mp(d2[v][0],mp(v,0))); 53 } 54 }else if(type==0){ 55 w*=-1; 56 if(d2[v][1]>d2[u][type]+w){ 57 d2[v][1]=d2[u][type]+w; 58 q.push(mp(d2[v][1],mp(v,1))); 59 } 60 } 61 } 62 } 63 } 64 void spfa(){} 65 int main(){ 66 int a,b,c,t; 67 ios::sync_with_stdio(false); 68 cin>>N>>M; 69 for(int i=1;i<=N;++i){ 70 for(int j=i;j<=N;++j){ 71 if(i==j)dis1[i][j]=dis2[i][j]=0; 72 else{ 73 dis1[i][j]=dis1[j][i]=inf; 74 dis2[i][j]=dis2[j][i]=inf; 75 } 76 } 77 } 78 while(M--){ 79 cin>>t>>a>>b>>c; 80 if(t==0){ 81 e1[a].push_back(mp(b,c)); 82 e1[b].push_back(mp(a,c)); 83 }else if(t==1){ 84 if(dis2[a][b]>c)dis2[a][b]=dis2[b][a]=c; 85 e2[a].push_back(mp(b,c)); 86 e2[b].push_back(mp(a,c)); 87 } 88 } 89 for(int i=1;i<=N;++i){ 90 dij(i,e2); 91 for(int j=1;j<=N;++j){ 92 if(i!=j && d[j]!=inf && d[j]<=1005){ 93 e1[i].push_back(mp(j,d[j]*d[j]*-1)); 94 } 95 } 96 } 97 dij2(); 98 cout<<min(d2[N][0],d2[N][1])<<endl; 99 return 0; 100 }