結果

問題 No.1078 I love Matrix Construction
コンテスト
ユーザー tamato
提出日時 2020-06-12 23:06:17
言語 PyPy3
(7.3.15)
結果
RE  
実行時間 -
コード長 3,020 bytes
コンパイル時間 209 ms
コンパイル使用メモリ 82,848 KB
実行使用メモリ 237,508 KB
最終ジャッジ日時 2024-06-24 05:49:40
合計ジャッジ時間 13,375 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 RE * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod = 1000000007
eps = 10**-9


def main():
    import sys
    from collections import deque
    input = sys.stdin.readline

    def fin():
        print(-1)
        exit()

    N = int(input())
    S = list(map(int, input().split()))
    T = list(map(int, input().split()))
    U = list(map(int, input().split()))

    def idx(i, j, flg):
        return i * N + j + 1 + flg * N*N

    adj = [[] for _ in range(N*N*2 + 1)]
    adj_rev = [[] for _ in range(N * N * 2 + 1)]
    for i in range(N):
        s = S[i] - 1
        t = T[i] - 1
        u = U[i]
        if u == 0:
            flg_s = 1
            flg_t = 1
        elif u == 1:
            flg_s = 0
            flg_t = 1
        elif u == 2:
            flg_s = 1
            flg_t = 0
        else:
            flg_s = 0
            flg_t = 0
        for j in range(N):
            a = idx(s, j, flg_s)
            b = idx(j, t, flg_t)
            if a > N*N:
                adj[a - N*N].append(b)
                adj_rev[b].append(a - N*N)
            else:
                adj[a + N*N].append(b)
                adj_rev[b].append(a + N*N)
            if b > N*N:
                adj[b - N*N].append(a)
                adj_rev[a].append(b - N*N)
            else:
                adj[b + N*N].append(a)
                adj_rev[a].append(b + N*N)

    # SCC
    seen = [0] * (N*N*2 + 1)
    order = []
    for v0 in range(1, N*N*2 + 1):
        if seen[v0]:
            continue
        st = deque()
        st.append(v0)
        seen[v0] = 1
        while st:
            v = st.pop()
            if v < 0:
                order.append(-v)
            else:
                st.append(-v)
                for u in adj[v]:
                    if not seen[u]:
                        seen[u] = 1
                        st.append(u)

    order.reverse()
    seen2 = [0] * (N*N*2 + 1)
    G = []
    G_idx = [-1] * (N*N*2 + 1)
    for v0 in order:
        if seen2[v0]:
            continue
        st = deque()
        st.append(v0)
        seen2[v0] = 1
        G.append([v0])
        G_idx[v0] = len(G) - 1
        while st:
            v = st.pop()
            for u in adj_rev[v]:
                if seen2[u]:
                    continue
                G[-1].append(u)
                G_idx[u] = len(G) - 1
                seen2[u] = 1
                st.append(u)
    """
    for k in range(len(G)):
        for v in G[k]:
            if v <= N*N:
                v_ = v + N*N
                if G_idx[v] == G_idx[v_]:
                    fin()
    """

    ans = [[0] * N for _ in range(N)]
    for i in range(N):
        for j in range(N):
            v = idx(i, j, 0)
            v_ = idx(i, j, 1)
            if G_idx[v] < G_idx[v_]:
                ans[i][j] = 1
    for i in range(N):
        print(*ans[i])

    # check
    for i in range(N):
        s = S[i] - 1
        t = T[i] - 1
        u = U[i]
        for j in range(N):
            if u == ans[s][j] + ans[j][t] * 2:
                assert False


if __name__ == '__main__':
    main()
0