20200428 T1 Arithmetic【Δ线段树(每个位置加的值不一样的线段树)】

题目描述：

题目分析：

（标题取成这样因为实在不知道叫什么好。。）

不难发现我们更新线段树时需要根据每个左端点对应的Max和Min进行加减，而不是像普通的线段树那样区间整体加上一个相同的值。
考虑我们的线段树需要支持什么操作，记每个位置的变化量为$\Delta$（即当区间整体+1时这个位置需要加的量）：

• 给区间$[l,r]$每个位置的$\Delta$加上一个值 $x$ (可正可负)
• 令区间$[l,r]$每个位置的答案加上其对应的$\Delta \ast k$，（此题中$k=1$，$[l,r]$可以只用$[1,n]$）
• 求区间$[l,r]$的答案的和

为此，我们的线段树需要记录这几个量：

• $su{m}_i$ ：表示该区间答案的和。
• $s{d}_i$：表示该区间$\Delta$的和。
• $k_i$：表示答案整体加$\Delta$次数的懒标记。
• $ta{g}_i$：表示修改$\Delta$的懒标记。
• $ad{d}_i$：（这是一个比较特别的标记，不同的人可能有不同的写法，在后面介绍）
• $le{n}_i$：该区间的长度，方便pushdown。

$pushup$就这样写：

inline void upd(int i){
	t[i].sum=t[lc].sum+t[rc].sum;
	t[i].sd=t[lc].sd+t[rc].sd;
}

问题主要在于$pushdown$，我们需要考虑标记优先级的问题，这里假设定义$k$的优先级高于$tag$（即先加答案，再修改$\Delta$）。

先忽略$add$标记，考虑标记$(k_i,ta{g}_i)$怎么传给$x$，首先肯定有：
$su{m}_x+=k_i\ast s{d}_x$
$k_x+=k_i$
$ta{g}_x+=ta{g}_i$
$s{d}_x+=ta{g}_i\ast le{n}_x$
那么这时候问题就来了！对于$x$的儿子$y$，当$x$ $pushdown$的时候，它本来应该先答案加上$k_x\ast s{d}_y$，再修改$s{d}_y$，然后当$i$下传给$x$再下传给$y$的时候，给答案加上$k_i\ast s{d}_y$，再修改$s{d}_y$，可以现在$i$提前传给了$x$，这就导致了传给$y$的时候，$k_i$本来应该是在$s{d}_y$修改之后再加的，但是由于$k$的优先级高于$tag$，所以$su{m}_y$就加少了！少了多少呢？就是少加给$\Delta$的$k_i\ast ta{g}_x$！这就是$add$所记录的信息。（此处可以稍作思考）

即$ad{d}_i$记录的是当标记$(k_{rt},ta{g}_{rt})$下传给$i$时，将会少加给$i$的儿子的$sum$的关于每个位置的$\Delta$的信息，即$k_{rt}\ast ta{g}_i$（说的不是特别清楚，结合代码理解一下）

重新梳理一下$pushdown$应该怎么写：

inline void pd(int i,int x){
	t[x].sum+=t[i].k*t[x].sd+t[i].add*t[x].len;//add记录的是每个位置对答案的额外贡献
	t[x].add+=t[i].k*t[x].tag+t[i].add;
	t[x].k+=t[i].k;
	t[x].tag+=t[i].tag;
	t[x].sd+=t[i].tag*t[x].len;
}

另外，对于题目中的$-(r-l)$，其实就是对每个$r$将$[L_r,r-1]$的$\Delta -1$，然后令$[L_r,r-1]$的$k+1$

Code：

