现在有一个 m m m个点 e e e条边的无向连通图,总共有 n n n天,你每天都选取一条路径使得 1 1 1到 m m m的花费最小,可是在 [ a i , b i ] [a_i,b_i] [ai,bi]天,某个节点 d i d_i di不能走,现在你要求这 n n n天总花费最少为多少。
我们先标记编号为 i i i的点是否能在 j j j天行得通,然后预处理出从第 i i i天到第 j j j天在不改变路径情况下的最短路。
接下来进行DP:
我们令 d p i dp_i dpi表示从第一天到第 i i i天的最少花费。那么我们枚举更改路径的时间 j j j,那么既有:
d p i = m i n ( d p i , d p j − 1 + K + c o s t i j ∗ ( j − i ) + 1 ) dp_i=min(dp_i,dp_{j-1}+K+cost_{ij}*(j-i)+1) dpi=min(dpi,dpj−1+K+costij∗(j−i)+1).
#include
using namespace std;
const int N = 300;
typedef long long ll;
const int INF = 0x3f3f3f3f;
ll MOD = 1e9+7;
ll unp[N][N];//
vector<pair<ll,ll> > G[N];
struct node{
ll id;
ll DIS;
};
bool operator<(node a,node b){
return a.DIS>b.DIS;
}
ll n,m,K,e;
ll dij(ll s,ll t){
bool vis[N];
bool uu[N];
int dis[N];
memset(dis,0x3f,sizeof(dis));
memset(vis,0,sizeof(vis));
memset(uu,0,sizeof(uu));
priority_queue<node> q;
q.push({1,0});dis[1]=0;
for(ll i=1;i<=m;i++){
for(ll j=s;j<=t;j++){
if(unp[i][j]){
uu[i]=1;break;//i此时不能用
}
}
}
while(!q.empty()){
node rt=q.top();q.pop();
ll u=rt.id;
if(vis[u]) continue;
if(uu[u]) continue;
vis[u]=1;
for(int j=0;j<(int)G[u].size();j++){
int v=G[u][j].first;int DD=G[u][j].second;
if(!vis[v]&&!uu[v]){
if(dis[v]>dis[u]+DD){
dis[v]=dis[u]+DD;
q.push({v,dis[v]});
}
}
}
}
return dis[m];
}
ll dis[N][N];
ll dp[N];
int main(){
cin>>n>>m>>K>>e;
ll u,v,w;
for(ll i=1;i<=e;i++){
cin>>u>>v>>w;
G[u].push_back({v,w});
G[v].push_back({u,w});
}
ll d;
cin>>d;
while(d--){
ll id,a,b;
cin>>id>>a>>b;
for(ll i=a;i<=b;i++) unp[id][i]=1; //第id个港口第i天不能使用
}
for(ll i=1;i<=n;i++){
for(ll j=1;j<=n;j++){
dis[i][j]=dij(i,j);
}
}
for(ll i=1;i<=n;i++){
dp[i]=dis[1][i]*i;
for(ll j=1;j<=i;j++){
dp[i]=min(dp[i],dp[j-1]+K+(i-j+1)*dis[j][i]);
}
}
cout<<dp[n]<<endl;
return 0;
}
``