QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#620244#7904. Rainbow SubarraySkyEyeControllerTL 815ms68004kbC++234.2kb2024-10-07 17:08:452024-10-07 17:08:47

Judging History

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

  • [2024-10-07 17:08:47]
  • 评测
  • 测评结果:TL
  • 用时:815ms
  • 内存:68004kb
  • [2024-10-07 17:08:45]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
#define DEBUG system("pause");
#define i64 long long
#define int long long
struct node
{
    int l, r;
    i64 sum;
    i64 cnt;
};
int tot = 0;
const int maxn = 5e5 + 9;
i64 bs[maxn]; // 离散化数组,下标离散化,内存为原数字
i64 a[maxn];
node tree[maxn << 5];
int root[maxn];
void init(int n)
{
    for (int i = 0; i <= n; i++)
    {
        root[i] = 0;
    }
    tot = 0;
    return;
}
int newnode()
{
    tot++;
    tree[tot] = {0, 0, 0, 0};
    return tot;
}
void update(int &rt, int pos, int cl, int cr)
{
    if (!rt)
    {
        rt = newnode();
    }
    tree[rt].cnt++;
    tree[rt].sum += bs[pos];
    if (cl == cr)
    {
        return;
    }
    int mid = (cl + cr) >> 1;
    if (pos <= mid)
        update(tree[rt].l, pos, cl, mid);
    else
        update(tree[rt].r, pos, mid + 1, cr);
}
void merge(int &rt1, int &rt2, int cl, int cr)
{
    if (!rt1 || !rt2)
    {
        rt1 |= rt2;
        return;
    }
    tree[rt1].cnt += tree[rt2].cnt;
    tree[rt1].sum += tree[rt2].sum;
    if (cl == cr)
    {
        return;
    }
    int mid = (cl + cr) >> 1;
    merge(tree[rt1].l, tree[rt2].l, cl, mid);
    merge(tree[rt1].r, tree[rt2].r, mid + 1, cr);
    return;
}
tuple<i64, i64, int> querykth(int rt1, int rt2, int k, int cl, int cr)
{
    if (cl == cr)
    {
        return {bs[cl], tree[rt2].sum - tree[rt1].sum, tree[rt2].cnt - tree[rt1].cnt};
    }
    int mid = (cl + cr) >> 1;
    int nowcnt = tree[tree[rt2].l].cnt - tree[tree[rt1].l].cnt;
    int nowsum = tree[tree[rt2].l].sum - tree[tree[rt1].l].sum;
    if (nowcnt >= k)
        return querykth(tree[rt1].l, tree[rt2].l, k, cl, mid);
    auto [kthrl, rsum, rcnt] = querykth(tree[rt1].r, tree[rt2].r, k - nowcnt, mid + 1, cr);
    return {kthrl, nowsum + rsum, nowcnt + rcnt};
}
int sz = 0;
void solve()
{
    int n, K;
    cin >> n >> K;
    init(n);
    for (int i = 1; i <= n; i++)
    {
        cin >> a[i];
        a[i] -= i;
        bs[i] = a[i];
    }
    sort(bs + 1, bs + 1 + n);
    sz = unique(bs + 1, bs + 1 + n) - bs - 1;
    for (int i = 1; i <= n; i++)
    {
        a[i] = lower_bound(bs + 1, bs + 1 + sz, a[i]) - bs;
    }

    /*cout << sz << endl;
    for (int i = 1; i <= n; i++)
        cout << a[i] << " ";
    cout << endl;
    system("pause");*/

    for (int i = 1; i <= n; i++)
    {
        update(root[i], a[i], 1, sz);
        merge(root[i], root[i - 1], 1, sz);
        // 主席树初始化
    }

    int ans = 0;

    for (int l = 1; l <= n; l++)
    {
        int L = l, R = n + 1;

        auto check = [&](int r) -> bool
        {
            int len = r - l + 1;
            int kth = (len + 1) / 2;
            i64 res = 0;
            auto [realkth, prevsum, prevcnt] = querykth(root[l - 1], root[r], kth, 1, sz);
            auto [_, sumall, cntall] = querykth(root[l - 1], root[r], len, 1, sz);

            // cout << "realkth: " << realkth << " prevsum: " << prevsum << " prevcnt: " << prevcnt << endl;
            // cout << "realkth: " << _ << " allsum: " << prevsum << " allcnt: " << prevcnt << endl;

            res = (sumall - prevsum) - (cntall - prevcnt) * realkth + prevcnt * realkth - prevsum;
            // cout << "res: " << res << endl;
            return res <= K;
        };

        while (R - L > 1)
        {
            int mids = (L + R) >> 1;

            // cout << "Check l:" << l << " r:" << mids << endl;

            /*for (int i = l; i <= mids; i++)
            {
                cout << bs[a[i]] << " ";
            }
            cout << endl;*/

            if (check(mids))
            {
                // cout << "ok" << endl;
                // system("pause");
                L = mids;
            }
            else
            {
                // cout << "fatal" << endl;
                // system("pause");
                R = mids;
            }
        }
        ans = max(ans, L - l + 1);
    }
    cout << ans << endl;
}
signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);
    int t;
    cin >> t;
    while (t--)
        solve();
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

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

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: 0
Accepted
time: 815ms
memory: 68004kb

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
7
8
7
7
1
7
6
2
4
3
1
6
7
7
3
4
3
9
3
8
6
6
3
1
6
3
1
2
4
6
4
6
4
1
4
7
1
6
3
5
6
6
1
7
5
3
1
6
4
5
3
2
2
6
2
3
10
1
4
3
2
4
5
1
7
5
5
5
8
5
3
6
3
5
5
8
5
4
5
2
1
5
2
3
3
4
8
1
3
1
2
2
8
3
1
6
8
1
8
4
5
6
6
8
4
8
3
2
8
4
5
6
2
6
2
4
1
5
4
5
3
2
4
1
2
1
4
5
8
3
7
3
3
3...

result:

ok 11102 lines

Test #3:

score: -100
Time Limit Exceeded

input:

1
500000 17244641009859
54748096 75475634 204928248 276927808 84875072 103158867 27937890 322595515 186026685 45468307 69240390 139887597 188586447 373764525 121365644 310156469 185188306 60350786 211308832 384695957 370562147 208427221 35937909 267590963 126478310 275357775 55361535 335993561 36696...

output:


result: