結果

問題 No.1078 I love Matrix Construction
コンテスト
ユーザー はまやんはまやん
提出日時 2020-06-12 23:55:00
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,631 bytes
コンパイル時間 2,750 ms
コンパイル使用メモリ 205,472 KB
最終ジャッジ日時 2025-01-11 03:16:22
ジャッジサーバーID
（参考情報） 		judge2 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 9 WA * 4 RE * 9
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<b;i++)
#define rrep(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,a) for(auto &i:a)
#define all(x) (x).begin(),(x).end()
//#pragma GCC optimize ("-O3")
using namespace std; void _main(); int main() { cin.tie(0); ios::sync_with_stdio(false); _main(); }
typedef long long ll; const int inf = INT_MAX / 2; const ll infl = 1LL << 60;
template<class T>bool chmax(T& a, const T& b) { if (a < b) { a = b; return 1; } return 0; }
template<class T>bool chmin(T& a, const T& b) { if (b < a) { a = b; return 1; } return 0; }
//---------------------------------------------------------------------------------------------------
#define ZERO(a) memset(a,0,sizeof(a))
class SCC {
public:
	static const int MV = 201010; // ここをしっかり変えること！！！
	vector<vector<int> > SC; int NV, GR[MV], rev[MV];
private:
	vector<int> E[MV], RE[MV], NUM; int vis[MV];
public:
	void init(int NV) { this->NV = NV; for (int i = 0; i<MV; i++) { E[i].clear(); RE[i].clear(); } }
	void add_edge(int x, int y) { E[x].push_back(y); RE[y].push_back(x); }
	void dfs(int cu) { vis[cu] = 1; for (int i = 0; i<E[cu].size(); i++) if (!vis[E[cu][i]]) dfs(E[cu][i]); NUM.push_back(cu); }
	void revdfs(int cu, int ind) {
		int i; vis[cu] = 1; GR[cu] = ind; SC[ind].push_back(cu); rev[cu] = ind;
		rep(i, 0, RE[cu].size()) if (!vis[RE[cu][i]]) revdfs(RE[cu][i], ind);
	}
	void scc() {
		int c = 0; SC.clear(); SC.resize(MV); NUM.clear();
		ZERO(vis); for (int i = 0; i<NV; i++) if (!vis[i]) dfs(i);
		ZERO(vis); for (int i = NUM.size() - 1; i >= 0; i--) if (!vis[NUM[i]]) {
			SC[c].clear(); revdfs(NUM[i], c); sort(SC[c].begin(), SC[c].end()); c++;
		}
		SC.resize(c);
	}
};
class TwoSat {
	int NV;
	SCC sc;
public:
	vector<int> val;
	void init(int NV) { this->NV = NV * 2; sc.init(NV * 2); val.resize(NV); }
	void add_edge(int x, int y) { // k+0:normal k+NV:inverse
		sc.add_edge((x + NV / 2) % NV, y);
		sc.add_edge((y + NV / 2) % NV, x);
	}
	bool sat() { // empty:false 
		sc.scc();
		for (int i = 0; i<NV / 2; i++) if (sc.GR[i] == sc.GR[i + NV / 2]) return false;
		for (int i = 0; i<NV / 2; i++) val[i] = sc.GR[i]>sc.GR[i + NV / 2];
		return true;
	}
};
/*---------------------------------------------------------------------------------------------------
　　　　　　　　　　　 ∧＿∧
　　　　　 ∧＿∧ 　（´<_｀ ）　 Welcome to My Coding Space!
　　　　 （ ´_ゝ`）　/　 ⌒i     @hamayanhamayan0
　　　　／　　　＼　 　  |　|
　　　 /　　 /￣￣￣￣/　　|
　 ＿_(__ﾆつ/　    ＿/ .| .|＿＿＿＿
　 　　　＼/＿＿＿＿/　（u　⊃
---------------------------------------------------------------------------------------------------*/









int N, S[501], T[501], U[501];
TwoSat ts;
//---------------------------------------------------------------------------------------------------
void _main() {
	cin >> N;
	rep(i, 0, N) cin >> S[i];
	rep(i, 0, N) cin >> T[i];
	rep(i, 0, N) cin >> U[i];

	ts.init(N * N);

	rep(i, 0, N) rep(j, 0, N) {
		int a = (S[i] - 1) * N + j;
		int b = j * N + (T[i] - 1);
		rep(x, 0, 2) rep(y, 0, 2) {
			if (x + y * 2 == U[i]) {
				int aa = a;
				if (x == 0) aa += N * N;

				int bb = b;
				if (y == 0) bb += N * N;

				ts.add_edge(aa, bb);
			}
		}
	}

	auto res = ts.sat();
	if (!res) {
		printf("-1\n");
		return;
	}

	rep(y, 0, N) {
		rep(x, 0, N) {
			if(x) printf(" ");
			if(ts.val[y * N + x]) printf("1");
			else printf("0");
		}
		printf("\n");
	}
}
0