CF1528D(n² dij+增边优化建图)

13 篇文章 1 订阅
5 篇文章 0 订阅

传送门

题意:给定有向图,n个点,m条边。第i条边在0时刻由a[i]->b[i]代价c[i] 。每过t时刻,第i条边变成由a[i]->(b[i]+t)%n ,代价仍是c[i]
对于每个出发点0~n-1,求出到其余点的最短路。
(2<=n<=600,n<=m<=n²)

分析:

  • 首先肯定是考虑跑dijkstra,这题的边数会达到n²,堆优化的dij时间的mlogn,所以得用普通版本的dij,时间复杂度是O(n²)。
  • 每个点在松弛过程中,由于我们可以在某个点停留,所以需要枚举停留时间和相邻点,这样松弛相邻边的时间就是O(n²),单次最短路变成了O(n³)。
  • 优化:我们可以对每个点连一条i -> (i+1)%n 代价为1的假边,表示i的上一个点可以停留1s然后走到 (i+1)%n ,我们保证第一步走的是真边,后面走真边假边都可以了。
  • 坑点:出发点不能一开始就更新为0,而是更新成INF,因为他会作为桥梁更新其他点的最短路。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
const int maxn = 6e2+7;
const int mx = 40;
const int mod = 1e9+7;
const ll inf = 34359738370;
const int INF = 2e9+7;
//给定无向图 m条边
//0 s的时候 边都由a[i] -> b[i] 花费时间 c[i] 
//每过一秒 边i a[i] -> (b[i]+1)%n
//对每个顶点求到其他点的最短路
//对于和出发点st 0s时刻相连的边 可以直接更新
//后面在松弛其他点的时候 可以等x秒在走
//如果我们暴力枚举时间+顶点 会达到N²时间复杂度 
//我们可以连一条假边:i (i+1)%n 1 
//表示i的上一个顶点出发的时候 等一秒  然后走到 (i+1)%n 
int g[maxn][maxn];
int n,m;
int dist[maxn];
bool vis[maxn];
void dj(int st) 
{
    for(int i=0;i<n;i++)
    {
        dist[i]=g[st][i];//第一步走真边
        vis[i]=0;
    }
    dist[st]=INF;//坑 st不能一开始就变成0 它还得作为桥梁更新其他点
    while(1) 
    {
        int index=-1;
        for(int i=0;i<n;i++) 
        {
            if(!vis[i] && (index==-1 || dist[i]<dist[index])) index=i;
        }
        if(index == -1) break;
        vis[index]=1;
        //假边 v -> (v+1)%n 1
        int to=(index+1)%n;
        if(!vis[to]) dist[to]=min(dist[to],dist[index]+1);
        for(int i=0;i<n;i++) 
        {
            if(g[index][i] == INF) continue;
            int to=(i+dist[index])%n;
            if(!vis[to] && dist[to] > dist[index]+g[index][i]) dist[to]=dist[index]+g[index][i];
        }
    }
    dist[st]=0;//st到st的实际最短路
    printf("%d",dist[0]);
    for(int i=1;i<n;i++) printf(" %d",dist[i]); 
    printf("\n");
}
int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++) 
    {
        for(int j=1;j<=n;j++) g[i-1][j-1]=INF;
    }
    for(int i=1;i<=m;i++) 
    {
        int a,b,c;
        scanf("%d %d %d",&a,&b,&c);
        g[a][b]=min(g[a][b],c);
    }
    for(int i=0;i<n;i++) dj(i);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值