#include <bits/stdc++.h>
#ifdef xay5421
#define D(...) fprintf(stderr,__VA_ARGS__)
#else
#define D(...) ((void)0)
#endif
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
using namespace std;

using LL=long long;
const long double PI=acos(-1);
const int INF=0X3F3F3F3F;
mt19937_64 rng(0);

const int N=10005;

int n;
double A,B;

LL dp[N];
vector<int>e[N];
int X[N][2],Y[N][2];
double mn[N],mx[N];

int dfn[N],idx;
void dfs(int u){
	dfn[++idx]=u;
	for(auto&x:e[u]){
		dfs(x);
	}
}

int main(){
#ifdef xay5421
	freopen("a.in","r",stdin);
#endif

	cin>>n;
	rep(i,1,n){
		int x;
		cin>>x;
		if(x){
			e[i].resize(x);
			rep(j,0,x-1){
				cin>>e[i][j];
			}
		}else{
			cin>>X[i][0]>>Y[i][0];
			X[i][1]=X[i][0];
			Y[i][1]=Y[i][0];
		}
	}
	
	dfs(1);
	
	LL ans=0;
	
	const int T=10000;
	rep(tc,1,T){
		A=cos(PI*2.*tc/T);
		B=sin(PI*2.*tc/T);
		/*A=(rng()%int(2e9))-1e9;
		B=(rng()%int(2e9))-1e9;*/
		per(_,n,1){
			int u=dfn[_];
			if(e[u].empty()){
				mn[u]=mx[u]=A*X[u][0]+B*Y[u][0];
			}else{
				double t=-1e18;
				int id=-1;
				double t_=1e18;
				int id_=-1;
				mx[u]=0;
				mn[u]=0;
				
				X[u][0]=X[u][1]=Y[u][0]=Y[u][1]=0;
				for(auto&x:e[u]){
					if(mn[x]+mx[x]>t){
						t=mn[x]+mx[x];
						id=x;
					}
					if(mn[x]+mx[x]<t_){
						t_=mn[x]+mx[x];
						id_=x;
					}
					mx[u]-=mn[x];
					mn[u]-=mx[x];
					X[u][0]-=X[x][1];
					Y[u][0]-=Y[x][1];
					
					X[u][1]-=X[x][1];
					Y[u][1]-=Y[x][1];
				}
				mn[u]+=t_;
				mx[u]+=t;
				
				X[u][0]+=X[id][1]+X[id][0];
				Y[u][0]+=Y[id][1]+Y[id][0];
				
				X[u][1]+=X[id_][1]+X[id_][0];
				Y[u][1]+=Y[id_][1]+Y[id_][0];
				
			}
		}
		ans=max(ans,LL(X[1][0])*X[1][0]+LL(Y[1][0])*Y[1][0]);
		ans=max(ans,LL(X[1][1])*X[1][1]+LL(Y[1][1])*Y[1][1]);
	}
	
	printf("%lld\n",ans);
	return 0;
}