AtCoder 4351前缀和+二分+求顺序对数

题意:

给一个长为N的序列,求出所有子区间中位数组成的新序列的中位数

思路:
首先我们需要找到这个题的可二分性:

  • 我们可以知道中位数的性质:
    在序列中,比中位数小的个数有差不多N/2个,呢么比中位数大的差不多也有N/2个,由此我们找到了可二分性,所以我们就可以二分答案
  • 得到序列中所有中位数比mid大的区间个数:
    ( 1 ) (1) 首先利用前缀和处理数组,将不小于mid的赋为1,其余赋为-1,这样我们就可以得到原序列中的前缀区间的值
    ( 2 ) (2) 利用树状数组求前缀和位置的所有顺序对数,这样就可以求出中位数大于mid的区间个数

参考代码:

#include <cstdio>
#include <algorithm>
#include <vector>
#include <iostream>
#include <map>
#include <queue>
#include <cstring>
#include <string>
#include <cstdlib>
#include <cmath>
#include <set>
using namespace std;
const int inf = (1 << 30);
typedef long long ll;
const ll maxn = 1e5 + 5;
vector<ll> vec;
ll n, sum[maxn * 10], c[maxn * 10], a[maxn];
ll lowbit(ll x)
{
    return x & (-x);
}
void modify(ll i)
{
    while (i <= 2 * maxn)
    {
        ++c[i];
        i += lowbit(i);
    }
}
ll getsum(ll i)
{
    ll res = 0;
    while (i > 0)
    {
        res += c[i];
        i -= lowbit(i);
    }
    return res;
}
bool check(ll mid)
{
    memset(c, 0, sizeof c);
    for (ll i = 1; i <= n; i++)
        sum[i] = sum[i - 1] + (a[i] >= mid ? 1 : -1);
    ll res = 0;
    for (ll i = 0; i <= n; i++)
    {
        res += getsum(sum[i] + maxn);
        modify(sum[i] + maxn);
    }
    if (res >= ((n + 1) * n / 4))
        return true;
    return false;
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    ll l = 0, r = 0;
    cin >> n;
    for (ll i = 1; i <= n; i++)
        cin >> a[i], r = max(r, a[i]);
    ll ans = 0;
    while (l <= r)
    {
        ll mid = l + r >> 1;
        if (check(mid))
            l = mid + 1, ans = mid;
        else
            r = mid - 1;
    }
    cout << ans << endl;
}

猜你喜欢

转载自blog.csdn.net/yangzijiangac/article/details/105929431