#613275 (C++17) No.1371 交換門松列・松

提出ソース

結果

問題	No.1371 交換門松列・松
ユーザー	rniya
提出日時	2021-02-05 17:36:44
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	101 ms / 4,000 ms
コード長	6,793 bytes
コンパイル時間	2,125 ms
コンパイル使用メモリ	221,236 KB
実行使用メモリ	15,924 KB
最終ジャッジ日時	2024-07-02 07:44:15
合計ジャッジ時間	5,062 ms
ジャッジサーバーID （参考情報）	judge3 / judge1

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 6,816 KB
testcase_01	AC	2 ms 6,944 KB
testcase_02	AC	2 ms 6,944 KB
testcase_03	AC	2 ms 6,940 KB
testcase_04	AC	2 ms 6,940 KB
testcase_05	AC	2 ms 6,940 KB
testcase_06	AC	2 ms 6,940 KB
testcase_07	AC	2 ms 6,940 KB
testcase_08	AC	1 ms 6,940 KB
testcase_09	AC	1 ms 6,940 KB
testcase_10	AC	1 ms 6,944 KB
testcase_11	AC	1 ms 6,944 KB
testcase_12	AC	2 ms 6,940 KB
testcase_13	AC	2 ms 6,940 KB
testcase_14	AC	2 ms 6,940 KB
testcase_15	AC	2 ms 6,944 KB
testcase_16	AC	1 ms 6,940 KB
testcase_17	AC	2 ms 6,944 KB
testcase_18	AC	60 ms 13,608 KB
testcase_19	AC	60 ms 13,692 KB
testcase_20	AC	73 ms 15,596 KB
testcase_21	AC	65 ms 15,684 KB
testcase_22	AC	62 ms 15,924 KB
testcase_23	AC	101 ms 14,184 KB
testcase_24	AC	99 ms 14,164 KB
testcase_25	AC	98 ms 13,808 KB
testcase_26	AC	94 ms 13,840 KB
testcase_27	AC	93 ms 13,796 KB
testcase_28	AC	96 ms 13,860 KB
testcase_29	AC	94 ms 13,872 KB
testcase_30	AC	97 ms 13,620 KB
testcase_31	AC	95 ms 13,772 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#define LOCAL
#include <bits/stdc++.h>
using namespace std;
#pragma region Macros
typedef long long ll;
typedef __int128_t i128;
typedef unsigned int uint;
typedef unsigned long long ull;
#define ALL(x) (x).begin(), (x).end()

template <typename T> istream& operator>>(istream& is, vector<T>& v) {
    for (T& x : v) is >> x;
    return is;
}
template <typename T> ostream& operator<<(ostream& os, const vector<T>& v) {
    for (int i = 0; i < (int)v.size(); i++) {
        os << v[i] << (i + 1 == (int)v.size() ? "" : " ");
    }
    return os;
}
template <typename T, typename U> ostream& operator<<(ostream& os, const pair<T, U>& p) {
    os << '(' << p.first << ',' << p.second << ')';
    return os;
}
template <typename T, typename U, typename V> ostream& operator<<(ostream& os, const tuple<T, U, V>& t) {
    os << '(' << get<0>(t) << ',' << get<1>(t) << ',' << get<2>(t) << ')';
    return os;
}
template <typename T, typename U, typename V, typename W> ostream& operator<<(ostream& os, const tuple<T, U, V, W>& t) {
    os << '(' << get<0>(t) << ',' << get<1>(t) << ',' << get<2>(t) << ',' << get<3>(t) << ')';
    return os;
}
template <typename T, typename U> ostream& operator<<(ostream& os, const map<T, U>& m) {
    os << '{';
    for (auto itr = m.begin(); itr != m.end();) {
        os << '(' << itr->first << ',' << itr->second << ')';
        if (++itr != m.end()) os << ',';
    }
    os << '}';
    return os;
}
template <typename T, typename U> ostream& operator<<(ostream& os, const unordered_map<T, U>& m) {
    os << '{';
    for (auto itr = m.begin(); itr != m.end();) {
        os << '(' << itr->first << ',' << itr->second << ')';
        if (++itr != m.end()) os << ',';
    }
    os << '}';
    return os;
}
template <typename T> ostream& operator<<(ostream& os, const set<T>& s) {
    os << '{';
    for (auto itr = s.begin(); itr != s.end();) {
        os << *itr;
        if (++itr != s.end()) os << ',';
    }
    os << '}';
    return os;
}
template <typename T> ostream& operator<<(ostream& os, const multiset<T>& s) {
    os << '{';
    for (auto itr = s.begin(); itr != s.end();) {
        os << *itr;
        if (++itr != s.end()) os << ',';
    }
    os << '}';
    return os;
}
template <typename T> ostream& operator<<(ostream& os, const unordered_set<T>& s) {
    os << '{';
    for (auto itr = s.begin(); itr != s.end();) {
        os << *itr;
        if (++itr != s.end()) os << ',';
    }
    os << '}';
    return os;
}
template <typename T> ostream& operator<<(ostream& os, const deque<T>& v) {
    for (int i = 0; i < (int)v.size(); i++) {
        os << v[i] << (i + 1 == (int)v.size() ? "" : " ");
    }
    return os;
}

