結果

問題 No.1640 簡単な色塗り
ユーザー 🍮かんプリン🍮かんプリン
提出日時 2021-08-07 22:58:39
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 380 ms / 2,000 ms
コード長 4,046 bytes
コンパイル時間 2,120 ms
コンパイル使用メモリ 188,792 KB
実行使用メモリ 57,412 KB
最終ジャッジ日時 2024-06-29 17:01:00
合計ジャッジ時間 18,320 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 224 ms
45,564 KB
testcase_05 AC 234 ms
47,168 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,940 KB
testcase_08 AC 1 ms
6,944 KB
testcase_09 AC 1 ms
6,940 KB
testcase_10 AC 184 ms
29,632 KB
testcase_11 AC 144 ms
24,852 KB
testcase_12 AC 127 ms
22,392 KB
testcase_13 AC 318 ms
43,608 KB
testcase_14 AC 309 ms
43,372 KB
testcase_15 AC 97 ms
17,280 KB
testcase_16 AC 108 ms
20,224 KB
testcase_17 AC 253 ms
36,664 KB
testcase_18 AC 22 ms
6,940 KB
testcase_19 AC 125 ms
22,260 KB
testcase_20 AC 173 ms
28,816 KB
testcase_21 AC 144 ms
24,068 KB
testcase_22 AC 15 ms
6,944 KB
testcase_23 AC 194 ms
30,752 KB
testcase_24 AC 50 ms
11,004 KB
testcase_25 AC 124 ms
22,632 KB
testcase_26 AC 212 ms
33,656 KB
testcase_27 AC 88 ms
16,512 KB
testcase_28 AC 290 ms
41,536 KB
testcase_29 AC 220 ms
33,396 KB
testcase_30 AC 24 ms
9,728 KB
testcase_31 AC 237 ms
51,572 KB
testcase_32 AC 186 ms
43,108 KB
testcase_33 AC 107 ms
30,308 KB
testcase_34 AC 166 ms
40,748 KB
testcase_35 AC 108 ms
28,744 KB
testcase_36 AC 25 ms
9,984 KB
testcase_37 AC 35 ms
13,184 KB
testcase_38 AC 218 ms
42,116 KB
testcase_39 AC 67 ms
21,516 KB
testcase_40 AC 70 ms
19,596 KB
testcase_41 AC 150 ms
35,060 KB
testcase_42 AC 74 ms
22,204 KB
testcase_43 AC 82 ms
27,392 KB
testcase_44 AC 77 ms
24,440 KB
testcase_45 AC 57 ms
19,960 KB
testcase_46 AC 28 ms
10,368 KB
testcase_47 AC 18 ms
8,448 KB
testcase_48 AC 210 ms
48,860 KB
testcase_49 AC 8 ms
6,940 KB
testcase_50 AC 2 ms
6,940 KB
testcase_51 AC 2 ms
6,944 KB
testcase_52 AC 380 ms
54,892 KB
testcase_53 AC 369 ms
57,412 KB
07_evil_01.txt AC 940 ms
93,648 KB
07_evil_02.txt AC 1,580 ms
139,500 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

/**
 *   @FileName	a.cpp
 *   @Author	kanpurin
 *   @Created	2021.08.07 22:58:36
**/

#include "bits/stdc++.h" 
using namespace std; 
typedef long long ll;



struct StronglyConnectedComponents {
private:
    vector<vector<int>> g, rg;
    vector<int> order; 
    vector<bool> visited;
    int V;
    void dfs(int v) {
        if (visited[v]) return;
        visited[v] = true;
        for (int u : g[v]) dfs(u);
        order.push_back(v);
    }
    void rdfs(int v, int c) {
        if (comp[v] != -1) return;
        comp[v] = c;
        for (int u : rg[v]) rdfs(u, c);
    }
public:
    vector<vector<int>> Graph; 
    vector<int> comp; 
    StronglyConnectedComponents(int v) : V(v) {
        g.resize(v);
        rg.resize(v);
        visited.resize(v);
        comp.resize(v);
    }
    
