K短路 (A*算法) [Usaco2008 Mar]牛跑步&[Sdoi2010]魔法猪学院
A*属于搜索的一种,启发式搜索,即:每次搜索时加一个估价函数
这个算法可以用来解决K短路问题,常用的估价函数是:已经走过的距离+期望上最短的距离
通常和Dijkstra一起解决K短路
BZOJ1598:牛跑步
求前K短路
因为A*算法我们每次用来向外拓展的是估价函数最小的点,那么,我们必定能够得到,第一个到达n的是最短路。(Dijkstra的贪心,可证)
那么,我们思考一下,第二个到达n的就是次短路!
由此观之:第K个到达的就是K短路
因此,A*算法可以用来解决K短路问题
附上代码:
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <iostream>
#include <map>
using namespace std;
#define N 10005
#define ll long long
struct node
{
int to,next,val;
}e[N*10],E[N*10];
int head[N],cnt,cnt1,n,m,K,vis[N],head1[N];
int dis[N];
void add(int x,int y,int z)
{
e[cnt].to=y;
e[cnt].next=head[x];
e[cnt].val=z;
head[x]=cnt++;
}
void add1(int x,int y,int z)
{
E[cnt1].to=y;
E[cnt1].next=head1[x];
E[cnt1].val=z;
head1[x]=cnt1++;
}
priority_queue<pair<int ,int > >q;
void Dijkstra()
{
memset(dis,0x3f,sizeof(dis));
dis[n]=0;
q.push(make_pair(0,n));
while(!q.empty())
{
int x=q.top().second;q.pop();
if(vis[x])continue;
vis[x]=1;
for(int i=head1[x];i!=-1;i=E[i].next)
{
int to1=E[i].to;
if(dis[to1]>dis[x]+E[i].val)
{
dis[to1]=dis[x]+E[i].val;
q.push(make_pair(-dis[to1],to1));
}
}
}
}
void Astar(int s)
{
q.push(make_pair(-dis[s],s));
while(!q.empty())
{
int d=q.top().first;int x=q.top().second;q.pop();
if(x==n)
{
K--;
printf("%d\n",-d);
if(!K)return ;
}
for(int i=head[x];i!=-1;i=e[i].next)
{
int to1=e[i].to;
q.push(make_pair(d+dis[x]-e[i].val-dis[to1],to1));
}
}
}
int main()
{
memset(head1,-1,sizeof(head1));
memset(head,-1,sizeof(head));
scanf("%d%d%d",&n,&m,&K);
for(int i=1;i<=m;i++)
{
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
add(y,x,z);
add1(x,y,z);
}
Dijkstra();
Astar(1);
while(K--)puts("-1");
return 0;
}
1975: [Sdoi2010]魔法猪学院
这个贪心贪的很显然,为了尽可能的多,自然是选择前K短路
那么,同样,我们用Dijkstra的贪心方法求n的前K短路
但是这种方法实际的时间复杂度不是很对,理论上可以过很多题。这道题BZOJ可以过,但是洛谷上过不去
附上不完美的代码:
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <iostream>
#include <map>
using namespace std;
#define N 5005
#define ll long long
struct node
{
int to,next;
double val;
}e[200005],E[200005];
int head[N],cnt,cnt1,n,m,vis[N],head1[N],ans;
double dis[N],K;
void add(int x,int y,double z)
{
e[cnt].to=y;
e[cnt].next=head[x];
e[cnt].val=z;
head[x]=cnt++;
}
void add1(int x,int y,double z)
{
E[cnt1].to=y;
E[cnt1].next=head1[x];
E[cnt1].val=z;
head1[x]=cnt1++;
}
priority_queue<pair<double ,int > >q;
void Dijkstra()
{
for(int i=0;i<N;i++)dis[i]=1e9;
dis[n]=0;
q.push(make_pair(0,n));
while(!q.empty())
{
int x=q.top().second;q.pop();
if(vis[x])continue;
vis[x]=1;
for(int i=head1[x];i!=-1;i=E[i].next)
{
int to1=E[i].to;
if(dis[to1]>dis[x]+E[i].val)
{
dis[to1]=dis[x]+E[i].val;
q.push(make_pair(-dis[to1],to1));
}
}
}
}
void Astar(int s)
{
q.push(make_pair(-dis[s],s));
while(!q.empty())
{
double d=q.top().first;int x=q.top().second;q.pop();
if(x==n)
{
if(K+d<0-1e-9)return ;
K+=d;
ans++;
}
for(int i=head[x];i!=-1;i=e[i].next)
{
int to1=e[i].to;
q.push(make_pair(d+dis[x]-e[i].val-dis[to1],to1));
}
}
}
int main()
{
memset(head1,-1,sizeof(head1));
memset(head,-1,sizeof(head));
scanf("%d%d%lf",&n,&m,&K);
for(int i=1;i<=m;i++)
{
int x,y;double z;
scanf("%d%d%lf",&x,&y,&z);
add(x,y,z);
add1(y,x,z);
}
Dijkstra();
Astar(1);
printf("%d\n",ans);
return 0;
}
正解是用可持久化可并堆+A*算法,挖坑待填
K短路 (A*算法) [Usaco2008 Mar]牛跑步&[Sdoi2010]魔法猪学院的更多相关文章
- bzoj 1598: [Usaco2008 Mar]牛跑步 [k短路 A*] [学习笔记]
1598: [Usaco2008 Mar]牛跑步 题意:k短路 ~~貌似A*的题目除了x数码就是k短路~~ \[ f(x) = g(x) + h(x) \] \(g(x)\)为到达当前状态实际代价,\ ...
- bzoj 1598: [Usaco2008 Mar]牛跑步 -- 第k短路,A*
1598: [Usaco2008 Mar]牛跑步 Time Limit: 10 Sec Memory Limit: 162 MB Description BESSIE准备用从牛棚跑到池塘的方法来锻炼 ...
- Bzoj 1598: [Usaco2008 Mar]牛跑步 dijkstra,堆,K短路,A*
1598: [Usaco2008 Mar]牛跑步 Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 427 Solved: 246[Submit][St ...
- BZOJ_1598_[Usaco2008 Mar]牛跑步_A*
BZOJ_1598_[Usaco2008 Mar]牛跑步_A* Description BESSIE准备用从牛棚跑到池塘的方法来锻炼. 但是因为她懒,她只准备沿着下坡的路跑到池塘, 然后走回牛棚. B ...
- bzoj1975: [Sdoi2010]魔法猪学院【k短路&A*算法】
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 2446 Solved: 770[Submit][Statu ...
- [BZOJ1975][SDOI2010]魔法猪学院(k短路,A*)
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 2748 Solved: 883[Submit][Statu ...
- Bzoj 1975: [Sdoi2010]魔法猪学院 dijkstra,堆,A*,K短路
1975: [Sdoi2010]魔法猪学院 Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 1357 Solved: 446[Submit][Statu ...
- [SDOI2010]魔法猪学院(A*,最短路)
[SDOI2010]魔法猪学院(luogu) Description 题目描述 iPig在假期来到了传说中的魔法猪学院,开始为期两个月的魔法猪训练.经过了一周理论知识和一周基本魔法的学习之后,iPig ...
- bzoj 1598: [Usaco2008 Mar]牛跑步【A*K短路】
A*K短路模板,详见https://blog.csdn.net/z_mendez/article/details/47057461 算法流程: 把有向图全建成反向边,跑一遍所有点到t的最短路记为dis ...
随机推荐
- Linux笔记2
touch 创建文件. echo 输出 >> 将输出写入到文件中 echo sss >> a.txt cat 查看文件内容 帮助命令 man 命令 man ...
- JAVA 综合面试题
JAVA 综合面试题 2007-08-12 目录 TOC \o "1-3" \h \z \u Java面试题整理 9 Java面向对象 9 1. super()与this()的区别 ...
- Oracle数据库date类型与Java中Date的联系与转化
以下是对Java中的日期对象与Oracle中的日期之间的区别与联系做点说明,以期对大家有所帮助.new Date():分配 Date 对象并初始化此对象,以表示分配它的时间(精确到毫秒),就是系统当前 ...
- C++堆栈问题
编写C++中的两个类 一个只能在栈中分配空间 一个只能在堆中分配. 解答: (1)代码如下 (2)堆栈分配内存的介绍 一.一个经过编译的C/C++的程序占用的内存分成以下几个部分:1.栈区(stack ...
- 使用FFmpeg捕获一帧摄像头图像
最近在研究FFmpeg,比较惊讶的是网上一大堆资料都是在说如何从已有的视频中截取一帧图像,却很少说到如何直接从摄像头中捕获一帧图像,其实我一直有个疑问,就是在Linux下,大家是用什么库来采集摄像头的 ...
- Ubuntu下定时任务和自启动任务的部署
1.定时任务的部署,最简单的方法是执行 crontab -e 然后在下面加上世间周期配置和要执行的命令,一般情况下,可以把要执行的任务用bash脚本封装一下,格式如下所示: minute hour ...
- Spring-mvc设置@RequestMapping标签更改返回头及@RequestMapping简述
1. 引子:设置返回头 2. 简述 3. value 4. method 5. consumes/produces 6. params 7. headers 1. 引子:设置返回头 返回JSON内容时 ...
- eclipse 安装svn和gradle
公司项目用的eclispe svn和gradle 所以需要配置 SVN教程:https://blog.csdn.net/jieshaowang1229/article/details/5159499 ...
- Spring Boot 发送邮件
需求 最近因为业务的变更,需要对老用户进行发送邮件处理.目前市面上也有很多代发邮件的接口,可以接入.由于量不是特别大,放弃了这个途径.改用我们自己通过 smtp 发送邮件来处理. 技术选择 Java ...
- 用Promise实现:带延时功能的链式调用
// 1) 调用方式 new People('whr').sleep(3).eat('apple').sleep(5).eat('durian'); // 2) 打印结果 'hello, whr' - ...