洛谷 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;
}