結果

問題 No.1078 I love Matrix Construction
コンテスト
ユーザー ningenMe
提出日時 2020-06-12 22:37:05
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 5,789 bytes
コンパイル時間 3,744 ms
コンパイル使用メモリ 209,140 KB
最終ジャッジ日時 2025-01-11 02:45:10
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 17 WA * 5
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#define ALL(obj) (obj).begin(),(obj).end()
#define SPEED cin.tie(0);ios::sync_with_stdio(false);

template<class T> using PQ = priority_queue<T>;
template<class T> using PQR = priority_queue<T,vector<T>,greater<T>>;

constexpr long long MOD = (long long)1e9 + 7;
constexpr long long MOD2 = 998244353;
constexpr long long HIGHINF = (long long)1e18;
constexpr long long LOWINF = (long long)1e15;
constexpr long double PI = 3.1415926535897932384626433L;

template <class T> vector<T> multivector(size_t N,T init){return vector<T>(N,init);}
template <class... T> auto multivector(size_t N,T... t){return vector<decltype(multivector(t...))>(N,multivector(t...));}
template <class T> void corner(bool flg, T hoge) {if (flg) {cout << hoge << endl; exit(0);}}
template <class T, class U>ostream &operator<<(ostream &o, const map<T, U>&obj) {o << "{"; for (auto &x : obj) o << " {" << x.first << " : " << x.second << "}" << ","; o << " }"; return o;}
template <class T>ostream &operator<<(ostream &o, const set<T>&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;}
template <class T>ostream &operator<<(ostream &o, const multiset<T>&obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr) o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;}
template <class T>ostream &operator<<(ostream &o, const vector<T>&obj) {o << "{"; for (int i = 0; i < (int)obj.size(); ++i)o << (i > 0 ? ", " : "") << obj[i]; o << "}"; return o;}
template <class T, class U>ostream &operator<<(ostream &o, const pair<T, U>&obj) {o << "{" << obj.first << ", " << obj.second << "}"; return o;}
template <template <class tmp>  class T, class U> ostream &operator<<(ostream &o, const T<U> &obj) {o << "{"; for (auto itr = obj.begin(); itr != obj.end(); ++itr)o << (itr != obj.begin() ? ", " : "") << *itr; o << "}"; return o;}
void print(void) {cout << endl;}
template <class Head> void print(Head&& head) {cout << head;print();}
template <class Head, class... Tail> void print(Head&& head, Tail&&... tail) {cout << head << " ";print(forward<Tail>(tail)...);}
template <class T> void chmax(T& a, const T b){a=max(a,b);}
template <class T> void chmin(T& a, const T b){a=min(a,b);}
void YN(bool flg) {cout << (flg ? "YES" : "NO") << endl;}
void Yn(bool flg) {cout << (flg ? "Yes" : "No") << endl;}
void yn(bool flg) {cout << (flg ? "yes" : "no") << endl;}

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

class StronglyConnectedComponents{
	int nodeNum;
	vector<vector<int>> edge,revEdge;

	vector<int> label,visited,order;

	void dfs(int curr){
		visited[curr] = 1;
		for(int next:edge[curr]) if(!visited[next]) dfs(next);
		order.push_back(curr);
	}

	void revDfs(int curr,int labelId){
		visited[curr] = 1;
		label[curr] = labelId;
		for(int next:revEdge[curr]) if(!visited[next]) revDfs(next,labelId);
	}

public:

	StronglyConnectedComponents(const int& nodeNum) : 
	nodeNum(nodeNum), edge(nodeNum), revEdge(nodeNum), label(nodeNum), visited(nodeNum) {
		// do nothing        
	}

	int operator[](int idx) {
		return label[idx];
	}

	void makeEdge(const int from,const int to) {
		edge[from].push_back(to);
		revEdge[to].push_back(from);
	}

	void solve(void) {
		for(int i = 0; i < nodeNum; ++i) visited[i] = 0;
		for(int i = 0; i < nodeNum; ++i) if(!visited[i]) dfs(i);
		for(int i = 0; i < nodeNum; ++i) visited[i] = 0;
		reverse(order.begin(),order.end());
		int labelId = 0;
		for(int i:order) if(!visited[i]) revDfs(i,labelId++);
	}

	void print(void) {
		for(auto labelId:label) cout << labelId << " ";
		cout << endl;
	}

};

int main(void){
	int N; cin >> N;
	//[0,N*N)       i番目は0
	//[N*N,2*N*N)   i番目は1
	StronglyConnectedComponents scc(2*N*N);
    vector<int> S(N),T(N),U(N);
	for(int i = 0; i < N; ++i) cin >> S[i],S[i]--;
	for(int i = 0; i < N; ++i) cin >> T[i],T[i]--;
	for(int i = 0; i < N; ++i) cin >> U[i];
    for(int i = 0; i < N; ++i) {
        for(int j = 0; j < N; ++j) {
            int a=S[i],b=j,c=j,d=T[i];
            int s=a*N+b,t=c*N+d;
            if(U[i]==0) {
                //ab=0&&cd=0がだめ
                scc.makeEdge(s,t+N*N);
                scc.makeEdge(t,s+N*N);
            }
            if(U[i]==1) {
                //ab=1&&cd=0がだめ
                scc.makeEdge(s+N*N,t+N*N);
                scc.makeEdge(t,s);
            }
            if(U[i]==2) {
                //ab=0&&cd=1がだめ
                scc.makeEdge(s,t);
                scc.makeEdge(t+N*N,s+N*N);
            }
            if(U[i]==3) {
                //ab=1&&cd=1がだめ
                scc.makeEdge(s+N*N,t);
                scc.makeEdge(t+N*N,s);
            }
        }
    }
    scc.solve();
	int flg = 1;
	for(int i = 0; i < N*N; ++i) if(scc[i]==scc[i+N*N]) flg = 0;
    corner(!flg,-1);

    vector<vector<int>> v(2*N*N);
    for(int i = 0; i < 2*N*N; ++i) {
        v[scc[i]].push_back(i);
    }

    vector<int> bl(2*N*N,-1);
    for(int i = 0; i < 2*N*N; ++i) {
        if(bl[scc[i]]!=-1) continue;
        queue<pair<int,int>> q;
        q.push({scc[i],1});
        bl[scc[i]]=1;
        while(q.size()){
            auto p = q.front();
            q.pop();
            int c = p.first;
            for(int j:v[c]) {
                int k;
                if(j<N*N ) k=j+N*N;
                if(N*N<=j) k=j-N*N;
                if(bl[scc[k]]!=-1) continue;
                bl[scc[k]]=p.second^1;
                q.push({scc[k],p.second^1});
            }
        }
    }
    for(int i = 0; i < N; ++i) {
        for(int j = 0; j < N; ++j) {
            int k = i*N+j;
            cout << bl[scc[k]] << " ";
        }
        cout << endl;
    }
	return 0;
}
0