jian

暴力
枚举

#include<cstdio>
using namespace std;
#define ll long long
int ans,n,s[1999999],a[199999],l,r;
int gcd(int x,int y){
    if(!y) return x;
    return gcd(y,x%y);
}
int main(){
    freopen("jian.in","r",stdin);
    freopen("jian.out","w",stdout);
    scanf("%d%d%d",&n,&l,&r);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        s[i]=s[i-1]+a[i];
    }
    for(int i=1;i<=n;i++)
    for(int j=i;j<=n;j++)
    {
         int w=(s[j]-s[i-1]);
         int lll=l*(j-i+1),rr=r*(j-i+1);
        if(w<=rr&&w>=lll) ans++;
    }
    int w=gcd(ans,n*(n+1)/2);
    if(w==ans) printf("1\n");
    else printf("%d/%d",ans/w,n*(n+1)/2/w);
}

正解
逆序对。。。。。。
要求 区间平均值>=l&&<=r 的个数

现在我们来求区间平均值在1~r的个数和1~l(不包括l)的个数 前减后即为所求
以求1~r为例
(a[i]+a[i+1]+……+a[i+k-1])/k<=r
(a[i]+a[i+1]+……+a[i+k-1])<=k*r
(a[i]+a[i+1]+……+a[i+k-1])-k*r<=0
(a[i]-r)+(a[i+1]-r)+……+(a[i+k-1]-r)<=0
令c[i]=a[i]-r得到一个新数组c
即求c数组区间和<=0的个数
令s为数组c的前缀和数组
c[i]+c[i+1]+…+c[i+k-1]<=0
s[i+k-1]-s[i]<=0
s[i+k-1]<=s[i]

i< i+k-1
s[i]>=s[i+k-1]
即求s数组逆序对数
这我们就可以用树状数组或归并排序来维护
树状数组需要离散化,因为有负的情况。。。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
ll n,s[1999999],s2[1999991],a[1999999],l,r,sum[1999999],b[1999999],b2[1999999];ll t1,t2;
ll gcd(ll x,ll y){
    if(!y) return x;
    return gcd(y,x%y);
}
void add(ll x){
    for(ll i=x;i<=n;i+=i&-i)
    sum[i]++;
}
ll ask(ll x){
    ll ans=0;
    for(ll i=x;i>=1;i-=i&-i)
    ans+=sum[i];
    return ans;
}
void add2(ll x){
    for(ll i=x+1;i<=n;i+=i&-i)
    sum[i]++;
}
ll ask2(ll x){
    ll ans=0;
    for(ll i=x;i>=1;i-=i&-i)
    ans+=sum[i];
    return ans;
}
ll ans=0,ans2=0;

int main(){
    freopen("jian.in","r",stdin);
    freopen("jian.out","w",stdout);
    scanf("%lld%lld%lld",&n,&l,&r);
    for(ll i=1;i<=n;i++){
        scanf("%lld",&a[i]);
        s[i]=s[i-1]+a[i]-r;
        s2[i]=s2[i-1]+a[i]-l;
        if(s[i]<=0) ans++;
        if(s2[i]<0) ans2++;
        b[i]=s[i],b2[i]=s2[i];
    }
    sort(b+1,b+n+1);sort(b2+1,b2+n+1);
     t1=unique(b+1,b+n+1)-b-1;
     t2=unique(b2+1,b2+n+1)-b2-1;
    for(ll i=1;i<=n;i++)
    {
        s[i]=lower_bound(b+1,b+t1+1,s[i])-b;

        s2[i]=lower_bound(b2+1,b2+t2+1,s2[i])-b2;
    }
    for(ll i=n;i>=1;i--){ans+=ask(s[i]);
        add(s[i]);

    }
    memset(sum,0,sizeof sum);
    for(ll i=n;i>=1;i--){ans2+=ask2(s2[i]);
        add2(s2[i]);

    }   
    ans=ans-ans2;

    ll w=gcd(ans,n*(n+1)/2);
    if(ans==0) printf("0\n");
    else if(n*(n+1)/2==ans) printf("1\n");
    else printf("%lld/%lld",ans/w,n*(n+1)/2/w);
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值