結果

問題 No.1078 I love Matrix Construction
ユーザー はまやんはまやんはまやんはまやん
提出日時 2020-06-13 00:00:50
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 344 ms / 2,000 ms
コード長 3,631 bytes
コンパイル時間 2,316 ms
コンパイル使用メモリ 214,092 KB
実行使用メモリ 79,940 KB
最終ジャッジ日時 2024-06-24 06:07:05
合計ジャッジ時間 7,930 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 19 ms
43,308 KB
testcase_01 AC 19 ms
44,144 KB
testcase_02 AC 46 ms
48,456 KB
testcase_03 AC 113 ms
56,528 KB
testcase_04 AC 162 ms
63,104 KB
testcase_05 AC 130 ms
58,916 KB
testcase_06 AC 46 ms
47,832 KB
testcase_07 AC 29 ms
44,948 KB
testcase_08 AC 141 ms
59,944 KB
testcase_09 AC 23 ms
44,676 KB
testcase_10 AC 344 ms
79,940 KB
testcase_11 AC 162 ms
64,252 KB
testcase_12 AC 280 ms
73,600 KB
testcase_13 AC 320 ms
77,200 KB
testcase_14 AC 215 ms
67,148 KB
testcase_15 AC 307 ms
75,784 KB
testcase_16 AC 26 ms
44,896 KB
testcase_17 AC 21 ms
42,412 KB
testcase_18 AC 39 ms
47,548 KB
testcase_19 AC 75 ms
51,592 KB
testcase_20 AC 71 ms
51,624 KB
testcase_21 AC 20 ms
42,692 KB
権限があれば一括ダウンロードができます

ソースコード

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 = 501010; // ここをしっかり変えること！！！
	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 == 1) aa += N * N;

				int bb = b;
				if (y == 1) 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