#!/usr/bin/python2

file = open("allcurves","r").read()
isoclass=[]
for f in file.split("\n"):
    a = f.split()
    if len(isoclass) == 0:
        isoclass = [a[0],a[1],[a[3]]]
    elif isoclass[0:2] == a[0:2]:
        isoclass[2].append(a[3])
    else:
        eqns = isoclass[2]
        cmd = 'eqns=%s; allbad=1; for(i=1,length(eqns),e=ellinit(eqns[i]); if(abs(e.disc)<=10^12 && abs(e.c4)<=1.44*10^12, allbad=0)); if(allbad,print("%s\t%s"));'%(eqns,isoclass[0],isoclass[1])
        cmd = cmd.replace("'","")
        print cmd
        isoclass=[a[0],a[1],[a[3]]]