BZOJ 1050: [HAOI2006]旅行comf|动点spfa

最新推荐文章于 2020-05-20 22:18:27 发布

ws_yzy

最新推荐文章于 2020-05-20 22:18:27 发布

阅读量705

点赞数

分类专栏：动点spfa 最短路文章标签：动点spfa

本文链接：https://blog.csdn.net/ws_yzy/article/details/50666797

版权

最短路同时被 2 个专栏收录

7 篇文章 0 订阅

订阅专栏

动点spfa

0 篇文章 0 订阅

订阅专栏

用Kruscal和并查集做似乎很显然，动点spfa似乎也是裸题。
按边权把边从大到小排序，枚举不断加边，将当前边的两端点入队spfa，不断更新mn数组（在当前边的情况下从起点到每一个点所经过的最大边的最小值）
因为从大到小枚举，所以当前枚举到的边肯定是最短的，然后要让最长的边最短，这就是动点spfa的精髓所在。
第一次写动点spfa之蜜汁错误之打反n，m毁一生！！！！！！
woc….调了半小时没找出n，m打反..手残毁一生

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<iostream>
#include<algorithm>
#define N 111111
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>'9'||c<'0'){if(f=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i;
}
struct W{int x,y,z;}a[N];
int head[N],nxt[N],lst[N],v[N],q[N];
int inq[N],mn[N];
int tot,n,m,S,T;
int l,r,A,B;
double ans=1e8;
bool cmp(W a,W b){return a.z>b.z;}
void insert(int x,int y,int z)
{
    lst[++tot]=y;
    nxt[tot]=head[x];
    v[tot]=z;
    head[x]=tot;
}
void spfa()
{
    for(;l!=r;l=(l+1)%N)
    {
        int x=q[l]; inq[x]=0;
        for(int i=head[x];i;i=nxt[i])
            if(mn[lst[i]]>max(mn[x],v[i]))
            {
                mn[lst[i]]=max(mn[x],v[i]);
                if(!inq[lst[i]])
                    inq[q[r++]=lst[i]]=1,
                    r%=N;
            }
    }
}
int main()
{
    n=sc(),m=sc();
    for(int i=1;i<=m;i++)
    {
        a[i].x=sc();
        a[i].y=sc();
        a[i].z=sc();
    }
    memset(mn,0x3f,sizeof(mn));
    sort(a+1,a+m+1,cmp);
    mn[S=sc()]=0;
    T=sc();
    for(int i=1;i<=m;i++)
    {
        int x=a[i].x,y=a[i].y,z=a[i].z;
        insert(x,y,z);
        insert(y,x,z);
        if(mn[x]<mn[y]) swap(x,y);
        if(mn[x]<=max(mn[y],z))  
            continue;
        mn[x]=max(mn[y],z);
        q[r++]=x; r%=N;
        spfa();
        if((double)mn[T]/z<ans)
        {
            ans=(double)mn[T]/z;
            A=mn[T];
            B=z;
        }
    }
    if(ans>1e6)puts("IMPOSSIBLE");
    else
    {
        if(A%B==0)cout<<A/B<<endl;
        else
            printf("%d/%d\n",A/__gcd(A,B),B/__gcd(A,B));
    }
    return 0;
}