【NOI2014】魔法森林 - 动态加边SPFA

题目描述

为了得到书法大家的真传，小 E 同学下定决心去拜访住在魔法森林中的隐 士。魔法森林可以被看成一个包含 n 个节点 m 条边的无向图，节点标号为 1,2,3,…,n，边标号为 1,2,3,…,m。初始时小 E 同学在 1 号节点，隐士则住在 n 号节点。小 E 需要通过这一片魔法森林，才能够拜访到隐士。

魔法森林中居住了一些妖怪。每当有人经过一条边的时候，这条边上的妖怪 就会对其发起攻击。幸运的是，在 1 号节点住着两种守护精灵：A 型守护精灵与 B 型守护精灵。小 E 可以借助它们的力量，达到自己的目的。

只要小 E 带上足够多的守护精灵，妖怪们就不会发起攻击了。具体来说，无 向图中的每一条边 ei 包含两个权值 ai 与 bi 。若身上携带的 A 型守护精灵个数不 少于 ai ，且 B 型守护精灵个数不少于 bi ，这条边上的妖怪就不会对通过这条边 的人发起攻击。当且仅当通过这片魔法森林的过程中没有任意一条边的妖怪向 小 E 发起攻击，他才能成功找到隐士。

由于携带守护精灵是一件非常麻烦的事，小 E 想要知道，要能够成功拜访到 隐士，最少需要携带守护精灵的总个数。守护精灵的总个数为 A 型守护精灵的 个数与 B 型守护精灵的个数之和。

输入输出格式

输入格式：

输入文件的第 1 行包含两个整数 n,m，表示无向图共有 n 个节点，m 条边。 接下来 m 行，第i+ 1 行包含 4 个正整数 Xi,Yi,ai,bi，描述第i条无向边。 其中Xi与 Yi为该边两个端点的标号，ai 与 bi 的含义如题所述。 注意数据中可能包含重边与自环。

输出格式：

输出一行一个整数：如果小 E 可以成功拜访到隐士，输出小 E 最少需要携 带的守护精灵的总个数；如果无论如何小 E 都无法拜访到隐士，输出“-1”（不 含引号）。 

 

思路

正解LCT调不出，SPFA的做法很简单。

把每条边按照a依次加入，每加入一条边

就跑一遍SPFA，不断更新答案。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 50000 + 10;
vector<pair<int,int> > edges[maxn];
int n,m,dis[maxn],ans = 9999999;
bool vis[maxn];
struct Node {
    int u,v,a,b;
    inline bool operator < (Node cmp) const {
        return a < cmp.a;
    }
}edge[maxn*10];
inline void add(int a,int b,int c) {
    edges[a].push_back(make_pair(b,c));
    edges[b].push_back(make_pair(a,c));
}
queue<int> Q;
void spfa(int num) {
    Q.push(edge[num].u);
    Q.push(edge[num].v);
    vis[edge[num].u] = true;
    vis[edge[num].v] = true;
    while (Q.size()) {
        int now = Q.front(); Q.pop();
        vis[now] = false;
        for (size_t i = 0;i < edges[now].size();i++) {
            int to = edges[now][i].first,w = edges[now][i].second;
            if (dis[to] > max(dis[now],w)) {
                dis[to] = max(dis[now],w);
                if (!vis[to]) {
                    Q.push(to);
                    vis[to] = true;
                }
            }
        }
    }
}
int main() {
    scanf("%d%d",&n,&m);
    for (int i = 1;i <= m;i++) scanf("%d%d%d%d",&edge[i].u,&edge[i].v,&edge[i].a,&edge[i].b);
    memset(dis,127/3,sizeof(dis));
    sort(edge+1,edge+m+1);
    Q.push(1);
    dis[1] = 0;
    vis[1] = true;
    for (int i = 1;i <= m;i++) {
        add(edge[i].u,edge[i].v,edge[i].b);
        spfa(i);
        ans = min(ans,dis[n]+edge[i].a);
    }
    if (ans == 9999999) printf("-1");
    else  printf("%d",ans);
    return 0;
}