Jzoj5422 天才绅士少女助手克里斯蒂娜

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define LL long long
#define N 1000010
#define M 20170927
using namespace std;
int n,m;
LL x
,y
,sx
,sy
,sxy
;
inline LL sqr(LL x){ return x*x%M; }
inline void ad(LL& x,LL y){ x=(x+y)%M; }
inline void add(LL* s,LL x,LL k){ for(;x<=n;x+=x&-x) ad(s[x],k); }
inline LL sum(LL* s,LL x,LL S=0){ for(;x;x^=x&-x) ad(S,s[x]); return S; }
int main(){
freopen("kurisu.in","r",stdin);
freopen("kurisu.out","w",stdout);
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i){
scanf("%lld%lld",x+i,y+i);
add(sx,i,sqr(x[i]));
add(sy,i,sqr(y[i]));
add(sxy,i,x[i]*y[i]);
}
for(LL o,l,r;m--;){
scanf("%lld%lld%lld",&o,&l,&r);
if(o==1){
scanf("%lld",&o);
add(sx,l,sqr(r)-sqr(x[l]));
add(sy,l,sqr(o)-sqr(y[l]));
add(sxy,l,o*r-x[l]*y[l]);
x[l]=r; y[l]=o;
} else {
LL s1=sum(sx,r)-sum(sx,l-1);
LL s2=sum(sy,r)-sum(sy,l-1);
LL s3=sum(sxy,r)-sum(sxy,l-1);
printf("%lld\n",(M+(s1*s2%M-sqr(s3))%M)%M);
}
}
}