POJ2449：Remmarguts' Date
——problem:求起点到终点的非严格最短路
——solution：dijkstra+A*
以下转载自：http://www.cppblog.com/linyangfei/archive/2008/07/26/47662.html
【
先说说启发式搜索吧。通常在解决问题的时候，我们需要用到搜索算法，由已知状态推出新的状态，然后检验新的状态是不是就是我们要求的最优解。检验完所有的状态实际上就相当于遍历了一张隐式图。遗憾的是，所有的状态组成的状态空间往往是成指数级别增长的，也就造成了遍历需要用到指数级别的时间，因此，纯粹的暴力搜索，时空效率都比较低。当然，我们在生活中遇到了类似于搜索的问题，我们并不会盲目地去搜寻每一种状态，我们会通过我们的思维，选择一条最接近于目标的路径或者是近似于最短的路径去完成搜索任务。当我们想要计算机去完成这样一项搜索任务的时候，就得让计算机像人一样能够区分尽量短的路径，以便高效地找到最优解。这时可以把计算机看作是一种智能体(agent)可以实现由初始状态向目标状态的转移。
有一种贪心策略，即每一步转移都由计算机选择当前的最优解生成新的状态，一直到达目标状态为止。这样做的时间效率虽然较高，但是贪心的策略只是用到了局部的最优解，并不能保证最后到达目标状态得到的是全局最优解。在能保证全局最优解的范围内，贪心算法还是很有用的。比如说我们熟知的Dijkstra算法求单源最短路。每次选择距离源节点最短距离的待扩展节点进行扩展，最后就能生成源节点到所有节点的最短路径。下面会讲到Dijkstra的扩展，当理解了这个算法之后，我想，你会对Dijkstra有更深入的理解。
这就是A*算法。定义初始状态S，目标状态t，g(s)是由初始状态转移到当前状态s所经过的路径长度，h*(s)是当前状态s距离目标状态t的实际长度，但是一般情况下我们是不知道h*(s)的值的，所以还要定义一个估价函数h(s)，是对h*(s)函数值的下界的估计，也就是有h(s)<=h*(s)，这样需要一个条件，使得由s1生成的每状态s2，都有h(s1)<=h(s2)，这是一个相容的估价函数。再定义f(s)=g(s)+h(s)为启发函数，因为h(s)是单调递增的，所以f(s)也是单调递增的。这样f(s)就估计出了由初始状态的总体代价。A*算法就通过构造这样一个启发函数，将所有的待扩展状态加入到队列里，每次从队列里选择f(s)值最小的状态进行扩展。由于启发函数的作用，使得计算机在进行状态转移的时候尽量避开了不可能产生最优解的分支，而选择相对较接近最优解的路径进行搜索，提高了搜索效率。
讲到这里，可能已经对A*算法的概念有点眉目了。下面我再来做一个比较，就用上面讲到的Dijkstra的例子。Dijkstra算法说的是每次选择距离源点最短距离的点进行扩展。当然可以看做事先将源点到所有节点距离的值保存在一个优先队列里，每次从优先队列里出队最短距离的点扩展，每个点的扩展涉及到要更新队列里所有待扩展节点的距离值，每个节点只能进队一次，就需要有一个表来记录每个节点的入队次数（就是算法中用到的标记数组）。将Dijkstra求最短路的方法扩展，这道题目要求的是两点间第k短路。类比于Dijkstra算法可以首先确定下面几个搜索策略：
1、用优先队列保存节点进行搜索。
2、放开每个节点的入队次数，求k短路，每个节点可以入队k次。
首先看第一个策略，在A*算法中用优先队列就是要用到启发函数f(s)确定状态在优先队列里面的优先级。其实Dijkstra用到的优先队列实际上就是估价函数值为0，启发函数f(s)=g(s)，即是选取到源点距离最近的点进行扩展。因为h(s)=0满足了估价函数相容这个条件。这题求k短路就不能单纯的使用h(s)=0这个估价函数。解决这道题的时候选取h(x)=dt(x), dt(x)是x节点到目标节点的最短距离。最短距离可以开始由Dijkstra直接求得。
再看第二个策略，控制每个节点的入队（或出队）次数为k次，可以找到第k短路径。可能这样想有点主观的套用，那么我就先来证明这样一个结论：
如果x是s到t的第k短路径上的一个节点，那么由这条路径s到x是s到x的第m短路径，则不可能有m>k。用反证法很容易得出：如果这条路径是s到x的第m短路径，如果m>k，那么经过x到t的路径就有m-1条比当前路径要短，不符合当前路径是s到t的第k短路径。
】
注意：
1、A*求的是非严格最短路，即一个点可以在路径中出现多次。
2、另外本题是有向边，因此在用DIJKSTRA的时候需要存下反向边。
View Code
1 #include<stdio.h>
2 #include<queue>
3 #include<memory.h>
4 using namespace std;
5 #define oo 1000000000
6 #define V 1005
7 #define E 200010
8 int i, j, ans, u, v, l, n, m, s, t, num, k;
9 int dist[V], head[V], nxt[E], ev[E], ew[E], vis[V], head_r[V], nxt_r[E],
10 ev_r[E], ew_r[E];
11 struct node
12 {
13 int u, l;
14 friend bool operator <(node a, node b)
15 {
16 return a.l > b.l;
17 }
18 node(int a = 0, int b = 0, int c = 0, int d = 0) :
19 u(a), l(b)
20 {
21 }
22 } x, y;
23 priority_queue<node> q;
24 void add_edge(int u, int v, int l)
25 {
26 nxt[++num] = head[u];
27 head[u] = num;
28 ev[num] = v;
29 ew[num] = l;
30 }
31 void add_r_edge(int u, int v, int l)
32 {
33 nxt_r[++num] = head_r[u];
34 head_r[u] = num;
35 ev_r[num] = v;
36 ew_r[num] = l;
37 }
38 void dijkstra()
39 {
40 int i;
41 q = priority_queue<node> ();
42 memset(vis, 0, sizeof(vis));
43 for (i = 1; i <= n; i++)
44 dist[i] = oo;
45 dist[t] = 0;
46 q.push(node(t, 0));
47 while (!q.empty())
48 {
49 x = q.top();
50 q.pop();
51 if (vis[x.u] == 0)
52 {
53 for (i = head_r[x.u]; i; i = nxt_r[i])
54 {
55 y.u = ev_r[i];
56 y.l = x.l + ew_r[i];
57 if (y.l < dist[y.u])
58 {
59 dist[y.u] = y.l;
60 q.push(y);
61 }
62 }
63 }
64 vis[x.u] = 1;
65 }
66 }
67 int A_star()
68 {
69 int i;
70 q = priority_queue<node> ();
71 memset(vis, 0, sizeof(vis));
72 if (dist[s] == oo)
73 return -1;
74 q.push(node(s, dist[s]));
75 while (!q.empty())
76 {
77 x = q.top();
78 q.pop();
79 vis[x.u]++;
80 if (x.u == t && vis[x.u] == k)
81 return x.l;
82 if (vis[x.u] <= k)
83 for (i = head[x.u]; i; i = nxt[i])
84 {
85 y.u = ev[i];
86 y.l = x.l - dist[x.u] + ew[i] + dist[y.u];
87 q.push(y);
88 }
89 }
90 return -1;
91 }
92 int main()
93 {
94 scanf("%d%d", &n, &m);
95 num = 0;
96 memset(head, 0, sizeof(head));
97 memset(nxt, 0, sizeof(nxt));
98 for (i = 0; i < m; i++)
99 {
100 scanf("%d%d%d", &u, &v, &l);
101 add_edge(u, v, l);
102 add_r_edge(v, u, l);
103 }
104 scanf("%d%d%d", &s, &t, &k);
105 if (s == t)
106 k++;
107 dijkstra();
108 ans = A_star();
109 printf("%d\n", ans);
110 return 0;
111 }
