結果

問題 No.1078 I love Matrix Construction
ユーザー beetbeet
提出日時 2020-06-12 21:40:46
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 344 ms / 2,000 ms
コード長 2,973 bytes
コンパイル時間 2,541 ms
コンパイル使用メモリ 218,708 KB
実行使用メモリ 95,456 KB
最終ジャッジ日時 2024-06-24 04:38:24
合計ジャッジ時間 8,071 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 42 ms
18,044 KB
testcase_03 AC 122 ms
38,940 KB
testcase_04 AC 167 ms
51,316 KB
testcase_05 AC 154 ms
43,844 KB
testcase_06 AC 35 ms
16,144 KB
testcase_07 AC 13 ms
8,280 KB
testcase_08 AC 138 ms
43,704 KB
testcase_09 AC 7 ms
6,944 KB
testcase_10 AC 344 ms
95,456 KB
testcase_11 AC 182 ms
54,052 KB
testcase_12 AC 292 ms
79,524 KB
testcase_13 AC 328 ms
89,360 KB
testcase_14 AC 229 ms
63,644 KB
testcase_15 AC 311 ms
84,404 KB
testcase_16 AC 12 ms
7,148 KB
testcase_17 AC 2 ms
6,940 KB
testcase_18 AC 25 ms
13,120 KB
testcase_19 AC 72 ms
26,064 KB
testcase_20 AC 67 ms
25,048 KB
testcase_21 AC 4 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}
template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}
using Int = long long;
const char newl = '\n';


struct SCC{
  vector< vector<int> > G,R,T,C;
  vector<int> vs,used,blg;
  SCC(){}
  SCC(int n):G(n),R(n),used(n),blg(n){}

  void add_edge(int u,int v){
    G[u].emplace_back(v);
    R[v].emplace_back(u);
  }

  void dfs(int v){
    used[v]=1;
    for(int u:G[v])
      if(!used[u]) dfs(u);
    vs.emplace_back(v);
  }

  void rdfs(int v,int k){
    used[v]=1;
    blg[v]=k;
    C[k].emplace_back(v);
    for(int u:R[v])
      if(!used[u]) rdfs(u,k);
  }

  int build(){
    int n=G.size();
    for(int v=0;v<n;v++)
      if(!used[v]) dfs(v);

    fill(used.begin(),used.end(),0);
    int k=0;
    for(int i=n-1;i>=0;i--){
      if(!used[vs[i]]){
        T.emplace_back();
        C.emplace_back();
        rdfs(vs[i],k++);
      }
    }

    for(int v=0;v<n;v++)
      for(int u:G[v])
        if(blg[v]!=blg[u])
          T[blg[v]].push_back(blg[u]);

    for(int i=0;i<k;i++){
      sort(T[i].begin(),T[i].end());
      T[i].erase(unique(T[i].begin(),T[i].end()),T[i].end());
    }
    return k;
  }

  int operator[](int k) const{return blg[k];}
};


struct TwoSat{
  int n;
  SCC scc;
  TwoSat(int n):n(n),scc(n*2){}
  int negate(int v){return (n+v)%(n*2);}
  void add_if(int u,int v){
    // u -> v <=> !v -> !u
    scc.add_edge(u,v);
    scc.add_edge(negate(v),negate(u));
  }
  void add_or(int u,int v){
    // u or v <=> !u -> v
    add_if(negate(u),v);
  }
  void add_nand(int u,int v){
    // u nand v <=> u -> !v
    add_if(u,negate(v));
  }
  void set_true(int v){
    //  v <=> !v ->  v
    scc.add_edge(negate(v),v);
  }
  void set_false(int v){
    // !v <=>  v -> !v
    scc.add_edge(v,negate(v));
  }
  vector<int> build(){
    scc.build();
    vector<int> res(n);
    for(int i=0;i<n;i++){
      if(scc[i]==scc[n+i]) return {};
      res[i]=scc[i]>scc[n+i];
    }
    return res;
  }
};


template<typename T> void drop(const T &x){cout<<x<<endl;exit(0);}

//INSERT ABOVE HERE
signed main(){
  cin.tie(0);
  ios::sync_with_stdio(0);

  int n;
  cin>>n;
  vector<int> ss(n),ts(n),us(n);
  for(int i=0;i<n;i++) cin>>ss[i],ss[i]--;
  for(int i=0;i<n;i++) cin>>ts[i],ts[i]--;
  for(int i=0;i<n;i++) cin>>us[i];

  auto idx=[&](int y,int x){return y*n+x;};
  TwoSat G(n*n);
  for(int i=0;i<n;i++){
    for(int j=0;j<n;j++){
      int u=idx(ss[i],j);
      int v=idx(j,ts[i]);
      for(int a=0;a<2;a++){
        for(int b=0;b<2;b++){
          if(a+b*2==us[i]){
            if(!a) u=G.negate(u);
            if(!b) v=G.negate(v);
            G.add_or(G.negate(u),G.negate(v));
          }
        }
      }
    }
  }

  auto res=G.build();
  if(res.empty()) drop(-1);

  for(int i=0;i<n;i++){
    for(int j=0;j<n;j++){
      if(j) cout<<' ';
      cout<<res[idx(i,j)];
    }
    cout<<newl;
  }
  return 0;
}
0