题目背景

zzs和zzy正在被寒假作业折磨，然而他们有答案可以抄啊。

题目描述

他们共有n项寒假作业。zzy给每项寒假作业都定义了一个疲劳值Ai，表示抄这个作业所要花的精力。zzs现在想要知道，有多少组连续的寒假作业的疲劳值的平均值不小于k？

简单地说，给定n个正整数A1,A2,A3,...,An，求出有多少个连续的子序列的平均值不小于k。

输入输出格式

输入格式：

第一行两个正整数，n和k。

第二行到第n+1行，每行一个正整数Ai。

输出格式：

一个非负整数。

输入输出样例

3 2
1
2
3

4

说明

样例解释：共有6个连续的子序列，分别是(1)、(2)、(3)、(1,2)、(2,3)、(1,2,3)，平均值分别为1、2、3、1.5、2.5、2，其中平均值不小于k的共有4个。

对于20%的数据，1<=n<=100；

对于50%的数据，1<=n<=5000；

对于100%的数据，1<=n<=100000；

对于100%的数据，1<=Ai<=10000，1<=k<=10000。

解题思路

　　先来简化一下题意：给定一个数列a[]，求有多少个子序列元素平均值大于等于k。

　　遍历所有子区间复杂度最小O(n^2)，可用前缀和O(1)得到任意区间和。

　　想个办法吧，这类题要么用奇奇怪怪的数据结构(可能有这样的数据结构?我不知道)，要么推公式。

　　那就推公式——

　　首先取出一段区间吧，设1<=i<=j<=n，然后这段区间（a[i]~a[j]）的平均值为

　　　　(a[i]+a[i+1]+…+a[j]) / (j-i+1)>=k

　　　　分母不太好看，乘到右边吧——

　　　　a[i]+a[i+1]+…+a[j] >= k*(j-i+1)

　　　　a[i]+a[i+1]+…+a[j]>= k+…+k //(j-i+1)个k

　　　　接下来有点关键啦

　　　　a[i]+a[i+1]+…+a[j] - k-…-k>=0

　　　　(a[i]-k) + (a[i+1]-k) +……+ (a[j]-k)>=0//哦？那么整齐，有意思

　　　　那我们另设一个数组b[]，使b[i]=a[i]-k吧

　　　　b[i]+b[i+1]+…+b[j]>=0

　　　　哦？b[]的一段区间和大于等于零？

　　记得前面想暴力的时候说过用前缀和能O(1)取得任意区间和吗？给b[]套上区间和吧

　　设s[i]=b[1]+b[2]+b[3]+…+b[i]，特别的，令s[0]=0，那么b[i]+b[i+1]+…+b[j] = s[j]-s[i-1]。

　　　　s[j] - s[i-1]>=0

　　　　移项一下

　　　　s[j]>=s[i-1]，哇，快了。

　　　　因为之前定义过i<=j，所以i-1<j，这是什么？

　　　　是的没错！逆序对（大雾，应该叫顺序对的）！

　　　　i-1<j且s[i-1]<=s[j]，对s数组求逆序对顺序对即可。

　　不会求逆序对的就先掌握这个姿势吧，也挺简单的。

　　最后说一下，我不知道为什么把求逆序对的程序类比着改一下，求出来“顺序对”就是对的，反正套上去就AC了，如果有那个dalao知道，能在评论区为本蒟蒻解惑吗？谢谢了。

-------------------------------------------------分割线--------------------------------------------

2019年03月05日20:06:16：当年我找这个解法是从百度第五页找出来的，当年洛谷上这题也只有十多个人AC。今天回来一看，博客园的相似推荐，一个平衡树，其他全是这个解法……感觉自己该抒发点什么，但啥也说不出来……剧场效应导致前浪被拍死在沙滩上？或是一代代OIer的努力降低了查找资料的难度，也极大增加了获奖的难度？（今年NOI免费了，感觉自己错亿）（心里真乱）

源代码

#include<stdio.h>int n,k;
int a[100010]={0},temp[100010]={0};
long long ans=0;void merge(int s1,int e1,int s2,int e2)
{int s=s1,e=e2,i=s1;while(s1<=e1&&s2<=e2){if(a[s1]<=a[s2]){ans+=e2-s2+1;//求逆序对在"a[s1]>a[s2]"里,求顺序对反过来temp[i++]=a[s1++];}elsetemp[i++]=a[s2++];}while(s1<=e1) temp[i++]=a[s1++];while(s2<=e2) temp[i++]=a[s2++];for(i=s;i<=e;i++) a[i]=temp[i];
}
void msort(int l,int r)
{if(l==r) return;int mid=(l+r)>>1;msort(l,mid);msort(mid+1,r);merge(l,mid,mid+1,r);
}int main()
{scanf("%d%d",&n,&k);for(int i=1,b;i<=n;i++){scanf("%d",&b);b-=k;a[i]=a[i-1]+b;//输入a[i]时就顺带处理出了s[i]，依然存在a[i]里
    }msort(0,n);printf("%lld",ans);return 0;
}