题目链接:
题目大意:
有N个城市和M条道路,构成一条无向图。
在每个城市里边都有一个加油站,不同的加油站的价格都不一样。
通过一条道路的油耗就是该道路的边权。
现在你需要回答不超过100个问题,在每个问题中,请计算出一条油箱容量为C的车子,从起点S到终点T至少要花多少钱。
题解:
将每个点
i
i
i 拆成
101
101
101 个点:
V
(
i
,
j
)
V(i,j)
V(i,j)表示车子在点
i
i
i 且当前油量为
j
j
j 升,从
V
(
i
,
j
)
V(i,j)
V(i,j) 到
V
(
i
,
j
+
1
)
V(i,j+1)
V(i,j+1) 建边,边权为
p
p
pi。
对于原图中的每一条双向边
(
u
,
v
,
w
)
(u,v,w)
(u,v,w),对于所有大于等于
w
w
w 的
j
j
j ,从
V
(
u
,
j
)
V(u,j)
V(u,j) 到
V
(
v
,
j
−
w
)
V(v,j-w)
V(v,j−w) 建边,从
V
(
v
,
j
)
V(v,j)
V(v,j) 到
V
(
u
,
j
−
w
)
V(u,j-w)
V(u,j−w) 建边,边权为
0
0
0 。
对于每一个询问
(
c
,
s
,
t
)
(c,s,t)
(c,s,t) :以
V
(
s
,
0
)
V(s,0)
V(s,0) 为起点跑最短路,答案为
m
i
n
(
d
i
s
[
V
(
t
,
j
)
]
)
min(dis[V(t,j)])
min(dis[V(t,j)]) (
j
≤
c
j\leq c
j≤c)。
代码:
#include <bits/stdc++.h>
#define LL long long
#define LD long double
#define ULL unsigned long long
#define UI unsigned int
#define PII pair<int,int>
#define MPII(x,y) pair<int,int>{x,y}
#define _for(i,j,k) for(int i=j;i<=k;i++)
#define for_(i,j,k) for(int i=j;i>=k;i--)
#define efor(i,u) for(int i=head[u];i;i=net[i])
#define lowbit(x) (x&-x)
#define ls(x) x<<1
#define rs(x) x<<1|1
#define inf 0x3fffffff
//#pragma comment(linker, "/STACK:10240000000,10240000000")
using namespace std;
const int maxn = 2e6 + 5;
const int M = 1e9 + 7;
inline int mad(int a,int b){return (a+=b)>=M?a-M:a;}
int n,m;
struct E{
int v,w;
}e[maxn];
int head[maxn],net[maxn],cnt;
void add(int u,int v,int w){
e[++cnt]=E{v,w};
net[cnt]=head[u];
head[u]=cnt;
}
int dis[maxn];
priority_queue<PII,vector<PII>,greater<PII> > Q;
void dj(int u,int c){
_for(i,1,n+c*n) dis[i]=inf;
dis[u]=0;
Q.push(PII{0,u});
int w;
while(!Q.empty()){
u=Q.top().second;
w=Q.top().first;
Q.pop();
if(w>dis[u]) continue;
for(int i=head[u];i;i=net[i]){
if(e[i].v<=n+n*c&&dis[e[i].v]>dis[u]+e[i].w){
dis[e[i].v]=dis[u]+e[i].w;
Q.push(PII{dis[e[i].v],e[i].v});
}
}
}
}
int main(){
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
cin>>n>>m;
int u,v,w,p,q,c,s,t;
_for(i,1,n){
cin>>p;
_for(j,1,100){
add(i+(j-1)*n,i+j*n,p);
}
}
_for(i,1,m){
cin>>u>>v>>w;u++;v++;
_for(j,w,100){
add(u+j*n,v+(j-w)*n,0);
add(v+j*n,u+(j-w)*n,0);
}
}
cin>>q;
while(q--){
cin>>c>>s>>t;s++;t++;
dj(s,c);
int an=inf;
_for(i,0,c) an=min(an,dis[t+i*n]);
if(an==inf) cout<<"impossible\n";
else cout<<an<<"\n";
}
return 0;
}