首页 > 其他 > 详细

P4062 [Code+#1]Yazid 的新生舞会

时间:2019-10-08 09:08:12      阅读:142      评论:0      收藏:0      [点我收藏+]

思路:分治

提交:2次

错因:数组开小

题解:

我们枚举一下众数\(x\)
\(s[n]=\sum_{i=1}^n [a[i]==x]\)
那么对于区间\((l,r]\),有\(s[r]-s[l]>\frac{r-l}{2}\)
\(2*s[r]-r>2*s[l]-l\)
考虑分治,我们求出所有过中点的区间\([l,r]\)的贡献。
如何求呢?首先观察一个性质,两个子区间的众数至少有1个是大区间的众数,反之亦然。
那么我们先求出子区间中的众数,作为大区间的可行众数。然后我们枚举每个可行众数,依次求出符合\(2*s[r]-r>2*s[l]-l\)的区间。
还有一个特别神的数组做法,可是我不会

代码:

#include<bits/stdc++.h>
#define ll long long
#define R register int
using namespace std;
namespace Luitaryi {
inline int g() { R x=0,f=1;
    register char ch; while(!isdigit(ch=getchar())) f=ch=='-'?-1:f;
    do x=x*10+(ch^48); while(isdigit(ch=getchar())); return x*f;
} const int N=500010;
int n,a[N],pos[N],num[N],cnt[N*2]; ll ans;
inline void solve(int l,int r) {
    if(l==r) {++ans; return ;} R md=l+r>>1,tot=0;
    solve(l,md),solve(md+1,r);
    for(R i=md;i>=l;--i) if(++cnt[a[i]]>(md-i+1)/2)
        if(!pos[a[i]]) num[pos[a[i]]=++tot]=a[i];//左边的众数 
    for(R i=l;i<=r;++i) cnt[a[i]]=0;
    for(R i=md+1;i<=r;++i) if(++cnt[a[i]]>(i-md)/2) 
        if(!pos[a[i]]) num[pos[a[i]]=++tot]=a[i];//右边的众数 
    for(R i=l;i<=r;++i) pos[a[i]]=0;
    for(R i=l;i<=r;++i) cnt[a[i]]=0;
    for(R i=1;i<=tot;++i) {//枚举可能的众数 
        R sum=r-l+1,mx=sum,mn=sum;//mn,mx记录桶的上下界,sum作为初值相当于偏移量 
        cnt[sum]=1; for(R j=l;j<md;++j) {//先处理出左边的桶
            a[j]==num[i]?++sum:--sum;
            mx=max(mx,sum),mn=min(mn,sum); ++cnt[sum];
        } a[md]==num[i]?++sum:--sum; 
        for(R j=mn;j<=mx;++j) cnt[j]+=cnt[j-1];//前缀和。
        for(R j=md+1;j<=r;++j) {//处理右边
            a[j]==num[i]?++sum:--sum;
            ans+=cnt[min(mx,sum-1)];    //求出顺序对个数
        } for(R j=mn;j<=mx;++j) cnt[j]=0;
    } 
}

inline void main() {
    n=g(); g(); for(R i=1;i<=n;++i) a[i]=g();
    solve(1,n); printf("%lld\n",ans);
}
} signed main() {Luitaryi::main(); return 0;}

2019.10.08
38

P4062 [Code+#1]Yazid 的新生舞会

原文:https://www.cnblogs.com/Jackpei/p/11633386.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!