#173454 (C++14) No.274 The Wall

• 問題: No.274 The Wall
• ユーザー: oyas
• 提出日時: 2017-05-14 09:50:42
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 425 ms / 2,000 ms
• コード長: 3,471 bytes
• コンパイル時間: 1,059 ms
• コンパイル使用メモリ: 82,892 KB
• 実行使用メモリ: 132,128 KB
• 最終ジャッジ日時: 2023-09-04 02:19:23
• 合計ジャッジ時間: 3,063 ms
• ジャッジサーバーID （参考情報）: judge15 / judge14

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 4,380 KB
testcase_01	AC	2 ms; 4,380 KB
testcase_02	AC	2 ms; 4,380 KB
testcase_03	AC	134 ms; 69,948 KB
testcase_04	AC	2 ms; 4,376 KB
testcase_05	AC	2 ms; 4,376 KB
testcase_06	AC	2 ms; 4,376 KB
testcase_07	AC	2 ms; 4,380 KB
testcase_08	AC	2 ms; 4,376 KB
testcase_09	AC	2 ms; 4,380 KB
testcase_10	AC	2 ms; 4,376 KB
testcase_11	AC	425 ms; 132,128 KB
testcase_12	AC	10 ms; 4,380 KB
testcase_13	AC	2 ms; 4,380 KB
testcase_14	AC	5 ms; 4,380 KB
testcase_15	AC	10 ms; 4,380 KB
testcase_16	AC	96 ms; 38,656 KB
testcase_17	AC	91 ms; 36,980 KB
testcase_18	AC	99 ms; 39,972 KB
testcase_19	AC	17 ms; 4,376 KB
testcase_20	AC	19 ms; 4,384 KB
testcase_21	AC	20 ms; 4,376 KB
testcase_22	AC	20 ms; 4,380 KB
testcase_23	AC	21 ms; 4,380 KB
testcase_24	AC	20 ms; 4,376 KB
testcase_25	AC	21 ms; 4,376 KB

ソースコード

// 2sat.cpp
#include <vector>
#include <iostream>

#define rep(i,n) for(int i=0; i<(n); i++)

using namespace std;

typedef vector<int> vi;

// scc.cppのコピー
// 強連結成分分解 (Strongly Connected Component)
struct SCCGraph: public vector<vi> {

	vector<int> order;    // 属する強連結成分の番号(トポロジカル順序になっている)
	vector<int> vs;   // 帰りがけ順
	vector<bool> used;    // すでに調べたか

	using vector::vector;	// 継承コンストラクタ
	SCCGraph(const vector<vi> &v): vector::vector(v) {}	// コピーコンストラクタを使用可能にする

	// a -> b の辺を追加する。push_back()でも同じ
	void add_edge(int a, int b){
		(*this)[a].push_back(b);
	}

	// scc()で使う
	void dfs(int n, int k, vector<vi> &v, vector<vi> &rv){
		used[n] = true;
		//order[n] = k;
		for(auto t: v[n]){
			rv[t].push_back(n);	// 逆辺のグラフを作成
			if( !used[t] ) dfs(t, k, v, rv);
		}
		vs.push_back(n);
	}
	void rdfs(int n, int k, vector<vi> &v){
		used[n] = true;
		order[n] = k;
		for(auto t: v[n]) if( !used[t] ) rdfs(t, k, v);
	}

	// 強連結成分分解を行う
	int scc(){
		int N = size();
		used.assign(N, false);
		order.resize(N);
		vs.clear();
		vs.reserve(2*N);
		vector<vi> rG(N), tmp(N);   // 辺を逆にしたグラフ用
		for(int n=0; n < N; n++){
			if( !used[n] ) dfs(n, n, (*this), rG);
		}
		used.assign(N, false);
		int k = 0;		// 強連結成分の番号
		for(int i = vs.size()-1; i >= 0; i--){
			if( !used[ vs[i] ] ) rdfs(vs[i], k++, rG);
		}
		return k;
	}

	// 属する強連結成分が同じかどうか判定
	bool find(int x, int y){
		return order[x] == order[y];
	}

};
// 関数としてsccを使う。戻り値は、各ノードが所属する強連結成分の番号
vector<int> scc(const vector<vi> &v){
	SCCGraph g(v);
	g.scc();
	return g.order;
}


// 2-SAT
struct TwoSAT {
	SCCGraph graph;
	int N;

	TwoSAT(int n): graph( (n+1)*2 ) {
		N = n;
	}

	// 辺( a -> b )を追加する。負の値は論理変数の否定(~a)を表す。
	void add_edge(int a, int b){
		if( a < 0 ) a = N - a;
		if( b < 0 ) b = N - b;
		graph.add_edge(a, b);
	}

	// 論理式( a V b ) を追加する。a,bは、論理変数の番号(0以外の整数)。否定は負の値で表す。
	void add_or(int a, int b){
		add_edge(-a, b);
		add_edge(-b, a);
	}

	// 判定を行う
	bool solve(){
		graph.scc();
		bool ret = true;
		for(int i=1; i <= N; i++){
			if( graph.order[i] == graph.order[N+i] ) ret = false;
		}
		return ret;
	}

	// 変数の値を取得する。
	bool value(int a){
		if( a < 0 ) a = N - a;
		return graph.order[a] > graph.order[(N+a-1)%(2*N)+1];
	}
};


// test
// yukicoder No.274 The Wall (http://yukicoder.me/problems/no/274)
int main(){
	int N, M;
	int L[2005], R[2005];
	int rL[2005], rR[2005];
	TwoSAT sat(2005);

	cin >> N >> M;
	rep(i,N){
		cin >> L[i] >> R[i];
		// 180度回転
		rR[i] = M - L[i] - 1;
		rL[i] = M - R[i] - 1;
	}
	auto fn = [&](int Li, int Lj, int Ri, int Rj, int a, int b){
		if( !(Rj < Li || Ri < Lj) ){
			// 重なる可能性があるものはadd
			sat.add_or( -a, -b );
		}
	};
	rep(i,N) rep(j,i){
		fn(  L[i],  L[j],  R[i],  R[j],  (i+1),  (j+1) );
		fn( rL[i],  L[j], rR[i],  R[j], -(i+1),  (j+1) );
		fn(  L[i], rL[j],  R[i], rR[j],  (i+1), -(j+1) );
		fn( rL[i], rL[j], rR[i], rR[j], -(i+1), -(j+1) );
	}

	cout << (sat.solve() ? "YES" : "NO") << endl;
	//rep(i,N) cout << sat.value(i+1) << endl;

	return 0;
}