結果

問題 No.1078 I love Matrix Construction
ユーザー knshnb
提出日時 2020-06-12 22:42:09
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 361 ms / 2,000 ms
コード長 3,320 bytes
コンパイル時間 3,246 ms
コンパイル使用メモリ 206,972 KB
最終ジャッジ日時 2025-01-11 02:46:44
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>  // clang-format off
using Int = long long;
#define REP_(i, a_, b_, a, b, ...) for (Int i = (a), lim##i = (b); i < lim##i; i++)
#define REP(i, ...) REP_(i, __VA_ARGS__, __VA_ARGS__, 0, __VA_ARGS__)
struct SetupIO { SetupIO() { std::cin.tie(nullptr), std::ios::sync_with_stdio(false), std::cout << std::fixed << std::setprecision(13); } } setup_io;
#ifndef dump
#define dump(...)
#endif  // clang-format on
/**
 *    author:  knshnb
 *    created: Fri Jun 12 22:02:09 JST 2020
 **/
/// @docs src/Graph/StronglyConnectedComponents.md
struct StronglyConnectedComponents {
    int n, size;  // sizeはbuild()後に強連結成分の数を格納
    std::vector<std::vector<int>> edge, redge;
    std::vector<int> belong_to;
    StronglyConnectedComponents(int n_) : n(n_), edge(n_), redge(n_) {}
    void add_edge(int u, int v) {
        assert(0 <= u && u < n && 0 <= v && v < n);
        edge[u].push_back(v);
        redge[v].push_back(u);
    }
    void build() {
        std::vector<int> ord;  // post-order
        // 正方向の辺でdfs、post-orderをふる
        std::vector<bool> visited(n, false);
        auto dfs1 = [&](auto f, int v) -> void {
            if (visited[v]) return;
            visited[v] = true;
            for (int s : edge[v]) {
                f(f, s);
            }
            ord.push_back(v);
        };
        for (int i = 0; i < n; i++) dfs1(dfs1, i);
        // post-order逆順にdfs、到達可能な同値類に分ける
        belong_to.assign(n, -1);
        int cur_group = 0;
        auto dfs2 = [&](auto f, int v) -> void {
            if (belong_to[v] != -1) return;
            belong_to[v] = cur_group;
            for (int s : redge[v]) {
                f(f, s);
            }
        };
        std::reverse(ord.begin(), ord.end());
        for (int i : ord) {
            if (belong_to[i] == -1) {
                dfs2(dfs2, i);
                cur_group++;
            }
        }
        size = cur_group;
    }
    bool has_loop() {
        build();
        return size < n;
    }
};
template <class T, class S> T make_vec(S x) { return x; }
template <class T, class... Ts> auto make_vec(size_t n, Ts... ts) {
    return std::vector<decltype(make_vec<T>(ts...))>(n, make_vec<T>(ts...));
}
signed main() {
    Int n;
    std::cin >> n;
    std::vector<Int> s(n), t(n), u(n);
    REP(i, n) std::cin >> s[i], s[i]--;
    REP(i, n) std::cin >> t[i], t[i]--;
    REP(i, n) std::cin >> u[i];
    StronglyConnectedComponents g(2 * n * n);
    auto enc = [&](Int i, Int j) { return n * i + j; };
    auto rev = [&](Int x) { return x < n * n ? x + n * n : x - n * n; };
    REP(i, n) {
        REP(j, n) {
            Int x = enc(s[i], j);
            if (u[i] & 1) x = rev(x);
            Int y = enc(j, t[i]);
            if (u[i] >> 1) y = rev(y);
            g.add_edge(rev(x), y);
            g.add_edge(rev(y), x);
        }
    }
    g.build();
    auto ans = make_vec<Int>(n, n, 0);
    REP(i, n) REP(j, n) {
        Int x = enc(i, j);
        if (g.belong_to[x] == g.belong_to[rev(x)]) {
            std::cout << -1 << std::endl;
            return 0;
        }
        ans[i][j] = g.belong_to[x] > g.belong_to[rev(x)];
    }
    REP(i, n) {
        REP(j, n) { std::cout << ans[i][j] << " "; }
        std::cout << "\n";
    }
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0