Axel and Marston in Bitland CodeForces - 782F (bitset优化)

题目链接

 

$dp[0/1][i][x][y]$表示起始边为0/1, 走$2^i$ 步, 是否能从$x$走到$y$

则有转移方程
$dp[z][i][x][y]\mid=dp[z][i-1][x][k]\&dp[z\wedge1][i-1][k][y]$

复杂度 $O(k_0n^3)$,  其中$k_0=log(1e18)$

这里可以用bitset优化第四维的递推, bitset底层相当于若干个64bit数, 可以优化64的常数
复杂度$O(\frac{k_0n^3}{\omega})$

 

#include <iostream>
#include <bitset>
#include <cstdio>
#define PER(i,a,n) for(int i=n;i>=a;--i)
#define REP(i,a,n) for(int i=a;i<=n;++i)
using namespace std;
typedef long long ll;

int n, m;
bitset<555> f[2][66][555], pre, now;

int main() {
    scanf("%d%d", &n, &m);
    REP(i,1,m) {
        int u, v, w;
        scanf("%d%d%d", &u, &v, &w);
        f[w][0][u][v] = 1;
    }
    REP(i,1,60) REP(z,0,1) REP(x,1,n) REP(y,1,n) {
        if (f[z][i-1][x][y]) f[z][i][x]|=f[z^1][i-1][y];
    }
    ll ans = 0, z = 0;
    pre[1] = 1;
    PER(i,0,60) {
        now.reset();
        REP(j,1,n) if (pre[j]) now|=f[z][i][j];
        if (now.count()) z^=1, ans^=1ll<<i, pre=now;
    }
    printf("%lld\n",ans>1e18?-1:ans);
}

 

转载于:https://www.cnblogs.com/uid001/p/10229692.html