結果

問題 No.1975 Zigzag Sequence
ユーザー monnu
提出日時 2022-06-11 07:38:13
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 103 ms / 2,000 ms
コード長 1,708 bytes
コンパイル時間 4,047 ms
コンパイル使用メモリ 237,444 KB
実行使用メモリ 8,972 KB
最終ジャッジ日時 2024-09-21 15:48:04
合計ジャッジ時間 7,287 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 33
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
using ll=long long;
using Graph=vector<vector<int>>;
#define INF 1000000000000000000
#define MOD 1000000007
#define MAX 500

class BIT{
  int n;
  vector<ll> a;
public:
  BIT(int n_):n(n_){
    a.resize(n+1,0);
  }
  void add(int i,ll x){
    while(i<=n){
      a[i]+=x;
      a[i]%=MOD;
      i+=i&(-i);
    }
  }
  ll sum(int i){
    ll ret=0;
    while(i>0){
      ret+=a[i];
      ret%=MOD;
      i-=i&(-i);
    }
    return ret;
  }
  ll sum(int l,int r){
    return (sum(r)+MOD-sum(l-1))%MOD;
  }
};

int main(){
  int N;
  cin>>N;
  vector<pair<int,int>> A(N);
  for(int i=0;i<N;i++){
    cin>>A[i].first;
    A[i].second=i;
  }
  sort(A.begin(),A.end());

  vector<ll> p2(N+1,1);
  for(int i=0;i<N;i++){
    p2[i+1]=2*p2[i]%MOD;
  }

  ll ans=0;
  BIT left_tree1(N);
  BIT right_tree1(N);
  int i=0;

  while(i<N){
    int a=A[i].first;
    vector<int> index;
    while(i<N&&A[i].first==a){
      index.push_back(A[i].second);
      i++;
    }
    for(int k:index){
      ans+=left_tree1.sum(k+1)*right_tree1.sum(k+1,N)%MOD;
      ans%=MOD;
    }
    for(int k:index){
      left_tree1.add(k+1,p2[k]);
      right_tree1.add(k+1,p2[N-(k+1)]);
    }
  }

  reverse(A.begin(),A.end());
  BIT left_tree2(N);
  BIT right_tree2(N);
  i=0;
  while(i<N){
    int a=A[i].first;
    vector<int> index;
    while(i<N&&A[i].first==a){
      index.push_back(A[i].second);
      i++;
    }
    for(int k:index){
      ans+=left_tree2.sum(k+1)*right_tree2.sum(k+1,N)%MOD;
      ans%=MOD;
    }
    for(int k:index){
      left_tree2.add(k+1,p2[k]);
      right_tree2.add(k+1,p2[N-(k+1)]);
    }
  }

  cout<<ans<<'\n';
}
0