传送门
题意:给定有向图,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;
}