【BZOJ3669】魔法森林(SPFA动态加边)

题目连接

题目描述

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

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

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

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

题意简述: 给你一张无向带权图，每条边有两种权值。求出一条从 1到n 的路径使得路径上的两种权值的最大值之和最小，并输出这个最小值

题解

首先此题不可二分,因为二分跨度大,可能导致漏掉最优解。
既然二分不行我们就枚举?
新套路 SPFA 动态加边。(看起来是个暴力，但每次重新跑SPFA时实际复杂度应该不高)
先把边按 a的值 从小到大排序,让后枚举最大的a值,再一条条把边加进来,以b为权值跑SPFA。
注意我们保证了a单调递增,那么每次加边后重新SPFA时 b 的值一定是递减的,因此不用重新设dis为INF。
另外加边的时候是把边的端点加到队列里,新加的边要起作用肯定是要经过它,那么只需加相连的点(也只有加相连的点才有用)。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<queue>
using namespace std;
inline int read()
{
    int x=0;char ch=getchar();int t=1;
    for(;ch>'9'||ch<'0';ch=getchar()) if(ch=='-') t=-1;
    for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+(ch-48);
    return x*t;
}
const int M=1e5+10;
const int N=5e4+10;
const int INF=1e9;
struct edge{
    int from,to,next;int a,b;
}ed[M],e[M<<1];
int head[N];
int dis[N];
int id[M];
bool g[N];
int n,m;
queue<int> Q;
inline bool cmp(int x,int y)
{
    return ed[x].a<ed[y].a;
}
inline void add(int x,int y,int a,int b,int& cnt,edge* ed,int op)
{
    ed[++cnt]=(edge){x,y,head[x],a,b};
    if(op==1) head[x]=cnt;
}
inline void SPFA()
{
    while(!Q.empty()){
        register int p=Q.front();Q.pop();
        for(register int v,i=head[p];i;i=e[i].next)
        {
            v=e[i].to;
            if(dis[v]<=max(dis[p],e[i].b)) continue;
            dis[v]=max(dis[p],e[i].b);
            if(!g[v]) Q.push(v),g[v]=1;
        }
        g[p]=0;
    }
}
int main()
{
    n=read();m=read();register int x,y,a,b;dis[1]=0;
    int cnt1=0,cnt2=0;
    for(register int i=2;i<=n;i++) dis[i]=INF;
    for(register int i=1;i<=m;i++) {x=read();y=read();a=read();b=read();add(x,y,a,b,cnt1,ed,0);id[i]=i;}
    sort(id+1,id+1+m,cmp);
    int h=1;
    int ans=INF;int T=0;
    while(h<=m){
        int sta=ed[id[h]].a;
        for(;h<=m;h++){
            if(ed[id[h]].a!=sta) break;
            edge ne=ed[id[h]];
            add(ne.from,ne.to,ne.a,ne.b,cnt2,e,1);
            add(ne.to,ne.from,ne.a,ne.b,cnt2,e,1);
            if(!g[ne.from]) Q.push(ne.from),g[ne.from]=1;
            if(!g[ne.to]) Q.push(ne.to),g[ne.to]=1;
        }
        SPFA();
        ans=min(ans,sta+dis[n]);
    }
    if(dis[n]==INF) printf("-1\n");
    else printf("%d\n",ans);
}