QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#643273#7904. Rainbow Subarraytkt0506TL 2ms7764kbC++142.6kb2024-10-15 20:16:332024-10-15 20:16:35

Judging History

你现在查看的是最新测评结果

  • [2024-10-15 20:16:35]
  • 评测
  • 测评结果:TL
  • 用时:2ms
  • 内存:7764kb
  • [2024-10-15 20:16:33]
  • 提交

answer

#include <bits/stdc++.h>
#pragma comment(linker, "/stack:200000000")
#pragma GCC optimize("Ofast,no-stack-protector")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
#pragma GCC optimize("unroll-loops")
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma G++ optimize(2)
#pragma G++ optimize(3)

#define int long long
#define pii pair<int, int>
#define all(a) a.begin(), a.end()
using namespace std;

int INF = 1e18; 
const int MOD = 998244353;
const int N = 5e5 + 10;
const int LOG = 20;
const int dx[4] = {0, 0, 1, -1};
const int dy[4] = {1, -1, 0, 0};

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

int BIT[N], cBIT[N], a[N];
int sz;

void add(int x, int val, int c){
    for(int i = x; i <= sz; i += i&-i){
        BIT[i] += val, cBIT[i] += c;
    }
}

int qry(int x){
    int res = 0;
    for(int i = x; i > 0; i -= i&-i) res += BIT[i];
    return res;
}
    
int cqry(int x){
    int res = 0;
    for(int i = x; i > 0; i -= i&-i) res += cBIT[i];
    return res;
}

vector<int> nums;
int n, k;

int find(int x){
    int res = lower_bound(all(nums), x) - nums.begin();
    if(nums[res] == x) res++;
    return res;
}

int cal(int l){
        int res = INF;

        for(int i = 1; i <= l-1; i++){
            add(find(a[i]), a[i], 1);
        }

        for(int i = l; i <= n; i++){
            add(find(a[i]), a[i], 1);

            int avg = qry(nums.size()) / l;
            for(int j = -1; j <= 1; j++){
                avg += j;

                int l_c = cqry(find(avg));
                int r_c = cqry(nums.size()) - l_c;
                int l_sum = qry(find(avg));
                int r_sum = qry(nums.size()) - l_sum;
                 
                res = min(res, l_c * avg - l_sum + r_sum - r_c * avg);

                avg -= j;
            }

            add(find(a[i - l + 1]), -a[i - l + 1], -1);
        }

        for(int i = n - l + 2; i <= n; i++) add(find(a[i]), -a[i], -1);

        return res;
    };

void solve(){
    cin >> n >> k;    
    nums.clear();
    nums.push_back(0);
    for(int i = 1; i <= n; i++){
        cin >> a[i];
        a[i] -= i;
        if(a[i] < 0) a[i] = -a[i];
        nums.push_back(a[i]);
    }

    sort(all(nums));
    nums.erase(unique(all(nums)), nums.end());
    sz = nums.size();
    int i = 1;
    for(int j = n-1; j > 0; j /= 2){
        while(i + j <= n && cal(i + j) <= k) i += j;
    }

    cout << i << '\n';
}

int32_t main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    
    int tt;
    cin >> tt;

    while(tt--){
        solve();
    }
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 2ms
memory: 7764kb

input:

5
7 5
7 2 5 5 4 11 7
6 0
100 3 4 5 99 100
5 6
1 1 1 1 1
5 50
100 200 300 400 500
1 100
3

output:

4
3
5
1
1

result:

ok 5 lines

Test #2:

score: -100
Time Limit Exceeded

input:

11102
2 167959139
336470888 134074578
5 642802746
273386884 79721198 396628655 3722503 471207868
6 202647942
268792718 46761498 443917727 16843338 125908043 191952768
2 717268783
150414369 193319712
6 519096230
356168102 262263554 174936674 407246545 274667941 279198849
9 527268921
421436316 3613460...

output:

1
4
3
2
6
5
7
2
4
1
4
1
1
3
2
2
6
8
7
7
1
7
6
2
4
3
1
4
7
7
3
4
3
9
3
8
6
6
2
1
6
3
1
2
4
6
4
5
4
1
4
7
1
6
3
5
6
6
1
7
5
3
1
6
3
4
3
2
2
6
2
3
10
1
4
3
2
4
5
1
7
5
5
4
6
5
3
6
2
5
5
8
5
4
5
2
1
5
2
3
3
4
7
1
3
1
2
2
6
3
1
6
8
1
8
4
5
5
6
7
4
8
3
2
8
4
5
5
2
6
2
4
1
5
4
5
3
2
4
1
2
1
4
4
8
3
7
3
3
3...

result: