思路:
首先预处理出cost数组,cost[i][j]表示从第i天到第j天满足这些天所有的条件能跑出来的最短路,如果没有路径,值为inf。
然后用
O(n2)
的动态规划求出答案。
状态转移方程:
f[i]=min{f[j]+cost[j+1][i]∗(i−j)+k}
因为一开始不需要转换路径,所以答案为
f[n]−k
代码:
#include <cstdio>
#include <iostream>
#include <queue>
using namespace std;
const int maxn = 1010;
const int inf = 1e9;
int n, m, k, e, p;
int head[25], next[maxn<<1], to[maxn<<1], dis[maxn<<1], cnt;
long long ok[25][110], cost[110][110], res[25], f[110];
bool can[25], vis[25];
void add(int a, int b, int c){
next[++cnt] = head[a];
head[a] = cnt;
to[cnt] = b;
dis[cnt] = c;
}
queue <int> q1;
int spfa(){
while(!q1.empty()) q1.pop();
for(int i = 0; i <= m; i ++) res[i] = inf, vis[i] = 0;
q1.push(1), res[1] = 0;
while(!q1.empty()){
int x = q1.front();
q1.pop(), vis[x] = 0;
for(int i = head[x]; i; i = next[i]){
if(!can[to[i]] && res[to[i]] > dis[i] + res[x]){
res[to[i]] = dis[i] + res[x];
if(vis[to[i]] == 0){
q1.push(to[i]);
vis[to[i]] = 1;
}
}
}
}
return res[m];
}
int main(){
scanf("%d%d%d%d", &n, &m, &k, &e);
for(int i = 1; i <= e; i ++){
int a, b, c;
scanf("%d%d%d", &a, &b, &c);
add(a, b, c), add(b, a, c);
}
scanf("%d", &p);
for(int i = 1; i <= p; i ++){
int a, b, c;
scanf("%d%d%d", &a, &b, &c);
ok[a][b] ++;
ok[a][c+1] --;
}
for(int i = 1; i <= m; i ++){
for(int j = 1; j <= n; j ++){
ok[i][j] += ok[i][j-1];
}
}
for(int i = 1; i <= n; i ++){
for(int j = 1; j <= n; j ++) cost[i][j] = inf;
for(int j = 1; j <= m; j ++) can[j] = false;
for(int j = i; j <= n; j ++){
for(int k = 1; k <= m; k ++)
if(ok[k][j] != 0) can[k] = true;
cost[i][j] = spfa();
if(cost[i][j] == inf) break;
}
}
for(int i = 1; i <= n; i ++) f[i] = inf;
for(int i = 1; i <= n; i ++){
for(int j = 0; j < i; j ++){
f[i] = min(f[i], f[j] + cost[j+1][i]*(i-j)+k);
}
}
printf("%d", f[n]-k);
return 0;
}