結果

問題 No.2046 Ans Mod? Mod Ans!
ユーザー t98slider
提出日時 2022-08-24 09:22:54
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 168 ms / 4,000 ms
コード長 1,267 bytes
コンパイル時間 1,784 ms
コンパイル使用メモリ 171,244 KB
実行使用メモリ 7,168 KB
最終ジャッジ日時 2024-12-16 08:35:47
合計ジャッジ時間 3,638 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;

template <class T> struct fenwick_tree {
    using U = T;

  public:
    fenwick_tree() : _n(0) {}
    fenwick_tree(int n) : _n(n), data(n) {}

    void add(int p, T x) {
        assert(0 <= p && p < _n);
        p++;
        while (p <= _n) {
            data[p - 1] += U(x);
            p += p & -p;
        }
    }

    T sum(int l, int r) {
        assert(0 <= l && l <= r && r <= _n);
        return sum(r) - sum(l);
    }

  private:
    int _n;
    std::vector<U> data;

    U sum(int r) {
        U s = 0;
        while (r > 0) {
            s += data[r - 1];
            r -= r & -r;
        }
        return s;
    }
};

int main(){
    int n, v, r = 200000;
    cin >> n;
    vector<int> tb(r + 1);
    fenwick_tree<ll> fw1(r + 1), fw2(r + 1);
    for(int i = 0; i < n; i++){
        cin >> v;
        tb[v]++;
        fw1.add(v, v);
        fw2.add(v, 1);
    }
    ll ans = 0;
    for(int i = 1; i <= r; i++){
        if(tb[i] == 0)continue;
        ans += fw2.sum(i + 1, r + 1) * tb[i] * i;
        for(int j = 1; i * j <= r; j++){
            ans -= fw1.sum(i * j, min(r + 1, i * (j + 1))) - fw2.sum(i * j, min(r + 1, i * (j + 1))) * i * j;
        }
    }
    cout << ans << '\n';
}
0