hdu6290 - 奢侈的旅行 - dijk+堆优化（详细）

最新推荐文章于 2024-09-05 20:21:40 发布

努力写题的tyf

最新推荐文章于 2024-09-05 20:21:40 发布

阅读量232

点赞数

分类专栏：最短路文章标签： dijk+堆优化数学

本文链接：https://blog.csdn.net/m0_37579232/article/details/90143145

版权

最短路专栏收录该内容

19 篇文章 0 订阅

订阅专栏

思路：

qwq，本菜菜第一次遇到卡spfa的题，据说要用dijk+堆优化.

我们可以发现，若从1走到n经过的是1 2……n点的话，距离是：

$log2(\frac{1+a_2}{1})+log2(\frac{1+a_2+a_3}{1+a_2})+...+log2(\frac{1+a_2+...+a_n}{1+a_2+...+a_{n-1}})$ =log2(1+a_2+...+a_n)

可以看出，到i点的距离只和ai相关，即把ai看做是边的权重

还有一点要注意的是：log会丢失精度，我们先不取log，最后再取

还有就是每次等级都提升ai，那么用u点更新v点时，判断条件是：

$log2(\frac{level_u+a[v]}{level_u})$ 与b的关系

可知levelu=d[u]，所以直接用d[u]求就可以啦

（写题的时候多想想，log的性质）

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cstring>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>
#include<functional>
using namespace std;
#define ll long long

typedef pair<ll,int>P;
const ll INF=1e17+10;
const int N=100005,mod=1e9+7;
int cnt,head[N],n;
ll d[N];

struct A{
    int to,nex,a,b;
}edge[2*N];

void addEdge(int from,int to,int a,int b){
    edge[cnt].to=to;
    edge[cnt].nex=head[from];
    edge[cnt].a=a;
    edge[cnt].b=b;
    head[from]=cnt++;
}

void dijk(){
    priority_queue<P,vector<P>,greater<P> >Q;
    for(int i=1;i<=n;i++){
        d[i]=INF;
    }
    d[1]=1;
    Q.push(P(d[1],1));
    while(!Q.empty()){
        P p=Q.top();Q.pop();
        int u=p.second;
        if(d[u]<p.first)continue;
        for(int i=head[u];i!=-1;i=edge[i].nex){
            int v=edge[i].to;
            if(log2((d[u]*1.0+edge[i].a*1.0)/(d[u]*1.0))<edge[i].b*1.0){
                continue;
            }
            if(d[v]>d[u]+edge[i].a){
                d[v]=d[u]+edge[i].a;
                Q.push(P(d[v],v));
            }
        }
    }
}

int main(){
    int t,m;
    scanf("%d",&t);
    while(t--){
        memset(head,-1,sizeof(head));
        cnt=0;
        scanf("%d%d",&n,&m);
        while(m--){
            int u,v,a,b;
            scanf("%d%d%d%d",&u,&v,&a,&b);
            addEdge(u,v,a,b);
        }
        dijk();
        if(d[n]==INF){
            printf("-1\n");
        }
        else printf("%lld\n",(ll)log2(d[n]));
    }
}