[BZOJ3538]坑爹的GPS

最新推荐文章于 2024-10-31 16:45:22 发布

XIN13261615514

最新推荐文章于 2024-10-31 16:45:22 发布

阅读量166

点赞数

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/FYH-SSGSS/p/7719041.html

版权

题目描述 Description

有一天，$FJ$ 买了一辆车，但是，他一手下载了两个$GPS$ 系统。好了现在麻烦的事情来了，$GPS$ 有一个功能大概大家也知道，如果$FJ$ 没有按照$GPS$ 内置地图的最短路走，$GPS$ 就会报错来骚扰你。现在$FJ$ 准备从他的农舍（在1这个点）开车到他的谷屋（$n$ 这个点）。$FJ$ 给了你两个$GPS$ 系统内置地图的信息，他想知道，他最少会听到多少次报错（如果$FJ$ 走的路同时不满足两个$GPS$ ，报错次数+2）

输入描述 Input Description

第一行：$n$ ,$k$ ；$n$ 表示有$FJ$ 的谷屋在哪，同时保证$GPS$ 内置地图里的点没有超过$n$ 的点。$K$ 表示$GPS$ 内置地图里的路有多少条，如果两个点没有连接则表明这不是一条通路。
接下来$k$ 行，每行4个数$X$ ,$Y$ ,$A$ ,$B$ 分别表示从$X$ 到$Y$ 在第一个$GPS$ 地图里的距离是$A$ ，在第二个$GPS$ 地图里的是$B$ 。注意由于地形的其他因素$GPS$ 给出的边是有向边。

输出描述 Output Description

一个值，表示$FJ$ 最少听到的报错次数。

样例输入 Sample Input

样例输出 Sample Output

数据范围及提示 Data Size & Hint

$FJ$ 选择的路线是1 2 4 5，但是$GPS$ 1认为的最短路是1到3，所以报错一次，对于剩下的2 4 5，两个$GPS$ 都不会报错。
$N<=10000$ ，至于路有多少条自己算吧。数据保证所有的距离都在$2^{31}-1 $ 以内。

之前的一些废话

离期中还有两周，离NOIP还有不到三周，下午开始停课计划，是时候拼了！

题解

对于第一个$GPS$ 求一遍最短路，对第二个也求一遍最短路，注意是以$n$ 为源点进行最短路，所以还要存一个反图。然后我们尝试建一个报错的地图。在原图上搜索，对于一条边，如果走这一条边不满足第一个$GPS$ 上的最短路这条路边权加一，同理判断是否第二个的$GPS$ ，建好图之后再跑一个最短路即可。

代码

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
#define mem(a,b) memset(a,b,sizeof(a))
inline int read()
{
    int x=0,f=1;char c=getchar();
    while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
    while(isdigit(c)){x=x*10+c-'0';c=getchar();}
    return x*f;
}
const int maxn=10010,maxm=50010;
struct Edge
{
    int u,v,w1,w2,next;
    Edge() {}
    Edge(int _1,int _2,int _3,int _4,int _5):u(_1),v(_2),w1(_3),w2(_4),next(_5) {}
}e[maxm<<1],rev[maxm<<1],E[maxm<<1];
int n,m,first[maxn],rfirst[maxn],First[maxn],c1=-1,c2=-1,ce=-1,a,b,c,d,dis1[maxn],dis2[maxn],dis3[maxn];
bool vis[maxn];
queue<int> Q;
void addEdge(int a,int b,int c,int d)
{
    e[++c1]=Edge(a,b,c,d,first[a]);first[a]=c1;
    rev[++c2]=Edge(b,a,c,d,rfirst[b]);rfirst[b]=c2;
}
void AddEdge(int a,int b,int c){E[++ce]=Edge(a,b,c,0,First[a]);First[a]=ce;}
void SPFA1()
{
    mem(vis,0);mem(dis1,42);
    while(Q.size())Q.pop();
    Q.push(n);vis[n]=1;dis1[n]=0;
    while(Q.size())
    {
        int now=Q.front();Q.pop();
        for(int i=rfirst[now];i!=-1;i=rev[i].next)
            if(dis1[now]+rev[i].w1<dis1[rev[i].v])
            {
                dis1[rev[i].v]=dis1[now]+rev[i].w1;
                if(!vis[rev[i].v])
                {
                    Q.push(rev[i].v);
                    vis[rev[i].v]=1;
                } 
            }
        vis[now]=0;
    }
    return;
}
void SPFA2()
{
    mem(vis,0);mem(dis2,42);
    while(Q.size())Q.pop();
    Q.push(n);vis[n]=1;dis2[n]=0;
    while(Q.size())
    {
        int now=Q.front();Q.pop();
        for(int i=rfirst[now];i!=-1;i=rev[i].next)
            if(dis2[now]+rev[i].w2<dis2[rev[i].v])
            {
                dis2[rev[i].v]=dis2[now]+rev[i].w2;
                if(!vis[rev[i].v])
                {
                    Q.push(rev[i].v);
                    vis[rev[i].v]=1;
                } 
            }
        vis[now]=0;
    }
    return;
}
void SPFA3()
{
    mem(vis,0);mem(dis3,42);
    while(Q.size())Q.pop();
    Q.push(1);vis[1]=1;dis3[1]=0;
    while(Q.size())
    {
        int now=Q.front();Q.pop();
        for(int i=First[now];i!=-1;i=E[i].next)
            if(dis3[now]+E[i].w1<dis3[E[i].v])
            {
                dis3[E[i].v]=dis3[now]+E[i].w1;
                if(!vis[E[i].v])
                {
                    Q.push(E[i].v);
                    vis[E[i].v]=1;
                } 
            }
        vis[now]=0;
    }
    return;
}
int main()
{
    mem(first,-1);mem(rfirst,-1);mem(First,-1);
    n=read();m=read();
    for(int i=1;i<=m;i++)a=read(),b=read(),c=read(),d=read(),addEdge(a,b,c,d);
    SPFA1();SPFA2();
    for(int point=1;point<=n;point++)
        for(int i=first[point];i!=-1;i=e[i].next)
        {
            int cnt=0;
            if(dis1[point]!=dis1[e[i].v]+e[i].w1)cnt++;
            if(dis2[point]!=dis2[e[i].v]+e[i].w2)cnt++;
            AddEdge(point,e[i].v,cnt);
        }
    SPFA3();
    printf("%d\n",dis3[n]);
    return 0;
}