    void add_vertex(int n = 1) {
        V += n;
        g.resize(V);
        rg.resize(V);
        visited.resize(V);
        comp.resize(V);
    }
    void add_edge(int from, int to) {
        g[from].push_back(to);
        rg[to].push_back(from);
    }
    void build() {
        Graph = vector<vector<int>>();
        order = vector<int>();
        visited.assign(V, false);
        comp.assign(V, -1);
        for (int i = 0; i < V; i++) dfs(i);
        reverse(order.begin(), order.end());
        int number = 0;
        for (int i = 0; i < V; i++) {
            if (comp[order[i]] == -1) {
                rdfs(order[i], number++);
            }
        }
        Graph.resize(number);
        for (int i = 0; i < V; i++) {
            for (int v : g[i]) {
                if (comp[i] == comp[v]) continue;
                Graph[comp[i]].push_back(comp[v]);
            }
        }
    }
};

struct TwoSAT {
private:
    int _n;
    std::vector<bool> _answer;
    StronglyConnectedComponents scc;
    void only_one_naive(const vector<int> &x, const vector<bool> &t) {
        assert(x.size() == t.size());
        int _m = x.size();
        for (int i = 0; i < _m-1; i++) {
            for (int j = i+1; j < _m; j++) {
                add_clause(x[i],!t[i],x[j],!t[j]);
            }
        }
    }
public:
    TwoSAT():_n(0),scc(0){}
    explicit TwoSAT(int n) : _n(n), _answer(n), scc(2*n) {}
    void add_clause(int i, bool f, int j, bool g) {
        assert(0 <= i && i < _n);
        assert(0 <= j && j < _n);
        scc.add_edge(2 * i + (f ? 0 : 1), 2 * j + (g ? 1 : 0));
        scc.add_edge(2 * j + (g ? 0 : 1), 2 * i + (f ? 1 : 0));
    }
    
    void only_one(const vector<int> &x, const vector<bool> &t) {
        assert(x.size() == t.size());
        int _m = x.size();
        if (_m <= 12) {
            only_one_naive(x,t);
            return;
        }
        int _nn = _n;
        
        _n += _m;
        
        scc.add_vertex(_m*2);
        _answer.resize(_n);
        for (int i = 0; i < _m; i++) {
            add_clause(_nn+i,false,x[i],!t[i]);
            
        }
        for (int i = 1; i < _m; i++) {
            add_clause(_nn+i,false,_nn+i-1,true);
            add_clause(x[i],!t[i],_nn+i-1,true);
            
            
        }
    }
    bool satisfiable() {
        scc.build();
        auto id = scc.comp;
        for (int i = 0; i < _n; i++) {
            if (id[2 * i] == id[2 * i + 1]) return false;
            _answer[i] = id[2 * i] < id[2 * i + 1];
        }
        return true;
    }
    std::vector<bool> answer() { return _answer; }
};
int main() {
    int n;cin >> n;
    TwoSAT ts(n);
    vector<int> a(n),b(n);
    vector<vector<int>> x(n);
    vector<vector<bool>> t(n);
    for (int i = 0; i < n; i++) {
        cin >> a[i] >> b[i];
        x[a[i]-1].push_back(i);
        x[b[i]-1].push_back(i);
        t[a[i]-1].push_back(true);
        t[b[i]-1].push_back(false);
    }
    for (int i = 0; i < n; i++) {
        ts.only_one(x[i],t[i]);
    }
    if (!ts.satisfiable()) {
        puts("No");
    }
    else {
        puts("Yes");
        auto ans = ts.answer();
        for (int i = 0; i < n; i++) {
            cout << (ans[i]?a[i]:b[i]) << endl;
        }
    }
    return 0;
}
0