結果
| 問題 |
No.1975 Zigzag Sequence
|
| コンテスト | |
| ユーザー |
monnu
|
| 提出日時 | 2022-06-11 00:21:00 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,692 bytes |
| コンパイル時間 | 4,283 ms |
| コンパイル使用メモリ | 236,296 KB |
| 実行使用メモリ | 8,968 KB |
| 最終ジャッジ日時 | 2024-09-21 07:29:14 |
| 合計ジャッジ時間 | 7,205 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 26 WA * 7 |
ソースコード
#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];
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';
}
monnu