結果

問題 No.1078 I love Matrix Construction
ユーザー yukudoyukudo
提出日時 2020-06-13 00:05:17
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 232 ms / 2,000 ms
コード長 2,450 bytes
コンパイル時間 2,158 ms
コンパイル使用メモリ 175,376 KB
実行使用メモリ 51,084 KB
最終ジャッジ日時 2024-06-24 06:07:21
合計ジャッジ時間 6,456 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 15 ms
29,940 KB
testcase_01 AC 14 ms
29,844 KB
testcase_02 AC 31 ms
33,060 KB
testcase_03 AC 76 ms
37,944 KB
testcase_04 AC 109 ms
40,840 KB
testcase_05 AC 98 ms
39,168 KB
testcase_06 AC 30 ms
32,960 KB
testcase_07 AC 20 ms
31,248 KB
testcase_08 AC 99 ms
39,448 KB
testcase_09 AC 17 ms
30,544 KB
testcase_10 AC 232 ms
51,084 KB
testcase_11 AC 117 ms
41,492 KB
testcase_12 AC 182 ms
47,372 KB
testcase_13 AC 212 ms
49,472 KB
testcase_14 AC 143 ms
43,928 KB
testcase_15 AC 200 ms
48,712 KB
testcase_16 AC 20 ms
30,828 KB
testcase_17 AC 14 ms
29,916 KB
testcase_18 AC 26 ms
32,060 KB
testcase_19 AC 63 ms
35,004 KB
testcase_20 AC 47 ms
34,848 KB
testcase_21 AC 15 ms
30,028 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
#define REP(i,n) for(int i=0,_n=(int)(n);i<_n;++i)
#define ALL(v) (v).begin(),(v).end()
#define CLR(t,v) memset(t,(v),sizeof(t))
template<class T1,class T2>ostream& operator<<(ostream& os,const pair<T1,T2>&a){return os<<"("<<a.first<<","<<a.second<< ")";}
template<class T>void pv(T a,T b){for(T i=a;i!=b;++i)cout<<(*i)<<" ";cout<<endl;}
template<class T>void chmin(T&a,const T&b){if(a>b)a=b;}
template<class T>void chmax(T&a,const T&b){if(a<b)a=b;}

ll nextLong() { ll x; scanf("%lld", &x); return x;}

const int MAX_V = 500 * 500 * 2 + 10;
vector<int> g[MAX_V];
vector<int> r[MAX_V];


const int MAX_N = 505;
int A[MAX_N][MAX_N];
int N;
int nodeId(int r, int c, int f) {
  return (r * N + c) * 2 + f;
}

int sccid[MAX_V];
int sccnum;
void dfs(int p, vector<int>& order) {
  if (sccid[p] != 0) return;
  sccid[p] = 1;
  for (int i : g[p]) dfs(i, order);
  order.push_back(p);
}
void rdfs(int p, int id) {
  if (sccid[p] != -1) return;
  sccid[p] = id;
  for (int i : r[p]) rdfs(i, id);
}
void scc() {
  vector<int> order;
  CLR(sccid, 0);
  REP(i, 2*N*N) dfs(i, order);
  CLR(sccid, -1);
  sccnum = 0;
  for (int i = (int)order.size() - 1; i >= 0; i--) if (sccid[order[i]] == -1) {
    rdfs(order[i], sccnum);
    sccnum++;
  }
}


int main2() {
  REP(i, MAX_V) g[i].clear();
  REP(i, MAX_V) r[i].clear();

  N = nextLong();

  vector<int> S(N), T(N), U(N);
  REP(i, N) S[i] = nextLong()-1;
  REP(i, N) T[i] = nextLong()-1;
  REP(i, N) U[i] = nextLong();


  REP(i, N) REP(j, N) {
    int a = nodeId(S[i], j, U[i] & 1);
    int b = nodeId(j, T[i], U[i] >> 1);
    g[b].push_back(a^1);
    g[a].push_back(b^1);
    r[a^1].push_back(b);
    r[b^1].push_back(a);
  }

  scc();

  // REP(i, N*N) cout << sccid[i * 2] << " "; cout << endl;
  // REP(i, N*N) cout << sccid[i * 2 + 1] << " "; cout << endl;

  bool ng = false;
  REP(r, N) REP(c, N) {
    if (sccid[nodeId(r, c, 0)] != -1 && sccid[nodeId(r, c, 0)] == sccid[nodeId(r, c, 1)]) {
      ng = true;
    }
  }
  if (ng) {
    cout << -1 << endl;
  } else {

    REP(i, N*N) {
      int r = i / N;
      int c = i % N;
      A[r][c] = (sccid[nodeId(r, c, 1)] > sccid[nodeId(r, c, 0)]) ? 1 : 0;
    }

    REP(r, N) {
      REP(c, N) {
        if (c) cout << ' ';
        cout << A[r][c];
      }
      cout << '\n';
    }
  }
  return 0;
}

int main() {

#ifdef LOCAL
  for (;!cin.eof();cin>>ws)
#endif
    main2();
  return 0;
}
0