[10.02]T1求和

题意

现有二维矩阵$d_{i,j}=i+j-1$，给两个坐标$(x1,y1),(x2,y2)$
求以两坐标为顶点的$d_{i,j}$之和，答案大无边，因此给出m，输出对m取模的答案
$1\le x2,y2\le 10^{18}，1\le x1\le x2,1\le y1\le y2$

思路

看到这道题我直接敲暴力打表了，

这个表看着非常有规律，但显然没有直接的公式能求出任意两点为顶点的矩阵的和，
所以我们可以先转化成对于以$(1,1),(x,y)$为顶点的矩阵和，
可以看到的是， 每一行都是一个等差数列。
第一行是首相$a_1=1$公差$d=1$的等差数列，这一行的求和就是$$\frac{(1+x)\ast x}{2}$$
第二行每一项都比第一行大了1，
那么第二行的求和就是$$\frac{(2+x+1)\ast x}{2}$$
然后我们把每个1都提出来就是$$\frac{(1+x)\ast x}{2}+x$$
第三行要每个都提一个2，第四行要提一个3，然后这也是一个等差数列，求和就是$$\frac{(y-1)\ast x\ast y}{2}$$
然后除了提出来的数，还有y个第一行的求和，就是$$\frac{(1+x)\ast x\ast y}{2}$$
两个加一起就是$$\frac{(x+y)\ast x\ast y}{2}$$

我们会求这个以后，就考虑对要求的矩阵进行容斥求和，先求出$(1,1),(x2,y2)$的和，然后减去两边多出来的矩阵$(1,1)(x2,y1-1)$和$(1,1)(x1-1,y2)$，再加上多减掉的矩阵$(1,1),(x1-1,y1-1)$

于是我们就可以得到70分，

数据大还真不是吹的，裸乘也会爆$2^{64}$,然后我们就用快（龟）速乘,防止乘法溢出就行了

#include<bits/stdc++.h>

#define MAXN 100010
#define ll long long 
#define rg register
#define INF 0x3f3f3f3f
#define gc() getchar()

using namespace std;

template <class T>
inline void read(T &s){
	T w = 1, ch = gc(); s = 0;
	while(!isdigit(ch)){if(ch == '-') w = -1; ch = gc();}
	while(isdigit(ch)){ s = s * 10 + ch - '0'; ch = gc();}
	s *= w; return ;
}

ll x1, yy, x2, y2, m;

ll mul(ll n, ll m, ll p){
	ll a = n, ans = 0;
	for(; m; m >>= 1, a = (a << 1) % p){
		if(m & 1) ans = (ans + a) % p;
	}
	return ans;
}

ll work(ll x, ll y, ll m){
	ll sum = x + y;
	if(~x & 1) x >>= 1;
	else if(~y & 1) y >>= 1;
	else sum >>= 1;
	return mul(mul(x, y, m), sum, m);
}

int main()
{
	//freopen("sum.in", "r", stdin);
	//freopen("sum.out", "w", stdout);
	
	read(x1), read(yy), read(x2), read(y2), read(m);
	ll first = work(x2, y2, m), sum = work(x1 - 1, yy- 1, m);
	ll left = work(x1 - 1, y2, m), right = work(x2, yy - 1, m);
	
	printf("%lld\n", (first + sum - left - right + 2 * m) % m);
	return 0;
}