2018HDU多校补题 HDU 6331（dp）

题意：

给定$n$为$50$的有向图，给出$1e6$的问题，问从$s$到$t$至少走$k$步的最短路

题解：

定义$dp[st][i][j]$表示恰好走st步从i到j的最短路。
那么我们输入的值则是$dp[1][u][v]$
之后就可以想到一个很容易想到的转移 类似于flord算法 dp[st][i][j]=dp[st-1][i][k]+dp[1][k][j]
这样我们发现时间复杂度为$m\ast n^3$算上T组基本会$TEL$，因此对于$m$进行根号分块，分别表示恰好走$st$步和每$st\ast 100$步的最短路，最后从大步数去更新小步数来使得数组定义变为至少走st步

代码：

/**
 *     author:     TelmaZzzz
 *     create:     2019-07-10-13.57.23
**/
#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <ctime>
#include <string>
#include <vector>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <random>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef double db;
#define PB push_back
#define MP make_pair
#define INF 12331313312000000
#define inf 1152921504606846976
#define pi 3.14159265358979323846
//#pragma comment(linker,"/STACK:10240000,10240000")
mt19937 rand_(time(0));
const int N=3e5+7,M=2e6;
const long long mod=1e9+7;
inline int read(){int ret=0;char ch=getchar();bool f=1;for(;!isdigit(ch);ch=getchar()) f^=!(ch^'-');for(;isdigit(ch);ch=getchar()) ret=(ret<<1)+(ret<<3)+ch-48;return f?ret:-ret;}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
ll ksm(ll a,ll b,ll mod){int ans=1;while(b){if(b&1) ans=(ans*a)%mod;a=(a*a)%mod;b>>=1;}return ans;}
ll inv2(ll a,ll mod){return ksm(a,mod-2,mod);}//逆元
//int head[N],NEXT[M],ver[M],tot;void link(int u,int v){ver[++tot]=v;NEXT[tot]=head[u];head[u]=tot;}



ll a[110][53][53],b[110][53][53];
int main(){
    //freopen("1.txt","r",stdin);
    //ios::sync_with_stdio(false);
    int t;
    scanf("%d",&t);
    int n,m,u,v,w;
    while(t--){
        scanf("%d%d",&n,&m);
        for(int k=0;k<110;k++){
            for(int i=1;i<=n;i++){
                for(int j=1;j<=n;j++){
                    a[k][i][j]=b[k][i][j]=INF;
                    if(i==j&&k==0) a[k][i][j]=b[k][i][j]=0;
                }
            }
        }
        for(int i=1;i<=m;i++){
            scanf("%d%d%d",&u,&v,&w);
            a[1][u][v]=min(a[1][u][v],1LL*w);
        }
        for(int p=2;p<=100;p++){
            for(int k=1;k<=n;k++){
                for(int i=1;i<=n;i++){
                    for(int j=1;j<=n;j++){
                        a[p][i][j]=min(a[p][i][j],a[p-1][i][k]+a[1][k][j]);
                    }
                }
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                b[1][i][j]=a[100][i][j];
            }
        }
        for(int p=2;p<=100;p++){
            for(int k=1;k<=n;k++){
                for(int i=1;i<=n;i++){
                    for(int j=1;j<=n;j++){
                        b[p][i][j]=min(b[p][i][j],b[p-1][i][k]+b[1][k][j]);
                    }
                }
            }
        }
        for(int p=100;p>=0;p--){
            for(int i=1;i<=n;i++){
                for(int j=1;j<=n;j++){
                    a[p][i][j]=min(a[p][i][j],a[p+1][i][j]);
                    b[p][i][j]=min(b[p+1][i][j],b[p][i][j]);
                }
            }
        }
        int q;
        scanf("%d",&q);
        while(q--){
            scanf("%d%d%d",&u,&v,&w);
            int A=w%100;
            int B=w/100;
            ll ans=INF;
            for(int i=1;i<=n;i++){
                ans=min(ans,b[B][u][i]+a[A][i][v]);
                ans=min(ans,b[B+1][u][i]+a[0][i][v]);
            }
            if(ans==INF) puts("-1");
            else printf("%lld\n",ans);
        }
    }
    //cout << "time: " << (long long)clock() * 1000 / CLOCKS_PER_SEC << " ms" << endl;
    return 0;
}