CODE[VS] 2273 USACO 扬帆远洋大战牧师妹酱 双端队列优化SPFA

戳我可穿越

题目描述 Description
扬帆远洋为了打败牧师妹酱,不惜一切代价通过各种各样的手段得到了传奇道具“れいじまいご”来释放能量使得全球充满扬帆远洋的威严气场和能量场。但是,“れいじまいご”释放的能量只能在某个很小的区域中存在,为了使得效果可以逸散的全球,扬帆远洋决定效仿心灵操控者——尤里,在世界各地修建了N个“れいじまいご·ベクトル変換”使得其能量场可以在全世界张开。

根据计算机的计算结果,这N个地点大多都是非常偏僻的城镇。因此,之间的道路并不定都是畅通的。这些地区之间通过R条道路连接,为了使得之间前往的方便,扬帆远洋TEL还帮助扬帆远洋产生了P个空间隧道,扬帆远洋称之为“時間と空間の扉”。每条道路i或者“時間と空間の扉”i连接地区Ai到Bi,花费为Ci。道路是双向的,可以从Ai到Bi,也可以从Bi到Ai,花费都是Ci。然而“時間と空間の扉”与之不同,只可以从Ai到Bi。事实上,由于牧师妹酱的墨镜党(看看我的头像就知道了)的不断干扰使得链接非常不稳定。牧师妹酱的墨镜党的干扰如下:如果有一条“時間と空間の扉”可以从Ai到Bi,那么保证不可能通过一些道路和“時間と空間の扉”从Bi回到Ai。所以有些地区可能无法到达,需要扬帆远洋对自己进行Teleport才可以。

扬帆远洋现在呆在S城镇,由于他需要把每个“れいじまいご·ベクトル変換”的装置调整完毕,因此扬帆远洋写了一个程序计算从S出发,对于每个地区而言计算出了最小的花费,如果不可能,则输出“Need
Teleport To Repair”。最后输出,一共有多少个地方是需要用Teleport的。

但是由于墨镜党不断的捣乱,这个程序现在无法工作,请你复述出来这个程序。

输入描述 Input Description
第1行有四个空格隔开的整数T、R、P和S。
第2至R+1行每行有3个空格隔开的正整数,表示一条道路的Ai、Bi和Ci。
第R+2至R+P+1行每行有3个空格隔开的正整数,表示一条“時間と空間の扉”的Ai、Bi和Ci。

输出描述 Output Description
第1到T行每行输出的是从S到达地区i的最小花费,如果不存在输出“Need Teleport
To Repair”(没有引号)。 最后一行为一共有多少个地方是需要用Teleport的。

样例输入 Sample Input
6 3 3 4
1 2 5
5 6 10
3 5 -100
4 6 -100
1 3 -10

样例输出 Sample Output
Need Teleport To Repair
Need Teleport To Repair
5
0
-95
-100 2

数据范围及提示 Data Size & Hint
1≤T≤25000,1≤R≤50000,1≤P≤50000,1≤S≤T
1≤Ai≤T,1≤Bi≤T,-10000≤Ci≤10000(对于“時間と空間の扉”来说花费才可能是负数)

牧师妹酱, 扬帆远洋, 扬帆远洋TEL均为真实存在的人物.

这里写题解:
题目共有25000个点 50000条边 在最差的状态下(图为一条链) 时间复杂度肯定是超时的 而且要注意 数组一定要够大 虽然内存限制非常小 但数组太小会导致RE

这是裸的SPFA

这是数组开小的锅

注意 因为有负权所以不能用dijkstra 本质为贪心的它看不到远处的负权

那么怎么办呢???
我们可以想一下 如果队首对答案的贡献更小 也就是说要加入的新点会使最短路更短 我们就先用要加入的点 后用队首点。
很明显 我们可以用一个双端队列来优化一下 这就是SPFA的SLF优化
伪代码:

if(!vis[es[i].to)
{
    vis[es[i].to] = 1;  
    if(d[es[i].to] < d[q.front] q.push_front(es[i].to);
    else q.push_back(es[i].to);
}

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 233333;
const int inf = 2147483647;
int tot;
int first[maxn],nxt[maxn],d[maxn];
struct edge{
    int from,to,cost;
}es[maxn];
void init()
{
    memset(first,-1,sizeof(first));
    tot=0;
}
void build(int ff,int tt,int dd)
{
    es[++tot]=(edge){ff,tt,dd};
    nxt[tot]=first[ff];
    first[ff]=tot;
}
bool vis[maxn];
deque<int>q;
void spfa(int s)
{
    d[s]=0;
    vis[s]=1;
    q.push_front(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop_front();
        vis[u]=0;
        for(int i=first[u];i!=-1;i=nxt[i])
        {
            int v=es[i].to;
            if(d[v]>d[u]+es[i].cost)
            {
                d[v]=d[u]+es[i].cost;
                if(!vis[v])
                {
                    vis[v]=1;
                    if(q.empty())
                        q.push_front(v);
                    else if(d[v]>d[q.front()])
                        q.push_back(v);
                    else
                        q.push_front(v);
                }
            }
        }
    }
}
int main()
{
    init();
    int T,E,P,S;
    scanf("%d%d%d%d",&T,&E,&P,&S);
    for(int i=1;i<=E;i++)
    {
        int ff,tt,dd;
        scanf("%d%d%d",&ff,&tt,&dd);
        build(ff,tt,dd);
        build(tt,ff,dd);
    }
    for(int i=1;i<=P;i++)
    {
        int ff,tt,dd;
        scanf("%d%d%d",&ff,&tt,&dd);
        build(ff,tt,dd);
    }
    for(int i = 1;i <= T;i ++)  d[i] = inf;
    int cnt=0;
    spfa(S);
    for(int i=1;i<=T;i++)
    {
        if(d[i]==inf)
        {
            puts("Need Teleport To Repair");
            cnt++;
        }
        else
            printf("%d\n",d[i]);
    }
    printf("%d\n",cnt);
    return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值