#include<bits/stdc++.h>
#define maxn 300005
#define LL long long
using namespace std;
char cb[1<<20],*cs,*ct;
#define getc() (cs==ct&&(ct=(cs=cb)+fread(cb,1,1<<20,stdin),cs==ct)?0:*cs++)
template<class T>inline void read(T &a){
	char c;while(!isdigit(c=getc()));
	for(a=c-'0';isdigit(c=getc());a=a*10+c-'0');
}
int n,d,m,a[maxn],A[maxn],An,B[maxn],Bn,last[10000005];
LL ans[maxn];
struct node{
	int l,r,id; bool operator < (const node &p)const{return r<p.r;}
}Q[maxn];
struct Seg{
	int len;
	LL sum,sd,tag,k,add;
	//sd: sum of delta, tag: whole delta of delta, k: times of delta affected. add: when k+=1, additional tag left for sum.
}t[maxn<<2];
#define lc (i<<1)
#define rc (i<<1|1)
inline void upd(int i){
	t[i].sum=t[lc].sum+t[rc].sum;
	t[i].sd=t[lc].sd+t[rc].sd;
}
inline void pd(int i,int x){
	t[x].sum+=t[i].k*t[x].sd+t[i].add*t[x].len;
	t[x].add+=t[i].k*t[x].tag+t[i].add;
	t[x].k+=t[i].k;
	t[x].tag+=t[i].tag;
	t[x].sd+=t[i].tag*t[x].len;
}
inline void down(int i){pd(i,lc),pd(i,rc),t[i].tag=t[i].k=t[i].add=0;}
void build(int i,int l,int r){
	t[i].len=r-l+1;
	if(l==r) return;
	int mid=(l+r)>>1;
	build(lc,l,mid),build(rc,mid+1,r);
}
void mdf(int i,int l,int r,int x,int y,LL d){
	if(x<=l&&r<=y) {t[i].sd+=d*t[i].len,t[i].tag+=d;return;}
	down(i);
	int mid=(l+r)>>1;
	if(x<=mid) mdf(lc,l,mid,x,y,d);
	if(y>mid) mdf(rc,mid+1,r,x,y,d);
	upd(i);
}
LL qry(int i,int l,int r,int x,int y){
	if(x<=l&&r<=y) return t[i].sum;
	down(i);
	int mid=(l+r)>>1; LL ret=0;
	if(x<=mid) ret+=qry(lc,l,mid,x,y);
	if(y>mid) ret+=qry(rc,mid+1,r,x,y);
	return ret;
}
void zero(int i,int l,int r,int x){
	if(l==r) return void(t[i].sd=0);
	down(i);
	int mid=(l+r)>>1;
	x<=mid?zero(lc,l,mid,x):zero(rc,mid+1,r,x);
	upd(i);
}
inline void accumulate(){
	t[1].sum+=t[1].sd,t[1].add+=t[1].tag,t[1].k++;
}
int main()
{
	freopen("arithmetic.in","r",stdin);
	freopen("arithmetic.out","w",stdout);
	read(n),read(d),read(m);
	if(!n||!m) return 0;
	for(int i=1;i<=n;i++) read(a[i]);
	for(int i=1;i<=m;i++) read(Q[i].l),read(Q[i].r),Q[i].id=i;
	sort(Q+1,Q+1+m);
	int L=1; build(1,1,n);
	for(int i=1,j=1;i<=n;i++){
		int pre=L;
		if(a[i]%d==a[i-1]%d) L=max(L,last[a[i]]+1); else L=i;
		last[a[i]]=i;
		while(pre<L) zero(1,1,n,pre++);//let delta be 0.
		
		while(A[An]>=L&&a[i]>a[A[An]]) mdf(1,1,n,max(A[An-1]+1,L),A[An],-a[A[An]]/d),--An;
		mdf(1,1,n,max(A[An]+1,L),i,a[i]/d),A[++An]=i;//max/d
	
		while(B[Bn]>=L&&a[i]<a[B[Bn]]) mdf(1,1,n,max(B[Bn-1]+1,L),B[Bn],a[B[Bn]]/d),--Bn;
		mdf(1,1,n,max(B[Bn]+1,L),i,-a[i]/d),B[++Bn]=i;//-min/d
		
		if(L<i) mdf(1,1,n,L,i-1,-1);//-(r-l)
		accumulate();
		
		for(;j<=m&&Q[j].r<=i;j++) ans[Q[j].id]=qry(1,1,n,Q[j].l,i);
	}
	for(int i=1;i<=m;i++) printf("%lld\n",ans[i]);
}