飞行计划(最短路有点转换)

飞行计划(最短路有点转换)

题目描述

已知地图上一共有n个机场，Alice要从1号机场飞到n号机场。

从i号机场飞行到j号机场需要消耗dist(i,j)个单位的油量，其中 dist(i,j)为i号机场和j号机场之间的欧几里得距离。

Alice可以选择在任意机场起降和加油。飞机油箱最多可以保存u单位的燃油。

考虑到可能遇到的极端天气，某些机场在某些时候不具备降落的条件。所以出于安全的考虑，当他飞到一个机场的时候，必须保证飞机剩下的油足以飞到离他最近的机场。

假设一次起飞和降落分别需要花费a秒和b秒的时间，加1单位的油需要花费c秒的时间，飞行1单位的距离需要花费d秒的时间，请问 Alice最早什么时候可以到达终点？

飞机初始的油量为0，也就是说，飞机需要先在起点加油再起飞。

输入

第一行6个整数n,a,b,c,d,u。
接下来n行，每行两个整数xi,yi表示第i个机场的坐标。

输出

对于每组数据输出一行一个数，表示最早的到达时间，只要你的答案和标准答案相对误差或者绝对误差在10 ^-5以内就算通过

数据保证有解

样例输入 Copy

4 1 1 1 1 3
1 0
2 0
3 0
4 0

 

样例输出 Copy

11

 

 

题解：粗体部分。

当从第一个机场起飞。花费为起飞时间+降落时间+加油时间+飞行时间。

飞行时间为两点的欧几里得距离。

加油时间由于粗体字可知，为（离目的地最近的那个点的距离+此时两个点的距离）*c;

当飞出第一个点外。我们还是需要多加油，但是前面的点我们已经足够飞到当前点，所以我们的代价位（离目的地最近的那个点的距离+此时两个点的距离-离起点最近的距离）*c;

 

注意 double

 

#pragma GCC optimize(3 , "Ofast" , "inline")

#include <bits/stdc++.h>

#define rep(i , a , b) for(register int i=(a);i<=(b);i++)
#define rop(i , a , b) for(register int i=(a);i<(b);i++)
#define per(i , a , b) for(register int i=(a);i>=(b);i--)
#define por(i , a , b) for(register int i=(a);i>(b);i--)

using namespace std;
int n,s;
typedef long long ll;
typedef pair<double ,double > pi;

const int N =1111;
const double eps = 1e-8;
double dist[N];
double dis[N];
pi p[N];
ll inf=ll(1e10);
bool inque[N];
struct edge{
    int v;
    double time;
};

struct youxian{
    int index;
    double dist;
    bool operator<(const youxian &tmp)
    const {
        return dist>tmp.dist;
    }
};
vector<edge>e[N];

void dil()
{
    int i;
    for(i=1;i<=n;i++)
    {
        dis[i]=inf;
        inque[i]=false;
    }
    s=1;
    dis[s]=0;
    priority_queue<youxian>Q;
    Q.push({s,0});
    while(!Q.empty()){
        int u=Q.top().index;
        Q.pop();
        if(inque[u]) continue;
        inque[u]=true;
        int size=int (e[u].size());
        for(i=0;i<size;i++){
            int v=e[u][i].v;
            double time=e[u][i].time;
            if(!inque[v]&&dis[v]>dis[u]+time+eps){
                dis[v]=dis[u]+time;
                Q.push({v,dis[v]});
            }
        }
    }
}


double dd(int x,int y) {
    return sqrt ((p[x].first-p[y].first)*(p[x].first-p[y].first)+(p[x].second-p[y].second)*(p[x].second-p[y].second));
}

int main()
{
    ios::sync_with_stdio (false);
    double a,b,c,d,u;
    cin>>n>>a>>b>>c>>d>>u;
    rep (i,1,n) {
        cin>>p[i].first>>p[i].second;
        dist[i]=inf;
    }
    rep (i,1,n) {
        rep (j,1,n) {
            if(i!=j) {
                dist[i]=min (dist[i],dd(i,j));
            }
        }
    }
    rep(i,2,n) {
        double di = dd (1,i);
        if(di+dist[i]>u) continue;
        e[1].push_back ({i,a+b+d*di+c*(di+dist[i])});
    }
    rep(i,2,n) {
        rep (j,2,n) {
            double di = dd (i,j);
            if(di+dist[j]>u) continue;
            e[i].push_back ({j,a+b+d*di+c*(di+dist[j]-dist[i])});
        }
    }
    dil ();
    printf ("%f\n",dis[n]);
    return 0;
}