HDU6181-求解次短路-A*或者dijkstra
(有任何问题欢迎留言或私聊 && 欢迎交流讨论哦
题意:传送门
原题目描述在最下面。
次短路裸题。
思路:
在dijstra的过程过维护两个数组变量:\(dis1[],dis2[]\)。分别表示到达\(i\)节点的最短距离和次短距离。
如果当前距离大于次短距离,则更新次短距离。如果更新后次短距离小于最短距离,就交换二者的值。到最后两个数组值就更新就完成了。
AC代码:
dijkstra
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<string>
#include<cmath>
#define test printf("***\n")
#define ka getchar();getchar()
#define ka1 getchar()
#define iis std::ios::sync_with_stdio(false)
using namespace std;
typedef long long LL;
const int N = 100005;
const int INF = 0x3f3f3f3f;
const LL mod = 10000;
const double eps = 1e-8;
struct lp{
int to;
LL w;
lp(int a,LL b){to=a;w=b;}
bool operator <(const lp &a)const {
if(w!=a.w) return w>a.w;
return to<a.to;
}
};
vector<lp>mp[N];
LL dis1[N],dis2[N];
int n,m;
void init(){
for(int i=0;i<=n;++i)mp[i].clear();
}
void dij(int s){
memset(dis1,0x3f,sizeof(dis1));
memset(dis2,0x3f,sizeof(dis2));
dis1[s]=0;
priority_queue<lp>Q;
Q.push(lp(s,0));
while(!Q.empty() ){
lp x=Q.top();Q.pop();
int u=x.to;
if(dis2[u]<x.w)continue;
for(int i=0;i<mp[u].size();++i){
lp y=mp[u][i];
if(dis2[y.to]>(x.w+y.w)){
dis2[y.to]=x.w+y.w;
Q.push(lp(y.to,dis2[y.to]));
}
if(dis2[y.to]<dis1[y.to]){
swap(dis2[y.to],dis1[y.to]);
}
}
}
printf("%lld\n",dis2[n]);
}
int main(){
int t;
scanf("%d",&t);
while(t--){
int a,b,c;
scanf("%d%d",&n,&m);
init();
for(int i=0;i<m;++i){
scanf("%d%d%d",&a,&b,&c);
mp[a].push_back(lp(b,c));
mp[b].push_back(lp(a,c));
}
dij(1);
}
return 0;
}
A*
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<assert.h>
#include<bitset>
#include<vector>
#include<queue>
#define lson rt<<1
#define rson rt<<1|1
#define lowbit(x) (x)&(-(x))
#define all(x) (x).begin(),(x).end()
#define mk make_pair
#define pb push_back
#define fi first
#define se second
using namespace std;
typedef long long LL;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int N = (int)1e5 +107;
int n, m, k, st, ed;
int vis[N],time[N];
LL dis[N];
struct lp{
LL f,g;
int v;
friend bool operator <(const lp &a,const lp &b){
if(a.f==b.f)return a.g>b.g;
return a.f>b.f;
}
}aa,bb;
struct lh{
int v,nex;
LL w;
}cw[200000+5],rev[200000+5];
int head[N],tot,headd[N],tum;
int q[2500005];
void add(int u,int v,LL w){
cw[++tot].v=v;cw[tot].nex=head[u];cw[tot].w=w;
head[u]=tot;
rev[++tum].v=u;rev[tum].nex=headd[v];rev[tum].w=w;
headd[v]=tum;
}
void spfa(){
for(int i = 1; i <= n; i++) dis[i] = INF;
memset(vis, 0, sizeof(vis));
int h = 0, t = 1;
q[0] = ed;
dis[ed] = 0;
while(h < t){
int u = q[h++];
vis[u] = 0;
for(int i = headd[u] ; ~i ; i = rev[i].nex){
int v = rev[i].v;
LL w = rev[i].w;
if(dis[v] > dis[u] + w){
dis[v] = dis[u] + w;
if(!vis[v]){
q[t++] = v;
vis[v] = 1;
}
}
}
}
}
LL Astar(){
if(dis[st]==INF)return -1;
memset(time,0,sizeof(time));
aa.v=st;aa.f=dis[st];aa.g=0;
priority_queue<lp>Q;
Q.push(aa);
while(!Q.empty()){
bb = Q.top();Q.pop();
int u = bb.v;
time[u]++;
if(time[u]==k&&u==ed)return bb.g;
if(time[u]>k)continue;
for(int i=head[u];~i;i=cw[i].nex){
int v = cw[i].v;
aa.v = v;
aa.g=bb.g+cw[i].w;
aa.f=aa.g+dis[v];
Q.push(aa);
}
}
return -1;
}
int main(){
int tim;
scanf("%d",&tim);
while(tim--){
scanf("%d%d", &n, &m);
tot=tum=-1;
memset(head,-1,sizeof(head));
memset(headd,-1,sizeof(headd));
for(int i = 0, u, v; i < m; ++i){
LL w;
scanf("%d%d%lld", &u, &v, &w);
add(u,v,w);add(v,u,w);
}
st=1;ed=n;
k=2;
spfa();
if(st == ed)k++;
printf("%lld\n", Astar());
}
return 0;
}
####原题目描述:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<string>
#include<cmath>
#define test printf("***\n")
#define ka getchar();getchar()
#define ka1 getchar()
#define iis std::ios::sync_with_stdio(false)
using namespace std;
typedef long long LL;
const int N = 100005;
const int INF = 0x3f3f3f3f;
const LL mod = 10000;
const double eps = 1e-8;
struct lp{
int to;
LL w;
lp(int a,LL b){to=a;w=b;}
bool operator <(const lp &a)const {
if(w!=a.w) return w>a.w;
return to<a.to;
}
};
vector<lp>mp[N];
LL dis1[N],dis2[N];
int n, m, pre[N];
void init(){
for(int i=0;i<=n;++i)mp[i].clear();
}
void dij(int s,LL *dis){
dis[s]=0;
priority_queue<lp>Q;
Q.push(lp(s,0));
while(!Q.empty() ){
lp x=Q.top();Q.pop();
int u=x.to;
for(int i=0;i<mp[u].size();++i){
lp y=mp[u][i];
if(dis[y.to]>(x.w+y.w)){
dis[y.to]=x.w+y.w;
pre[y.to]=u;
Q.push(lp(y.to,dis[y.to]));
}
}
}
}
int main(){
int t;
scanf("%d",&t);
while(t--){
int a,b,c;
scanf("%d%d",&n,&m);
init();
for(int i=0;i<m;++i){
scanf("%d%d%d",&a,&b,&c);
mp[a].push_back(lp(b,c));
mp[b].push_back(lp(a,c));
}
memset(dis1,0x3f,sizeof(dis1));
memset(dis2,0x3f,sizeof(dis2));
dij(n,dis2);
memset(pre,0,sizeof(pre));
dij(1,dis1);
LL ans = 1e18;
for(int i=1;i<=n;++i){
for(int j=0;j<mp[i].size();++j){
int v= mp[i][j].to;
LL tmp =dis1[i]+dis2[v]+mp[i][j].w;
if(pre[v]!=i){
if(tmp<ans){
ans=tmp;
}
}
}
}
printf("%lld\n", ans);
}
return 0;
}
HDU6181-求解次短路-A*或者dijkstra的更多相关文章
- ACM-最短路(SPFA,Dijkstra,Floyd)之最短路——hdu2544
***************************************转载请注明出处:http://blog.csdn.net/lttree************************** ...
- Bellman-ford算法、SPFA算法求解最短路模板
Bellman-ford 算法适用于含有负权边的最短路求解,复杂度是O( VE ),其原理是依次对每条边进行松弛操作,重复这个操作E-1次后则一定得到最短路,如果还能继续松弛,则有负环.这是因为最长的 ...
- ACM - 最短路 - AcWing 849 Dijkstra求最短路 I
AcWing 849 Dijkstra求最短路 I 题解 以此题为例介绍一下图论中的最短路算法.先让我们考虑以下问题: 给定一个 \(n\) 个点 \(m\) 条边的有向图(无向图),图中可能存在重边 ...
- 最短路模板(Dijkstra & Dijkstra算法+堆优化 & bellman_ford & 单源最短路SPFA)
关于几个的区别和联系:http://www.cnblogs.com/zswbky/p/5432353.html d.每组的第一行是三个整数T,S和D,表示有T条路,和草儿家相邻的城市的有S个(草儿家到 ...
- BZOJ-1922 大陆争霸 多限制、分层图最短路 (堆+dijkstra)
1922: [Sdoi2010]大陆争霸 Time Limit: 10 Sec Memory Limit: 64 MB Submit: 1154 Solved: 478 [Submit][Status ...
- POJ 1797 Heavy Transportation 最短路变形(dijkstra算法)
题目:click here 题意: 有n个城市,m条道路,在每条道路上有一个承载量,现在要求从1到n城市最大承载量,而最大承载量就是从城市1到城市n所有通路上的最大承载量.分析: 其实这个求最大边可以 ...
- 最短路算法之Dijkstra算法通俗解释
Dijkstra算法 说明:求解从起点到任意点的最短距离,注意该算法应用于没有负边的图. 来,看图. 用邻接矩阵表示 int[][] m = { {0, 0, 0, 0, 0, 0}, {0, 0, ...
- 【最短路算法】Dijkstra+heap和SPFA的区别
单源最短路问题(SSSP)常用的算法有Dijkstra,Bellman-Ford,这两个算法进行优化,就有了Dijkstra+heap.SPFA(Shortest Path Faster Algori ...
- CCPC-Wannafly Winter Camp Day1 Div1 - 爬爬爬山 - [最短路][堆优化dijkstra]
题目链接:https://zhixincode.com/contest/3/problem/F?problem_id=39 样例输入 1 4 5 1 1 2 3 4 1 2 1 1 3 1 1 4 ...
- 基础最短路(模板 dijkstra)
Description 某省自从实行了很多年的畅通工程计划后,终于修建了很多路.不过路多了也不好,每次要从一个城镇到另一个城镇时,都有许多种道路方案可以选择,而某些方案要比另一些方案行走的距离要短很多 ...
随机推荐
- leetcood学习笔记-35-二分法
题目: 第一次提交; class Solution: def searchInsert(self, nums: List[int], target: int) -> int: for i in ...
- Delphi ResourceString的用法
在Delphi编程的那段“古老”的日子里(就是在版本4之前),在程序中使用字符串有两个基本的方法.你可以使用字符串将它们嵌入到源程序中,例如: MessageDlg( 'Leave your stin ...
- centos修改、保存文件的详细步骤
[一]修改文件 如果是使用普通用户登录的,需要先切换到管理员用户,打开终端,输入:su,接着按提示输入密码即可:然后使用命令进入需要修改文件的所在目录,常用的几个命令如下: ① cd + 目录名 ② ...
- Chrome 调试跨域问题解决方案之插件篇
跨域,就是A域名下的js,想请求B域名下的接口数据.跨域,只存在于浏览器端.App和小程序不存在跨域问题.跨域,分浏览器策略和服务器策略. 如果服务器配置了允许跨域,那就没有跨域问题 如果uni-ap ...
- (转)OpenFire源码学习之十七:HTTP Service插件
转:http://blog.csdn.net/huwenfeng_2011/article/details/43457645 HTTP Service插件 这里的http接口插件是神马? Openfi ...
- 2018ACM-ICPC EC-Final 现场赛I题 Misunderstanding...Missing 倒着DP
目录 Catalog Solution: (有任何问题欢迎留言或私聊 && 欢迎交流讨论哦 Catalog @ Problem: 很多acm群里有题面PDF了,我就不赘述了.简单说一 ...
- qemu的动态翻译机制
qemu的作者在QEMU, a Fast and Portable Dynamic Translator一文提到了qemu的动态翻译机制, 大致可以总结为如下过程: 目标代码中的一条指令 | |--( ...
- PAT_A1080#Graduate Admission
Source: PAT A1080 Graduate Admission (30 分) Description: It is said that in 2011, there are about 10 ...
- 【React-Native】---Android环境配置
一.前言 本文主要内容是在Window系统下配置Android APP的开发环境,要配置RN的Android开发环境需要3个依赖 1.JDK(版本必须是 1.8) 2.Node(版本必须高于8.3) ...
- 前端(十一)—— JavaScript基础:JS存在的位置、JS变量定义、调试方式、JS数据类型、运算符
JS存在的位置.JS变量定义.调试方式.JS数据类型.运算符 一.JS语言介绍 1.概念 浏览器脚本语言 可以编写运行在浏览器上的代码程序 属于解释性.弱语言类型编程语言 2.组成 ES语法:ECMA ...