結果

問題 No.1975 Zigzag Sequence
ユーザー SSRS
提出日時 2022-06-10 21:35:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 87 ms / 2,000 ms
コード長 1,541 bytes
コンパイル時間 1,799 ms
コンパイル使用メモリ 177,280 KB
実行使用メモリ 9,600 KB
最終ジャッジ日時 2024-09-21 06:02:28
合計ジャッジ時間 5,946 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 33
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 1000000007;
template <typename T>
struct binary_indexed_tree{
  int N;
  vector<T> BIT;
  binary_indexed_tree(int N): N(N), BIT(N + 1, 0){
  }
  void add(int i, T x){
    i++;
    while (i <= N){
      BIT[i] += x;
      BIT[i] %= MOD;
      i += i & -i;
    }
  }
  T sum(int i){
    T ans = 0;
    while (i > 0){
      ans += BIT[i];
      ans %= MOD;
      i -= i & -i;
    }
    return ans;
  }
  T sum(int L, int R){
    return (sum(R) - sum(L) + MOD) % MOD;
  }
};
int main(){
  int N;
  cin >> N;
  vector<int> A(N);
  for (int i = 0; i < N; i++){
    cin >> A[i];
  }
  vector<int> B = A;
  sort(B.begin(), B.end());
  B.erase(unique(B.begin(), B.end()), B.end());
  int cnt = B.size();
  for (int i = 0; i < N; i++){
    A[i] = lower_bound(B.begin(), B.end(), A[i]) - B.begin();
  }
  vector<long long> pow2(N);
  pow2[0] = 1;
  for (int i = 0; i < N - 1; i++){
    pow2[i + 1] = pow2[i] * 2 % MOD;
  }
  vector<long long> L1(N), L2(N), R1(N), R2(N);
  binary_indexed_tree<long long> BIT1(cnt);
  for (int i = 0; i < N; i++){
    L1[i] = BIT1.sum(0, A[i]);
    L2[i] = BIT1.sum(A[i] + 1, cnt);
    BIT1.add(A[i], pow2[i]);
  }
  binary_indexed_tree<long long> BIT2(cnt);
  for (int i = N - 1; i >= 0; i--){
    R1[i] = BIT2.sum(0, A[i]);
    R2[i] = BIT2.sum(A[i] + 1, cnt);
    BIT2.add(A[i], pow2[N - 1 - i]);
  }
  long long ans = 0;
  for (int i = 0; i < N; i++){
    ans += L1[i] * R1[i];
    ans += L2[i] * R2[i];
    ans %= MOD;
  }
  cout << ans << endl;
}
0