HDU6290 奢侈的旅行（堆优化的迪杰斯特拉算法）

最新推荐文章于 2024-03-26 10:27:41 发布

对你说的对

最新推荐文章于 2024-03-26 10:27:41 发布

阅读量269

点赞数 1

分类专栏：数据结构

本文链接：https://blog.csdn.net/jinghui_7/article/details/83986652

版权

数据结构专栏收录该内容

23 篇文章 0 订阅

订阅专栏

这是一道关于寻宝游戏和升级养成类游戏的算法题，玩家小Q需要从1号城镇到达n号城镇，途中经过的道路需要支付积分并提升等级。每条道路都有最低消费限制，若积分不足则无法通行。问题要求找到最小积分的路径，输出积分的整数部分。解决方案是转化成本问题，用堆优化的迪杰斯特拉算法求解，判断并选择能通过的最低消费路径。

摘要由CSDN通过智能技术生成

题目链接

Time Limit: 14000/7000 MS (Java/Others) Memory Limit: 512000/512000 K (Java/Others)
Total Submission(s): 1643 Accepted Submission(s): 336

Problem Description

高玩小Q不仅喜欢玩寻宝游戏，还喜欢一款升级养成类游戏。在这个游戏的世界地图中一共有n个城镇，编号依次为1到n。

这些城镇之间有m条单向道路，第i 条单项道路包含四个参数ui,vi,ai,bi，表示一条从ui号城镇出发，在vi号城镇结束的单向道路，因为是单向道路，这不意味着小Q可以从vi沿着该道路走到ui。小Q的初始等级level为1，每当试图经过一条道路时，需要支付cost=log2（（level+ai）/level）点积分，并且经过该道路后，小Q的等级会提升ai级，到达level+ai级。但是每条道路都会在一定意义上歧视低消费玩家，准确地说，如果该次所需积分cost<bi，那么小Q不能经过该次道路，也不能提升相应的等级。

注意：本游戏中等级为正整数，但是积分可以是任意实数。

小Q位于1号城镇，等级为1，现在为了做任务要到n号城镇去。这将会是一次奢侈的旅行，请写一个程序帮助小Q找到需要支付的总积分最少的一条路线，或判断这是不可能的。

Input

第一行包含一个正整数T(1≤T≤30)，表示测试数据的组数。

每组数据第一行包含两个整数n,m(2≤n≤100000,1≤m≤200000)，表示城镇数和道路数。

接下来m行，每行四个整数ui,vi,ai,bi(1≤ui,vi≤n,ui≠vi,0≤ai≤109,0≤bi≤60)，分别表示每条单向道路。

Output

对于每组数据，输出一行一个整数，即最少所需的总积分的整数部分，如：4.9999输出4，1.0输出1。若不存在合法路线请输出−1。

Sample Input

3 3

1 2 3 2

2 3 1 6

1 3 5 0

Sample Output

思路：

将求花费的公式转化为 $1+a_{i}/level=2^{cost}$ ,将所有的log去掉，转化为指数形式便于比较。有公式可以看出，level越大，cost越低，该条路不能通过的可能性越大。

判断该条路能否通过， $2^{cost}>=2^{b_{i}}$

对于要求的最后的总花费，假设最后选了k条路，在这k条路上的等级提升分别为 $a_{1},a_{2},a_{3},......,a_{k}$ ,总花费为 $\log_{2}(1+a_{1})+\log_{2}((1+a_{1}+a_{2}+a_{3})/(1+a_{1}+a_{2}))+\log_{2}$ $((1+a_{1}+a_{2}+a_{3}+a_{4})/(1+a_{1}+a_{2}+a_{3}))$ $+......+\log_{2}((1+a_{1}+a_{2}+......+a_{k})/(1+a_{1}+a_{2}+......+a_{k-1}))$ 经过化简得最后花费的积分为 $\log_{2}(1+a_{1}+a_{2}+a_{3}+......+a_{k})$ ,去掉log，得 $1+a_{1}+a_{2}+a_{3}+......+a_{k}$ ，求的是 $2^{cost}$ ，同样是等级越低越好。

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<ll,ll>Node;

const int maxn=300010;
const ll INF=0x3f3f3f3f3f3f3f3f;

ll Pow[65],dis[maxn],head[maxn],n,m,cnt;

struct Edge
{
    ll to,next;
    ll a,b;
}edge[maxn];
void init()
{
    for(int i=0;i<65;i++)
    {
        if(i==0)
            Pow[i]=1;
        else
            Pow[i]=Pow[i-1]<<1;
    }
    memset(head,-1,sizeof(head));
    cnt=0;
}
void add_edge(ll from,ll to,ll a,ll b)
{
    edge[cnt].to=to;
    edge[cnt].a=a;
    edge[cnt].b=b;
    edge[cnt].next=head[from];
    head[from]=cnt++;
}

void dijk(ll x)
{
    memset(dis,INF,sizeof(dis));
    priority_queue<Node,vector<Node>,greater<Node> >q;
    dis[x]=1;//到达1点花费的积分
    q.push(Node(1,x));
    while(!q.empty())
    {
        Node now=q.top();
        q.pop();
        if(dis[now.second]<now.first) continue;
        ll id=now.second;
        for(int i=head[id];i!=-1;i=edge[i].next)
        {
            ll to=edge[i].to;
            ll aa=edge[i].a;
            ll bb=edge[i].b;
            if(aa/now.first>=bb-1&&now.first+aa<dis[to])//添加条件，积分>=bb时才可以通过这条路
            {//判断是否连通
                dis[to]=now.first+aa;
                q.push(Node(dis[to],to));
            }
        }
    }
    if(dis[n]==INF)
        puts("-1");
    else
        printf("%lld\n",(ll)log2(dis[n]));
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        init();
        for(int i=0;i<m;i++)
        {
            ll a,b,c,d;
            scanf("%lld%lld%lld%lld",&a,&b,&c,&d);
            add_edge(a,b,c,Pow[d]);
        }
        dijk(1);
    }
    return 0;
}