P3106 [USACO14OPEN]GPS的决斗(最短路)
化简:够简的了.....但是!翻译绝对有锅。
这个最短路是从n到每个点的单源最短路,也就是最短路径树。
那么,思路就很明确了。建两个图,然后跑两边SPFA,记录下最短路径。
然后,对于两点之间的边,如果最短路不经过它,那么最终图边权+1;
然后在最终图上(边权为0,1,2)跑一遍SPFA即可。
一开始我想复杂了,在想怎么记录路径,怎么重构图.balabala。
然后发现,怎么才能让两点不在最短路径上呢?
SPFA的松弛操作,依据是三角不等式。于是,如果两点之间的最短路的距离如果不等于边权(也就是最短路径不过它俩之间的边)那么就它就是一条会报警的边。
$$看来SPFA最重要的是三角不等式$$
于是,只要暴力跑三遍SPFA即可。
用尽浑身解数,信仰SPFA,各种常数优化,读入挂,我还是没能跑到最优解的第一面....在第二页前几个徘徊.....
代码:(这么长的图论题也不多见,但是好像逛公园就那么长)
#include<iostream>
#include<cstdio>
#include<queue>
#define rg register
using namespace std;
const int maxn=1e6+;
int n,m;
inline int read()
{
int x=,f=;char s=getchar();
while(s>''||s<''){if(s=='-')f=-;s=getchar();}
while(s<=''&&s>=''){x=x*+s-'';s=getchar();}
return x*f;
}
struct edge
{
int to,next,dis;
}e1[maxn],e2[maxn],e[maxn];
int head[maxn],head1[maxn],head2[maxn];
int cnt1,cnt2,cnt;
inline void addedge1(int from,int to,int dis)
{
e1[++cnt1].next=head1[from];
e1[cnt1].to=to;
e1[cnt1].dis=dis;
head1[from]=cnt1;
}
inline void addedge2(int from,int to,int dis)
{
e2[++cnt2].next=head2[from];
e2[cnt2].to=to;
e2[cnt2].dis=dis;
head2[from]=cnt2;
}
inline void addedge(int from,int to,int dis)
{
e[++cnt].next=head[from];
e[cnt].to=to;
e[cnt].dis=dis;
head[from]=cnt;
} int dis1[maxn],vis1[maxn],pre1[maxn];
struct cmp1
{
bool operator () (int a,int b)
{
return dis1[a]>dis1[b];
}
};
inline void spfa1()
{
priority_queue < int , vector < int > , cmp1 > q;
for(rg int i=;i<=n;i++)
{
dis1[i]=;
vis1[i]=;
}
q.push(n);
dis1[n]=;
vis1[n]=;
while(!q.empty())
{
int u=q.top();
q.pop();
vis1[u]=;
for(rg int i=head1[u];i;i=e1[i].next)
{
int v=e1[i].to;
if(dis1[v]>dis1[u]+e1[i].dis)
{
dis1[v]=dis1[u]+e1[i].dis;
if(vis1[v]==)
{
vis1[v]=;
q.push(v);
}
}
}
}
}
int dis2[maxn],vis2[maxn],pre2[maxn];
struct cmp2
{
bool operator () (int a,int b)
{
return dis2[a]>dis2[b];
}
};
inline void spfa2()
{
priority_queue < int , vector < int > , cmp2 > q;
for(rg int i=;i<=n;i++)
{
dis2[i]=;
vis2[i]=;
}
q.push(n);
dis2[n]=;
vis2[n]=;
while(!q.empty())
{
int u=q.top();
q.pop();
vis2[u]=;
for(rg int i=head2[u];i;i=e2[i].next)
{
int v=e2[i].to;
if(dis2[v]>dis2[u]+e2[i].dis)
{
dis2[v]=dis2[u]+e2[i].dis;
if(vis2[v]==)
{
vis2[v]=;
q.push(v);
}
}
}
}
}
int dis[maxn],vis[maxn];
struct cmp
{
bool operator () (int a,int b)
{
return dis[a]>dis[b];
}
};
inline void rebuild()
{
for(rg int i=;i<=n;i++)
{
for(int j=head1[i];j;j=e1[j].next)
{
int v1=e1[j].to;
int v2=e2[j].to;
//cout<<i<<' '<<v1<<endl;
int d=;
if(dis1[v1]-dis1[i]!=e1[j].dis)
d++;
if(dis2[v2]-dis2[i]!=e2[j].dis)
d++;
addedge(v1,i,d);
}
}
}
inline void spfa()
{
priority_queue < int , vector < int > , cmp > q;
for(rg int i=;i<=n;i++)
{
dis[i]=;
vis[i]=;
}
q.push();
dis[]=;
vis[]=;
while(!q.empty())
{
int u=q.top();
q.pop();
vis[u]=;
for(rg int i=head[u];i;i=e[i].next)
{
int v=e[i].to;
if(dis[v]>dis[u]+e[i].dis)
{
dis[v]=dis[u]+e[i].dis;
if(vis[v]==)
{
vis[v]=;
q.push(v);
}
}
}
}
}
int main()
{
n=read();
m=read();
for(rg int i=;i<=m;i++)
{
int a=read(),b=read(),c=read(),d=read();
//scanf("%d%d%d%d",&a,&b,&c,&d);
addedge1(b,a,c);
addedge2(b,a,d);
//addedge(a,b,0);
}
spfa1();
spfa2();
rebuild();
spfa();
printf("%d",dis[n]);
return ;
}
(完)
P3106 [USACO14OPEN]GPS的决斗(最短路)的更多相关文章
- Luogu P3106 [USACO14OPEN]GPS的决斗Dueling GPS's(最短路)
P3106 [USACO14OPEN]GPS的决斗Dueling GPS's 题意 题目描述 Farmer John has recently purchased a new car online, ...
- BZOJ 3538 == 洛谷 P3106 [USACO14OPEN]GPS的决斗Dueling GPS's
P3106 [USACO14OPEN]GPS的决斗Dueling GPS's 题目描述 Farmer John has recently purchased a new car online, but ...
- [USACO14OPEN]GPS的决斗Dueling GPS's
题目概况 题目描述 给你一个\(N\)个点的有向图,可能有重边. 有两个\(GPS\)定位系统,分别认为经过边\(i\)的时间为\(P_i\),和\(Q_i\). 每走一条边的时候,如果一个系统认为走 ...
- 洛谷 3106 [USACO14OPEN]GPS的决斗Dueling GPS's 3720 [AHOI2017初中组]guide
[题解] 这两道题是完全一样的. 思路其实很简单,对于两种边权分别建反向图跑dijkstra. 如果某条边在某一种边权的图中不是最短路上的边,就把它的cnt加上1.(这样每条边的cnt是0或1或2,代 ...
- 2018.07.22 洛谷P3106 GPS的决斗Dueling GPS's(最短路)
传送门 图论模拟题. 这题直接写3个(可以压成一个)spfa" role="presentation" style="position: relative;&q ...
- USACO Dueling GPS's
洛谷 P3106 [USACO14OPEN]GPS的决斗Dueling GPS's 洛谷传送门 JDOJ 2424: USACO 2014 Open Silver 2.Dueling GPSs JDO ...
- CSP-S 2019图论总结
CSP-S 2019图论总结 一.最短路问题 模板 Floyd算法 void floyd() { memset(map,0x3f,sizeof(map)); for(int i=1;i<=n;i ...
- [USACO14OPEN] Dueling GPS's[最短路建模]
题目描述 Farmer John has recently purchased a new car online, but in his haste he accidentally clicked t ...
- [BZOJ3538]坑爹的GPS
题目描述 Description 有一天,\(FJ\) 买了一辆车,但是,他一手下载了两个\(GPS\) 系统.好了现在麻烦的事情来了,\(GPS\) 有一个功能大概大家也知道,如果\(FJ\) 没有 ...
随机推荐
- C# WinForm 跨线程访问控件(实用简洁写法)
C# WinForm 跨线程访问控件(实用简洁写法) 1.<C# WinForm 跨线程访问控件(实用简洁写法)> 2.<基于.NET环境,C#语言 实现 TCP NAT ...
- springmvc中将servlet api对象作为处理方法的入参使用
在springmvc中,控制器不依赖任何servlet api对象,也可以将servlet api对象作为处理方法的入参使用,非常方便,比如需要使用HttpSession对象,那么就可以直接将Http ...
- 关于jQery中$.Callbacks()的理解
$.Callbacks()主要使用了回调,而说到回调又不得不说javascript的事件循环机制了. 所以想了解回调最好先看看js运行机制. $.Callbacks()可以理解为创建一个回调队列 va ...
- 基于mosquitto的MQTT客户端实现C语言
在对MQTT的学习过程中 一下的内容对我提供了帮助 https://www.runoob.com/w3cnote/mqtt-intro.html 对MQTT的入门级介绍 很基础讲解了什么是MQTT h ...
- C语言-查找一个元素在数组中的位置
#include<stdio.h> #include <stdlib.h> #include <time.h> int search(int key, int a[ ...
- unittest中的方法调用时报错ValueError: no such test method in <class 'mytestcase.MyTestCase'>: runTest
调用unittest中的方法时报错: ValueError: no such test method in <class 'mytestcase.MyTestCase'>: runTest ...
- Windows系统调用中API从3环到0环(下)
Windows内核分析索引目录:https://www.cnblogs.com/onetrainee/p/11675224.html Windows系统调用中API从3环到0环(下) 如果对API在 ...
- Centos7升级Git版本
centos 升级 Git 版本 问题描述 centos7 系统默认的 git 安装版本是 1.8,但是在项目构建中发现 git 版本过低,于是用源码编译的方式进行升级.同时该文章也适用于安装新的 g ...
- 致远OA_0day批量植Cknife马一步到位
最近各位师傅都在刷这个嘛,原本的exp是上传一个test123456.jsp的命令执行的马子,不过我在试的时候发现替换成C刀一句话出错,原因未知,并且test123456.jsp如果存在的话用原来ex ...
- [HNOI2007] 理想正方形 二维ST表
题目描述 有一个a*b的整数组成的矩阵,现请你从中找出一个n*n的正方形区域,使得该区域所有数中的最大值和最小值的差最小. 输入输出格式 输入格式: 第一行为3个整数,分别表示a,b,n的值 第二行至 ...