#include<bits/stdc++.h>
using namespace std;
using ll=long long;
template<typename T>
ostream& operator << (ostream &out,const vector<T>&x){
	if(x.empty())return out<<"[]";
	out<<'['<<x[0];
	for(int len=x.size(),i=1;i<len;i++)out<<','<<x[i];
	return out<<']';
}
template<typename T>
vector<T> ary(const T *a,int l,int r){
	return vector<T>{a+l,a+1+r};
}
template<typename T>
void debug(T x){
	cerr<<x<<'\n';
}
template<typename T,typename ...S>
void debug(T x,S ...y){
	cerr<<x<<' ',debug(y...);
}
using vec=complex<ll>;
ll dot(const vec &a,const vec &b){
	return real(conj(a)*b);
}
const int N=1e4+10;
int n;
vector<int>to[N];
vec cur,a[N];
bool cmp(vec a,vec b){
	return dot(a,cur)<dot(b,cur);
}
pair<vec,vec> dfs(int u){
	if(to[u].empty())return {a[u],a[u]};
	vector<vec>p;
	vec sl,sr;
	for(int v:to[u]){
		vec l,r;
		tie(l,r)=dfs(v);
		p.push_back(l+r);
		sl+=l,sr+=r;
	}
	vec l=p.back()-sr,r=p.back()-sl;
	for(vec x:p)l=min(l,x-sr,cmp),r=min(r,x-sl,cmp);
	return {l,r};
}
ll ans;
pair<vec,vec> solve(vec cur){
	::cur=cur;
	auto res=dfs(1);
	ans=max({ans,norm(res.first),norm(res.second)});
	return res;
}
void find(vec l,vec r){
	vec p=(r-l)*vec(0,1),mid=solve(p).second;
	if(dot(p,mid)>max(dot(p,l),dot(p,r))){
		find(l,mid),find(mid,r);
	}
}
int main(){
	scanf("%d",&n);
	for(int i=1,k,x,y;i<=n;i++){
		scanf("%d",&k);
		if(!k){
			scanf("%d%d",&x,&y);
			a[i]=vec(x,y);
		}else{
			for(;k--;){
				scanf("%d",&x);
				to[i].push_back(x);
			}
		}
	}
	vec l,r;
	tie(l,r)=solve(vec(0,1));
	find(l,r);
	cout<<ans<<endl;
	return 0;
}