BZOJ#1003物流运输
思路:
动态规划,f[i]表示以i结尾的最短时间,f[i]=min(f[j]+c[j+1][i]*(i-j)+k)
c[i][j]表示从时间i到j的最短路
代码:
#include <bits/stdc++.h>
#define int long long
int _ = 0, Case = 1;
using namespace std;
#define all(v) begin(v),end(v)
#define nline '\n'
const int N = 210;
int f[N];
int c[N][N];
struct T {
int a, b, w;
} edge[N];
using PII = pair<int, int> ;
vector<PII> h[N];
int Time[N][N];
int n, m, e, k;
int dist[N], vis[N];
int dijkstra(){
for(int i=1;i<=m;i++) dist[i]=1e8,vis[i]=false;
dist[1]=0;
priority_queue<PII,vector<PII>,greater<PII>> q;
q.push({0,1});
while(q.size()){
auto [x,y]=q.top();
q.pop();
if(vis[y]) continue;
vis[y]=true;
for(auto [id,w]:h[y]){
if(dist[id]>x+w){
dist[id]=x+w;
q.push({dist[id],id});
}
}
}
return dist[m];
}
void init() {
for (int len = 1; len <= n; len++) {
for (int l = 1; l + len - 1 <= n; l++) {
for (int i = 1; i <= m; i++) {
h[i].clear();
}
int r = l + len - 1;
for (int k = 1; k <= e; k++) {
auto[a, b, w] = edge[k];
bool ok = true;
for(int i=l;i<=r;i++){
if(Time[a][i] or Time[b][i]) ok=false;
}
if (ok) {
h[a].push_back({b, w});
h[b].push_back({a, w});
}
}
c[l][r] = dijkstra();
}
}
}
void solve(int Case) {
cin >> n >> m >> k >> e;
for (int i = 1; i <= e; i++) {
auto&[a, b, w] = edge[i];
cin >> a >> b >> w;
}
int d;
cin >> d;
for (int i = 1; i <= d; i++) {
int p, a, b;
cin >> p >> a >> b;
for(int j=a;j<=b;j++){
Time[p][j]=1;
}
}
init();
memset(f, 0x3f, sizeof f);
f[0] = 0;
for (int i = 1; i <= n; i++) {
for (int j = 0; j < i; j++) {
f[i] = min(f[i], f[j] + c[j + 1][i] * (i - j) + k);
}
}
cout << f[n] - k << nline;
}
signed main() {
ios::sync_with_stdio(false); cin.tie(nullptr);
// cin >> _; for (Case = 1; Case <= _; Case++)
solve(Case);
return 0;
}