結果

問題 No.1698 Face to Face
ユーザー eSeFeSeF
提出日時 2021-09-11 18:21:22
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 429 ms / 5,000 ms
コード長 5,908 bytes
コンパイル時間 5,384 ms
コンパイル使用メモリ 270,700 KB
実行使用メモリ 76,388 KB
最終ジャッジ日時 2023-09-26 15:45:00
合計ジャッジ時間 17,169 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 1 ms
4,388 KB
testcase_03 AC 232 ms
76,100 KB
testcase_04 AC 255 ms
76,176 KB
testcase_05 AC 354 ms
76,276 KB
testcase_06 AC 356 ms
76,036 KB
testcase_07 AC 360 ms
75,564 KB
testcase_08 AC 364 ms
75,776 KB
testcase_09 AC 363 ms
75,840 KB
testcase_10 AC 284 ms
76,020 KB
testcase_11 AC 380 ms
76,188 KB
testcase_12 AC 401 ms
76,324 KB
testcase_13 AC 2 ms
4,380 KB
testcase_14 AC 1 ms
4,380 KB
testcase_15 AC 429 ms
76,376 KB
testcase_16 AC 426 ms
76,388 KB
testcase_17 AC 425 ms
76,280 KB
testcase_18 AC 364 ms
67,276 KB
testcase_19 AC 349 ms
66,424 KB
testcase_20 AC 344 ms
63,500 KB
testcase_21 AC 170 ms
35,624 KB
testcase_22 AC 125 ms
27,844 KB
testcase_23 AC 29 ms
9,512 KB
testcase_24 AC 203 ms
40,708 KB
testcase_25 AC 388 ms
70,868 KB
testcase_26 AC 17 ms
6,736 KB
testcase_27 AC 87 ms
20,584 KB
testcase_28 AC 96 ms
22,412 KB
testcase_29 AC 66 ms
16,752 KB
testcase_30 AC 236 ms
46,288 KB
testcase_31 AC 171 ms
35,660 KB
testcase_32 AC 285 ms
55,428 KB
testcase_33 AC 292 ms
55,792 KB
testcase_34 AC 391 ms
72,072 KB
testcase_35 AC 1 ms
4,380 KB
testcase_36 AC 2 ms
4,380 KB
testcase_37 AC 2 ms
4,384 KB
testcase_38 AC 2 ms
4,380 KB
testcase_39 AC 2 ms
4,380 KB
testcase_40 AC 2 ms
4,380 KB
testcase_41 AC 1 ms
4,380 KB
testcase_42 AC 2 ms
4,384 KB
testcase_43 AC 1 ms
4,384 KB
testcase_44 AC 2 ms
4,384 KB
testcase_45 AC 302 ms
76,344 KB
testcase_46 AC 291 ms
76,264 KB
testcase_47 AC 376 ms
76,336 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
class DSU_IDX{
    private:
        int n;
        vector<int> par_size, num_pair;

    public:
        int score;
        DSU_IDX():n(0),score(0){}
        DSU_IDX(int size):n(size),score(0),par_size(size,-1),num_pair(size,0){}
        int root(int v){
            if(par_size[v] < 0)return v;
            while(par_size[par_size[v]] >= 0){
                int ps = par_size[v];
                par_size[v] = par_size[ps];
                v = ps;
            }
            return par_size[v];
        }

        // 連結の際にスコアを更新する
        void unite(int u,int v){
            int ru = root(u),rv = root(v);
            if(ru == rv)return;

            score -= sub_score(ru) + sub_score(rv);
            if(par_size[ru] > par_size[rv]){
                par_size[rv] += par_size[ru];
                par_size[ru] = rv;
                num_pair[rv] += num_pair[ru];
                score += sub_score(rv);
            }
            else{
                par_size[ru] += par_size[rv];
                par_size[rv] = ru;
                num_pair[ru] += num_pair[rv];
                score += sub_score(ru);
            }
        }
        // v を含む連結成分の表す区間に配置可能なペアを増加
        // スコアを更新する
        void add_pair(int v){
            int rv = root(v);
            score -= sub_score(rv);
            num_pair[rv]++;
            score += sub_score(rv);
        }

        bool same(int u,int v){
            return root(u) == root(v);
        }

        int size(int v){
            return -par_size[root(v)];
        }

        int sub_score(int v){
            return min(-par_size[root(v)],num_pair[root(v)]);
        }


};
class DSU_AB{
    private:
        int n;
        vector<int> par_size,L,R;