void debug_out() { cerr << '\n'; }
template <class Head, class... Tail> void debug_out(Head&& head, Tail&&... tail) {
    cerr << head;
    if (sizeof...(Tail) > 0) cerr << ", ";
    debug_out(move(tail)...);
}
#ifdef LOCAL
#define debug(...)                                                                   \
    cerr << " ";                                                                     \
    cerr << #__VA_ARGS__ << " :[" << __LINE__ << ":" << __FUNCTION__ << "]" << '\n'; \
    cerr << " ";                                                                     \
    debug_out(__VA_ARGS__)
#else
#define debug(...) 42
#endif

template <typename T> T gcd(T x, T y) { return y != 0 ? gcd(y, x % y) : x; }
template <typename T> T lcm(T x, T y) { return x / gcd(x, y) * y; }

template <class T1, class T2> inline bool chmin(T1& a, T2 b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template <class T1, class T2> inline bool chmax(T1& a, T2 b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
#pragma endregion

/**
 * @brief Binary Indexed Tree
 * @docs docs/datastructure/BinaryIndexedTree.md
 */
template <typename T> class BinaryIndexedTree {
    T sum(int i) {
        T res = T();
        for (; i > 0; i -= (i & -i)) res += dat[i];
        return res;
    }

public:
    int n;
    vector<T> dat;
    BinaryIndexedTree(int n_) : n(n_ + 1), dat(n + 1, 0) {}
    void add(int i, const T& x) {
        for (++i; i <= n; i += (i & -i)) dat[i] += x;
    }
    T query(int l, int r) { return sum(r) - sum(l); }
    int lower_bound(T x) const {
        if (x <= 0) return 0;
        int pos = 0, k = 1;
        while (k < n) k <<= 1;
        for (; k > 0; k >>= 1) {
            if (pos + k <= n && dat[pos + k] < x) {
                x -= dat[pos + k];
                pos += k;
            }
        }
        return pos;
    }
    int upper_bound(T x) const { return lower_bound(x + 1); }
    T operator[](int i) { return query(i, i + 1); }
};

const int INF = 1e9;
const long long IINF = 1e18;
const int dx[4] = {1, 0, -1, 0}, dy[4] = {0, 1, 0, -1};
const char dir[4] = {'D', 'R', 'U', 'L'};
const long long MOD = 1000000007;
// const long long MOD = 998244353;

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    int N;
    cin >> N;
    vector<int> A(N);
    cin >> A;

    vector<int> group(N);  // 1:両隣より大きい 0:両隣より小さい
    for (int i = 0; i < N; i++) group[i] = ((i & 1) ^ (A[0] > A[1]));

    // 入れ替えは1及び0同士でしか発生し得ない
    // 隣接での交換は許されえない(元々門松列)
    // 各場所での有効範囲を求めて、お互いにその範囲に入っているような2数の数え上げ
    // 自分同士が入れ替え可能の判定になってしまうことに注意

    vector<int> Min(N), Max(N);
    for (int i = 0; i < N; i++) {
        if (group[i]) {
            Min[i] = max((i - 1 >= 0 ? A[i - 1] : 1), (i + 1 < N ? A[i + 1] : 1));
            Max[i] = N + 1;
        } else {
            Min[i] = 1;
            Max[i] = min((i - 1 >= 0 ? A[i - 1] : N + 1), (i + 1 < N ? A[i + 1] : N + 1));
        }
    }

    vector<tuple<int, int, int>> u, v;
    for (int i = 0; i < N; i++) {
        u.emplace_back(A[i], Min[i], Max[i]);
        v.emplace_back(Min[i], A[i], Max[i]);
    }

    sort(ALL(u));
    sort(ALL(v));

    priority_queue<pair<int, int>, vector<pair<int, int>>, greater<pair<int, int>>> pq;
    BinaryIndexedTree<int> BIT(N + 1);
    ll ans = 0;

    auto PUSH = [&](int i) {
        int val, l, r;
        tie(l, val, r) = v[i];
        pq.emplace(r, val);
        BIT.add(val, 1);
    };
    auto POP = [&]() {
        auto p = pq.top();
        pq.pop();
        int val = p.second;
        BIT.add(val, -1);
    };

    for (int i = 0, j = 0; i < N; i++) {
        int val, l, r;
        tie(val, l, r) = u[i];
        while (j < N && get<0>(v[j]) <= val) PUSH(j++);
        while (!pq.empty() && pq.top().first <= val) POP();
        ans += BIT.query(l, r) - 1;
    }

    cout << ans / 2 << '\n';
    return 0;
}

yukicoder

結果

テストケース

ソースコード