結果
問題 |
No.3305 Shift Sort
|
ユーザー |
![]() |
提出日時 | 2025-10-05 15:30:13 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 794 ms / 2,000 ms |
コード長 | 4,833 bytes |
コンパイル時間 | 5,891 ms |
コンパイル使用メモリ | 334,484 KB |
実行使用メモリ | 58,952 KB |
最終ジャッジ日時 | 2025-10-05 15:30:35 |
合計ジャッジ時間 | 21,085 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#include <bits/stdc++.h> using namespace std; #include <atcoder/all> using namespace atcoder; template<typename T> inline bool chmax(T &a, T b) { return ((a < b) ? (a = b, true) : (false)); } template<typename T> inline bool chmin(T &a, T b) { return ((a > b) ? (a = b, true) : (false)); } #define rep(i, n) for (long long i = 0; i < (long long)(n); i++) #define rep2(i, m ,n) for (int i = (m); i < (long long)(n); i++) #define REP(i, n) for (long long i = 1; i < (long long)(n); i++) typedef long long ll; #define updiv(N,X) (N + X - 1) / X #define l(n) n.begin(),n.end() #define YesNo(Q) Q==1?cout<<"Yes":cout<<"No" using P = pair<int, int>; using mint = modint; const int MOD = 998244353LL; const ll INF = 999999999999LL; vector<long long> fact, fact_inv, inv; /* init_nCk :二項係数のための前処理 計算量:O(n) */ template <typename T> void input(vector<T> &v){ rep(i,v.size()){cin>>v[i];} return; } void init_nCk(int SIZE) { fact.resize(SIZE + 5); fact_inv.resize(SIZE + 5); inv.resize(SIZE + 5); fact[0] = fact[1] = 1; fact_inv[0] = fact_inv[1] = 1; inv[1] = 1; for (int i = 2; i < SIZE + 5; i++) { fact[i] = fact[i - 1] * i % MOD; inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD; fact_inv[i] = fact_inv[i - 1] * inv[i] % MOD; } } /* nCk :MODでの二項係数を求める(前処理 int_nCk が必要) 計算量:O(1) */ long long nCk(int n, int k) { assert(!(n < k)); assert(!(n < 0 || k < 0)); return fact[n] * (fact_inv[k] * fact_inv[n - k] % MOD) % MOD; } long long modpow(long long a, long long n, long long mod) { long long res = 1; while (n > 0) { if (n & 1) res = res * a % mod; a = a * a % mod; n >>= 1; } return res; } ll POW(ll a,ll n){ long long res = 1; while (n > 0) { if (n & 1) res = res * a; a = a * a; n >>= 1; } return res; } mint pro(ll a){ if(a<=1){return (mint)(0);} else{ mint x = (a-1)/2; mint y = a/2; mint ret = (x*x+x+(y-1)*y); ret /= 2; return ret; } } mint pro2(ll a,ll b){ mint r1 = (b-1+(b-a)); r1 *= a; r1 /= 2; if((a-1)%2==0){ r1 -= (a+1)/2; } else{ r1 -= a/2; } r1 /= 2; return r1; } ll up=1; bool ch(ll t){ return up>=t; //upは、それ自身 } void test(){ ll n,k;cin>>n>>k; mint ret = pro(POW(10LL,k)); //mint ret = 1; ret *= pro(n/POW(10LL,k)); mint ret2 = (n/POW(10LL,k))/2; ret2 *= pro2(n%POW(10LL,k),POW(10LL,k)); cout << (ret2+ret).val() << endl; } ll e(){ return -1; } ll mx(ll a,ll b){ return max(a,b); } template <typename T> struct merge_sort_tree { int N; std::vector<std::vector<T>> x; merge_sort_tree() = default; merge_sort_tree(const std::vector<T> &vec) : N(vec.size()) { x.resize(N * 2); for (int i = 0; i < N; i++) x[N + i] = {vec[i]}; for (int i = N - 1; i; i--) { std::merge(x[i * 2].begin(), x[i * 2].end(), x[i * 2 + 1].begin(), x[i * 2 + 1].end(), std::back_inserter(x[i])); } } int cntLess(int l, int r, T query) const { l += N, r += N; int ret = 0; while (l < r) { if (l & 1) ret += std::lower_bound(x[l].begin(), x[l].end(), query) - x[l].begin(), l++; if (r & 1) r--, ret += std::lower_bound(x[r].begin(), x[r].end(), query) - x[r].begin(); l >>= 1, r >>= 1; } return ret; } int cntLesseq(int l, int r, T query) const { l += N, r += N; int ret = 0; while (l < r) { if (l & 1) ret += std::upper_bound(x[l].begin(), x[l].end(), query) - x[l].begin(), l++; if (r & 1) r--, ret += std::upper_bound(x[r].begin(), x[r].end(), query) - x[r].begin(); l >>= 1, r >>= 1; } return ret; } int cntMore(int begin, int end, T query) const { int tot = std::max(0, std::min(end, N) - std::max(begin, 0)); return tot - cntLesseq(begin, end, query); } int cntMoreeq(int begin, int end, T query) const { int tot = std::max(0, std::min(end, N) - std::max(begin, 0)); return tot - cntLess(begin, end, query); } template <class OStream> friend OStream &operator<<(OStream &os, const merge_sort_tree &clt) { os << '['; for (int i = 0; i < clt.N; i++) os << clt.x[clt.N + i][0] << ','; return os << ']'; } }; int main(){ int n;cin>>n;int q;cin>>q; segtree<ll,mx,e> seg(n); vector<ll> v(n);rep(i,n){cin>>v[i];seg.set(i,v[i]);} vector<ll> a(n); rep(i,n){up = v[i];a[i] = seg.min_left<ch>(i+1)-1;} merge_sort_tree<ll> mg(a); while(q--){ ll x,y;cin>>x>>y;x--; cout << mg.cntMoreeq(x, y, x) << endl; } }