CoderForces-1005E2 Median on Segments

传送门：CoderForces-1005E2

思路分析

计算中位数为m的区间个数，可以通过计算大于等于m的答案-大于等于m+1的答案
怎么计算中位数大于等于x的区间个数呢，如果一个区间的中位数是大于等于x的，那么这个区间内大于等于x-小于x一定大于等于1，所以可以利用前缀和差分的性质计算出答案。

样例输入

15 2
1 2 3 1 2 3 1 2 3 1 2 3 1 2 3

样例输出

97

AC代码

#include <bits/stdc++.h>
#define ls x<<1
#define rs x<<1|1
#define fi first
#define se second
#define ll long long
#define pb push_back
#define mp make_pair
#define fun function
#define sz(x) (x).size()
#define lowbit(x) (x)&(-x)
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))
#define int long long
namespace FastIO {
#define BUF_SIZE 100000
#define OUT_SIZE 100000
    bool IOerror=0;
    inline char nc() {
        static char buf[BUF_SIZE],*p1=buf+BUF_SIZE,*pend=buf+BUF_SIZE;
        if(p1==pend) {
            p1=buf;
            pend=buf+fread(buf,1,BUF_SIZE,stdin);
            if(pend==p1) {
                IOerror=1;
                return -1;
            }
        }
        return *p1++;
    }
    inline bool blank(char ch) {
        return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';
    }
    template<class T> inline bool read(T &x) {
        bool sign=0;
        char ch=nc();
        x=0;
        for(; blank(ch); ch=nc());
        if(IOerror)return false;
        if(ch=='-')sign=1,ch=nc();
        for(; ch>='0'&&ch<='9'; ch=nc())x=x*10+ch-'0';
        if(sign)x=-x;
        return true;
    }
    template<class T,class... U>bool read(T& h,U&... t) {
        return read(h)&&read(t...);
    }
#undef OUT_SIZE
#undef BUF_SIZE
};
using namespace std;
using namespace FastIO;

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

const int INF = 0x3f3f3f3f;
const int N = 2e5+10;

int a[N],b[N<<1],n,m; 

void add(int x){
    for(int i=x;i<N*2;i+=lowbit(i)) b[i]++;
}

int get(int x){
    int ans=0;
    for(int i=x;i;i-=lowbit(i)) ans+=b[i];
    return ans;
} 

ll solve(int x){
    mem(b,0);
    int sum=n+1;
    ll ans=0;
    add(sum);
    for(int i=1;i<=n;i++){
        if(a[i]>=x) sum++;
        else sum--;
        ans+=get(sum-1);
        add(sum);
    }
    return ans;
}

signed main() {

#ifdef xiaofan
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
#endif

    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>a[i];
    cout<<solve(m)-solve(m+1)<<endl;



    return 0;
}