import sys
input = sys.stdin.readline
t=int(input())
mod=998244353
for _ in range(t):
  n=int(input())
  xy=[]
  for i in range(n):
    x,y=map(int,input().split())
    xy.append((x,y))
  xy.reverse()
  tmp1=0
  tmp2=0
  tmp3=0
  tmp4=0
  tmp5=0
  ans=0
  for i in range(n):
    x,y=xy[i]
    val=x**2+y**2
    inv=pow(val,mod-2,mod)
    ans+=val*tmp2+inv*tmp1
    co=x*x*inv*tmp3+y*y*inv*tmp4+2*x*y*inv*tmp5
    ans-=4*co
    ans%=mod
    tmp1+=val
    tmp1%=mod
    tmp2+=inv
    tmp2%=mod
    tmp3+=x*x*inv
    tmp3%=mod
    tmp4+=y*y*inv
    tmp4%=mod
    tmp5+=x*y*inv
    tmp5%=mod
  ans+=((n)*(n-1))
  print(ans%mod)