bzoj4773

矩阵乘法

。。。爆零了。。。

想到Floyd，却不知道怎么限制点数。。。

其实我们只要给Floyd加一维，dp[i][j][k]表示当前走过了i个点，从j到k的最短距离，然后这样可以倍增，最后看是否有i->i的距离<0

做dp或最短路之类的题的时候，如果限制条件较多，可以考虑加维度

#include<bits/stdc++.h>
using namespace std;
const int N = ;
int n, m, ans;
int Log[N];
struct matrix {
    int a[N][N];
    matrix() { memset(a, 0x3f3f, sizeof(a)); }
    friend matrix operator * (const matrix &a, const matrix &b) {
        matrix ret;
        for(int k = ; k <= n; ++k)
            for(int i = ; i <= n; ++i)
                for(int j = ; j <= n; ++j) ret.a[i][j] = min(ret.a[i][j], a.a[i][k] + b.a[k][j]);
        return ret;
    }
} A[N], tmp;
bool judge(const matrix &a)
{
    for(int i = ; i <= n; ++i) if(a.a[i][i] < ) return true;
    return false;
}
int main()
{
    scanf("%d%d", &n, &m);
    for(int i = ; i <= n; ++i) A[].a[i][i] = tmp.a[i][i] = ;
    for(int i = , u, v, w; i <= m; ++i) scanf("%d%d%d", &u, &v, &w), A[].a[u][v] = w;
    Log[] = ;
    for(int i = ; i <= n; ++i) Log[i] = Log[i >> ] + ;
    for(int i = ; i <= Log[n]; ++i) A[i] = A[i - ] * A[i - ];
    for(int i = Log[n]; i >= ; --i) if(!judge(A[i] * tmp)) tmp = tmp * A[i], ans +=  << i;
    tmp = tmp * A[];
    printf("%d\n", judge(tmp) ? ans +  : );
    return ;
}