import numpy as np

f=open('PopIII_haloes_0038_type5.txt','r')

fo = open('Overlaped_haloes_list.txt','w')
line = f.readline()
radius = []
cx = []
cy = []
cz = []
i=0
while line !='':
   radius.append(float(line.split()[2]))
   cx.append(float(line.split()[3]))
   cy.append(float(line.split()[4]))
   cz.append(float(line.split()[5].replace(":","")))
   line = f.readline()
   pn = line.split()[0] 
   for j in xrange(int(pn)):
     line=f.readline()
   for j in xrange(i):
     if ((cx[i]-cx[j])**2+(cy[i]-cy[j])**2+(cz[i]-cz[j])**2)**0.5<radius[j]:
          fo.write(str(i)+"\n")
   f.readline();f.readline()
   line = f.readline()
   i=i+1

f.close()
fo.close()