洛谷 P2387 魔法森林

洛谷 P2387 魔法森林

---

题目

题目描述

为了得到书法大家的真传，小 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”（不 含引号）。

输入输出样例

输入样例#1：

4 5 
1 2 19 1 
2 3 8 12 
2 4 12 15 
1 3 17 8 
3 4 1 17 

输出样例#1：

32

输入样例#2：

3 1 
1 2 1 1 

输出样例#2：

-1

说明

解释1

如果小 E 走路径 1→2→4，需要携带 19+15=34 个守护精灵； 如果小 E 走路径 1→3→4，需要携带 17+17=34 个守护精灵； 如果小 E 走路径 1→2→3→4，需要携带 19+17=36 个守护精灵； 如果小 E 走路径 1→3→2→4，需要携带 17+15=32 个守护精灵。 综上所述，小 E 最少需要携带 32 个守护精灵。

解释2

小 E 无法从 1 号节点到达 3 号节点，故输出-1。

---

题解

按照b的从小到大排序，然后每加一条边，做一次spfa（在上一次的结果上做，不要傻乎乎的从头开始做），很显然答案就是d[n]+b（d[n]就是这次修改以后到n所需A型守护精灵的数量，b就是该次加边时，所对应的B型守护精灵的数量）（无论这条边是否对spfa中修改有影响，因为如果没影响的话，上一次修改的答案一定比这次小）

---

代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;

int n,m,tot,now,ans;
int lnk[50005],a[100005],b[100005],x[100005],y[100005],s[100005],d[50005];
bool vis[50005];
struct edge
{
    int nxt,y,v;
} e[200005];

int readln()
{
    int x=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while ('0'<=ch&&ch<='9') x=x*10+ch-48,ch=getchar();
    return x;
}

void qsort(int l,int r)
{
    int i=l,j=r,mid=b[rand()%(r-l+1)+l],t;
    do {
        while (b[i]<mid) i++;
        while (b[j]>mid) j--;
        if (i<=j) {
            t=x[i];x[i]=x[j];x[j]=t;
            t=y[i];y[i]=y[j];y[j]=t;
            t=a[i];a[i]=a[j];a[j]=t;
            t=b[i];b[i]=b[j];b[j]=t;
            i++;j--;
        }
    } while (i<=j);
    if (i<r) qsort(i,r);
    if (l<j) qsort(l,j);
}

void add(int x,int y,int v)
{
    tot++;e[tot].nxt=lnk[x];lnk[x]=tot;e[tot].y=y;e[tot].v=v;
}

int max(int x,int y){return x>y?x:y;}

int min(int x,int y){return x<y?x:y;}

int spfa(int x,int y)
{
    s[++now]=x;s[++now]=y;
    while (now!=0)
    {
        int a=s[now--];
        vis[a]=false;
        for (int i=lnk[a];i;i=e[i].nxt)
        {
            int b=e[i].y,w=max(d[a],e[i].v);
            if (d[b]>w) {
                d[b]=w;
                if (!vis[b]) {
                    vis[b]=true;s[++now]=b;
                }
            }
        }
    }
    return d[n];
}

int main()
{
    n=readln();m=readln();
    for (int i=1;i<=m;i++) x[i]=readln(),y[i]=readln(),a[i]=readln(),b[i]=readln();
    qsort(1,m);
    memset(d,63,sizeof(d));
    ans=1e9;s[++now]=1;d[1]=0;
    for (int i=1;i<=m;i++)
    {
        add(x[i],y[i],a[i]);
        add(y[i],x[i],a[i]);
        ans=min(ans,b[i]+spfa(x[i],y[i]));
    }
    if (ans==1e9) printf("-1");else printf("%d",ans);
    return 0;
}