#include <bits/stdc++.h>
#define ll long long int
using
namespace
std;
void
cumm(vector<ll>& x, vector<ll>& y,
vector<ll>& cummx, vector<ll>& cummy,
vector<ll>& cummx2, vector<ll>& cummy2, ll n)
{
for
(
int
i = 1; i <= n; i++) {
cummx[i] = cummx[i - 1] + x[i];
cummy[i] = cummy[i - 1] + y[i];
cummx2[i] = cummx2[i - 1] + x[i] * x[i];
cummy2[i] = cummy2[i - 1] + y[i] * y[i];
}
}
int
calHammeredDistance(
int
n, vector<ll>& x, vector<ll>& y)
{
vector<ll> cummx(n + 1, 0), cummy(n + 1, 0);
vector<ll> cummx2(n + 1, 0), cummy2(n + 1, 0);
cumm(x, y, cummx, cummy, cummx2, cummy2, n);
ll hdx = 0, hdy = 0;
for
(
int
i = 1; i <= n; i++) {
hdx += (i - 1) * x[i] * x[i] + cummx2[i - 1]
- 2 * x[i] * cummx[i - 1];
hdy += (i - 1) * y[i] * y[i] + cummy2[i - 1]
- 2 * y[i] * cummy[i - 1];
}
ll total = hdx + hdy;
return
total;
}
int
main()
{
int
n = 3;
vector<ll> x = {0, 0, 1, 0};
vector<ll> y = {1, 0, 0, 0};
cout << calHammeredDistance(n, x, y);
return
0;
}