2018.10.22【COGS2633】数列操作e（线段树）传送门解析：

传送门

---

解析：

输出又一次阻挡了我一A的步伐。。。

思路：

线段树维护区间加二次函数，区间求和。

题目中要求的加操作是在第$i$个位置加上(i−l+1)2×v(i-l+1)^2times v(i−l+1)2×v，显然我们需要展开成二次函数一般形式，这样才能够令它满足加法性质，便于合并标记。

展开后以$i$为主元就是i2∗v+i∗(2∗v∗(1−pos))+(pos−1)2∗vi^2*v+i*(2*v*(1-pos))+(pos-1)^2*vi2∗v+i∗(2∗v∗(1−pos))+(pos−1)2∗v，显然这三个多项式的合并都可以直接加法，并且我们预处理平方前缀和可以做到$O(1)$计算出者三个的区间和，具体实现就在代码里面。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const
#define ull unsigned ll

inline int getint(){
	re int num;
	re char c;
	while(!isdigit(c=gc()));num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return num;
}

cs int N=100005;

ull sum[N<<2],A[N<<2],B[N<<2],C[N<<2];
ull square[N];

inline void pushnow(int k,int l,int r,cs ull &a,cs ull &b,cs ull &c){
	sum[k]+=(r-l+1)*a+(1ll*r-l+1)*(l+r)/2*b+(square[r]-square[l-1])*c;
	A[k]+=a;B[k]+=b;C[k]+=c;
}

inline void pushdown(int k,int l,int r){
	if(A[k]||B[k]||C[k]){
		int mid=(l+r)>>1;
		pushnow(k<<1,l,mid,A[k],B[k],C[k]);
		pushnow(k<<1|1,mid+1,r,A[k],B[k],C[k]);
		A[k]=B[k]=C[k]=0;
	}
}

inline void pushup(int k){
	sum[k]=sum[k<<1]+sum[k<<1|1];
}

void modify(int k,int l,int r,cs int &ql,cs int &qr,cs ull &pos,cs ull &v){
	if(ql<=l&&r<=qr)return pushnow(k,l,r,(pos-1)*(pos-1)*v,v*2*(1-pos),v);
	pushdown(k,l,r);
	int mid=(l+r)>>1;
	if(ql<=mid)modify(k<<1,l,mid,ql,qr,pos,v);
	if(qr>mid)modify(k<<1|1,mid+1,r,ql,qr,pos,v);
	pushup(k);
}

ull query(int k,int l,int r,cs int &ql,cs int &qr){
	if(ql<=l&&r<=qr)return sum[k];
	int mid=(l+r)>>1;
	pushdown(k,l,r);
	if(ql>mid)return query(k<<1|1,mid+1,r,ql,qr);
	if(qr<=mid)return query(k<<1,l,mid,ql,qr);
	return query(k<<1,l,mid,ql,qr)+query(k<<1|1,mid+1,r,ql,qr);
}

ull ans;
int n,m;
signed main(){
	freopen("rneaty.in","r",stdin);
	freopen("rneaty.out","w",stdout);
	n=getint();
	m=getint();
	for(int re i=1;i<=n;++i)square[i]=square[i-1]+(ull)i*i;
	while(m--){
		int op=getint(),l=getint(),r=getint();
		switch(op){
			case 1:{
				ull x=getint();
				modify(1,1,n,l,r,l,x);
				break;
			}
			case 2:{
				ans^=query(1,1,n,l,r);
				break;
			}
		}
	}
	printf("%llu",ans);
	return 0;
}