题面
传送门:洛谷
Solution
这题的思想很巧妙. 首先,我们可以考虑一下最暴力的做法,对每个时刻的所有点都求一遍单元最短路 因为最多只有200个时刻,时间复杂度为$O(n^3log(n+m)))$ (堆优化的迪杰斯特拉) 显然对于$n=200$,并过不了 我们可有进一步分析 这一题,我们堆优化的迪杰斯特拉慢在每加入一个点,我们每一次都得对全图彻彻底底做一轮松弛 那换个角度考虑,如果我只松弛经过新加入的点的点对呢? 没错,就得用Floyd了. 因为Floyd本质就是一个DP,给了我们极大的魔改的空间 考虑到Floyd最外层循环就是枚举加入的点,我们就可以只枚举里面那两层枚举点对的循环. 也就是说我们只用考虑它有可能松弛到的点. 当然,在此之前,我们得先把这个点有关的边先连回去 然后先用两层循环(枚举中转点和起始点)来松弛终点为加入点的路径 接下来用刚刚说的两层循环来松弛经过新加入点路径就好 时间复杂度$O(n^3)$ 然后就OjbK了 具体请看代码
Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
   | 
 
  #include<iostream> #include<cstdio> #include<cstring> using namespace std; const int N=200+10; long long read() {     long long x=0,f=1; char c=getchar();     while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}     while(isdigit(c)){x=x*10+c-'0';c=getchar();}     return x*f; } int n,m,T[N],dis[N][N],e[N][N]; int main() {     n=read(),m=read();     memset(T,0x3f,sizeof T);     memset(dis,0x3f,sizeof dis);     memset(e,0x3f,sizeof e);     for(int i=0;i<n;i++)         T[i]=read();     for(int i=1;i<=m;i++)     {         int a=read(),b=read(),temp=read();         e[a][b]=e[b][a]=temp;     }
      for(int i=0;i<n;i++)         e[i][i]=dis[i][i]=0;     int Q=read(),to=0;     for(int i=1;i<=Q;i++)     {         int x=read(),y=read(),t=read();         while(T[to]<=t)         {             for(int j=0;T[j]<=t;j++)                 dis[to][j]=dis[j][to]=min(dis[to][j],e[to][j]);             for(int j=0;T[j]<=t;j++)                 for(int k=0;T[k]<=t;k++)                     dis[to][k]=dis[k][to]=min(dis[k][to],dis[k][j]+dis[j][to]);             for(int j=0;T[j]<=t;j++)                 for(int k=0;T[k]<=t;k++)                     dis[j][k]=min(dis[j][k],dis[j][to]+dis[to][k]);             to++;         }         if(dis[x][y]==0x3f3f3f3f)             printf("-1\n");         else             printf("%d\n",dis[x][y]);     }     return 0; }
 
 
 
  |