結果

問題 No.404 部分門松列
ユーザー anta
提出日時 2016-07-21 14:48:10
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 176 ms / 2,000 ms
コード長 3,000 bytes
コンパイル時間 1,639 ms
コンパイル使用メモリ 177,320 KB
実行使用メモリ 8,192 KB
最終ジャッジ日時 2024-06-28 16:53:06
合計ジャッジ時間 7,979 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include "bits/stdc++.h"
using namespace std;
#define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i))
#define rer(i,l,u) for(int (i)=(int)(l);(i)<=(int)(u);++(i))
#define reu(i,l,u) for(int (i)=(int)(l);(i)<(int)(u);++(i))
static const int INF = 0x3f3f3f3f; static const long long INFL = 0x3f3f3f3f3f3f3f3fLL;
typedef vector<int> vi; typedef pair<int, int> pii; typedef vector<pair<int, int> > vpii; typedef long long ll;
template<typename T, typename U> static void amin(T &x, U y) { if(y < x) x = y; }
template<typename T, typename U> static void amax(T &x, U y) { if(x < y) x = y; }
#ifdef NDEBUG
#error assert is disabled!
#endif
template<typename T>
T readNatural(T lo, T up) {
    assert(0 <= lo && lo <= up);
    T x = 0;
    while(1) {
        int d = getchar();
        if(!('0' <= d && d <= '9')) {
            ungetc(d, stdin);
            break;
        }
        d -= '0';
        assert(d <= up && x <= (up - d) / 10);
        x = x * 10 + d;
    }
    assert(lo <= x && x <= up);
    return x;
}
void readSpace() { int c = getchar(); assert(c == ' '); }
static bool read_eof = false;
void readEOL() {
    int c = getchar();
    if(c == EOF) {
        assert(!read_eof);
        read_eof = true;
    } else {
        assert(c == '\n');
    }
}
void readEOF() {
    assert(!read_eof);
    int c = getchar();
    assert(c == EOF);
    read_eof = true;
}
struct FenwickTree {
    typedef int T;
    vector<T> v;
    void init(int n) { v.assign(n, 0); }
    void add(int i, T x) {
        for(; i < (int)v.size(); i |= i + 1) v[i] += x;
    }
    T sum(int i) const {    //[0, i)
        T r = 0;
        for(-- i; i >= 0; i = (i & (i + 1)) - 1) r += v[i];
        return r;
    }
    T sum(int left, int right) const {  //[left, right)
        return sum(right) - sum(left);
    }
};
int main() {
    int N = readNatural(1, 200000);
    readEOL();
    vector<int> A(N);
    for(int i = 0; i < N; ++ i) {
        if(i != 0) readSpace();
        A[i] = readNatural(1, (int)1e9);
    }
    readEOL();
    vector<int> values = A;
    sort(values.begin(), values.end());
    values.erase(unique(values.begin(), values.end()), values.end());
    int X = (int)values.size();
    for(auto &x : A)
        x = (int)(lower_bound(values.begin(), values.end(), x) - values.begin());
    vector<long long> ansSum(X + 1);
    FenwickTree ftL, ftR;
    ftL.init(X);
    ftR.init(X);
    vector<int> numL(X), numR(X);
    ll sumProd = 0;
    for(int x : A) {
        ftR.add(x, 1);
        ++ numR[x];
    }
    for(int x : A) {
        ftR.add(x, -1);
        -- numR[x];
        sumProd -= numL[x];
        ll cnt = 0;
        cnt += (ll)ftL.sum(x) * ftR.sum(x);
        cnt += (ll)ftL.sum(x + 1, X) * ftR.sum(x + 1, X);
        cnt -= sumProd - (ll)numL[x] * numR[x];
        ansSum[x + 1] += cnt;
        ftL.add(x, 1);
        ++ numL[x];
        sumProd += numR[x];
    }
    rep(x, X)
        ansSum[x + 1] += ansSum[x];
    int Q = readNatural(1, 200000);
    readEOL();
    rep(ii, Q) {
        int L = readNatural(1, (int)1e9);
        readSpace();
        int R = readNatural(L, (int)1e9);
        readEOL();
        int l = (int)(lower_bound(values.begin(), values.end(), L) - values.begin());
        int r = (int)(upper_bound(values.begin(), values.end(), R) - values.begin());
        long long ans = ansSum[r] - ansSum[l];
        printf("%lld\n", ans);
    }
    readEOF();
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0