結果

問題 No.1400 すごろくで世界旅行
ユーザー yosswi414yosswi414
提出日時 2021-02-19 02:52:55
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 653 ms / 3,153 ms
コード長 3,867 bytes
コンパイル時間 2,386 ms
コンパイル使用メモリ 210,672 KB
実行使用メモリ 4,352 KB
最終ジャッジ日時 2023-10-13 10:17:11
合計ジャッジ時間 4,783 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,352 KB
testcase_01 AC 2 ms
4,352 KB
testcase_02 AC 2 ms
4,352 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 2 ms
4,352 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 4 ms
4,352 KB
testcase_07 AC 6 ms
4,348 KB
testcase_08 AC 2 ms
4,352 KB
testcase_09 AC 2 ms
4,348 KB
testcase_10 AC 3 ms
4,348 KB
testcase_11 AC 2 ms
4,348 KB
testcase_12 AC 16 ms
4,348 KB
testcase_13 AC 193 ms
4,348 KB
testcase_14 AC 653 ms
4,352 KB
testcase_15 AC 78 ms
4,348 KB
testcase_16 AC 77 ms
4,348 KB
testcase_17 AC 81 ms
4,348 KB
testcase_18 AC 79 ms
4,352 KB
testcase_19 AC 2 ms
4,352 KB
testcase_20 AC 2 ms
4,352 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// head of mylib.h

#include <bits/stdc++.h>

using ll = long long;
using ull = unsigned long long;
using pi = std::pair<ll, ll>;

#define chmin(x, a) ((x) = std::min({(x), (a)})
#define chmax(x, a) ((x) = std::max({(x), (a)})

template <class T>
class V : public std::vector<T> {
    typedef std::vector<T> vector;

  public:
    using vector::begin;
    using vector::end;
    using vector::operator[];
    using vector::back;
    using vector::clear;
    using vector::const_iterator;
    using vector::empty;
    using vector::erase;
    using vector::front;
    using vector::insert;
    using vector::iterator;
    using vector::pop_back;
    using vector::push_back;
    using vector::rbegin;
    using vector::rend;
    using vector::resize;
    using vector::size;
    using vector::vector;

    ll sum() const {
        ll x = 0;
        std::cout << "sum()" << std::endl;
        return x;
    }
    ll max() const {
        ll x = std::numeric_limits<ll>::min();
        std::cout << "max()" << std::endl;
        return x;
    }
    ll min() const {
        ll x = std::numeric_limits<ll>::max();
        std::cout << "min()" << std::endl;
        return x;
    }
};

template <class T>
using VV = V<V<T>>;

template <class T>
std::ostream& operator<<(std::ostream& os, const V<T>& v) {
    for (auto itr = v.begin(); itr != v.end();++itr)
        os << (itr != v.begin() ? " " : "") << *itr;
    os << std::endl;
    return os;
}
template <class T>
std::istream& operator>>(std::istream& is, V<T>& v) {
    for (auto& i : v) is >> i;
    return is;
}

void fin(const char* str){
    std::cout << str << std::endl;
    exit(0);
}

// end of mylib.h

using namespace std;

#define MAX_LEN 2000
using bs = bitset<MAX_LEN>;

// O(V^3 / w^2)
V<bs> prod(const V<bs>& A, const V<bs>& B, bool counted = false) {
    int v = A.size();
    V<bs> X(v);
    for (int i = 0; i < v; ++i) {
        for (int j = A[i]._Find_first(); j < v; j = A[i]._Find_next(j)) {
            X[i] |= B[j];
        }
    }
    return X;
};

// O(V^3 log(VD) / w^2)
// V^3 log(VD) / w^2 = 1.38e+8
int doubling() {
    ll v, d;
    cin >> v >> d;
    d = min(d, v * 2ll);
    V<bs> E(v);
    // O(V^2)
    for (auto& i : E) {
        string str;
        cin >> str;
        reverse(str.begin(), str.end());
        i = bs(str);
    }
    const int pwrdp = 12; // 2^12 > 2V > 2^11
    VV<bs> dp(pwrdp, V<bs>(v));
    dp[0] = E;

    // O(V^3 log(V) / w^2)
    for (int i = 1; i < pwrdp; ++i) {
        dp[i] = prod(dp[i - 1], dp[i - 1]);
    }

    V<bs> Rchbl(v);
    // O(V)
    for (int i = 0; i < v; ++i) Rchbl[i].set(i);
    // O(V^3 log(D) / w^2)
    for (int di = 0, dd = d; dd > 0; ++di, dd /= 2) {
        if (dd & 1) Rchbl = prod(Rchbl, dp[di]);
    }
    // O(V^2 / w)
    for (int i = 0; i < v; ++i)
        if ((ll)Rchbl[i].count() < v) return cout << "No" << endl, 0;
    cout << "Yes" << endl;
    return 0;
}

// O(V^2 min(V,D) / w)
// V^2 min(V,D) / w ~ V^3 / w = 1.25e+8
int normal() {
    ll v, d;
    cin >> v >> d;
    d = min(d, v * 2);
    V<bs> E(v);

    // O(V^2 / w)
    for (auto& i : E) {
        string str;
        cin >> str;
        reverse(str.begin(), str.end());
        i = bs(str);
    }

    // O(V^2 min(V,D)/w)
    for (int i = 0; i < v; ++i) {
        bs reachable = E[i], back1, back2;
        back1[i] = true;
        // already consumed 1 turn for reachable
        for (int j = 0; j < d - 1; ++j) {
            bs diff = reachable ^ back2;  // back1 ^ back2
            if (diff.none()) break;
            back2 = back1;
            for (int k = diff._Find_first(); k < v; k = diff._Find_next(k))
                back1 |= E[k];
            swap(reachable, back1);
        }
        if ((ll)reachable.count() < v) return cout << "No" << endl, 0;
    }
    cout << "Yes" << endl;
    return 0;
}

int main(){
    normal();
    // doubling();
}
0