Description
农夫Byteasar买了一片n亩的土地,他要在这上面种草。
他在每一亩土地上都种植了一种独一无二的草,其中,第i亩土地的草每天会长高a[i]厘米。
Byteasar一共会进行m次收割,其中第i次收割在第d[i]天,并把所有高度大于等于b[i]的部分全部割去。Byteasar想知道,每次收割得到的草的高度总和是多少,你能帮帮他吗?
Input
第一行包含两个正整数n,m(1<=n,m<=500000),分别表示亩数和收割次数。
第二行包含n个正整数,其中第i个数为ai,依次表示每亩种植的草的生长能力。
接下来m行,每行包含两个正整数d[i],bi,依次描述每次收割。
数据保证d[1]< d[2]< …< d[m],并且任何时刻没有任何一亩草的高度超过10^12。
Output
输出m行,每行一个整数,依次回答每次收割能得到的草的高度总和。
Sample Input
4 4
1 2 4 3
1 1
2 2
3 0
4 4
Sample Output
6
6
18
0
HINT
第1天,草的高度分别为1,2,4,3,收割后变为1,1,1,1。
第2天,草的高度分别为2,3,5,4,收割后变为2,2,2,2。
第3天,草的高度分别为3,4,6,5,收割后变为0,0,0,0。
第4天,草的高度分别为1,2,4,3,收割后变为1,2,4,3。
Source
By Claris
注意到所有位置草长度的大小关系总可以看成是不变的,而收割是无视区间下标问题的
那么我们可以先排序然后再构线段树.
每次二分出最小的那个长度大于等于b的位置就行了.
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 500010
#define GET (ch>='0'&&ch<='9')
#define lchild rt<<1,l,mid
#define rchild rt<<1|1,mid+1,r
#define ln rt<<1
#define rn rt<<1|1
#define LL long long
using namespace std;
int n,m,pos;
LL sum[MAXN],a[MAXN],d,b,last;
struct seg
{
int l,r;
LL maxn,sum,flag,t;
seg() { flag=-1; }
}tree[MAXN<<2];
inline void in(LL &x)
{
char ch=getchar();x=0;
while (!GET) ch=getchar();
while (GET) x=x*10+ch-'0',ch=getchar();
}
inline void update(int rt,LL t)
{
tree[rt].maxn+=t*a[tree[rt].r];
tree[rt].sum+=(sum[tree[rt].r]-sum[tree[rt].l-1])*t;
tree[rt].t+=t;
}
inline void push_down(int rt)
{
if (tree[rt].flag!=-1)
{
tree[ln].flag=tree[rn].flag=tree[rt].flag;tree[ln].t=tree[rn].t=0;
tree[ln].maxn=tree[rn].maxn=tree[rt].flag;
tree[ln].sum=(tree[ln].r-tree[ln].l+1)*tree[rt].flag;
tree[rn].sum=(tree[rn].r-tree[rn].l+1)*tree[rt].flag;
tree[rt].flag=-1;
}
if (tree[rt].t) update(ln,tree[rt].t),update(rn,tree[rt].t),tree[rt].t=0;
}
inline void push_up(int rt)
{
tree[rt].sum=tree[ln].sum+tree[rn].sum;
tree[rt].maxn=tree[rn].maxn;
}
int find(int rt)
{
int L=tree[rt].l,R=tree[rt].r;
if (L==R) return L;
push_down(rt);
if (tree[ln].maxn>=b) return find(ln); else return find(rn);
}
LL query(int rt)
{
int L=tree[rt].l,R=tree[rt].r;LL ret=0;
if (R<pos) return 0;
if (L>=pos)
{
ret=tree[rt].sum;tree[rt].sum=(R-L+1)*b;
tree[rt].t=0;tree[rt].maxn=b;tree[rt].flag=b;
return ret;
}
push_down(rt);ret=query(ln)+query(rn);push_up(rt);
return ret;
}
void build(int rt=1,int l=1,int r=n)
{
tree[rt].l=l;tree[rt].r=r;
if (l==r) return;
int mid=(l+r)>>1;build(lchild);build(rchild);
}
int main()
{
scanf("%d%d",&n,&m);build();
for (int i=1;i<=n;i++) in(a[i]);sort(a+1,a+n+1);
for (int i=1;i<=n;i++) sum[i]=sum[i-1]+a[i];
for (;m;m--)
{
in(d);in(b);
update(1,d-last);last=d;
if (tree[1].maxn<b) { puts("0");continue; }
pos=find(1);printf("%lld\n",query(1)-b*(n-pos+1));
}
}