結果
問題 | No.1975 Zigzag Sequence |
ユーザー |
![]() |
提出日時 | 2022-06-10 22:45:25 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 212 ms / 2,000 ms |
コード長 | 1,756 bytes |
コンパイル時間 | 2,152 ms |
コンパイル使用メモリ | 185,140 KB |
実行使用メモリ | 11,904 KB |
最終ジャッジ日時 | 2024-09-21 06:39:56 |
合計ジャッジ時間 | 7,116 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 33 |
ソースコード
#include <bits/stdc++.h>using namespace std;typedef long long int ll;typedef pair<ll,ll> P;typedef vector<ll> VI;typedef vector<VI> VVI;#define REP(i,n) for(ll i=0;i<(n);i++)#define ALL(v) v.begin(),v.end()template<typename T> bool chmax(T &x, const T &y) {return (x<y)?(x=y,true):false;};template<typename T> bool chmin(T &x, const T &y) {return (x>y)?(x=y,true):false;};constexpr ll MOD=1000000007;constexpr ll INF=2e18;struct BIT{int N;VI dat;BIT(int N):N(N){dat.resize(N+1,0);}ll sum(int k){k++;if(k<=0)return 0;ll res=0;for(;k>0;k-=k&-k)res=(res+dat[k])%MOD;return res;}ll sum(int l, int r){return (sum(r-1)-sum(l-1)+MOD)%MOD;}void add(int k, ll x){k++;for(;k<=N;k+=k&-k)dat[k]=(dat[k]+x)%MOD;}};ll power(ll x, ll y){x%=MOD;ll ret=1;while(y){if(y&1) ret=ret*x%MOD;x=x*x%MOD;y>>=1;}return ret;}int main(){int n; cin >> n;VI a(n); REP(i,n) cin >> a[i];VI sa=a; sort(ALL(sa));map<int,int> idx;REP(i,n) idx[sa[i]]=i;BIT btx(n), btl(n), btr(n);ll ans=0;REP(i,n){int j=idx[a[i]];btx.add(j,btx.sum(j,j+1));btx.add(j,btl.sum(j,j+1));btx.add(j,btr.sum(j,j+1));btx.add(j,1);ll lsl=btl.sum(0,j), lsr=btl.sum(j+1,n);ans=(ans+lsl*power(2,n-1-i)%MOD)%MOD;btr.add(j,lsl+btx.sum(0,j));btl.add(j,lsr);ll rsr=btr.sum(j+1,n), rsl=btr.sum(0,j);ans=(ans+rsr*power(2,n-1-i)%MOD)%MOD;btl.add(j,rsr+btx.sum(j+1,n));btr.add(j,rsl);}cout << ans << endl;return 0;}