【NOIP2017提高A组集训10.25】天才绅士少女助手克里斯蒂娜

题面

Description

Input

第一行两个整数n;m 表示电子个数和询问个数.
接下来n 行, 每行两个整数x; y 表示vi.
接下来m 行, 每行形如1 p x y 或2 l r, 分别表示两种操作.

Output

对于每个操作2, 输出一行一个整数表示飘升系数对20170927 取模的值.

Sample Input

9 5
13052925 5757314
9968857 11135327
13860145 3869873
6912189 3461377
2911603 7061332
6334922 7708411
5505379 5915686
6806727 588727
7603043 15687404
2 1 6
1 7 2602783 18398476
1 8 8636316 19923037
2 2 7
2 2 4

Sample Output

18529202
963126
19167545

Data Constraint

思路

$\sum |v_i\ast v_j{|}^2=\sum x^2\ast \sum y^2-(\sum xy{)}^2$
用树状数组维护即可。

Code

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#define p 20170927
int n,m;
long long c[1000005][3],a[1000005],b[1000005];
int lowbit(int x)
{
    
    
	return x&(x^(x-1));
}
void change(int k,int delta,int q)
{
    
    
	while(k<=n)
	{
    
    
		c[k][q]+=delta;
		c[k][q]%=p;
		k+=lowbit(k);
	}
}
long long get(int k,int q)
{
    
    
	long long s=0;
	while(k>0)
	{
    
    
		s+=c[k][q];
		s%=p;
		k-=lowbit(k);
	}
	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++)
	{
    
    
		long long x,y;
		scanf("%lld%lld",&x,&y);
		a[i]=x,b[i]=y;
		change(i,x*x%p,0);
		change(i,y*y%p,1);
		change(i,x*y%p,2);
	}
	for(int i=1;i<=m;i++)
	{
    
    
		int bz;
		scanf("%d",&bz);
		if(bz==1)
		{
    
    
			int k;
			long long x,y;
			scanf("%d%lld%lld",&k,&x,&y);
			change(k,(x*x%p-a[k]*a[k]%p)%p,0);
			change(k,(y*y%p-b[k]*b[k]%p)%p,1);
			change(k,(x*y%p-a[k]*b[k]%p)%p,2);
			a[k]=x,b[k]=y;
		}
		else
		{
    
    
			int l,r;
			scanf("%d%d",&l,&r);
			long long s1,s2,s3;
			s1=(get(r,0)-get(l-1,0)+p)%p;
			s2=(get(r,1)-get(l-1,1)+p)%p;
			s3=(get(r,2)-get(l-1,2)+p)%p;
			printf("%lld\n",((s1*s2%p)-(s3*s3)%p+5*p)%p);
		}
	}
}