結果

問題 No.3590 I Love Inversions
コンテスト
ユーザー V_Melville
提出日時 2026-07-17 22:19:33
言語 C++23
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 2,128 ms / 5,000 ms
+ 31µs
コード長 1,294 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,042 ms
コンパイル使用メモリ 335,192 KB
実行使用メモリ 5,888 KB
平均クエリ数 12104.46
最終ジャッジ日時 2026-07-17 22:19:58
合計ジャッジ時間 24,601 ms
ジャッジサーバーID
(参考情報)
judge3_0 / judge2_0
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code

#include <bits/stdc++.h>
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rep1(i, n) for (int i = 1; i <= (n); ++i)

using namespace std;
using ll = long long;

template<typename T>
struct BIT {
  int n;
  vector<T> d;
  BIT(int n=0):n(n),d(n+1) {}
  void add(int i, T x=1) {
    for (i++; i <= n; i += i&-i) {
      d[i] += x;
    }
  }
  T sum(int i) {
    T x = 0;
    for (i++; i; i -= i&-i) {
      x += d[i];
    }
    return x;
  }
  T sum(int l, int r) {
    return sum(r-1) - sum(l-1);
  }
  int kth(ll k) {
      int b = 1, i = 0;
      while (b*2 <= n) b <<= 1;
      while (b) {
          if (i+b <= n and d[i+b] < k) {
              k -= d[i+b];
              i += b;
          }
          b >>= 1;
      }
      return i;
  }
};

int main() {
    int n;
    cin >> n;
    
    vector<ll> inv(n+1);
    rep1(i, n) {
        cout << "? " << 1 << " " << i << endl;
        ll x;
        cin >> x;
        inv[i] = x;
    }
    
    vector<ll> g(n+1);
    rep1(i, n) g[i] = inv[i]-inv[i-1];
    
    BIT<ll> t(n);
    rep(i, n) t.add(i);
    
    vector<int> p(n+1);
    for (int i = n; i >= 1; --i) {
        int ni = t.kth(i-g[i]);
        p[i] = ni+1;
        t.add(ni, -1);
    }

    cout << "! ";
    rep1(i, n) cout << p[i] << " ";
    cout << endl;
    
    return 0;
}
0