先枚举任意两天,SPFA求出这两天之间的公共最短路,记为cost[i][j],这个可以把这两天之间的所有不能走到的点标记一下,然后求最短路.
求出了这个东西,就发现这个长得很像DP.设f[i]表示到i这天的最小费用.
f[i]=max(f[i],f[j]+cost[j+1][i]+k),f[i]初值为cost[1][i].
1 #include<set> 2 #include<map> 3 #include<queue> 4 #include<stack> 5 #include<ctime> 6 #include<cmath> 7 #include<string> 8 #include<vector> 9 #include<cstdio> 10 #include<cstdlib> 11 #include<cstring> 12 #include<iostream> 13 #include<algorithm> 14 #define maxn 22 15 #define inf 1999999999 16 using namespace std; 17 struct data{ 18 int nex,to,w; 19 }e[maxn*maxn*2]; 20 int head[maxn],edge=0,cost[110][110],bj[22][110],bj1[maxn],n,dis[maxn],f[110]; 21 bool vis[maxn]; 22 void add(int from,int to,int w){ 23 e[++edge].nex=head[from]; 24 e[edge].to=to; 25 e[edge].w=w; 26 head[from]=edge; 27 } 28 queue<int>q; 29 void SPFA(int s,int t){ 30 for(int i=1;i<=20;i++) dis[i]=inf,bj1[i]=0; 31 for(int i=s;i<=t;i++) 32 for(int j=1;j<=n;j++) 33 if(bj[j][i]) bj1[j]=1; 34 q.push(1);dis[1]=0; 35 while(!q.empty()){ 36 int u=q.front(); 37 q.pop(); 38 vis[u]=0; 39 for(int i=head[u];i;i=e[i].nex){ 40 int v=e[i].to; 41 if(bj1[v]) continue; 42 if(dis[v]>dis[u]+e[i].w){ 43 dis[v]=dis[u]+e[i].w; 44 if(!vis[v]) q.push(v),vis[v]=1; 45 } 46 } 47 } 48 if(dis[n]!=inf) 49 cost[s][t]=dis[n]*(t-s+1); 50 else cost[s][t]=inf; 51 } 52 int main() 53 { 54 freopen("!.in","r",stdin); 55 freopen("!.out","w",stdout); 56 int day,k,m,x,y,z; 57 scanf("%d%d%d%d",&day,&n,&k,&m); 58 for(int i=1;i<=m;i++) 59 scanf("%d%d%d",&x,&y,&z),add(x,y,z),add(y,x,z); 60 int d;scanf("%d",&d); 61 for(int i=1;i<=d;i++){ 62 scanf("%d%d%d",&x,&y,&z); 63 for(int i=y;i<=z;i++) 64 bj[x][i]=1; 65 } 66 for(int i=1;i<=day;i++) 67 for(int j=i;j<=day;j++) 68 SPFA(i,j);//,printf("%d %d %d ",i,j,cost[i][j]); 69 for(int i=1;i<=day;i++) f[i]=inf; 70 for(int i=1;i<=day;i++){ 71 f[i]=cost[1][i]; 72 for(int j=0;j<i;j++) 73 if(cost[j+1][i]<inf) 74 f[i]=min(f[i],f[j]+cost[j+1][i]+k); 75 } 76 printf("%d",f[day]); 77 return 0; 78 }