結果
| 問題 | No.1733 Sum of Sorted Subarrays | 
| コンテスト | |
| ユーザー |  daddy | 
| 提出日時 | 2022-12-09 03:12:30 | 
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 533 ms / 3,000 ms | 
| コード長 | 1,552 bytes | 
| コンパイル時間 | 1,277 ms | 
| コンパイル使用メモリ | 95,892 KB | 
| 実行使用メモリ | 15,752 KB | 
| 最終ジャッジ日時 | 2024-10-14 18:32:31 | 
| 合計ジャッジ時間 | 11,248 ms | 
| ジャッジサーバーID (参考情報) | judge2 / judge3 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 24 | 
ソースコード
#include <algorithm>
#include <atcoder/lazysegtree>
#include <atcoder/modint>
#include <functional>
#include <iostream>
#include <utility>
#include <vector>
using namespace std;
using mint = atcoder::modint998244353;
mint op(mint a, mint b) {
    return a + b;
}
mint e() {
    return mint(0);
}
mint mapping(mint f, mint x) {
    return f * x;
}
mint composition(mint f, mint g) {
    return f * g;
}
mint id() {
    return mint(1);
}
vector<size_t> sorted_idx(const vector<size_t> A, bool rev) {
    auto n = A.size();
    vector<pair<size_t, long long>> tmp(n);
    for (size_t i = 0; i < n; i++) tmp[i] = {A[i], rev ? -i : i};
    stable_sort(begin(tmp), end(tmp));
    vector<size_t> ret(n);
    for (size_t i = 0; i < n; i++) ret[i] = abs(tmp[i].second);
    return ret;
}
vector<mint> cumulative_sum(const vector<size_t> &A, bool rev = false) {
    auto n = A.size();
    atcoder::lazy_segtree<mint, op, e, mint, mapping, composition, id> lst(
        vector<mint>(n, 1));
    vector<mint> ret(n);
    auto &&tmp = sorted_idx(A, rev);
    for (auto &i : tmp) {
        ret[i] = lst.prod(0, i + 1) / lst.get(i);
        lst.apply(0, i + 1, mint(2));
    }
    return ret;
}
int main() {
    size_t n;
    cin >> n;
    vector<size_t> A(n);
    for (size_t i = 0; i < n; i++) cin >> A[i];
    auto &&x = cumulative_sum(A);
    reverse(begin(A), end(A));
    auto &&y = cumulative_sum(A, true);
    reverse(begin(x), end(x));
    mint ans;
    for (size_t i = 0; i < n; i++) { ans += x[i] * y[i] * A[i]; }
    cout << ans.val() << endl;
}
            
            
            
        