POJ 2449 求第K短路
第一道第K短路的题目 QAQ
拿裸的DIJKSTRA + 不断扩展的A* 给2000MS过了
题意:大意是 有N个station 要求从s点到t点 的第k短路 (不过我看题意说的好像是从t到s 可能是出题人写错了)
从这题中还真的学到了很多
1.第k短路的算法 A* 还有用边表实现dij
(注:以下部份资料来源于网上)
所谓A*就是启发是搜索 说白了就是给搜索一个顺序使得搜索更加合理减少无谓的搜索. 如何来确定搜索的顺序?..也就是用一个值来表示 这个值为f[n]..每次搜索取f[x]最小的拓展 那么这个f[n]=h[n]+g[n]
其中f(n) 是节点n的估价函数,g(n)是在状态空间中从初始节点到n节点的实际代价,h(n)是从n到目标节点最佳路径的估计代价。在这里主要是h(n)体现了搜索的启发信息,因为g(n)是已知的。如果说详细 点,g(n)代表了搜索的广度的优先趋势。但是当h(n) >> g(n)时,可以省略g(n),而提高效率。
A*算法的估价函数可表示为:
f’(n) = g’(n) + h’(n)
这里,f’(n)是估价函数,g’(n)是起点到终点的最短路径值,h’(n)是n到目标的最短路经的启发值。由 于这个f’(n)其实是无法预先知道的,所以我们用前面的估价函数f(n)做近似。g(n)代替g’(n),但 g(n)>=g’(n) 才可(大多数情况下都是满足的,可以不用考虑),h(n)代替h’(n),但h(n)<=h’(n)才可(这一点特别的重 要)。可以证明应用这样的估价函数是可以找到最短路径的,也就是可采纳的。我们说应用这种估价函数的 最好优先算法就是
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <iostream>
#include <stack>
#include <queue>
#include <algorithm> #define ll long long
using namespace std;
const int INF = 0x3f3f3f3f;
const int MAXN = ; struct vertex{
int sum, h, pos;
bool operator < (vertex a) const{
return a.sum + a.h < sum + h;
}
};
struct sc{
int u, v, w, next;
}line1[MAXN*MAXN],line2[MAXN*MAXN]; int link1[MAXN],link2[MAXN],h[MAXN],times[];
int n, m, s, e, k;
bool vis[MAXN];
priority_queue <vertex> que; void init(){
memset(link1, , sizeof(link1));
memset(link2, , sizeof(link2));
memset(vis,,sizeof(vis));
memset(h,0x3f,sizeof(h));//h函数初始值为最大
while (!que.empty()) que.pop();
memset(times,,sizeof(times));
} void djikstra(){
int i,k,p;
h[e] = ;
for (p = ; p <= n; ++p){
k = ;
for (i = ; i <= n; ++i)
if (!vis[i] && (!k||h[i]<h[k])) k = i;
vis[k] = true;
k = link2[k];
while (k){
if (h[line2[k].v] > h[line2[k].u] + line2[k].w)
h[line2[k].v] = h[line2[k].u] + line2[k].w;
k = line2[k].next;
}
}
} int Astar(){
int t;
vertex g,temp;
g.pos = s;
g.sum = ;
g.h = h[s];
que.push(g); if (s==e) ++k;
while (!que.empty()){
g = que.top();//每次取估价函数值最小的节点
que.pop();
++times[g.pos];
if (times[g.pos] == k && g.pos == e) return g.sum + g.h;
if (times[g.pos] > k) continue; t = link1[g.pos];
while (t){//扩展,并把其加入优先队列即openlist
temp.sum = g.sum + line1[t].w;
temp.h = h[line1[t].v];
temp.pos = line1[t].v;
que.push(temp);
t = line1[t].next;
}
}
return -;
} int main(){
int i, j;
while(cin >> n >> m){
init();
for (i = ; i <= m; ++i){
cin >> line1[i].u >> line1[i].v >> line1[i].w;
line1[i].next = link1[line1[i].u];//记录与节点u有直接边的节点
link1[line1[i].u] = i; line2[i].u = line1[i].v;
line2[i].v = line1[i].u;
line2[i].w = line1[i].w;
line2[i].next = link2[line2[i].u];
link2[line2[i].u] = i;
}
cin >> s >> e >> k;
djikstra();
cout << Astar() << endl;
}
return ;
}
POJ 2449 求第K短路的更多相关文章
- POJ 2449 Remmarguts' Date (K短路 A*算法)
题目链接 Description "Good man never makes girls wait or breaks an appointment!" said the mand ...
- poj 2449 Remmarguts' Date K短路+A*
题目链接:http://poj.org/problem?id=2449 "Good man never makes girls wait or breaks an appointment!& ...
- POJ 2449 Remmarguts' Date --K短路
题意就是要求第K短的路的长度(S->T). 对于K短路,朴素想法是bfs,使用优先队列从源点s进行bfs,当第K次遍历到T的时候,就是K短路的长度. 但是这种方法效率太低,会扩展出很多状态,所以 ...
- POJ 2449 Dijstra + A* K短路
这题一开始的思路应该是直接从源点进行BFS搜索K短路. 但这样的复杂度在点数和K的值增大后将会变得很大. 而A*算法则构造一个h(x),在进行BFS时,每次都抛出最小的h(x)从而使汇点的出队速度加快 ...
- poj 2449(A*求第K短路)
题目链接:http://poj.org/problem?id=2449 思路:我们可以定义g[x]为源点到当前点的距离,h[x]为当前点到目标节点的最短距离,显然有h[x]<=h*[x](h*[ ...
- A*算法的认识与求第K短路模板
现在来了解A*算法是什么 现在来解决A*求K短路问题 在一个有权图中,从起点到终点最短的路径成为最短路,第2短的路成为次短路,第3短的路成为第3短路,依此类推,第k短的路成为第k短路.那么,第k短路怎 ...
- POJ 2499 A*求第K短路
DES就是给你一个图.然后给你起点和终点.问你从起点到终点的第K短路. 第一次接触A*算法. 题目链接:Remmarguts' Date 转载:http://blog.csdn.net/mbxc816 ...
- POJ 2449Remmarguts' Date 第K短路
Remmarguts' Date Time Limit: 4000MS Memory Limit: 65536K Total Submissions: 29625 Accepted: 8034 ...
- aStar算法求第k短路
A*的概念主意在于估计函数,f(n)=g(n)+h(n),f(n)是估计函数,g(n)是n节点的当前代价,h(n)是n节点的估计代价:而实际中,存在最优的估计函数f'(n)=g'(n)+h'(n),那 ...
随机推荐
- Android UiAutomator 自动化测试编译运行---新手2
1.首先打开eclipse创建java项目
- (Problem 41)Pandigital prime
We shall say that an n-digit number is pandigital if it makes use of all the digits 1 to n exactly o ...
- AngularJS Factory Service Provider
先看看http://www.cnblogs.com/mbydzyr/p/3460501.html http://www.oschina.net/translate/angularjs-factory- ...
- (asp.net MVC学习)System.Web.Mvc.UrlHelper的学习与使用
上一次学习了HtmlHelper帮助类,这次我们学习一下UrlHelper帮 助类,看类名也都知道这个类是用来帮我们生成URL在ASP.NET MVC应用程序中.让我们来看看该类给我们带来了哪些方便的 ...
- poj 3630 Phone List(字典树)
题目链接: http://poj.org/problem?id=3630 思路分析: 求在字符串中是否存在某个字符串为另一字符串的前缀: 即对于某个字符串而言,其是否为某个字符串的前缀,或存在某个其先 ...
- vs2010经常使用快捷键
调试快捷键 F6: 生成解决方式 Ctrl+F6: 生成当前项目 F7: 查看代码 Shift+F7: 查看窗口设计器 F5: 启动调试 Ctrl+F5: 開始运行(不调试) Shift+F5: 停止 ...
- TCP传输连接建立与释放详解
一直以来有许多读者朋友对TCP的传输连接建立和释放过程不是很理解,而这又是几乎网络认证中必考的知识点,包括软考.CCNA\CCNP.H3CNA\H3CNE等,为此再把笔者年度巨作,广受好评的——< ...
- openstack之cinder
一.cinder概述: 1.F版之前,并没有cinder,对应的组件为nova-volume:以Rest API的形式提供服务 2.cinder目标: 减少nova的复杂性,降低nova的负载,支持多 ...
- for语句及switch case用法示例
for(初始条件:循环条件:状态改变){ 循环体:执行代码} 示例一: 示例二: 示例三: 示例四: 示例五: 示例六: 示例七: switch case: 示例一: 示例二: 示例三:
- js得到分页栏
自己写的,感觉返回html代码蠢蠢的,但是新手并不知道怎么写更好的,感觉这样子也蛮简单.记录下来,以后来越改越好. //获得分页栏.注意indexSize为奇数,这样也比较好看 //totalNum: ...