#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
typedef long long ll ;
#define sz(x) (ll)(x.size())
const ll MOD = (ll)(1e9),MX=(int)(1e5+5);
typedef pair<ll,ll> ii;
int par[MX];
ll ans,cnt[MX],comp,cl[MX];
vector<pair<ll,ii> > vec[MX];
vector<ii> vecD[2];
ii solve(ll node,ll p){
ii h={0,cnt[node]};
for(auto &it:vec[node])
if(it.fi!=p){
if(it.se.fi==-1)
it.se=solve(it.fi,node);
h.fi+=it.se.fi+it.se.se;
h.se+=it.se.se;
}
return h;
}
ll findSetXY(ll i){
if(par[i]==i)return i;
return par[i]=findSetXY(par[i]);
}
bool UnionXY(ll a,ll b){
ll x=findSetXY(a),y=findSetXY(b);
if(x!=y){par[x]=y;return 1;}
return 0;
}
void work(ll idx){
comp=0;
for (int i = 0; i < sz(vecD[idx]); ++i)
par[i]=i,cnt[i]=0,vec[i].clear();
for (int i = 0; i < sz(vecD[idx]);)
{
ll j=i;
while(j<sz(vecD[idx]) && vecD[idx][j].fi==vecD[idx][i].fi && (vecD[idx][i].se+(j-i))==vecD[idx][j].se){
cl[j]=comp;
ll k=lower_bound(vecD[idx].begin(),vecD[idx].end(),make_pair(vecD[idx][j].fi-1,vecD[idx][j].se))-vecD[idx].begin();
if(k!=sz(vecD[idx]) && vecD[idx][k]==make_pair(vecD[idx][j].fi-1,vecD[idx][j].se)){
ll a=cl[j],b=cl[k];
if(UnionXY(a,b)){vec[a].push_back({b,{-1,-1}}),vec[b].push_back({a,{-1,-1}});}
}
cnt[cl[j]]++,j++;
}
comp++,i=j;
}
for (int i = 0; i < comp; ++i)
ans+=solve(i,i).fi*cnt[i];
}
int DistanceSum(int N, int *X, int *Y) {
vecD[0].resize(N),vecD[1].resize(N);
for (ll i = 0; i < N; ++i)
vecD[0][i]={X[i],Y[i]},vecD[1][i]={Y[i],X[i]};
for (ll i = 0; i < 2; ++i)
sort(vecD[i].begin(),vecD[i].end());
work(0),work(1);
return (ans>>1LL)%MOD;
}
