結果
| 問題 |
No.1921 Range LthKth Query
|
| ユーザー |
りあん
|
| 提出日時 | 2022-05-05 14:30:30 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 499 ms / 4,000 ms |
| コード長 | 3,833 bytes |
| コンパイル時間 | 3,920 ms |
| コンパイル使用メモリ | 263,496 KB |
| 実行使用メモリ | 101,504 KB |
| 最終ジャッジ日時 | 2024-07-04 15:56:42 |
| 合計ジャッジ時間 | 18,958 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 23 |
ソースコード
// #pragma GCC target("avx2")
#pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
using namespace std;
using P = pair<int, int>;
const int M = 1000000007;
template<typename T> class BIT {
private:
int n,m; vector<vector<T> > bit;
public:
// (i, j) に val を加算する
void add(int i, int j, T val){
for(int i_ = i+1; i_ < n; i_ += i_ & -i_)
for(int j_ = j+1; j_ < m; j_ += j_ & -j_)
bit[i_][j_] += val;
}
// [0,i]×[0,j]の範囲の和を求める
T sum(int i, int j){
T s = 0;
for(int i_ = i+1; i_ > 0; i_ -= i_ & -i_)
for(int j_ = j+1; j_ > 0; j_ -= j_ & -j_)
s += bit[i_][j_];
return s;
}
// [lx, rx)×[ly, ry)の範囲の和を求める
T sum(int lx, int rx, int ly, int ry){
return sum(rx-1, ry-1) - sum(lx-1, ry-1) - sum(rx-1, ly-1) + sum(lx-1, ly-1);
}
BIT(int sz1, int sz2){
n = sz1 + 1, m = sz2 + 1;
bit.resize(n, vector<T>(m, 0));
}
BIT(vector<vector<T> >& v){
n = (int)v.size() + 1, m = (int)v[0].size() + 1;
bit.resize(n, vector<T>(m, 0));
for(int i = 0; i < n - 1; i++)
for(int j = 0; j < m - 1; j++)
add(i, j, v[i][j]);
}
void print(){
for(int i = 0; i < n; i++){
for(int j = 0; j < m; j++){
cout<< sum(i-1, i, j-1, j) << " ";
}
cout << "\n";
}
}
void print_sum(){
for(int i = 0; i < n; i++){
for(int j = 0; j < m; j++){
cout<< sum(i-1, j-1) << " ";
}
cout << "\n";
}
}
};
int main() {
cin.tie(0);
ios::sync_with_stdio(0);
int n, K, L;
cin >> n >> K >> L;
vector<int> a(n);
for (int i = 0; i < n; ++i) {
cin >> a[i];
}
vector<vector<int>> kth(n + 1, vector<int>(n + 1, M));
for (int i = 0; i < n + 1; ++i) {
for (int l = 0, r = 0, c = 0; l < n; ++l) {
while (r < n && c < K) {
if (a[r] <= i) ++c;
++r;
}
if (c == K) {
kth[l][r] = min(kth[l][r], i);
}
if (a[l] <= i) --c;
}
}
for (int d = 1; d <= n; ++d) {
for (int i = 0; i + d <= n; ++i) {
kth[i][i + d] = min(kth[i][i + d], kth[i + 1][i + d]);
kth[i][i + d] = min(kth[i][i + d], kth[i][i + d - 1]);
}
}
vector<vector<int>> cnth(n + 1, vector<int>(n + 1));
vector<vector<int>> cntw(n + 1, vector<int>(n + 1));
for (int i = 0; i < n + 1; ++i) {
for (int j = 0; j < n + 1; ++j) {
if (kth[i][j] != M) {
++cnth[i][kth[i][j]];
++cntw[j][kth[i][j]];
}
}
}
for (int i = 0; i < n + 1; ++i) {
for (int j = 0; j < n; ++j) {
cnth[i][j + 1] += cnth[i][j];
cntw[i][j + 1] += cntw[i][j];
}
}
vector<vector<int>> lth(n + 1, vector<int>(n + 1, M));
for (int i = 0; i < n + 1; ++i) {
for (int l = 0, r = 0, c = 0; l < n; ++l) {
while (r < n && c < L) {
++r;
c += max(cntw[r][i] - l, 0);
}
if (c >= L) {
lth[l][r] = min(lth[l][r], i);
}
c -= max(cnth[l][i] - (n - r), 0);
}
}
for (int d = 1; d <= n; ++d) {
for (int i = 0; i + d <= n; ++i) {
lth[i][i + d] = min(lth[i][i + d], lth[i + 1][i + d]);
lth[i][i + d] = min(lth[i][i + d], lth[i][i + d - 1]);
}
}
int q;
cin >> q;
for (int _ = 0; _ < q; ++_) {
int l, r;
cin >> l >> r;
--l;
cout << lth[l][r] << '\n';
}
return 0;
}
りあん