    public:
        DSU_AB():n(0){}
        DSU_AB(int size):n(size),par_size(size,-1),L(size),R(size){
            for(int i=0;i<size;i++)L[i] = R[i] = i;
        }
        int root(int v){
            if(par_size[v] < 0)return v;
            while(par_size[par_size[v]] >= 0){
                int ps = par_size[v];
                par_size[v] = par_size[ps];
                v = ps;
            }
            return par_size[v];
        }

        void unite(int u,int v){
            int ru = root(u),rv = root(v);
            if(ru == rv)return;
            if(par_size[ru] > par_size[rv]){
                par_size[rv] += par_size[ru];
                par_size[ru] = rv;
                L[rv] = min(L[rv],L[ru]);
                R[rv] = max(R[rv],R[ru]);
            }
            else{
                par_size[ru] += par_size[rv];
                par_size[rv] = ru;
                L[ru] = min(L[rv],L[ru]);
                R[ru] = max(R[rv],R[ru]);
            }
        }

        bool same(int u,int v){
            return root(u) == root(v);
        }

        int size(int v){
            return -par_size[root(v)];
        }

        int left(int v){
            return L[root(v)];
        }
        int right(int v){
            return R[root(v)];
        }
};
int main(){
    int n;
    cin >> n;
    vector<int>a(n),b(n),z(n);
    for(int i=0;i<n;i++)cin >> a[i];
    for(int i=0;i<n;i++)cin >> b[i];
    for(int i=0;i<n;i++)cin >> z[i];

    vector<int>ida(n),idb(n);
    for(int i=0;i<n;i++){
        ida[--a[i]] = i;
        idb[--b[i]] = i;
        --z[i];
    }

    //並列二分探索

    vector<int>pos_left(n),l(n),r(n);
    for(int i=0;i<n;i++){
        l[i] = -1;r[i] = n;
    }
    vector<queue<int>>mid(n,queue<int>());
    // mid[i] : k = i で判定を行うペア番号の集合
    for(int _ = 0;_ < 18;_++){
        for(int i=0;i<n;i++){
            if(r[i]-l[i]>1){
                mid[(l[i]+r[i])/2].push(i);
            }
        }
        DSU_AB ufa(n),ufb(n);
        // i を含む連結成分には、a で i が動ける区間 [l, r] を持たせる
        // 同じ連結成分の要素は動ける区間も同じ

        for(int k=0;k<n;k++){
            int ia = ida[k],ib = idb[k];
            if(ia>0&&a[ia-1]<=k)ufa.unite(ia-1,ia);
            if(ia+1<n&&a[ia+1]<=k)ufa.unite(ia,ia+1);
            if(ib>0&&b[ib-1]<=k)ufb.unite(ib-1,ib);
            if(ib+1<n&&b[ib+1]<=k)ufb.unite(ib,ib+1);

            while(!mid[k].empty()){
                int j = mid[k].front();
                mid[k].pop();

                int la = ufa.left(ida[j]);
                int ra = ufa.right(ida[j]);
                int lb = ufb.left(idb[z[j]]);
                int rb = ufb.right(idb[z[j]]);
                // この k の値において、A の j と B の z[j] が同じ位置に置けるか、
                // 置ける場合はその位置の左端はどこか

                if(ra < lb || rb < la){
                    l[j] = k;
                }
                else{
                    r[j] = k;
                    pos_left[j] = max(la,lb);
                }
            }
        }
    }
    for(int i=0;i<n;i++)mid[r[i]].push(i);


    vector<int>less_k(n);
    // less_k[i] = (a_i <= k) + (b_i <= k)
    DSU_IDX ufi(n);
    // インデックス [1,n] を管理する

    for(int k=0;k<n;k++){
        int ia = ida[k], ib = idb[k];
        less_k[ia]++;less_k[ib]++;
        if(less_k[ia]==2){
            if(ia>0&&less_k[ia-1]==2)ufi.unite(ia-1,ia);
            if(ia+1<n&&less_k[ia+1]==2)ufi.unite(ia,ia+1);
        }
        if(less_k[ib]==2){
            if(ib>0&&less_k[ib-1]==2)ufi.unite(ib-1,ib);
            if(ib+1<n&&less_k[ib+1]==2)ufi.unite(ib,ib+1);
        }

        while(!mid[k].empty()){
            ufi.add_pair(pos_left[mid[k].front()]);
            // pos_left[i] を含む連結成分が表す区間で、(i, z[i]) を同じ位置に配置可能になる
            mid[k].pop();
        }

        cout << ufi.score << "\n";
    }


    return 0;
}
0