P1266 速度限制 ( dij分层 + 图上dp

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
using VI = vector<int>;
using ll = long long;
using PII = pair <int , int>;
const int mod = 998244353;
int n,m,d;
int idx = 0 ;
struct edge{
     int to,from;
     int next;
     int v,l;
}g[200010];
int head[200010];
void add(int a , int b , int v , int l){
    idx++;
    g[idx].from = a;
    g[idx].to = b;
    g[idx].v = v;
    g[idx].l = l;
    g[idx].next = head[a];
    head[a] = idx;
}
double dis[200][510];//达到某个点的时间 , 速度为j
int vis[200][510];
PII fa[200][510];
double calc(int x, int y){
    return x * 1.0 / (y * 1.0);
}
int finalV;


void dij(int s){
    //时间,速度,点
    priority_queue<pair<double,PII>> pq;
    memset(dis , 126  ,sizeof dis);
    //printf("%lf\n" , dis[2][25]);
    dis[1][70] = 0;
    fa[1][70] = {0,0};
    pq.push({0,{70,1}});
    while(pq.size()){
        auto x = pq.top().second;
        pq.pop();
        int v = x.first;
        int now = x.second;
        if(vis[now][v]) continue;
        vis[now][v] = 1;
        for(int i = head[now] ; i ; i = g[i].next){
            int to = g[i].to;
            int l = g[i].l;
            if(g[i].v == 0){
                if(dis[to][v] > dis[now][v] + calc(l,v)){
                    dis[to][v] = dis[now][v] + calc(l,v);
                    fa[to][v] = {now , v};
                    pq.push({-dis[to][v] ,{v,to}});
                }
            }else{
                if(dis[to][g[i].v] > dis[now][v] + calc(l ,g[i].v)){
                    dis[to][g[i].v] = dis[now][v] + calc(l ,g[i].v);
                    fa[to][g[i].v] = {now , v};
                    pq.push({-dis[to][g[i].v],{g[i].v , to}});
                }
            }

        }
    }
}

int output(int p , int v){
    if(p == 0) return 0;
    output(fa[p][v].first , fa[p][v].second);
    printf("%d " , p - 1);
    return 0;
}




int main(){
    cin>>n>>m>>d;
    d++;
    for(int i = 1 ; i <= m ; i++){
        int a,b,v,l;
        cin>>a>>b>>v>>l;
        a++,b++;
        add(a,b,v,l);
    }

    dij(1);
    int t = -1;
    for(int i = 0 ; i <= 510 ; i++){
        if(t == -1 || dis[d][t] > dis[d][i]){
            t = i;
        }
    }
    //cout<<d<<"" <<t<<"\n";
    //printf("%lf\n", dis[d][t]);
    //printf("%d\n", fa[d][t].first);
    //cout<<dis[d][t];
    output(d , t);



}

其实并不是真的建立分层图,只是有分层图的思想,或者是dp思想在里面

dp[i][j] // 到达 i 点速度是 j 的花费

v == 0 dp[i][x] = min( dp[u][x] + calc(len[u][i] / x) )

v != 0 dp[i][v] = min( dp[u][x] + calc(len[u][i] / v) )

相关推荐
SweetCode9 分钟前
裴蜀定理:整数解的奥秘
数据结构·python·线性代数·算法·机器学习
ゞ 正在缓冲99%…23 分钟前
leetcode76.最小覆盖子串
java·算法·leetcode·字符串·双指针·滑动窗口
xuanjiong23 分钟前
纯个人整理,蓝桥杯使用的算法模板day2(0-1背包问题),手打个人理解注释,超全面,且均已验证成功(附带详细手写“模拟流程图”,全网首个
算法·蓝桥杯·动态规划
惊鸿.Jh42 分钟前
【滑动窗口】3254. 长度为 K 的子数组的能量值 I
数据结构·算法·leetcode
明灯L43 分钟前
《函数基础与内存机制深度剖析:从 return 语句到各类经典编程题详解》
经验分享·python·算法·链表·经典例题
碳基学AI1 小时前
哈尔滨工业大学DeepSeek公开课:探索大模型原理、技术与应用从GPT到DeepSeek|附视频与讲义免费下载方法
大数据·人工智能·python·gpt·算法·语言模型·集成学习
补三补四1 小时前
机器学习-聚类分析算法
人工智能·深度学习·算法·机器学习
独好紫罗兰1 小时前
洛谷题单3-P5718 【深基4.例2】找最小值-python-流程图重构
开发语言·python·算法
正脉科工 CAE仿真1 小时前
基于ANSYS 概率设计和APDL编程的结构可靠性设计分析
人工智能·python·算法
Dovis(誓平步青云)2 小时前
【数据结构】排序算法(中篇)·处理大数据的精妙
c语言·数据结构·算法·排序算法·学习方法