hdu 4179 Difficult Routes (SP)

Problem - 4179

　　坑了我一个晚上的SP题。

　　题意是，给出若干空间中的点，给出其中某些点之间是有直线线段路径相连的。要求求出一条从s开始到t结束的路径，它的难度是d。难度的计算是空间线段两点的高度差乘以100再除以投影到xOy平面上线段的长度。难度是d的路径的定义是路径中，经过的线段的难度最大是d。

　　其实比较容易可以想到的是枚举难度最大的线段的，同时构建难度小于等于给定值的正图和反图。正图是用来求出起点s到各个顶点的最短距离，因为是有向图，所以求出t到各个顶点的最短距离是要用到反图来求的。然后就将两条路径拼接到被枚举的线段上。

　　中间错的好多的是，写之前还想到要用反图求t到各个顶点的最短距离，结果写的时候就忘记了，最后多得队友ly的提醒！然后就是比较SB的一个行为了。明明就是直接拼接上去就好了，结果我想太多，没有考虑到我的图是有向图，搞多了两个比较，于是就狠狠的被坑了一个晚上了！！！_(:з」∠)_

代码如下：

 #include <cstdio>
 #include <cmath>
 #include <iostream>
 #include <algorithm>
 #include <cstring>

 using namespace std;

 const int N = ;
 const int M = ;
 int eh[][N], ec[];
 struct Edge {
     int id, nx, df;
     double l;
     Edge() {}
     Edge(int id, int nx, int df, double l) : id(id), nx(nx), df(df), l(l) {}
 } edge[][M << ];

 inline void init() {
     memset(eh, -, sizeof(eh));
     memset(ec, , sizeof(ec));
 }

 inline void addedge(int u, int v, int df, double l, int id) {
     edge[id][ec[id]] = Edge(v, eh[id][u], df, l);
     eh[id][u] = ec[id]++;
 }

 int x[N], y[N], z[N];
 template<class T> inline T sqr(T a) { return a * a;}
 int st[M << ], ed[M << ];

 inline int getdf(int a, int b) {
     if (z[a] >= z[b]) return ;
     double len = sqrt(sqr((double) x[a] - x[b]) + sqr((double) y[a] - y[b]));
     return (int) floor(100.0 * (z[b] - z[a]) / len);
 }

 inline double getdis(int a , int b) { return sqrt(sqr((double) x[a] - x[b]) + sqr((double) y[a] - y[b]) + sqr((double) z[a] - z[b]));}
 const double FINF = 1e50;
 double dis[][N];
 int q[M << ];
 bool vis[N];

 void spfa(int s, int id) {
     for (int i = ; i < N; i++) dis[id][i] = FINF;
     int qh, qt;
     qh = qt = ;
     dis[id][s] = ;
     q[qt++] = s;
     vis[s] = true;
     while (qh < qt) {
         int cur = q[qh++];
         vis[cur] = false;
         for (int i = eh[id][cur]; ~i; i = edge[id][i].nx) {
             if (dis[id][edge[id][i].id] > dis[id][cur] + edge[id][i].l) {
                 dis[id][edge[id][i].id] = dis[id][cur] + edge[id][i].l;
                 if (!vis[edge[id][i].id]) q[qt++] = edge[id][i].id, vis[edge[id][i].id] = true;
             }
         }
     }
 }

 int main() {
 //    freopen("in", "r", stdin);
     int n, m;
     while (~scanf("%d%d", &n, &m) && (n || m)) {
         for (int i = ; i <= n; i++) scanf("%d%d%d", &x[i], &y[i], &z[i]);
         for (int i = ; i < m; i++) scanf("%d%d", &st[i], &ed[i]);
         int A, B, D;
         scanf("%d%d%d", &A, &B, &D);
         init();
         for (int i = ; i < m; i++) {
             double tmp = getdis(st[i], ed[i]);
             int df = getdf(st[i], ed[i]);
             if (df <= D) {
                 addedge(st[i], ed[i], df, tmp, );
                 addedge(ed[i], st[i], df, tmp, );
             }
             df = getdf(ed[i], st[i]);
             if (df <= D) {
                 addedge(ed[i], st[i], df, tmp, );
                 addedge(st[i], ed[i], df, tmp, );
             }
         }
         spfa(A, );
         spfa(B, );
 //        for (int i = 1; i <= n; i++) cout << dis[0][i] << ' '; cout << endl;
 //        for (int i = 1; i <= n; i++) cout << dis[1][i] << ' '; cout << endl;
         double ans = FINF;
         for (int i = ; i <= n; i++) {
             for (int t = eh[][i]; ~t; t = edge[][t].nx) {
                 if (edge[][t].df == D) ans = min(ans, dis[][i] + edge[][t].l + dis[][edge[][t].id]);
             }
         }
         if (ans < FINF) printf("%.1f\n", ans);
         else puts("None");
     }
     return ;
 }

　　再搞搞几何就得开图论来练了！不然对于这种SB题目都毫无自己做错的意识，相当蛋疼的感觉。。_(:з」∠)_

——written by Lyon