python3 - <<'PY'
p=[2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71]
N=1048576
n=[0]*N
s=[0]*N
q=[0]*N
d=[0]*N
for m in range(1,N):
 l=m&(-m)
 k=l.bit_length()-1
 u=m^l
 n[m]=n[u]+1
 s[m]=s[u]+p[k]
 q[m]=q[u]+p[k]*p[k]
 d[m]=n[m]*q[m]-s[m]*s[m]
ps=[(2,3),(2,5),(2,7),(2,11),(3,5),(3,7)]
for a,b in ps:
 g=b-a
 mask=(1<p.index(a))|(1<p.index(b))
 found=0
 for m in range(1,N):
  if (m&mask)==mask and m!=mask:
   if 16*d[m]==g**4*n[m]*n[m]:
    found=1
    print('FOUND',a,b,m)
    break
 print('PAIR',a,b,found)
PY