【数学式子推导】梦熊[csp-noip十连测第六套] T2 距离

感谢 Jayun 的技术支持

原题式子：
$$\sum _{i=1}^n\sum _{j=i+1}^n((x_i-x_j{)}^2+(y_i-y_j{)}^2)$$
先只考虑 $x$。
$$\sum _{i=1}^n\sum _{j=i+1}^n(x_i-x_j{)}^2$$
化简：
$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=i+1}^n(x_i-x_j{)}^2\\ & =\sum _{i=1}^n\sum _{j=i+1}^n(x_i^2+x_j^2-2x_i{x}_j)\\ & =\sum _{i=1}^n\left((n-i)x_i^2+\sum _{j=i+1}^n(x_j^2-2x_i{x}_j)\right)\\ & =\sum _{i=1}^n\left((n-i)x_i^2+(i-1{)}^{[1]}{x}_i^2-2\sum _{j=i+1}^n{x}_i{x}_j\right)\\ & =\sum _{i=1}^n\left((n-1)x_i^2-2x_i\sum _{j=i+1}^n{x}_j\right)\end{array}$$

$(4)[1]:i<j$ 时，每个 $i$ 计算都加一遍 $x_j^2$，总共加 $i-1$ 次。

用前缀和处理后面一个 sigma 即可，即
$$\sum _{i=1}^n\left((n-1)x_i^2-2x_i(su{m}_{x_n}-su{m}_{x_i})\right)$$

$y$ 同理，复杂度 $O(n)$。

开 __int128，附带读写。

#include<bits/stdc++.h>
//#pragma GCC optimize(2)
using namespace std;
typedef long long ll;
ll n;
__int128 x[100002],y[100002],h[100002],ns,ans;
inline ll rd()
{
	ll rt=0;int fs=1;char c=getchar();
	while(c<'0'||c>'9') {if(c=='-') fs=-1;c=getchar();}
	while('0'<=c&&c<='9') {rt=rt*10-48+c;c=getchar();}
	return rt*fs;
}
void print(__int128 pt)
{
	if(pt<0)
	{
		putchar('-');
		pt=-pt;
	}
	if(pt>9) print(pt/10);
	putchar(char(pt%10+48));
}
int main()
{
//	freopen("dis.in","r",stdin);
//	freopen("dis.out","w",stdout);
	n=rd();
	for(int i=1;i<=n;++i)
	{
		x[i]=rd(),y[i]=rd();
		h[i]=x[i]+h[i-1];
	}
	ns=h[n];
	for(int i=1;i<=n;++i)
		ans+=(n-1)*x[i]*x[i]-2*x[i]*(ns-h[i]);
	for(int i=1;i<=n;++i)
		h[i]=y[i]+h[i-1];
	ns=h[n];
	for(int i=1;i<=n;++i)
		ans+=(n-1)*y[i]*y[i]-2*y[i]*(ns-h[i]);
	print(ans);
	return 0;//freopen?
}
/*
*/
//-static-libgcc -std=c++14 -Wall -Wl,--stack=134217728