Mod-Jon_Code_9-10-11

257 days ago by Ben.LeVeque

%cython from libc.stdint cimport uint8_t, int32_t, uint32_t, int64_t, uint64_t from psage.libs.smalljac.wrapper1 import elliptic_curve_ap from sage.rings.integer cimport Integer from sage.schemes.elliptic_curves.constructor import EllipticCurve from psage.modform.hilbert.sqrt5.sqrt5 import F as K cdef extern from "pari/pari.h": uint32_t Fl_sqrt(uint32_t, uint32_t) ctypedef uint32_t pr_cur[2] ctypedef uint32_t uint_cur[4] ctypedef int32_t int_cur[10] def gal_con(alpha): return alpha[0]+alpha[1]*(1-K.gen(0)) cdef class SievedSearch: cdef uint32_t a11,a12,a31,a32,A1,A2,B1,B2,p,numPrimes,L cdef uint32_t primes[5],lens[5],liftCoef[2] cdef public uint32_t r1,r2 cdef int64_t a21,a22,a41,a42,a61,a62 cdef int_cur var_lifts[2304] cdef int_cur *cur_lift cdef uint32_t q,v,numFirst,numSecond cdef mpz_t s,t,u cdef pr_cur *firstp, *secondp def __init__(self, *args): for i,p in enumerate(args): self.primes[i] = p self.numPrimes = len(args) self.p = 1 for i in range(self.numPrimes): self.p *= self.primes[i] self.q = self.p+1u>>1u self.v = (self.p+1u)/3u if self.p%3u==2u else ((self.p<<1u)+1u)/3u cpdef get_primes(self,uint32_t p): cdef uint32_t g = Fl_sqrt(5u, p)+1u cdef uint32_t w = (<uint64_t>(p+1u)>>1u)*(<uint64_t>g)%p g = p+1u-w if w > g: self.r1=w self.r2=g else: self.r1=g self.r2=w cdef uint_cur *split_prime_combine(self, uint32_t i): cdef uint32_t p = self.primes[i] cdef uint32_t r1,r2,diff_inv self.get_primes(p) r1 = self.r1 r2 = self.r2 diff_inv = self.inv_mod(r1-r2,p) cdef uint_cur *ret ret = <uint_cur *>sage_malloc(sizeof(uint_cur)*self.lens[i]) cdef uint32_t *E, *F, *lift cdef pr_cur *Eptr = self.firstp cdef pr_cur *Fptr cdef uint_cur *liftprt liftptr = ret cdef uint32_t j,k for j in range(self.numFirst): E = Eptr[0] Fptr = self.secondp for k in range(self.numSecond): F = Fptr[0] lift = liftptr[0] self.custom_crt(E[0],F[0],r1,diff_inv,p) lift[0] = self.liftCoef[0] lift[1] = self.liftCoef[1] self.custom_crt(E[1],F[1],r1,diff_inv,p) lift[2] = self.liftCoef[0] lift[3] = self.liftCoef[1] liftptr += 1 Fptr += 1 Eptr += 1 return ret #-------------------------------------------------------------------------------# cdef uint_cur *integral_combine(self, uint_cur *first_list, uint_cur *second_list, uint32_t i): cdef uint32_t oldL = self.L self.L *= self.lens[i] cdef uint_cur *ret = <uint_cur *>sage_malloc(sizeof(uint_cur)*self.L) cdef uint32_t j,k,oldpMod,newp if not i: for j in range(self.L): ret[j][0] = second_list[j][0] ret[j][1] = second_list[j][1] ret[j][2] = second_list[j][2] ret[j][3] = second_list[j][3] return ret cdef uint32_t oldp = 1u cdef uint_cur *Eptr,*Fptr,*liftptr cdef uint32_t *E,*F,*lift for j in range(i): oldp *= self.primes[j] newp = oldp*self.primes[i] oldpMod = self.inv_mod(oldp,self.primes[i]) oldpMod *= oldp Eptr = first_list liftptr = ret for j in range(oldL): E = Eptr[0] Fptr = second_list for k in range(self.lens[i]): F = Fptr[0] lift = liftptr[0] lift[0] = self.integral_crt(E[0],F[0],oldpMod,newp) lift[1] = self.integral_crt(E[1],F[1],oldpMod,newp) lift[2] = self.integral_crt(E[2],F[2],oldpMod,newp) lift[3] = self.integral_crt(E[3],F[3],oldpMod,newp) liftptr += 1 Fptr += 1 Eptr += 1 return ret cdef uint32_t integral_crt(self, uint32_t s, uint32_t t, uint32_t pMod, uint32_t newp): return ((newp+t-s)*pMod+s)%newp # cdef uint_cur *integral_combine(self, uint32_t first_prime, uint_cur *first_prime_list, uint_cur *second_prime_list, uint32_t second_prime, uint32_t len_first_prime, uint32_t len_second_prime): # from sage.all import crt # cdef uint32_t n = first_prime*second_prime # cdef uint_cur *ret # ret = <uint_cur *>sage_malloc(sizeof(uint_cur)*len_first_prime*len_second_prime) # cdef uint32_t *E, *F # cdef uint_cur *Eptr = first_prime_list # cdef uint_cur *Fptr # cdef uint_cur *lift = ret # cdef uint32_t i,j # for i in range(len_first_prime): # E = Eptr[0] # Fptr = second_prime_list # for j in range(len_second_prime): # F = Fptr[0] # lift[0][0] = crt([E[0],F[0]],[first_prime,second_prime]) # lift[0][1] = crt([E[1],F[1]],[first_prime,second_prime]) # lift[0][2] = crt([E[2],F[2]],[first_prime,second_prime]) # lift[0][3] = crt([E[3],F[3]],[first_prime,second_prime]) # lift+=1 # Fptr+=1 # Eptr+=1 # return ret #-------------------------------------------------------------------------------# cdef custom_crt(self, uint32_t s, uint32_t t, uint32_t r1, uint32_t diff_inv, uint32_t p): self.liftCoef[1] = (p+s-t)*diff_inv%p self.liftCoef[0] = (p+s-self.liftCoef[1]*r1%p)%p cdef uint32_t inv_mod(self, uint32_t x, uint32_t p): cdef uint32_t t,u1,u2,v1,v2 u1=x v1=1 v2=(p-p/x)%p u2=v2*x%p while u2>1: t=u1/u2 v1,v2=v2,(p+v1-t*v2%p)%p u1,u2=u2,(p+u1-t*u2%p)%p return v2 cdef pr_cur *curves_in_GF(self, uint32_t which, int32_t ap, uint32_t k): cdef uint32_t p = self.primes[k] cdef pr_cur *curves = <pr_cur *>sage_malloc(sizeof(pr_cur)<<4) cdef uint32_t cur_len = 16u cdef uint32_t ret = 0u cdef uint32_t i,j for i in range(p): for j in range(p): if 64u*i*i*i + 432u*j*j: if elliptic_curve_ap(i, j, p) == ap: if ret == cur_len: cur_len <<= 1u curves = <pr_cur *>sage_realloc(curves, sizeof(pr_cur)*cur_len) curves[ret][0] = i curves[ret][1] = j ret += 1u if which == 1u: self.numFirst = ret else: self.numSecond = ret return curves cdef do_clift(self): self.cur_lift[0][0] = self.a11 self.cur_lift[0][1] = self.a12 self.cur_lift[0][2] = self.a21 self.cur_lift[0][3] = self.a22 self.cur_lift[0][4] = self.a31 self.cur_lift[0][5] = self.a32 self.cur_lift[0][6] = self.a41 self.cur_lift[0][7] = self.a42 self.cur_lift[0][8] = self.a61 self.cur_lift[0][9] = self.a62 self.cur_lift += 1 cdef do_a4lift(self): cdef int64_t old1,old2 old1 = self.a61 old2 = self.a62 self.do_clift()#0 self.a61-=self.p self.do_clift()#1 self.a62-=self.p self.do_clift()#3 self.a61+=self.p self.do_clift()#2 self.a61 = old1 self.a62 = old2 ######################################################################################################### cdef do_t(self): cdef uint64_t s1 = self.q if self.a11 else 0u cdef uint64_t s2 = self.q if self.a12 else 0u cdef uint64_t t1 = self.q if self.a31 else 0u cdef uint64_t t2 = self.q if self.a32 else 0u cdef uint64_t r1,r2 cdef uint64_t xsqr1 = s1*s1+s2*s2 cdef uint64_t xsqr2 = s2*((s1<<1u)+s2) r1 = self.v*(self.p+xsqr1+self.a21)%self.p r2 = self.v*(xsqr2+self.a22)%self.p cdef uint64_t rpow1 = (r1*r1+r2*r2)%self.p cdef uint64_t rpow2 = r2*((r1<<1u)+r2)%self.p cdef uint64_t temp1 = s2*t2%self.p cdef uint64_t temp2 = (s1*t2+s2*t1+temp1)%self.p temp1 += s1*t1 self.a41 = (self.p+self.A1+3*rpow1-(temp1<<1u)%self.p)%self.p self.a42 = (self.p+self.A2+3*rpow2-(temp2<<1u)%self.p)%self.p temp1 = rpow2*r2 rpow1, rpow2 = temp1+rpow1*r1, temp1+rpow1*r2+rpow2*r1 temp1 = self.A2*r2 temp2 = self.A1*r2+self.A2*r1+temp1 temp1 += self.A1*r1 xsqr1 = (t1*t1+t2*t2)%self.p xsqr2 = t2*((t1<<1u)+t2)%self.p self.a61 = (self.p+self.B1+temp1+rpow1-xsqr1)%self.p self.a62 = (self.p+self.B2+temp2+rpow2-xsqr2)%self.p self.do_a4lift()#0 self.a41-=self.p self.do_a4lift()#1 self.a42-=self.p self.do_a4lift()#3 self.a41+=self.p self.do_a4lift()#2 cdef do_r(self): self.a31 = 0u self.a32 = 0u self.do_t() self.a31 = 1u self.do_t() self.a32 = 1u self.do_t() self.a31 = 0u self.do_t() cdef do_s(self): self.a21 = 0 self.a22 = 0 self.do_r() self.a21 = 1 self.do_r() self.a22 = 1 self.do_r() self.a22 = -1 self.do_r() self.a21 = 0 self.do_r() self.a22 = 1 self.do_r() self.a21 = -1 self.do_r() self.a22 = -1 self.do_r() self.a22 = 0 self.do_r() ################################################################## ################################################################## ################################################################## def __call__(self, cond, ap_list, *args): if len(args) != 2u*self.numPrimes: raise ValueError cdef uint32_t p = self.p cdef uint32_t i cdef uint_cur **big_list_newp = <uint_cur **>sage_malloc(sizeof(uint_cur *)*self.numPrimes) for i in range(self.numPrimes): self.firstp = self.curves_in_GF(1, args[2u*i],i) self.secondp = self.curves_in_GF(2, args[2u*i+1u],i) self.lens[i] = self.numFirst*self.numSecond big_list_newp[i] = self.split_prime_combine(i) sage_free(self.firstp) sage_free(self.secondp) cdef uint_cur *bigger_list_new, *temp self.L = 1u for i in range(self.numPrimes): temp = self.integral_combine(bigger_list_new, big_list_newp[i],i) if i: sage_free(bigger_list_new) sage_free(big_list_newp[i]) bigger_list_new = temp #sage_free(big_list_newp) cdef uint32_t L = self.L #print('Length of big_list_new%s:'%p,L) cdef uint32_t current = 1u cdef uint_cur *e = bigger_list_new gal_cond = gal_con(K(cond)) cN = Integer(K(cond).norm().abs()) cdef uint32_t N = cN tempfact = cN.factor() cdef int_cur *endLift = self.var_lifts+2304 cdef uint32_t Nfactlen = len(tempfact) cdef uint32_t *Nfact, *Nfactend a = K.gen(0) Nfact = <uint32_t *>sage_malloc(sizeof(uint32_t)*Nfactlen) for i in range(Nfactlen): Nfact[i] = tempfact[i][0]**tempfact[i][1] Nfactend = Nfact+Nfactlen mpz_init(self.s) mpz_init(self.t) mpz_init(self.u) for i in range(L): #print(current_inc,'of',L,'---',p) current += 1u self.A1 = e[0][0] self.A2 = e[0][1] self.B1 = e[0][2] self.B2 = e[0][3] self.cur_lift = self.var_lifts ########## self.a11 = 0u self.a12 = 0u self.do_s() self.a11 = 1u self.do_s() self.a12 = 1u self.do_s() self.a11 = 0u self.do_s() self.cur_lift = self.var_lifts while self.cur_lift != endLift: if self.check_disc(self.cur_lift[0],N,Nfact,Nfactend): E = EllipticCurve(K,[self.cur_lift[0][2*j]+a*self.cur_lift[0][2*j+1] for j in range(5)]) c = E.conductor() if c == cond: #print 1,E.a_invariants(), 'tentative' if self.check_aps_against(E,ap_list): #print E.a_invariants(), 'a_ps match' mpz_clear(self.t) mpz_clear(self.s) mpz_clear(self.u) sage_free(Nfact) sage_free(big_list_newp) sage_free(bigger_list_new) return E elif c == gal_cond: print E E = EllipticCurve(K,[gal_con(K(t)) for t in E.ainvs()]) #print E.a_invariants(), 'tentative' if self.check_aps_against(E,ap_list): #print E.a_invariants(), 'a_ps match' mpz_clear(self.t) mpz_clear(self.s) mpz_clear(self.u) sage_free(Nfact) sage_free(big_list_newp) sage_free(bigger_list_new) return E self.cur_lift+=1 e += 1 #if not (i+1u)%10000u: # print(i+1u,'of',L,'---',p) mpz_clear(self.t) mpz_clear(self.s) mpz_clear(self.u) sage_free(Nfact) sage_free(big_list_newp) sage_free(bigger_list_new) return None cdef bint check_disc(self, int_cur ainv, uint32_t N, uint32_t *fact, uint32_t *factEnd): cdef int32_t b21, b22, b41, b42, b61, b62, b81, b82 cdef int64_t temp b21 = ainv[0]+ainv[1] b22 = ainv[1]*(ainv[0]+b21)+(ainv[3]<<2) b21 += ainv[2]<<2 b81 = b22*ainv[9] b82 = b21*ainv[9]+b22*ainv[8]+b81 b81 += b21*ainv[8] b41 = ainv[1]*ainv[5] b42 = b41+ainv[0]*ainv[5]+ainv[1]*ainv[4] b41 += ainv[0]*ainv[4] temp = ainv[7]*b42 b82 -= temp + ainv[6]*b42 + ainv[7]*b41 b81 -= temp + ainv[6]*b41 b41 += ainv[6]<<1 b42 += ainv[7]<<1 b61 = ainv[4]+ainv[5] b62 = ainv[5]*(ainv[4]+b61) temp = ainv[3]*b62 b82 += temp + ainv[2]*b62 + ainv[3]*b61 b81 += temp + ainv[2]*b61 b61 += ainv[8]<<2 b62 += ainv[9]<<2 temp = ainv[7]*ainv[7] b81 -= temp + ainv[6]*ainv[6] b82 -= temp + ((ainv[6]*ainv[7])<<1) cdef int64_t tmp1,tmp2 tmp1 = b22*b42 tmp2 = tmp1+b22*b41+b21*b42-3*b62 tmp1 += b21*b41-3*b61 temp = b62*tmp2 tmp2 = 9*(temp+b62*tmp1+b61*tmp2) tmp1 = 9*(temp+tmp1*b61) mpz_set_si(self.s,tmp1) mpz_set_si(self.t,tmp2) tmp1 = b41*b41 tmp2 = b42*b42 temp = (b42+3*b41)*tmp2 tmp2 = temp+(tmp2+3*tmp1)*b42<<3 tmp1 = temp+tmp1*b41<<3 if tmp1 >= 0: mpz_sub_ui(self.s,self.s,tmp1) else: mpz_add_ui(self.s,self.s,-tmp1) if tmp2 >= 0: mpz_sub_ui(self.t,self.t,tmp2) else: mpz_add_ui(self.t,self.t,-tmp2) temp = b21*b21+b22*b22 b22 = ((b21<<1)+b22)*b22 b21 = temp tmp1 = b81 tmp2 = b82 temp = b22*tmp2 tmp2 = temp+tmp1*b22+tmp2*b21 tmp1 = temp+b21*tmp1 if tmp1 >= 0: mpz_sub_ui(self.s,self.s,tmp1) else: mpz_add_ui(self.s,self.s,-tmp1) if tmp2 >= 0: mpz_sub_ui(self.t,self.t,tmp2) else: mpz_add_ui(self.t,self.t,-tmp2) if not mpz_sgn(self.s) and not mpz_sgn(self.t): return 0 mpz_sub(self.u,self.s,self.t) mpz_mul(self.u,self.u,self.t) mpz_addmul(self.u,self.s,self.s) if mpz_mod_ui(self.t,self.u,N): return 0 while fact != factEnd: while not mpz_mod_ui(self.t,self.u,fact[0]): mpz_divexact_ui(self.u,self.u,fact[0]) fact += 1 if mpz_sgn(self.u) < 0: mpz_neg(self.u,self.u) return mpz_root(self.u, self.u, 12ul) def ap(self,E,p): return E.change_ring(p.residue_field()).trace_of_frobenius() def check_aps_against(self, E, aps): for i,p in enumerate([11,19,29,31,41,59,61,71,79,89]): if not isinstance(aps[2*i], str): r1,r2 = K.primes_above(p) if self.ap(E,r1) != aps[2*i] or self.ap(E,r2) != aps[2*i+1]: return False return True 
       
from psage.ellcurve.minmodel.sqrt5 import canonical_model from psage.modform.hilbert.sqrt5.sqrt5 import F as K a = K.gen(0) def ap(E,p): return E.change_ring(p.residue_field()).trace_of_frobenius() def aplist(E): aplist=[] for p in [2,3,5,7,11,19,29,31,41,59,61,71,79,89]: for i in K.primes_above(p): try: aplist.append(ap(E,i)) except ArithmeticError: aplist.append('?') return aplist 
       
from psage.ellcurve.minmodel.sqrt5 import canonical_model do = {} primes = [11,19,29,31,41,59,61,71,79,89] for i,p in enumerate(primes): do[p] = SievedSearch(p) for p2 in primes[i+1:]: do[p*p2] = SievedSearch(p,p2) def curve_finder(cond, aps): for i in range(10): pr = primes[i] ap1 = aps[2*i] if not isinstance(ap1, str): ap2 = aps[2*i+1] p1,p2 = K.primes_above(pr) if p1.integral_basis()[1][0] > p2.integral_basis()[1][0]: ap1,ap2=ap2,ap1 # E = do[pr](cond,aps,ap1,ap2) # if E is not None: # return canonical_model(E.global_minimal_model()) for j in range(i+1,10): pr2 = primes[j] ap3 = aps[2*j] if not isinstance(ap3, str): ap4 = aps[2*j+1] print 'Current search:', pr*pr2 p1,p2 = K.primes_above(pr2) if p1.integral_basis()[1][0] > p2.integral_basis()[1][0]: ap3,ap4=ap4,ap3 E = do[pr*pr2](cond,aps,ap1,ap2,ap3,ap4) if E is not None: return canonical_model(E.global_minimal_model()) return None 
       
%time #f1 = open('/Users/bleveque/Sage/unmatched_728.txt') missing_curves = ["1856 -8*a + 48 ['?',4,-4,-2,4,-2,-4,0,'?','?',4,4,10,8,12,-10,4,-6,10,0,-14,16,6,14]"] for s in missing_curves: S = s.split('[') aps = eval('['+S[1]) aps = aps[4:] print S[0] t = S[0].split(' ') cond = eval(''.join(t[1:])) E = curve_finder(cond,aps) if E is None: print 'none' else: print E.conductor(), E.a_invariants() 
       
1856 -8*a + 48 
Current search: 209
Current search: 341
Current search: 451
Current search: 649
Current search: 671
Current search: 781
Current search: 869
Current search: 979
Current search: 589
Current search: 779
Current search: 1121
Current search: 1159
Current search: 1349
1856 -8*a + 48 
Current search: 209
Current search: 341
Current search: 451
Current search: 649
Current search: 671
Current search: 781
Current search: 869
Current search: 979
Current search: 589
Current search: 779
Current search: 1121
Current search: 1159
Current search: 1349

 

#############################################################
ainv = (a + 1, -a + 1, a + 1, -62*a - 86, -411*a - 311) ######## Needed Functions ########## def ap(E,p): return E.change_ring(K.residue_field(p)).trace_of_frobenius() def ap_list(E): aplist=[] for p in [2,3,5,7,11,19,29,31,41,59,61,71,79,89]: for i in K.primes_above(p): try: aplist.append(ap(E,i)) except: aplist.append('?') return aplist ch = var('ch') K.<a>=NumberField(ch^2-ch-1) E = EllipticCurve(K,list(ainv)) S = str(E.conductor().norm())+' '+str(E.conductor().gens_reduced()[0])+' '+str(ap_list(E)).replace(' ','')+str(list(E.a_invariants())) print S 
       
1805 -a + 43
[1,2,'?',2,0,0,'?',-4,6,6,8,4,6,-6,-12,12,-14,10,0,-12,8,16,-6,-6][a +
1, -a + 1, a + 1, -62*a - 86, -411*a - 311]
1805 -a + 43 [1,2,'?',2,0,0,'?',-4,6,6,8,4,6,-6,-12,12,-14,10,0,-12,8,16,-6,-6][a + 1, -a + 1, a + 1, -62*a - 86, -411*a - 311]
29*31 
       
899
899
RR(3129/60) 
       
52.1500000000000
52.1500000000000
451 19*a - 10 Current search: 19 Current search: 551 Fractional ideal (19*a - 10) (a + 1, a - 1, 0, -193*a - 127, -1754*a - 1070) 549 21*a - 12 Current search: 11 Current search: 209 Fractional ideal (21*a - 12) (a, a - 1, a, 36*a - 65, -155*a + 226) 781 25*a - 13 Current search: 19 Current search: 551 Fractional ideal (25*a - 13) (a + 1, -a - 1, a + 1, -144*a - 95, 1048*a + 654) 836 26*a - 16 Current search: 29 Current search: 899 Fractional ideal (-26*a + 16) (1, a + 1, a, 47*a - 71, 173*a - 218) 836 2*a - 30 Current search: 29 Current search: 899 Fractional ideal (2*a - 30) (a, a - 1, a + 1, -49*a + 37, 137*a - 108) 
       
E1 = EllipticCurve(K,[a+1,a-1,0,-193*a-127,-1754*a-1070]) 
       
E1.conductor().norm() 
       
451
451
aplist(E1) 
       
[3, 1, -1, 4, 0, '?', -6, 1, -6, 2, 4, 0, '?', 0, -1, 2, -3, 6, -6, -12,
11, 14, -3, -16]
[3, 1, -1, 4, 0, '?', -6, 1, -6, 2, 4, 0, '?', 0, -1, 2, -3, 6, -6, -12, 11, 14, -3, -16]
E2 = EllipticCurve(K,[a, a - 1, a, 36*a - 65, -155*a + 226]) 
       
E2.conductor().norm() 
       
549
549
aplist(E2) 
       
[3, '?', -1, 3, -3, 3, -2, 7, -6, 1, 3, -5, -2, 6, 6, -14, -11, '?', -8,
-5, 14, 0, -1, -10]
[3, '?', -1, 3, -3, 3, -2, 7, -6, 1, 3, -5, -2, 6, 6, -14, -11, '?', -8, -5, 14, 0, -1, -10]
E3 = EllipticCurve(K,[a + 1, -a - 1, a + 1, -144*a - 95, 1048*a+ 654]) 
       
E3.conductor().norm() 
       
781
781
aplist(E3) 
       
[1, 2, 2, -2, '?', -4, 8, 4, -2, 6, -4, -8, 10, -2, -12, 0, 2, -2, -12,
'?', 4, 8, -6, -6]
[1, 2, 2, -2, '?', -4, 8, 4, -2, 6, -4, -8, 10, -2, -12, 0, 2, -2, -12, '?', 4, 8, -6, -6]
E4 = EllipticCurve(K,[1, a + 1, a, 47*a - 71, 173*a - 218]) 
       
E4.conductor().norm() 
       
836
836
aplist(E4) 
       
['?', -4, 2, -4, 4, '?', '?', -2, 4, 8, 10, 0, -4, -4, -14, -12, -2, -2,
-4, 0, 4, -16, 14, -10]
['?', -4, 2, -4, 4, '?', '?', -2, 4, 8, 10, 0, -4, -4, -14, -12, -2, -2, -4, 0, 4, -16, 14, -10]
E5 = EllipticCurve(K,[a, a - 1, a + 1, -49*a + 37, 137*a - 108]) 
       
E5.conductor().norm() 
       
836
836
aplist(E5) 
       
['?', 2, 1, 0, 4, '?', -2, '?', -2, 4, 7, 4, -5, 9, 0, -11, -1, -2, 6,
-7, -2, 8, 7, -8]
['?', 2, 1, 0, 4, '?', -2, '?', -2, 4, 7, 4, -5, 9, 0, -11, -1, -2, 6, -7, -2, 8, 7, -8]
E=EllipticCurve([a,1]) 
       
G = E.change_ring(K.residue_field(K.primes_above(11)[0])) 
       
F = E.change_ring(K.residue_field(K.primes_above(11)[0])) 
       
F.is_isomorphic(G) 
       
False
False
       
Elliptic Curve defined by y^2 = x^3 + 4*x + 1 over Residue field of
Fractional ideal (3*a - 2)
Elliptic Curve defined by y^2 = x^3 + 4*x + 1 over Residue field of Fractional ideal (3*a - 2)
       
Elliptic Curve defined by y^2 = x^3 + 8*x + 1 over Residue field of
Fractional ideal (3*a - 2)
Elliptic Curve defined by y^2 = x^3 + 8*x + 1 over Residue field of Fractional ideal (3*a - 2)
F.change_weierstrass_model([1,5,-4,2]) 
       
Elliptic Curve defined by y^2 = x^3 + 8*x + 1 over Residue field of
Fractional ideal (3*a - 1)
Elliptic Curve defined by y^2 = x^3 + 8*x + 1 over Residue field of Fractional ideal (3*a - 1)
F = E.change_weierstrass_model([1,3*a+3,5*a,8*a+4]); F.ainvs() 
       
(11*a, -22*a - 22, 22*a + 11, -109*a - 66, -11*a + 1)
(11*a, -22*a - 22, 22*a + 11, -109*a - 66, -11*a + 1)
-109%11 
       
1
1
F[4][0]%11 
       
8
8
f = open('/Users/bleveque/Sage/unmatched_728.txt') i = 0 for s in f.readlines(): S = s.split('[') aps = eval('['+S[1]) aps = aps[4:] print S[0] t = S[0].split(' ') cond = eval(''.join(t[1:])) E = curve_finder(cond,aps) g = open('/Users/bleveque/Sage/modjon728/output%s.txt'%i,'w') if E is None: g.write(s) g.write(' unfound\n') else: g.write(s[:-1]) g.write(' '+str(list(E.ainvs()))+'\n') g.close() i += 1 
       
179 -12*a + 5 
Traceback (click to the left of this block for traceback)
...
NameError: name 'a' is not defined
179 -12*a + 5 
Traceback (most recent call last):        print S[0]
  File "", line 1, in <module>
    
  File "/private/var/folders/+3/+3+bRf7bGhaMk-2vvrWEvU+++TI/-Tmp-/tmp9S8hAv/___code___.py", line 5, in <module>
    exec compile(u"for s in f.readlines():\n    S = s.split('[')\n    aps = eval('['+S[_sage_const_1 ])\n    aps = aps[_sage_const_4 :]\n    print S[_sage_const_0 ]\n    t = S[_sage_const_0 ].split(' ')\n    cond = eval(''.join(t[_sage_const_1 :]))\n    E = curve_finder(cond,aps)\n    g = open('/Users/bleveque/Sage/modjon728/output%s.txt'%i,'w')\n    if E is None:\n        g.write(s)\n        g.write(' unfound\\n')\n    else:\n        g.write(s[:-_sage_const_1 ])\n        g.write(' '+str(list(E.ainvs()))+'\\n')\n    g.close()\n    i += _sage_const_1 " + '\n', '', 'single')
  File "", line 7, in <module>
    
  File "<string>", line 1, in <module>
NameError: name 'a' is not defined
225%29 
       
22
22
E = EllipticCurve(K,[1-a, a-1 - 1, 1, 15*(1-a) - 26, 35*(1-a) - 57]) 
       
canonical_model(E) 
       
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + x^2 +
(-15*a-11)*x + (-50*a-33) over Number Field in a with defining
polynomial x^2 - x - 1
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + x^2 + (-15*a-11)*x + (-50*a-33) over Number Field in a with defining polynomial x^2 - x - 1
E.conductor() 
       
Fractional ideal (-12*a + 5)
Fractional ideal (-12*a + 5)
E.change_ring(K.residue_field(K.primes_above(29)[0])).trace_of_frobenius() 
       
0
0
%time N=205 cond=13*a - 9 aps=[-2,0,-4,4,-8,-6,10,-6,'?','?',12,4,6,2,-2,4,-16,-4,-6,-8] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i] ap2 = aps[2*i+1] if not isinstance(ap1, str): p1,p2 = K.primes_above(pr) if p1.integral_basis()[0][1] > p2.integral_basis()[0][1]: ap1,ap2=ap2,ap1 if Dopr(cond,ap1,ap2,aps): break 
       
current prime =  11
current prime =  19
(-a + 1, -a, 0, -4, 5*a - 2) tentative
(-a + 2, -a + 1, 0, -3*a + 5, 0) tentative
(-a + 1, -a, 0, 1, 0) tentative
current prime =  29
(a, a, 1, -8*a - 16, 25*a + 2) tentative
(1, a - 1, -a + 2, a, 3*a) tentative
current prime =  31
current prime =  41
current prime =  59
(10000L, 'of', 40368L, '---', 59L)
(20000L, 'of', 40368L, '---', 59L)
(-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4) tentative
(30000L, 'of', 40368L, '---', 59L)
(40000L, 'of', 40368L, '---', 59L)
current prime =  61
(10000L, 'of', 43200L, '---', 61L)
(20000L, 'of', 43200L, '---', 61L)
(30000L, 'of', 43200L, '---', 61L)
(-a + 1, 0, -a + 1, a - 1, a) tentative
(40000L, 'of', 43200L, '---', 61L)
(-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4) tentative
current prime =  71
(10000L, 'of', 19600L, '---', 71L)
(a, a, 1, -8*a - 16, 25*a + 2) tentative
current prime =  79
(10000L, 'of', 56784L, '---', 79L)
(20000L, 'of', 56784L, '---', 79L)
(30000L, 'of', 56784L, '---', 79L)
(-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4) tentative
(40000L, 'of', 56784L, '---', 79L)
(50000L, 'of', 56784L, '---', 79L)
current prime =  89
(10000L, 'of', 108416L, '---', 89L)
(20000L, 'of', 108416L, '---', 89L)
(30000L, 'of', 108416L, '---', 89L)
(40000L, 'of', 108416L, '---', 89L)
(50000L, 'of', 108416L, '---', 89L)
(60000L, 'of', 108416L, '---', 89L)
(70000L, 'of', 108416L, '---', 89L)
(80000L, 'of', 108416L, '---', 89L)
(90000L, 'of', 108416L, '---', 89L)
(100000L, 'of', 108416L, '---', 89L)
CPU time: 209.05 s,  Wall time: 211.04 s
current prime =  11
current prime =  19
(-a + 1, -a, 0, -4, 5*a - 2) tentative
(-a + 2, -a + 1, 0, -3*a + 5, 0) tentative
(-a + 1, -a, 0, 1, 0) tentative
current prime =  29
(a, a, 1, -8*a - 16, 25*a + 2) tentative
(1, a - 1, -a + 2, a, 3*a) tentative
current prime =  31
current prime =  41
current prime =  59
(10000L, 'of', 40368L, '---', 59L)
(20000L, 'of', 40368L, '---', 59L)
(-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4) tentative
(30000L, 'of', 40368L, '---', 59L)
(40000L, 'of', 40368L, '---', 59L)
current prime =  61
(10000L, 'of', 43200L, '---', 61L)
(20000L, 'of', 43200L, '---', 61L)
(30000L, 'of', 43200L, '---', 61L)
(-a + 1, 0, -a + 1, a - 1, a) tentative
(40000L, 'of', 43200L, '---', 61L)
(-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4) tentative
current prime =  71
(10000L, 'of', 19600L, '---', 71L)
(a, a, 1, -8*a - 16, 25*a + 2) tentative
current prime =  79
(10000L, 'of', 56784L, '---', 79L)
(20000L, 'of', 56784L, '---', 79L)
(30000L, 'of', 56784L, '---', 79L)
(-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4) tentative
(40000L, 'of', 56784L, '---', 79L)
(50000L, 'of', 56784L, '---', 79L)
current prime =  89
(10000L, 'of', 108416L, '---', 89L)
(20000L, 'of', 108416L, '---', 89L)
(30000L, 'of', 108416L, '---', 89L)
(40000L, 'of', 108416L, '---', 89L)
(50000L, 'of', 108416L, '---', 89L)
(60000L, 'of', 108416L, '---', 89L)
(70000L, 'of', 108416L, '---', 89L)
(80000L, 'of', 108416L, '---', 89L)
(90000L, 'of', 108416L, '---', 89L)
(100000L, 'of', 108416L, '---', 89L)
CPU time: 209.05 s,  Wall time: 211.04 s
E = EllipticCurve(K,[-a + 2, -a + 2, 0, -11*a - 12, 40*a - 4]) lst = [] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): for p in K.primes_above(pr): try: lst.append(ap(E,p)) except ArithmeticError: lst.append('?') print lst 
       
[0, -2, 4, -4, -6, -8, -6, 10, '?', 10, 4, 12, 2, 6, 4, -2, -4, -16, -8,
-6]
[0, -2, 4, -4, -6, -8, -6, 10, '?', 10, 4, 12, 2, 6, 4, -2, -4, -16, -8, -6]
-9+13*(1-a) 
       
Fractional ideal (-13*a + 4)
Fractional ideal (-13*a + 4)
E.conductor() 
       
Fractional ideal (13*a - 9)
Fractional ideal (13*a - 9)
aps=[-2,0,-4,4,-8,-6,10,-6,'?','?',12,4,6,2,-2,4,-16,-4,-6,-8] print lst print aps 
       
[0, -2, 4, -4, -6, -8, -6, 10, '?', 10, 4, 12, 2, 6, 4, -2, -4, -16, -8,
-6]
[-2, 0, -4, 4, -8, -6, 10, -6, '?', '?', 12, 4, 6, 2, -2, 4, -16, -4,
-6, -8]
[0, -2, 4, -4, -6, -8, -6, 10, '?', 10, 4, 12, 2, 6, 4, -2, -4, -16, -8, -6]
[-2, 0, -4, 4, -8, -6, 10, -6, '?', '?', 12, 4, 6, 2, -2, 4, -16, -4, -6, -8]
N=599 cond=a + 24 aps=[-4,6,6,-4,-2,2,2,0,-8,0,12,6,6,-8,-14,-8,2,-16,6,-2] for i,pr in enumerate([11,19,29,31,41,59]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if pr==11: ap1,ap2=ap2,ap1 if Dopr(cond,ap1,ap2,aps): break 
       
current prime =  11
current prime =  19
current prime =  29
current prime =  31
(-a + 1, -a, -a + 1, 5*a - 23, -24*a + 23) tentative
(-a + 1, -a, -a + 1, 5*a - 23, -24*a + 23) a_ps match
current prime =  11
current prime =  19
current prime =  29
current prime =  31
(-a + 1, -a, -a + 1, 5*a - 23, -24*a + 23) tentative
(-a + 1, -a, -a + 1, 5*a - 23, -24*a + 23) a_ps match
E11 = EllipticCurve(K,[a+1,a,a,5*a-22,-42*a+28]) print E11.conductor().norm(), E11.conductor() ap(E11,K.ideal(3*a-2)) 
       
599 Fractional ideal (a + 24)
-4
599 Fractional ideal (a + 24)
-4
K.<a> = NumberField(x**2-x-1) 
       
N=284 cond=2*a + 16 aps=[-2,4,2,6,0,4,-8,-8,-10,-6,0,2,8,-10,'?','?',4,12,-6,-10] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  61
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  71
current prime =  79
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  61
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  71
current prime =  79
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-9)*x + (3*a-16) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
E12 = EllipticCurve(K,[a+1,a+1,1,5*a-9,3*a-16]) print E12.conductor().norm(), E12.conductor() ap(E12,K.ideal(3*a-2)) 
       
284 Fractional ideal (2*a + 16)
-2
284 Fractional ideal (2*a + 16)
-2
N=380 cond=-4*a - 18 aps=[4,2,'?','?',-2,4,0,4,2,-12,0,-4,6,-2,-8,0,-16,-6,-14,6] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
current prime =  29
singular curve
singular curve
singular curve
current prime =  31
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  41
current prime =  59
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  61
singular curve
current prime =  71
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (a+1)*x^2 +
(-12*a+6)*x + (-76*a+64) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
current prime =  89
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (a+1)*x^2 +
(-12*a+6)*x + (-76*a+64) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
current prime =  11
current prime =  19
current prime =  29
singular curve
singular curve
singular curve
current prime =  31
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  41
current prime =  59
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  61
singular curve
current prime =  71
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (a+1)*x^2 + (-12*a+6)*x + (-76*a+64) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  89
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (a+1)*x^2 + (-12*a+6)*x + (-76*a+64) over Number Field in a with defining polynomial x^2 - x - 1 tentative
E22=EllipticCurve(K,[a+1, a + 1, 0,-12*a+6,-76*a+64]) #this curve is probably Gal. conj to a curve we want print E22.conductor() ap(E22,K.primes_above(29)[0]) 
       
Fractional ideal (-4*a - 18)
4
Fractional ideal (-4*a - 18)
4
N=395 cond=2*a + 19 aps=[-2,2,0,0,-4,6,0,-4,2,-6,12,4,-14,-10,4,-6,'?','?',0,0] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  19
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
current prime =  29
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
current prime =  41
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  59
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  61
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  71
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  79
current prime =  89
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  11
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  19
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
current prime =  29
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
current prime =  41
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  59
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  61
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  71
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  79
current prime =  89
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-6*a-23)*x + (2*a-32) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (2*a-7)*x + (3*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-11*a-8)*x + (-25*a-16) over Number Field in a with defining polynomial x^2 - x - 1 tentative
E23=EllipticCurve(K,[a,0,0,-6*a-23,2*a-32]) #one of these may work E24=EllipticCurve(K,[1,0,0,2*a-7,3*a-8]) E25=EllipticCurve(K,[a,0,0,-11*a-8,-25*a-16]) print E23.conductor(), E24.conductor(), E25.conductor() print ap(E23,K.primes_above(89)[1]), ap(E24,K.primes_above(89)[1]), ap(E25,K.primes_above(89)[1]) 
       
Fractional ideal (2*a + 19) Fractional ideal (2*a + 19) Fractional ideal
(2*a + 19)
0 0 0
Fractional ideal (2*a + 19) Fractional ideal (2*a + 19) Fractional ideal (2*a + 19)
0 0 0
N=396 cond=18*a - 12 aps=['?','?',-7,7,-7,0,6,-1,2,-5,4,11,-12,-5,12,5,-14,-14,-11,10] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  41
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  89
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x
+ (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  41
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  89
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-a-1)*x^2 + (-a-3)*x + (2*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
E13 = EllipticCurve(K,[a+1,-a-1,0,-a-3,2*a-9]) print E13.conductor().norm(), E13.conductor() ap(E13,K.primes_above(19)[1]) 
       
396 Fractional ideal (18*a - 12)
7
396 Fractional ideal (18*a - 12)
7
N=404 cond=18*a - 10 aps=[-1,-2,0,-2,3,0,5,-1,-6,-1,14,-9,-10,7,-8,-4,-8,-10,0,17] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
singular curve
singular curve
current prime =  19
singular curve
singular curve
singular curve
current prime =  29
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  31
singular curve
current prime =  41
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x +
(13*a+7) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x +
(13*a+7) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x +
(13*a+7) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x +
(13*a+7) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x +
(13*a+7) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  11
singular curve
singular curve
singular curve
current prime =  19
singular curve
singular curve
singular curve
current prime =  29
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  31
singular curve
current prime =  41
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x + (13*a+7) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x + (13*a+7) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x + (13*a+7) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x + (13*a+7) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (48*a+28)*x + (13*a+7) over Number Field in a with defining polynomial x^2 - x - 1 tentative
E27 = EllipticCurve(K,[a+1,0,a,48*a+28,13*a+7]) #gal-conj print E27.conductor().norm(), E27.conductor() ap(E27,K.primes_above(19)[0]) 
       
404 Fractional ideal (18*a - 10)
-2
404 Fractional ideal (18*a - 10)
-2
N=404 cond=18*a - 8 aps=[-1,6,4,2,9,-4,-5,-11,-2,3,6,-5,-6,1,-16,-12,-8,6,12,9] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  41
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
singular curve
singular curve
singular curve
current prime =  59
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  61
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x +
(-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x +
(-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 +
(-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2
- x - 1 a_ps match
current prime =  89
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x +
(-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x +
(-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  11
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  41
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
singular curve
singular curve
current prime =  59
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  61
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x + (-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x + (-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + x^2 + (-a-6)*x + (-3*a-9) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  89
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x + (-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + (-1)*x^2 + (-20*a-13)*x + (-50*a-31) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
E18 = EllipticCurve(K,[1,-1,0,-20*a-13,-50*a-31]) print E18.conductor().norm(), E18.conductor() ap(E18,K.primes_above(19)[0]) 
       
404 Fractional ideal (18*a - 8)
4
404 Fractional ideal (18*a - 8)
4
N=445 cond=-19*a + 7 aps=[0,0,-8,4,-10,-2,0,4,10,2,0,4,14,-2,-16,8,-12,-4,'?','?'] for i,pr in enumerate([11,19,29,31,41,59]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  19
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-1)*x^2 +
(-7*a-7)*x + (-12*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + a*x^2 + (-6*a-4)*x
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 = x^3 + x^2 + (-7*a-20)*x + (3*a+20) over
Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-1)*x^2 +
(-2*a-2)*x + (-a) over Number Field in a with defining polynomial x^2 -
x - 1 tentative
current prime =  31
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  41
singular curve
current prime =  59
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + x^2 +
(-23*a-16)*x + (52*a+30) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + x^2 +
(-23*a-16)*x + (52*a+30) over Number Field in a with defining polynomial
x^2 - x - 1 a_ps match
current prime =  11
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  19
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-1)*x^2 + (-7*a-7)*x + (-12*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + a*x^2 + (-6*a-4)*x over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 = x^3 + x^2 + (-7*a-20)*x + (3*a+20) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-1)*x^2 + (-2*a-2)*x + (-a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  31
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  41
singular curve
current prime =  59
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + x^2 + (-23*a-16)*x + (52*a+30) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + x^2 + (-23*a-16)*x + (52*a+30) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
E19 = EllipticCurve(K,[a,1,a+1,-23*a-16,52*a+30]) print E19.conductor().norm(), E19.conductor() ap(E19,K.primes_above(19)[0]) 
       
445 Fractional ideal (-19*a + 7)
-8
445 Fractional ideal (-19*a + 7)
-8
N=484 cond=22+0*a aps=['?','?',7,7,-3,-3,-5,-5,-8,-8,-12,-12,5,5,-7,-7,4,4,7,7] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+1] ap2 = aps[2*i] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  29
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  59
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  79
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  89
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x +
(-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  29
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  59
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  79
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
current prime =  89
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + x^2 + (-2*a-2)*x + (-16*a-12) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
E32=EllipticCurve(K,[a+1,1,0,-2*a-2,-16*a-12]) print E32.conductor() print ap(E32,K.primes_above(19)[1]) 
       
Fractional ideal (22)
7
Fractional ideal (22)
7
N=495 cond=-3*a + 24 aps=[3,'?','?',-2,'?','?',0,4,2,-2,-8,4,2,-10,4,0,-10,-6,0,-8,4,-16,-10,14] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (-a)*x^2 +
(-5*a-3)*x over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-21*a-7)*x + (-21*a-2) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-4)*x + (a+1) over Number Field in a with defining polynomial x^2 -
x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(a+1)*x over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
Elliptic Curve defined by y^2 = x^3 + a*x^2 + (-12*a-4)*x + (-12*a-4)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  41
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (-a)*x^2 + (-a-3)*x + (-a-3)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (-a)*x^2 +
(-5*a-3)*x over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
current prime =  71
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (a-1)*x^2 +
(-17*a-10)*x + (-63*a-39) over Number Field in a with defining
polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (a-1)*x^2 +
(-17*a-10)*x + (-63*a-39) over Number Field in a with defining
polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (a-1)*x^2 +
(-17*a-10)*x + (-63*a-39) over Number Field in a with defining
polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 +
(5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
current prime =  11
current prime =  19
singular curve
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (-a)*x^2 + (-5*a-3)*x over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-21*a-7)*x + (-21*a-2) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-4)*x + (a+1) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (a+1)*x over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 = x^3 + a*x^2 + (-12*a-4)*x + (-12*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  41
singular curve
Elliptic Curve defined by y^2 + x*y = x^3 + (-a)*x^2 + (-a-3)*x + (-a-3) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (-a)*x^2 + (-5*a-3)*x over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
current prime =  71
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (a-1)*x^2 + (-17*a-10)*x + (-63*a-39) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (a-1)*x^2 + (-17*a-10)*x + (-63*a-39) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (10*a-12)*x + (15*a-23) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (a-1)*x^2 + (-17*a-10)*x + (-63*a-39) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + (a+1)*x^2 + (5*a-17)*x + (-5*a-29) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
E41 = EllipticCurve(K,[1,-a,a+1,-5*a-3,0]) #42,43,45, or 49 could be gal conj to our curve E42 = EllipticCurve(K,[a+1,a+1,1,5*a-17,-5*a-29]) E43 = EllipticCurve(K,[a+1,a+1,1,10*a-12,15*a-23]) E44 = EllipticCurve(K,[a+1,-a+1,a+1,-21*a-7,-21*a-2]) E45 = EllipticCurve(K,[1,-a,0,-a-3,-a-3]) E46 = EllipticCurve(K,[0,a,0,-12*a-4,-12*a-4]) E47 = EllipticCurve(K,[a+1,a+1,a,a+1,0]) E48 = EllipticCurve(K,[a+1,a+1,a,-4*a-4,a+1]) E49 = EllipticCurve(K,[a,a-1,a+1,-17*a-10,-63*a-39]) print E41.conductor(), E42.conductor(), E43.conductor(),E44.conductor(), E45.conductor(), E46.conductor(),E47.conductor(), E48.conductor(), E49.conductor() for curve in [E41, E42, E43, E44, E45, E46, E47, E48, E49]: print ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(89)[0]),ap(curve,K.primes_above(89)[1]) 
       
Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional
ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a +
24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional
ideal (-3*a + 24) Fractional ideal (-3*a + 24)
-4 -8 -6 6 -6 2
4 0 -2 2 14 -10
4 0 -2 2 14 -10
-4 -8 -6 6 -6 2
4 0 -2 2 14 -10
-4 -8 -6 6 -6 2
-4 -8 -6 6 -6 2
-4 -8 -6 6 -6 2
4 0 -2 2 14 -10
Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24) Fractional ideal (-3*a + 24)
-4 -8 -6 6 -6 2
4 0 -2 2 14 -10
4 0 -2 2 14 -10
-4 -8 -6 6 -6 2
4 0 -2 2 14 -10
-4 -8 -6 6 -6 2
-4 -8 -6 6 -6 2
-4 -8 -6 6 -6 2
4 0 -2 2 14 -10
N=505 cond=a + 22 aps=[0,1,'?',9,5,-6,7,8,-6,0,2,-8,0,-9,5,-14,-5,0,-1,13,-6,-8,-1,-18] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
current prime =  29
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  31
singular curve
(0, -a - 1, 0, 9*a + 3, a - 8) tentative
(0, a, 1, -4*a - 2, 4*a + 2) tentative
current prime =  41
current prime =  59
current prime =  61
singular curve
current prime =  71
current prime =  79
singular curve
current prime =  89
current prime =  11
current prime =  19
current prime =  29
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  31
singular curve
(0, -a - 1, 0, 9*a + 3, a - 8) tentative
(0, a, 1, -4*a - 2, 4*a + 2) tentative
current prime =  41
current prime =  59
current prime =  61
singular curve
current prime =  71
current prime =  79
singular curve
current prime =  89
E505_1 = EllipticCurve(K,[0,-a-1,0,9*a+3,a-8]) #probably not... E505_2 = EllipticCurve(K,[0,a,1,-4*a-2,4*a+2]) for curve in [E505_1, E505_2]: print curve.conductor(), ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(31)[0]),ap(curve,K.primes_above(31)[1]) 
       
Fractional ideal (a + 22) 4 -4 -6 6 -4 4
Fractional ideal (a + 22) -5 0 10 0 2 -8
Fractional ideal (a + 22) 4 -4 -6 6 -4 4
Fractional ideal (a + 22) -5 0 10 0 2 -8
N=539 cond=21*a - 14 aps=[-3,-1,-1,'?','?','?',-6,-3,-1,-4,-4,-4,5,8,11,-7,13,-14,4,-2,5,-10,-9,6] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  61
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  79
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
current prime =  89
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 +
(-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
singular curve
current prime =  11
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  61
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  79
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  89
Elliptic Curve defined by y^2 + x*y + (a+1)*y = x^3 + (a+1)*x^2 + (-28*a-18)*x + (-13*a-8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-4*a-6)*x + (-2*a) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
E539_1 = EllipticCurve(K,[1,a+1,a+1,-28*a-18,-13*a-8]) #these work E539_2 = EllipticCurve(K,[a+1,-a+1,a+1,-4*a-6,-2*a]) for curve in [E539_1, E539_2]: print curve.conductor(), ap(curve,K.primes_above(2)[0]),ap(curve,K.primes_above(3)[0]), ap(curve,K.primes_above(5)[0]), ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(31)[0]),ap(curve,K.primes_above(31)[1]),ap(curve,K.primes_above(41)[0]),ap(curve,K.primes_above(41)[1]),ap(curve,K.primes_above(59)[0]),ap(curve,K.primes_above(59)[1]),ap(curve,K.primes_above(61)[0]),ap(curve,K.primes_above(61)[1]),ap(curve,K.primes_above(71)[0]),ap(curve,K.primes_above(71)[1]),ap(curve,K.primes_above(79)[0]),ap(curve,K.primes_above(79)[1]),ap(curve,K.primes_above(89)[0]),ap(curve,K.primes_above(89)[1]) 
       
Fractional ideal (21*a - 14) -3 -1 -1 -6 -3 -1 -4 -4 -4 5 8 11 -7 13 -14
4 -2 5 -10 -9 6
Fractional ideal (21*a - 14) -3 -1 -1 -6 -3 -1 -4 -4 -4 5 8 11 -7 13 -14
4 -2 5 -10 -9 6
Fractional ideal (21*a - 14) -3 -1 -1 -6 -3 -1 -4 -4 -4 5 8 11 -7 13 -14 4 -2 5 -10 -9 6
Fractional ideal (21*a - 14) -3 -1 -1 -6 -3 -1 -4 -4 -4 5 8 11 -7 13 -14 4 -2 5 -10 -9 6
N=539 cond=21*a - 7 aps=[-3,2,2,'?','?','?',6,6,2,2,-10,8,-10,8,2,8,-14,4,-8,-8,-4,-4,0,12] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  41
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  59
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
singular curve
current prime =  79
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x +
(-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1
tentative
singular curve
current prime =  11
current prime =  19
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  41
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  59
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  61
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  79
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  89
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + y = x^3 + a*x^2 + (a-11)*x + (-22*a+8) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
E539_3_1 = EllipticCurve(K,[a+1,a,1,a-11,-22*a+8]) #works for curve in [E539_3_1]: print curve.conductor(), ap(curve,K.primes_above(2)[0]),ap(curve,K.primes_above(3)[0]), ap(curve,K.primes_above(5)[0]), ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(31)[0]),ap(curve,K.primes_above(31)[1]),ap(curve,K.primes_above(41)[0]),ap(curve,K.primes_above(41)[1]),ap(curve,K.primes_above(59)[0]),ap(curve,K.primes_above(59)[1]),ap(curve,K.primes_above(61)[0]),ap(curve,K.primes_above(61)[1]),ap(curve,K.primes_above(71)[0]),ap(curve,K.primes_above(71)[1]),ap(curve,K.primes_above(79)[0]),ap(curve,K.primes_above(79)[1]),ap(curve,K.primes_above(89)[0]),ap(curve,K.primes_above(89)[1]) 
       
Fractional ideal (21*a - 7) -3 2 2 6 6 2 2 8 -10 8 -10 8 2 4 -14 -8 -8
-4 -4 12 0
Fractional ideal (21*a - 7) -3 2 2 6 6 2 2 8 -10 8 -10 8 2 4 -14 -8 -8 -4 -4 12 0
N=545 cond=21*a - 8 aps=[1,2,'?',2,2,0,0,6,2,2,-4,-8,-2,2,0,6,2,2,-10,-2,-8,-14,2,-6] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
singular curve
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 = x^3 + (a+1)*x^2 + (9*a-37)*x + (23*a-12)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-1)*x^2 + (8*a-13)*x
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-1)*x^2 +
(-32*a+52)*x + (31*a-50) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + a*x^2 + (a-2)*x over Number
Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + a*x^2 + (-4*a+8)*x + (3*a-2)
over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-a+2)*x + (18*a+12) over Number Field in a with defining polynomial x^2
- x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (-a+1)*x^2 +
(-6*a-3)*x + (5*a+3) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
singular curve
current prime =  61
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(11*a+1)*x + (73*a+36) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
current prime =  89
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2
- x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 +
(11*a+1)*x + (73*a+36) over Number Field in a with defining polynomial
x^2 - x - 1 tentative
current prime =  11
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  19
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  29
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  31
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  41
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  59
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 = x^3 + (a+1)*x^2 + (9*a-37)*x + (23*a-12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-1)*x^2 + (8*a-13)*x over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + (a+1)*x*y = x^3 + (-1)*x^2 + (-32*a+52)*x + (31*a-50) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + a*x^2 + (a-2)*x over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + x*y = x^3 + a*x^2 + (-4*a+8)*x + (3*a-2) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-a+2)*x + (18*a+12) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y + (a+1)*y = x^3 + (-a+1)*x^2 + (-6*a-3)*x + (5*a+3) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
current prime =  61
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
singular curve
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  71
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
current prime =  79
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (11*a+1)*x + (73*a+36) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
current prime =  89
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (-4*a-9)*x + (6*a-4) over Number Field in a with defining polynomial x^2 - x - 1 tentative
singular curve
Elliptic Curve defined by y^2 + (a+1)*x*y + a*y = x^3 + (a+1)*x^2 + (11*a+1)*x + (73*a+36) over Number Field in a with defining polynomial x^2 - x - 1 tentative
E545_1 = EllipticCurve(K,[a,-a+1,a+1,-6*a-3,5*a+3]) #E545_2 might be a valid Gal. conjugate E545_2 = EllipticCurve(K,[a+1,a+1,a,-4*a-9,6*a-4]) E545_3 = EllipticCurve(K,[a,-a+1,a+1,-a+2,18*a+12]) E545_4 = EllipticCurve(K,[a+1,a+1,a,11*a+1,73*a+36]) E545_5 = EllipticCurve(K,[1,a,0,-4*a+8,3*a-2]) E545_6 = EllipticCurve(K,[1,a,0,a-2,0]) E545_7 = EllipticCurve(K,[a+1,-1,0,8*a-13,0]) E545_8 = EllipticCurve(K,[0,a+1,0,9*a-37,23*a-12]) for curve in [E545_1,E545_2,E545_3,E545_4,E545_5,E545_6,E545_7,E545_8]: try: print curve.conductor(), ap(curve,K.primes_above(2)[0]),ap(curve,K.primes_above(3)[0]), ap(curve,K.primes_above(7)[0]), ap(curve,K.primes_above(11)[0]),ap(curve,K.primes_above(11)[1]), ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(31)[0]),ap(curve,K.primes_above(31)[1]),ap(curve,K.primes_above(41)[0]),ap(curve,K.primes_above(41)[1]),ap(curve,K.primes_above(59)[0]),ap(curve,K.primes_above(59)[1]),ap(curve,K.primes_above(61)[0]),ap(curve,K.primes_above(61)[1]),ap(curve,K.primes_above(71)[0]),ap(curve,K.primes_above(71)[1]),ap(curve,K.primes_above(79)[0]),ap(curve,K.primes_above(79)[1]),ap(curve,K.primes_above(89)[0]),ap(curve,K.primes_above(89)[1]) except ArithmeticError: pass 
       
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10
-6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) 1 2 2 0 2 6 0 2 2 -8 -4 2 -2 6 0 2 2 -2 -10
-14 -8 -6 2
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10
-6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) 1 2 2 0 2 6 0 2 2 -8 -4 2 -2 6 0 2 2 -2 -10
-14 -8 -6 2
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10
-6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10
-6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10
-6 6 -10 8 -6 -6
Fractional ideal (21*a - 8)
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10 -6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) 1 2 2 0 2 6 0 2 2 -8 -4 2 -2 6 0 2 2 -2 -10 -14 -8 -6 2
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10 -6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) 1 2 2 0 2 6 0 2 2 -8 -4 2 -2 6 0 2 2 -2 -10 -14 -8 -6 2
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10 -6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10 -6 6 -10 8 -6 -6
Fractional ideal (21*a - 8) -1 -2 -10 0 6 2 8 -6 -6 8 -4 6 6 6 0 -10 -10 -6 6 -10 8 -6 -6
Fractional ideal (21*a - 8)
aps=[1,2,2,2,0,0,6,2,2,-4,-8,-2,2,0,6,2,2,-10,-2,-8,-14,2,-6] 
       
N=549 cond=21*a - 12 aps=[3,'?',-1,3,3,-3,7,-2,1,-6,-5,3,6,-2,-14,6,'?','?',-5,-8,0,14,-10,-1] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
singular curve
current prime =  19
singular curve
current prime =  29
singular curve
current prime =  31
current prime =  41
singular curve
singular curve
singular curve
singular curve
current prime =  59
current prime =  61
current prime =  71
singular curve
singular curve
current prime =  79
singular curve
current prime =  89
singular curve
current prime =  11
singular curve
singular curve
current prime =  19
singular curve
current prime =  29
singular curve
current prime =  31
current prime =  41
singular curve
singular curve
singular curve
singular curve
current prime =  59
current prime =  61
current prime =  71
singular curve
singular curve
current prime =  79
singular curve
current prime =  89
singular curve
N=580 cond=22*a - 6 aps=['?',3,'?',-4,-4,1,7,2,'?','?',3,-10,-5,-12,-14,6,13,10,-6,2,-7,-14,-18,7] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
singular curve
current prime =  29
current prime =  31
singular curve
current prime =  41
singular curve
current prime =  59
current prime =  61
singular curve
current prime =  71
singular curve
singular curve
current prime =  79
current prime =  89
singular curve
current prime =  11
current prime =  19
singular curve
current prime =  29
current prime =  31
singular curve
current prime =  41
singular curve
current prime =  59
current prime =  61
singular curve
current prime =  71
singular curve
singular curve
current prime =  79
current prime =  89
singular curve
N=599 cond=a - 25 aps=[-1,2,2,11,5,-3,0,2,1,2,-1,-3,-5,6,3,-12,-12,1,7,-2,5,11,-15,-2] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
singular curve
singular curve
current prime =  19
singular curve
singular curve
singular curve
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
singular curve
singular curve
current prime =  29
singular curve
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  31
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  41
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  59
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  61
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
singular curve
current prime =  71
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  79
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  89
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  11
singular curve
singular curve
singular curve
current prime =  19
singular curve
singular curve
singular curve
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
singular curve
singular curve
current prime =  29
singular curve
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  31
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  41
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  59
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  61
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
singular curve
current prime =  71
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  79
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
current prime =  89
(a, 1, 0, -7*a - 5, -16*a - 11) tentative
E599_2_1 = EllipticCurve(K,[a, 1, 0, -7*a - 5, -16*a - 11]) #could be gal-conj. to our curve for curve in [E599_2_1]: try: print curve.conductor(), ap(curve,K.primes_above(2)[0]),ap(curve,K.primes_above(3)[0]), ap(curve,K.primes_above(5)[0]),ap(curve,K.primes_above(7)[0]), ap(curve,K.primes_above(11)[0]),ap(curve,K.primes_above(11)[1]), ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(31)[0]),ap(curve,K.primes_above(31)[1]),ap(curve,K.primes_above(41)[0]),ap(curve,K.primes_above(41)[1]),ap(curve,K.primes_above(59)[0]),ap(curve,K.primes_above(59)[1]),ap(curve,K.primes_above(61)[0]),ap(curve,K.primes_above(61)[1]),ap(curve,K.primes_above(71)[0]),ap(curve,K.primes_above(71)[1]),ap(curve,K.primes_above(79)[0]),ap(curve,K.primes_above(79)[1]),ap(curve,K.primes_above(89)[0]),ap(curve,K.primes_above(89)[1]) except ArithmeticError: pass 
       
Fractional ideal (a - 25) -1 2 2 11 -3 5 2 0 2 1 -3 -1 6 -5 -12 3 1 -12
-2 7 11 5 -2 -15
Fractional ideal (a - 25) -1 2 2 11 -3 5 2 0 2 1 -3 -1 6 -5 -12 3 1 -12 -2 7 11 5 -2 -15
####aps=[-1,2,2,11,5,-3,0,2,1,2,-1,-3,-5,6,3,-12,-12,1,7,-2,5,11,-15,-2] 
       
N=620 cond=2*a - 26 aps=['?',4,'?',-9,-3,1,5,6,7,-2,'?','?',-4,8,3,6,-14,-5,-5,14,10,-14,-9,10] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
current prime =  19
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
(1, 0, 1, a - 2, -2*a + 3) tentative
singular curve
current prime =  29
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
current prime =  31
current prime =  41
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
current prime =  59
singular curve
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
singular curve
singular curve
current prime =  61
singular curve
singular curve
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
singular curve
current prime =  71
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
current prime =  79
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
current prime =  89
singular curve
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
current prime =  11
current prime =  19
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
(1, 0, 1, a - 2, -2*a + 3) tentative
singular curve
current prime =  29
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
current prime =  31
current prime =  41
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
current prime =  59
singular curve
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
singular curve
singular curve
current prime =  61
singular curve
singular curve
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
singular curve
current prime =  71
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
current prime =  79
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
current prime =  89
singular curve
(a + 1, -1, 1, 5*a - 12, 7*a - 14) tentative
singular curve
E620_1 = EllipticCurve(K,[a + 1, -1, 1, 5*a - 12, 7*a - 14]) #E620_1 might be gal-conj to our curve E620_2 = EllipticCurve(K,[1, 0, 1, a - 2, -2*a + 3]) for curve in [E620_1,E620_2]: try: print curve.conductor(),ap(curve,K.primes_above(3)[0]),ap(curve,K.primes_above(7)[0]), ap(curve,K.primes_above(11)[0]),ap(curve,K.primes_above(11)[1]), ap(curve,K.primes_above(19)[0]),ap(curve,K.primes_above(19)[1]),ap(curve,K.primes_above(29)[0]),ap(curve,K.primes_above(29)[1]),ap(curve,K.primes_above(41)[0]),ap(curve,K.primes_above(41)[1]),ap(curve,K.primes_above(59)[0]),ap(curve,K.primes_above(59)[1]),ap(curve,K.primes_above(61)[0]),ap(curve,K.primes_above(61)[1]),ap(curve,K.primes_above(71)[0]),ap(curve,K.primes_above(71)[1]),ap(curve,K.primes_above(79)[0]),ap(curve,K.primes_above(79)[1]),ap(curve,K.primes_above(89)[0]),ap(curve,K.primes_above(89)[1]) except ArithmeticError: pass 
       
Fractional ideal (2*a - 26) 4 -9 1 -3 6 5 -2 7 8 -4 6 3 -5 -14 14 -5 -14
10 10 -9
Fractional ideal (2*a - 26) 4 5 -3 3 2 -7 6 9 -12 0 -6 -3 -1 -10 -6 9
-10 -10 -6 15
Fractional ideal (2*a - 26) 4 -9 1 -3 6 5 -2 7 8 -4 6 3 -5 -14 14 -5 -14 10 10 -9
Fractional ideal (2*a - 26) 4 5 -3 3 2 -7 6 9 -12 0 -6 -3 -1 -10 -6 9 -10 -10 -6 15
##aps=[4,-9,-3,1,5,6,7,-2,-4,8,3,6,-14,-5,-5,14,10,-14,-9,10] 
       
N=639 cond=3*a + 24 aps=[2,'?',1,2,3,0,-6,0,6,3,-2,-3,2,-9,10,-5,-6,-5,'?','?',10,4,10,-2] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
current prime =  19
current prime =  29
singular curve
current prime =  31
current prime =  41
singular curve
current prime =  59
current prime =  61
current prime =  71
current prime =  79
singular curve
singular curve
current prime =  89
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  11
singular curve
current prime =  19
current prime =  29
singular curve
current prime =  31
current prime =  41
singular curve
current prime =  59
current prime =  61
current prime =  71
current prime =  79
singular curve
singular curve
current prime =  89
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
E639_1_1 = EllipticCurve(K,[a + 1, a, a, 13*a - 25, 30*a - 49]) #probably not aplist1=[] aplist2=[] for i,curve in enumerate([E639_1_1]): aplist = eval('aplist%s'%(i+1)) #try: for pr in [2,3,5,7,11,19,29,31,41,59,61,71,79,89]: try: aplist.append(ap(curve,K.primes_above(pr)[0])) if len(K.primes_above(pr))>1: aplist.append(ap(curve,K.primes_above(pr)[1])) except ArithmeticError: aplist.append('?') if len(K.primes_above(pr))>1: aplist.append('?') print E639_1_1.conductor(), aplist1 
       
Fractional ideal (3*a + 24) [-1, '?', 4, -4, 0, 0, 6, 0, -6, 6, -2, 0,
-10, 0, 4, -2, -12, -2, -2, '?', '?', -14, 16, -2, 10]
Fractional ideal (3*a + 24) [-1, '?', 4, -4, 0, 0, 6, 0, -6, 6, -2, 0, -10, 0, 4, -2, -12, -2, -2, '?', '?', -14, 16, -2, 10]
##aps=[2,'?',1,2,3,0,-6,0,6,3,-2,-3,2,-9,10,-5,-6,-5,'?','?',10,4,10,-2] 
       
N=639 cond=3*a + 24 aps=[-1,'?',4,-4,0,0,6,0,-6,6,-2,0,-10,0,4,-2,-12,-2,'?','?',-14,16,-2,10] for i,pr in enumerate([11,19,29,31,41,59,61,71,79,89]): print 'current prime = ', pr Dopr = SievedSearch(pr) ap1 = aps[2*i+5] ap2 = aps[2*i+4] if type(ap1)!=str: if pr==11: ap1,ap2=ap2,ap1 Dopr(cond,ap1,ap2,aps) Dopr(cond,ap2,ap1,aps) 
       
current prime =  11
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  19
current prime =  29
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  31
singular curve
singular curve
singular curve
current prime =  41
singular curve
singular curve
singular curve
current prime =  59
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  61
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  71
current prime =  79
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  89
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
singular curve
current prime =  11
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  19
current prime =  29
singular curve
singular curve
singular curve
singular curve
singular curve
singular curve
current prime =  31
singular curve
singular curve
singular curve
current prime =  41
singular curve
singular curve
singular curve
current prime =  59
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  61
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  71
current prime =  79
singular curve
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
current prime =  89
(a + 1, a, a, 13*a - 25, 30*a - 49) tentative
singular curve
E639_1_1 = EllipticCurve(K,[a + 1, a, a, 13*a - 25, 30*a - 49]) #probably not aplist1=[] aplist2=[] for i,curve in enumerate([E639_1_1]): aplist = eval('aplist%s'%(i+1)) #try: for pr in [2,3,5,7,11,19,29,31,41,59,61,71,79,89]: try: aplist.append((ap(curve,K.primes_above(pr)[0]),pr)) if len(K.primes_above(pr))>1: aplist.append((ap(curve,K.primes_above(pr)[1]),pr)) except ArithmeticError: aplist.append(('?',pr)) if len(K.primes_above(pr))>1: aplist.append(('?',pr)) print E639_1_1.conductor(), aplist1 
       
Fractional ideal (3*a + 24) [(-1, 2), ('?', 3), (4, 5), (-4, 7), (0,
11), (0, 11), (6, 19), (0, 19), (-6, 29), (6, 29), (-2, 31), (0, 31),
(-10, 41), (0, 41), (4, 59), (-2, 59), (-12, 61), (-2, 61), (-2, 71),
('?', 71), ('?', 71), (-14, 79), (16, 79), (-2, 89), (10, 89)]
Fractional ideal (3*a + 24) [(-1, 2), ('?', 3), (4, 5), (-4, 7), (0, 11), (0, 11), (6, 19), (0, 19), (-6, 29), (6, 29), (-2, 31), (0, 31), (-10, 41), (0, 41), (4, 59), (-2, 59), (-12, 61), (-2, 61), (-2, 71), ('?', 71), ('?', 71), (-14, 79), (16, 79), (-2, 89), (10, 89)]
aps=[-1,'?',4,-4,0,0,6,0,-6,6,-2,0,-10,0,4,-2,-12,-2,'?','?',-14,16,-2,10] 
       
 
       
 
       
 
       
 
       

###############################################################################################################

Do = SievedSearch(19) 
       
E2=EllipticCurve(K,[1, -a + 1, a + 1, -11*a - 6, -18*a - 11]) E2.conductor().norm() 
       
Traceback (click to the left of this block for traceback)
...
ArithmeticError: Invariants [1, -a + 1, a + 1, -11*a - 6, -18*a - 11]
define a singular curve.
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_5.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("RTI9RWxsaXB0aWNDdXJ2ZShLLFsxLCAtYSArIDEsIGEgKyAxLCAtMTEqYSAtIDYsIC0xOCphIC0gMTFdKQpFMi5jb25kdWN0b3IoKS5ub3JtKCk="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/private/var/folders/+3/+3+bRf7bGhaMk-2vvrWEvU+++TI/-Tmp-/tmpUXECyo/___code___.py", line 3, in <module>
    E2=EllipticCurve(K,[_sage_const_1 , -a + _sage_const_1 , a + _sage_const_1 , -_sage_const_11 *a - _sage_const_6 , -_sage_const_18 *a - _sage_const_11 ])
  File "/Applications/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/constructor.py", line 281, in EllipticCurve
    return ell_number_field.EllipticCurve_number_field(x, y)
  File "/Applications/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/ell_number_field.py", line 163, in __init__
    EllipticCurve_field.__init__(self, [field(x) for x in ainvs])
  File "/Applications/sage/local/lib/python2.6/site-packages/sage/schemes/elliptic_curves/ell_generic.py", line 164, in __init__
    "Invariants %s define a singular curve."%ainvs
ArithmeticError: Invariants [1, -a + 1, a + 1, -11*a - 6, -18*a - 11] define a singular curve.
Do29 = SievedSearch(29) 
       
Do29(-12*a+5,7,5,[-5,3,1,-4,7,5,2,1,-6,-3,8,-4,-10,-11,12,10,12,0,4,4]) 
       
WARNING: Output truncated!  
full_output.txt



(1L, 'of', 392L, '---', 29L)
(2L, 'of', 392L, '---', 29L)
(3L, 'of', 392L, '---', 29L)
(4L, 'of', 392L, '---', 29L)
(5L, 'of', 392L, '---', 29L)
(6L, 'of', 392L, '---', 29L)
(7L, 'of', 392L, '---', 29L)
(8L, 'of', 392L, '---', 29L)
(9L, 'of', 392L, '---', 29L)
(10L, 'of', 392L, '---', 29L)
(11L, 'of', 392L, '---', 29L)
(12L, 'of', 392L, '---', 29L)
(13L, 'of', 392L, '---', 29L)
(14L, 'of', 392L, '---', 29L)
(15L, 'of', 392L, '---', 29L)
(16L, 'of', 392L, '---', 29L)
(17L, 'of', 392L, '---', 29L)
(18L, 'of', 392L, '---', 29L)
(19L, 'of', 392L, '---', 29L)
(20L, 'of', 392L, '---', 29L)
(21L, 'of', 392L, '---', 29L)
(22L, 'of', 392L, '---', 29L)
(23L, 'of', 392L, '---', 29L)
(24L, 'of', 392L, '---', 29L)
(25L, 'of', 392L, '---', 29L)
(26L, 'of', 392L, '---', 29L)
(27L, 'of', 392L, '---', 29L)
(28L, 'of', 392L, '---', 29L)
(29L, 'of', 392L, '---', 29L)
(30L, 'of', 392L, '---', 29L)
(31L, 'of', 392L, '---', 29L)
(32L, 'of', 392L, '---', 29L)
(33L, 'of', 392L, '---', 29L)
(34L, 'of', 392L, '---', 29L)
(35L, 'of', 392L, '---', 29L)
(36L, 'of', 392L, '---', 29L)
(37L, 'of', 392L, '---', 29L)
(38L, 'of', 392L, '---', 29L)
(39L, 'of', 392L, '---', 29L)
(40L, 'of', 392L, '---', 29L)
(41L, 'of', 392L, '---', 29L)
(42L, 'of', 392L, '---', 29L)
(43L, 'of', 392L, '---', 29L)
(44L, 'of', 392L, '---', 29L)
(45L, 'of', 392L, '---', 29L)
(46L, 'of', 392L, '---', 29L)
(47L, 'of', 392L, '---', 29L)
(48L, 'of', 392L, '---', 29L)
(49L, 'of', 392L, '---', 29L)
(50L, 'of', 392L, '---', 29L)
(51L, 'of', 392L, '---', 29L)
(52L, 'of', 392L, '---', 29L)
(53L, 'of', 392L, '---', 29L)
(54L, 'of', 392L, '---', 29L)
(55L, 'of', 392L, '---', 29L)
(56L, 'of', 392L, '---', 29L)
(57L, 'of', 392L, '---', 29L)
(58L, 'of', 392L, '---', 29L)
(59L, 'of', 392L, '---', 29L)

...

(333L, 'of', 392L, '---', 29L)
(334L, 'of', 392L, '---', 29L)
(335L, 'of', 392L, '---', 29L)
(336L, 'of', 392L, '---', 29L)
(337L, 'of', 392L, '---', 29L)
(338L, 'of', 392L, '---', 29L)
(339L, 'of', 392L, '---', 29L)
(340L, 'of', 392L, '---', 29L)
(341L, 'of', 392L, '---', 29L)
(342L, 'of', 392L, '---', 29L)
(343L, 'of', 392L, '---', 29L)
(344L, 'of', 392L, '---', 29L)
(345L, 'of', 392L, '---', 29L)
(346L, 'of', 392L, '---', 29L)
(347L, 'of', 392L, '---', 29L)
(348L, 'of', 392L, '---', 29L)
(349L, 'of', 392L, '---', 29L)
(350L, 'of', 392L, '---', 29L)
(351L, 'of', 392L, '---', 29L)
(352L, 'of', 392L, '---', 29L)
(353L, 'of', 392L, '---', 29L)
(354L, 'of', 392L, '---', 29L)
(355L, 'of', 392L, '---', 29L)
(356L, 'of', 392L, '---', 29L)
(357L, 'of', 392L, '---', 29L)
(358L, 'of', 392L, '---', 29L)
(359L, 'of', 392L, '---', 29L)
(360L, 'of', 392L, '---', 29L)
(361L, 'of', 392L, '---', 29L)
(362L, 'of', 392L, '---', 29L)
(363L, 'of', 392L, '---', 29L)
(364L, 'of', 392L, '---', 29L)
(365L, 'of', 392L, '---', 29L)
(366L, 'of', 392L, '---', 29L)
(367L, 'of', 392L, '---', 29L)
(368L, 'of', 392L, '---', 29L)
(369L, 'of', 392L, '---', 29L)
(370L, 'of', 392L, '---', 29L)
(371L, 'of', 392L, '---', 29L)
(372L, 'of', 392L, '---', 29L)
(373L, 'of', 392L, '---', 29L)
(374L, 'of', 392L, '---', 29L)
(375L, 'of', 392L, '---', 29L)
(376L, 'of', 392L, '---', 29L)
(377L, 'of', 392L, '---', 29L)
(378L, 'of', 392L, '---', 29L)
(379L, 'of', 392L, '---', 29L)
(380L, 'of', 392L, '---', 29L)
(381L, 'of', 392L, '---', 29L)
(382L, 'of', 392L, '---', 29L)
(383L, 'of', 392L, '---', 29L)
(384L, 'of', 392L, '---', 29L)
(385L, 'of', 392L, '---', 29L)
(386L, 'of', 392L, '---', 29L)
(387L, 'of', 392L, '---', 29L)
(388L, 'of', 392L, '---', 29L)
(389L, 'of', 392L, '---', 29L)
(390L, 'of', 392L, '---', 29L)
(391L, 'of', 392L, '---', 29L)
(392L, 'of', 392L, '---', 29L)
WARNING: Output truncated!  
full_output.txt



(1L, 'of', 392L, '---', 29L)
(2L, 'of', 392L, '---', 29L)
(3L, 'of', 392L, '---', 29L)
(4L, 'of', 392L, '---', 29L)
(5L, 'of', 392L, '---', 29L)
(6L, 'of', 392L, '---', 29L)
(7L, 'of', 392L, '---', 29L)
(8L, 'of', 392L, '---', 29L)
(9L, 'of', 392L, '---', 29L)
(10L, 'of', 392L, '---', 29L)
(11L, 'of', 392L, '---', 29L)
(12L, 'of', 392L, '---', 29L)
(13L, 'of', 392L, '---', 29L)
(14L, 'of', 392L, '---', 29L)
(15L, 'of', 392L, '---', 29L)
(16L, 'of', 392L, '---', 29L)
(17L, 'of', 392L, '---', 29L)
(18L, 'of', 392L, '---', 29L)
(19L, 'of', 392L, '---', 29L)
(20L, 'of', 392L, '---', 29L)
(21L, 'of', 392L, '---', 29L)
(22L, 'of', 392L, '---', 29L)
(23L, 'of', 392L, '---', 29L)
(24L, 'of', 392L, '---', 29L)
(25L, 'of', 392L, '---', 29L)
(26L, 'of', 392L, '---', 29L)
(27L, 'of', 392L, '---', 29L)
(28L, 'of', 392L, '---', 29L)
(29L, 'of', 392L, '---', 29L)
(30L, 'of', 392L, '---', 29L)
(31L, 'of', 392L, '---', 29L)
(32L, 'of', 392L, '---', 29L)
(33L, 'of', 392L, '---', 29L)
(34L, 'of', 392L, '---', 29L)
(35L, 'of', 392L, '---', 29L)
(36L, 'of', 392L, '---', 29L)
(37L, 'of', 392L, '---', 29L)
(38L, 'of', 392L, '---', 29L)
(39L, 'of', 392L, '---', 29L)
(40L, 'of', 392L, '---', 29L)
(41L, 'of', 392L, '---', 29L)
(42L, 'of', 392L, '---', 29L)
(43L, 'of', 392L, '---', 29L)
(44L, 'of', 392L, '---', 29L)
(45L, 'of', 392L, '---', 29L)
(46L, 'of', 392L, '---', 29L)
(47L, 'of', 392L, '---', 29L)
(48L, 'of', 392L, '---', 29L)
(49L, 'of', 392L, '---', 29L)
(50L, 'of', 392L, '---', 29L)
(51L, 'of', 392L, '---', 29L)
(52L, 'of', 392L, '---', 29L)
(53L, 'of', 392L, '---', 29L)
(54L, 'of', 392L, '---', 29L)
(55L, 'of', 392L, '---', 29L)
(56L, 'of', 392L, '---', 29L)
(57L, 'of', 392L, '---', 29L)
(58L, 'of', 392L, '---', 29L)
(59L, 'of', 392L, '---', 29L)

...

(333L, 'of', 392L, '---', 29L)
(334L, 'of', 392L, '---', 29L)
(335L, 'of', 392L, '---', 29L)
(336L, 'of', 392L, '---', 29L)
(337L, 'of', 392L, '---', 29L)
(338L, 'of', 392L, '---', 29L)
(339L, 'of', 392L, '---', 29L)
(340L, 'of', 392L, '---', 29L)
(341L, 'of', 392L, '---', 29L)
(342L, 'of', 392L, '---', 29L)
(343L, 'of', 392L, '---', 29L)
(344L, 'of', 392L, '---', 29L)
(345L, 'of', 392L, '---', 29L)
(346L, 'of', 392L, '---', 29L)
(347L, 'of', 392L, '---', 29L)
(348L, 'of', 392L, '---', 29L)
(349L, 'of', 392L, '---', 29L)
(350L, 'of', 392L, '---', 29L)
(351L, 'of', 392L, '---', 29L)
(352L, 'of', 392L, '---', 29L)
(353L, 'of', 392L, '---', 29L)
(354L, 'of', 392L, '---', 29L)
(355L, 'of', 392L, '---', 29L)
(356L, 'of', 392L, '---', 29L)
(357L, 'of', 392L, '---', 29L)
(358L, 'of', 392L, '---', 29L)
(359L, 'of', 392L, '---', 29L)
(360L, 'of', 392L, '---', 29L)
(361L, 'of', 392L, '---', 29L)
(362L, 'of', 392L, '---', 29L)
(363L, 'of', 392L, '---', 29L)
(364L, 'of', 392L, '---', 29L)
(365L, 'of', 392L, '---', 29L)
(366L, 'of', 392L, '---', 29L)
(367L, 'of', 392L, '---', 29L)
(368L, 'of', 392L, '---', 29L)
(369L, 'of', 392L, '---', 29L)
(370L, 'of', 392L, '---', 29L)
(371L, 'of', 392L, '---', 29L)
(372L, 'of', 392L, '---', 29L)
(373L, 'of', 392L, '---', 29L)
(374L, 'of', 392L, '---', 29L)
(375L, 'of', 392L, '---', 29L)
(376L, 'of', 392L, '---', 29L)
(377L, 'of', 392L, '---', 29L)
(378L, 'of', 392L, '---', 29L)
(379L, 'of', 392L, '---', 29L)
(380L, 'of', 392L, '---', 29L)
(381L, 'of', 392L, '---', 29L)
(382L, 'of', 392L, '---', 29L)
(383L, 'of', 392L, '---', 29L)
(384L, 'of', 392L, '---', 29L)
(385L, 'of', 392L, '---', 29L)
(386L, 'of', 392L, '---', 29L)
(387L, 'of', 392L, '---', 29L)
(388L, 'of', 392L, '---', 29L)
(389L, 'of', 392L, '---', 29L)
(390L, 'of', 392L, '---', 29L)
(391L, 'of', 392L, '---', 29L)
(392L, 'of', 392L, '---', 29L)
E1 = EllipticCurve(K,[a+1,1,a,-15*a-11,-50*a-33]) 
       
E1.is_isogenous(F1) 
       
False
False
F1 = canonical_model(EllipticCurve(K,[a+1,1,a+1,-5*a-8,-14*a-15]).global_minimal_model()) 
       
is_galois_conjugate_upto_isogeny(E1,F1) 
       
True
True
def is_galois_conjugate_upto_isogeny(E,F): #eainv = E.global_minimal_model().a_invariants() fainv = [K(t) for t in F.global_minimal_model().a_invariants()] fainv_conj = [] for i in range(5): fainv_conj.append(fainv[i][0]+fainv[i][1]*(1-a)) FC = EllipticCurve(K,fainv_conj).global_minimal_model() #if eainv == FC.a_invariants(): # return True if E.is_isogenous(FC): return True else: return False 
       
E = EllipticCurve(K,[a,-1,0,-6*a-1,-7*a-1]) 
       
F = EllipticCurve(K,[a,-1,0,13*a+18,12*a+18]) 
       
F.change_weierstrass_model([1,5+11*a,9*a,11*a+4]).a_invariants() 
       
(19*a, -57*a - 76, 38*a + 19, 156*a + 93, 4919*a + 3019)
(19*a, -57*a - 76, 38*a + 19, 156*a + 93, 4919*a + 3019)
G=EllipticCurve(K,[0,0,0,4*a+17, 17*a + 17]) 
       
-26%19 
       
12
12
mod(13*90,19) 
       
11
11
timeit('Do(-12*a + 5,-5,3,[-5,3,1,-4,7,5,2,1,-6,-3,8,-4,-10,-11,12,10,12,0,4,4])') 
       
5 loops, best of 3: 41.1 ms per loop
5 loops, best of 3: 41.1 ms per loop
Do(-6*a-29,2,-5,['?','?',-5,2,1,-1,8,-4,-2,0,-3,0,-5,-10,8,3,4,4,'?','?']) 
       
WARNING: Output truncated!  
full_output.txt



0 0 0 0 0 0 17 4 17 17
0 0 0 0 0 0 17 4 -2 17
0 0 0 0 0 0 17 4 -2 -2
0 0 0 0 0 0 17 4 17 -2
0 0 0 0 0 0 -2 4 17 17
0 0 0 0 0 0 -2 4 -2 17
0 0 0 0 0 0 -2 4 -2 -2
0 0 0 0 0 0 -2 4 17 -2
0 0 0 0 0 0 -2 -15 17 17
0 0 0 0 0 0 -2 -15 -2 17
0 0 0 0 0 0 -2 -15 -2 -2
0 0 0 0 0 0 -2 -15 17 -2
0 0 0 0 0 0 17 -15 17 17
0 0 0 0 0 0 17 -15 -2 17
0 0 0 0 0 0 17 -15 -2 -2
0 0 0 0 0 0 17 -15 17 -2
0 0 0 0 1 0 17 4 12 17
0 0 0 0 1 0 17 4 -7 17
0 0 0 0 1 0 17 4 -7 -2
0 0 0 0 1 0 17 4 12 -2
0 0 0 0 1 0 -2 4 12 17
0 0 0 0 1 0 -2 4 -7 17
0 0 0 0 1 0 -2 4 -7 -2
0 0 0 0 1 0 -2 4 12 -2
0 0 0 0 1 0 -2 -15 12 17
0 0 0 0 1 0 -2 -15 -7 17
0 0 0 0 1 0 -2 -15 -7 -2
0 0 0 0 1 0 -2 -15 12 -2
0 0 0 0 1 0 17 -15 12 17
0 0 0 0 1 0 17 -15 -7 17
0 0 0 0 1 0 17 -15 -7 -2
0 0 0 0 1 0 17 -15 12 -2
0 0 0 0 1 1 17 4 7 2
0 0 0 0 1 1 17 4 -12 2
0 0 0 0 1 1 17 4 -12 -17
0 0 0 0 1 1 17 4 7 -17
0 0 0 0 1 1 -2 4 7 2
0 0 0 0 1 1 -2 4 -12 2
0 0 0 0 1 1 -2 4 -12 -17
0 0 0 0 1 1 -2 4 7 -17
0 0 0 0 1 1 -2 -15 7 2
0 0 0 0 1 1 -2 -15 -12 2
0 0 0 0 1 1 -2 -15 -12 -17
0 0 0 0 1 1 -2 -15 7 -17
0 0 0 0 1 1 17 -15 7 2
0 0 0 0 1 1 17 -15 -12 2
0 0 0 0 1 1 17 -15 -12 -17
0 0 0 0 1 1 17 -15 7 -17
0 0 0 0 0 1 17 4 12 12
0 0 0 0 0 1 17 4 -7 12
0 0 0 0 0 1 17 4 -7 -7
0 0 0 0 0 1 17 4 12 -7
0 0 0 0 0 1 -2 4 12 12
0 0 0 0 0 1 -2 4 -7 12
0 0 0 0 0 1 -2 4 -7 -7
0 0 0 0 0 1 -2 4 12 -7
0 0 0 0 0 1 -2 -15 12 12
0 0 0 0 0 1 -2 -15 -7 12
0 0 0 0 0 1 -2 -15 -7 -7

...

0 1 -1 -1 1 0 5 10 -15 13
0 1 -1 -1 1 0 5 10 -15 -6
0 1 -1 -1 1 0 5 10 4 -6
0 1 -1 -1 1 0 -14 10 4 13
0 1 -1 -1 1 0 -14 10 -15 13
0 1 -1 -1 1 0 -14 10 -15 -6
0 1 -1 -1 1 0 -14 10 4 -6
0 1 -1 -1 1 0 -14 -9 4 13
0 1 -1 -1 1 0 -14 -9 -15 13
0 1 -1 -1 1 0 -14 -9 -15 -6
0 1 -1 -1 1 0 -14 -9 4 -6
0 1 -1 -1 1 0 5 -9 4 13
0 1 -1 -1 1 0 5 -9 -15 13
0 1 -1 -1 1 0 5 -9 -15 -6
0 1 -1 -1 1 0 5 -9 4 -6
0 1 -1 -1 1 1 5 0 18 17
0 1 -1 -1 1 1 5 0 -1 17
0 1 -1 -1 1 1 5 0 -1 -2
0 1 -1 -1 1 1 5 0 18 -2
0 1 -1 -1 1 1 -14 0 18 17
0 1 -1 -1 1 1 -14 0 -1 17
0 1 -1 -1 1 1 -14 0 -1 -2
0 1 -1 -1 1 1 -14 0 18 -2
0 1 -1 -1 1 1 -14 -19 18 17
0 1 -1 -1 1 1 -14 -19 -1 17
0 1 -1 -1 1 1 -14 -19 -1 -2
0 1 -1 -1 1 1 -14 -19 18 -2
0 1 -1 -1 1 1 5 -19 18 17
0 1 -1 -1 1 1 5 -19 -1 17
0 1 -1 -1 1 1 5 -19 -1 -2
0 1 -1 -1 1 1 5 -19 18 -2
0 1 -1 -1 0 1 15 10 4 8
0 1 -1 -1 0 1 15 10 -15 8
0 1 -1 -1 0 1 15 10 -15 -11
0 1 -1 -1 0 1 15 10 4 -11
0 1 -1 -1 0 1 -4 10 4 8
0 1 -1 -1 0 1 -4 10 -15 8
0 1 -1 -1 0 1 -4 10 -15 -11
0 1 -1 -1 0 1 -4 10 4 -11
0 1 -1 -1 0 1 -4 -9 4 8
0 1 -1 -1 0 1 -4 -9 -15 8
0 1 -1 -1 0 1 -4 -9 -15 -11
0 1 -1 -1 0 1 -4 -9 4 -11
0 1 -1 -1 0 1 15 -9 4 8
0 1 -1 -1 0 1 15 -9 -15 8
0 1 -1 -1 0 1 15 -9 -15 -11
0 1 -1 -1 0 1 15 -9 4 -11
0 1 -1 0 0 0 18 13 18 12
0 1 -1 0 0 0 18 13 -1 12
0 1 -1 0 0 0 18 13 -1 -7
0 1 -1 0 0 0 18 13 18 -7
0 1 -1 0 0 0 -1 13 18 12
0 1 -1 0 0 0 -1 13 -1 12
0 1 -1 0 0 0 -1 13 -1 -7
0 1 -1 0 0 0 -1 13 18 -7
0 1 -1 0 0 0 -1 -6 18 12
0 1 -1 0 0 0 -1 -6 -1 12
0 1 -1 0 0 0 -1 -6 -1 -7
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-1)*x^2 + (-6*a-1)*x +
(-7*a-1) over Number Field in a with defining polynomial x^2 - x - 1
tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-1)*x^2 + (-6*a-1)*x +
(-7*a-1) over Number Field in a with defining polynomial x^2 - x - 1
a_ps match
WARNING: Output truncated!  
full_output.txt



0 0 0 0 0 0 17 4 17 17
0 0 0 0 0 0 17 4 -2 17
0 0 0 0 0 0 17 4 -2 -2
0 0 0 0 0 0 17 4 17 -2
0 0 0 0 0 0 -2 4 17 17
0 0 0 0 0 0 -2 4 -2 17
0 0 0 0 0 0 -2 4 -2 -2
0 0 0 0 0 0 -2 4 17 -2
0 0 0 0 0 0 -2 -15 17 17
0 0 0 0 0 0 -2 -15 -2 17
0 0 0 0 0 0 -2 -15 -2 -2
0 0 0 0 0 0 -2 -15 17 -2
0 0 0 0 0 0 17 -15 17 17
0 0 0 0 0 0 17 -15 -2 17
0 0 0 0 0 0 17 -15 -2 -2
0 0 0 0 0 0 17 -15 17 -2
0 0 0 0 1 0 17 4 12 17
0 0 0 0 1 0 17 4 -7 17
0 0 0 0 1 0 17 4 -7 -2
0 0 0 0 1 0 17 4 12 -2
0 0 0 0 1 0 -2 4 12 17
0 0 0 0 1 0 -2 4 -7 17
0 0 0 0 1 0 -2 4 -7 -2
0 0 0 0 1 0 -2 4 12 -2
0 0 0 0 1 0 -2 -15 12 17
0 0 0 0 1 0 -2 -15 -7 17
0 0 0 0 1 0 -2 -15 -7 -2
0 0 0 0 1 0 -2 -15 12 -2
0 0 0 0 1 0 17 -15 12 17
0 0 0 0 1 0 17 -15 -7 17
0 0 0 0 1 0 17 -15 -7 -2
0 0 0 0 1 0 17 -15 12 -2
0 0 0 0 1 1 17 4 7 2
0 0 0 0 1 1 17 4 -12 2
0 0 0 0 1 1 17 4 -12 -17
0 0 0 0 1 1 17 4 7 -17
0 0 0 0 1 1 -2 4 7 2
0 0 0 0 1 1 -2 4 -12 2
0 0 0 0 1 1 -2 4 -12 -17
0 0 0 0 1 1 -2 4 7 -17
0 0 0 0 1 1 -2 -15 7 2
0 0 0 0 1 1 -2 -15 -12 2
0 0 0 0 1 1 -2 -15 -12 -17
0 0 0 0 1 1 -2 -15 7 -17
0 0 0 0 1 1 17 -15 7 2
0 0 0 0 1 1 17 -15 -12 2
0 0 0 0 1 1 17 -15 -12 -17
0 0 0 0 1 1 17 -15 7 -17
0 0 0 0 0 1 17 4 12 12
0 0 0 0 0 1 17 4 -7 12
0 0 0 0 0 1 17 4 -7 -7
0 0 0 0 0 1 17 4 12 -7
0 0 0 0 0 1 -2 4 12 12
0 0 0 0 0 1 -2 4 -7 12
0 0 0 0 0 1 -2 4 -7 -7
0 0 0 0 0 1 -2 4 12 -7
0 0 0 0 0 1 -2 -15 12 12
0 0 0 0 0 1 -2 -15 -7 12
0 0 0 0 0 1 -2 -15 -7 -7

...

0 1 -1 -1 1 0 5 10 -15 13
0 1 -1 -1 1 0 5 10 -15 -6
0 1 -1 -1 1 0 5 10 4 -6
0 1 -1 -1 1 0 -14 10 4 13
0 1 -1 -1 1 0 -14 10 -15 13
0 1 -1 -1 1 0 -14 10 -15 -6
0 1 -1 -1 1 0 -14 10 4 -6
0 1 -1 -1 1 0 -14 -9 4 13
0 1 -1 -1 1 0 -14 -9 -15 13
0 1 -1 -1 1 0 -14 -9 -15 -6
0 1 -1 -1 1 0 -14 -9 4 -6
0 1 -1 -1 1 0 5 -9 4 13
0 1 -1 -1 1 0 5 -9 -15 13
0 1 -1 -1 1 0 5 -9 -15 -6
0 1 -1 -1 1 0 5 -9 4 -6
0 1 -1 -1 1 1 5 0 18 17
0 1 -1 -1 1 1 5 0 -1 17
0 1 -1 -1 1 1 5 0 -1 -2
0 1 -1 -1 1 1 5 0 18 -2
0 1 -1 -1 1 1 -14 0 18 17
0 1 -1 -1 1 1 -14 0 -1 17
0 1 -1 -1 1 1 -14 0 -1 -2
0 1 -1 -1 1 1 -14 0 18 -2
0 1 -1 -1 1 1 -14 -19 18 17
0 1 -1 -1 1 1 -14 -19 -1 17
0 1 -1 -1 1 1 -14 -19 -1 -2
0 1 -1 -1 1 1 -14 -19 18 -2
0 1 -1 -1 1 1 5 -19 18 17
0 1 -1 -1 1 1 5 -19 -1 17
0 1 -1 -1 1 1 5 -19 -1 -2
0 1 -1 -1 1 1 5 -19 18 -2
0 1 -1 -1 0 1 15 10 4 8
0 1 -1 -1 0 1 15 10 -15 8
0 1 -1 -1 0 1 15 10 -15 -11
0 1 -1 -1 0 1 15 10 4 -11
0 1 -1 -1 0 1 -4 10 4 8
0 1 -1 -1 0 1 -4 10 -15 8
0 1 -1 -1 0 1 -4 10 -15 -11
0 1 -1 -1 0 1 -4 10 4 -11
0 1 -1 -1 0 1 -4 -9 4 8
0 1 -1 -1 0 1 -4 -9 -15 8
0 1 -1 -1 0 1 -4 -9 -15 -11
0 1 -1 -1 0 1 -4 -9 4 -11
0 1 -1 -1 0 1 15 -9 4 8
0 1 -1 -1 0 1 15 -9 -15 8
0 1 -1 -1 0 1 15 -9 -15 -11
0 1 -1 -1 0 1 15 -9 4 -11
0 1 -1 0 0 0 18 13 18 12
0 1 -1 0 0 0 18 13 -1 12
0 1 -1 0 0 0 18 13 -1 -7
0 1 -1 0 0 0 18 13 18 -7
0 1 -1 0 0 0 -1 13 18 12
0 1 -1 0 0 0 -1 13 -1 12
0 1 -1 0 0 0 -1 13 -1 -7
0 1 -1 0 0 0 -1 13 18 -7
0 1 -1 0 0 0 -1 -6 18 12
0 1 -1 0 0 0 -1 -6 -1 12
0 1 -1 0 0 0 -1 -6 -1 -7
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-1)*x^2 + (-6*a-1)*x + (-7*a-1) over Number Field in a with defining polynomial x^2 - x - 1 tentative
Elliptic Curve defined by y^2 + a*x*y = x^3 + (-1)*x^2 + (-6*a-1)*x + (-7*a-1) over Number Field in a with defining polynomial x^2 - x - 1 a_ps match
E=EllipticCurve(K,[5*a,25+22*a]) E.change_ring(K.residue_field(K.primes_above(29)[1])).trace_of_frobenius() 
       
5
5
len(E.points()) 
       
Traceback (click to the left of this block for traceback)
...
AttributeError: 'EllipticCurve_number_field' object has no attribute
'points'
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_81.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("bGVuKEUucG9pbnRzKCkp"),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/private/var/folders/+3/+3+bRf7bGhaMk-2vvrWEvU+++TI/-Tmp-/tmpCVW6kQ/___code___.py", line 2, in <module>
    exec compile(u'len(E.points())
  File "", line 1, in <module>
    
  File "parent.pyx", line 705, in sage.structure.parent.Parent.__getattr__ (sage/structure/parent.c:5410)
  File "parent.pyx", line 172, in sage.structure.parent.raise_attribute_error (sage/structure/parent.c:2638)
AttributeError: 'EllipticCurve_number_field' object has no attribute 'points'
do(-12*a+5 
       
a = K.gen(0) do(-12*a+5,5,7,[-5,3,1,-4,7,5,2,1,-6,-3,8,-4,-10,-11,12,10,12,0,4,4],29) 
       
<type 'long'>
0
0
(0L, 0L, 201581136L)
<type 'long'>
0
0
(0L, 0L, 201581136L)
(18+9*6)%29 
       
14
14
from psage.modform.hilbert.sqrt5.sqrt5 import F as K 
       
a = K.gen(0) do(-12*a+5,5,7,[-5,3,1,-4,7,5,2,1,-6,-3,8,-4,-10,-11,12,10,12,0,4,4],29) 
       
(0L, 0L, 201581008L)
(0L, 0L, 201581008L)
%cython from libc.stdint cimport int32_t, uint32_t cdef mpz_t x,y,z mpz_init(x) mpz_init(y) mpz_init(z) mpz_set_ui(x,2) mpz_set_ui(y,1213) mpz_mul_ui(z,x,-4) mpz_sub(z,z,y) cdef int32_t c c = mpz_get_si(z) print(c) mpz_clear(x) mpz_clear(y) mpz_clear(z) 
2^32 
       
4294967296
4294967296
%cython from libc.stdint cimport int32_t, uint32_t cdef uint32_t quick_disc_norm(uint32_t a_11, uint32_t a_12, uint32_t a_21, uint32_t a_22, uint32_t a_31, uint32_t a_32, uint32_t a_41, uint32_t a_42, uint32_t a_61, uint32_t a_62): cdef mpz_t A,B,C,D,E,F,G,H,I,J, b21,b22,b41,b42,b61,b62,b81,b82, temp1,temp2,temp3, sim1,sim2,sim3,sim4,sim5,sim6, disc1,disc2 mpz_init(A) mpz_init(B) mpz_init(C) mpz_init(D) mpz_init(E) mpz_init(F) mpz_init(G) mpz_init(H) mpz_init(I) mpz_init(J) mpz_init(b21) mpz_init(b22) mpz_init(b41) mpz_init(b42) mpz_init(b61) mpz_init(b62) mpz_init(b81) mpz_init(b82) mpz_init(sim1) mpz_init(sim2) mpz_init(sim3) mpz_init(sim4) mpz_init(sim5) mpz_init(sim6) mpz_init(temp1) mpz_init(temp2) mpz_init(temp3) mpz_init(disc1) mpz_init(disc2) mpz_set_ui(A,a_11) mpz_set_ui(B,a_12) mpz_set_ui(C,a_21) mpz_set_ui(D,a_22) mpz_set_ui(E,a_31) mpz_set_ui(F,a_32) mpz_set_ui(G,a_41) mpz_set_ui(H,a_42) mpz_set_ui(I,a_61) mpz_set_ui(J,a_62) #b21# mpz_mul(temp1,A,A) mpz_mul(temp2,B,B) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,C,4u) mpz_add(b21,temp1,temp2) #b22# mpz_mul(temp1,A,B) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,B,B) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,D,4) mpz_add(b22,temp1,temp2) #b41# mpz_mul(temp1,A,E) mpz_mul(temp2,B,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,G,2u) mpz_add(b41,temp1,temp2) #b42# mpz_mul(temp1,A,F) mpz_mul(temp2,B,E) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,B,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,H,2u) mpz_add(b42,temp1,temp2) #b61# mpz_mul(temp1,E,E) mpz_mul(temp2,F,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,I,4u) mpz_add(b61,temp1,temp2) #b62# mpz_mul(temp1,E,F) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,F,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,J,4) mpz_add(b62,temp1,temp2) #b8similarities# mpz_mul(temp1,A,F) mpz_mul(sim1,temp1,H) mpz_mul(temp1,B,E) mpz_mul(sim2,temp1,H) mpz_mul(temp1,B,F) mpz_mul(sim3,temp1,G) mpz_mul(temp1,B,F) mpz_mul(sim4,temp1,H) mpz_mul(temp1,C,F) mpz_mul(sim5,temp1,F) mpz_mul(temp1,D,F) mpz_mul(sim6,temp1,F) #b81# mpz_mul(temp1,b21,I) mpz_mul(temp2,b22,J) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,A,E) mpz_mul(temp2,temp2,G) mpz_sub(temp1,temp1,temp2) mpz_sub(temp1,temp1,sim1) mpz_sub(temp1,temp1,sim2) mpz_sub(temp1,temp1,sim3) mpz_sub(temp1,temp1,sim4) mpz_mul(temp2,C,E) mpz_mul(temp2,temp2,E) mpz_add(temp1,temp1,temp2) mpz_sub(temp1,temp1,sim5) mpz_mul_ui(temp2,E,2u) mpz_mul(temp2,temp2,F) mpz_mul(temp2,temp2,D) mpz_add(temp1,temp1,temp2) mpz_add(temp1,temp1,sim6) mpz_mul(temp2,G,G) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,H,H) mpz_sub(b81,temp1,temp2) #b82# mpz_mul(temp1,b21,J) mpz_mul(temp2,b22,I) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,b22,J) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,A,E) mpz_mul(temp2,temp2,H) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,A,F) mpz_mul(temp2,temp2,G) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,B,E) mpz_mul(temp2,temp2,G) mpz_sub(temp1,temp1,temp2) mpz_sub(temp1,temp1,sim1) mpz_sub(temp1,temp1,sim2) mpz_sub(temp1,temp1,sim4) mpz_sub(temp1,temp1,sim3) mpz_mul(temp2,C,E) mpz_mul(temp2,temp2,F) mpz_mul_ui(temp2,temp2,2u) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,C,F) mpz_mul(temp2,temp2,F) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,D,E) mpz_mul(temp2,temp2,E) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,D,E) mpz_mul(temp2,temp2,F) mpz_mul_ui(temp2,temp2,2u) mpz_add(temp1,temp1,temp2) mpz_add(temp1,temp1,sim6) mpz_add(temp1,temp1,sim6) mpz_mul(temp2,G,H) mpz_mul_ui(temp2,temp2,2u) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,H,H) mpz_mul_ui(temp2,temp2,2u) mpz_sub(b82,temp1,temp2) mpz_clear(A) mpz_clear(B) mpz_clear(C) mpz_clear(D) mpz_clear(E) mpz_clear(F) mpz_clear(G) mpz_clear(H) mpz_clear(I) mpz_clear(J) #disc similarities# mpz_mul(temp1,b21,b41) mpz_mul(temp2,b22,b42) mpz_add(sim1,temp1,temp2) mpz_add(temp1,b41,b42) mpz_mul(temp1,temp1,b22) mpz_mul(temp2,b21,b42) mpz_add(sim2,temp1,temp2) mpz_mul(temp1,b21,b21) mpz_mul(temp2,b22,b22) mpz_add(sim3,temp1,temp2) mpz_mul(temp1,b21,b22) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,b22,b22) mpz_add(sim4,temp1,temp2) #disc1# mpz_mul(temp1,sim1,b61) mpz_mul(temp2,sim2,b62) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp1,temp1,9u) #this finishes +9b2b4b6 mpz_mul(temp2,sim3,b81) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,sim4,b82) mpz_sub(temp1,temp1,temp2) #this finishes -b2^2b8 mpz_mul(temp2,b41,b41) mpz_mul(temp3,b42,b42) mpz_mul_ui(temp3,temp3,3u) mpz_add(temp2,temp2,temp3) mpz_mul(temp2,temp2,b41) mpz_mul(temp3,b42,b42) mpz_mul(temp3,temp3,b42) mpz_add(temp2,temp2,temp3) mpz_mul_ui(temp2,temp2,8u) mpz_sub(temp1,temp1,temp2) #this finishes -8b4^3 mpz_mul(temp2,b61,b61) mpz_mul(temp3,b62,b62) mpz_add(temp2,temp2,temp3) mpz_mul_ui(temp2,temp2,27u) mpz_sub(disc1,temp1,temp2) #this gives disc1 #disc2# mpz_mul(temp1,sim1,b62) mpz_add(temp2,b61,b62) mpz_mul(temp2,temp2,sim2) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp1,temp1,9u) #this finishes +9b2b4b6 mpz_mul(temp2,sim3,b82) mpz_sub(temp1,temp1,temp2) mpz_add(temp2,b81,b82) mpz_mul(temp2,temp2,sim4) mpz_sub(temp1,temp1,temp2) #this finishes -b2^2b8 mpz_add(temp2,b41,b42) mpz_mul(temp2,temp2,b41) mpz_mul_ui(temp2,temp2,3u) mpz_mul(temp3,b42,b42) mpz_mul_ui(temp3,temp3,2u) mpz_add(temp2,temp2,temp3) mpz_mul(temp2,temp2,b42) mpz_mul_ui(temp2,temp2,8u) mpz_sub(temp1,temp1,temp2) #this finishes -8b4^3 mpz_mul(temp2,b61,b62) mpz_mul_ui(temp2,temp2,2u) mpz_mul(temp3,b62,b62) mpz_add(temp2,temp2,temp3) mpz_mul_ui(temp2,temp2,27u) mpz_sub(disc2,temp1,temp2) #this gives disc2 mpz_clear(b21) mpz_clear(b22) mpz_clear(b41) mpz_clear(b42) mpz_clear(b61) mpz_clear(b62) mpz_clear(b81) mpz_clear(b82) mpz_clear(sim1) mpz_clear(sim2) mpz_clear(sim3) mpz_clear(sim4) mpz_clear(sim5) mpz_clear(sim6) #norm# mpz_mul(temp1,disc1,disc1) mpz_mul(temp2,disc1,disc2) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,disc2,disc2) mpz_sub(temp1,temp1,temp2) ret = mpz_get_ui(temp1) mpz_clear(temp1) mpz_clear(temp2) mpz_clear(temp3) return ret def qn(): cdef uint32_t A,B,C,D,E,F,G,H,I,J A=2 B=3 C=0 D=0 E=0 F=0 G=-1 H=1 I=3 J=8 return quick_disc_norm(A,B,C,D,E,F,G,H,I,J) 
qn() 
       
python(24680) malloc: *** error for object 0x10c211f80: pointer being
freed was not allocated
*** set a breakpoint in malloc_error_break to debug
python(24680) malloc: *** error for object 0x10c211f80: pointer being freed was not allocated
*** set a breakpoint in malloc_error_break to debug
       
Traceback (click to the left of this block for traceback)
...
NameError: name 'E' is not defined
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_2.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("RQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/private/var/folders/+3/+3+bRf7bGhaMk-2vvrWEvU+++TI/-Tmp-/tmpOwPB72/___code___.py", line 2, in <module>
    exec compile(u'E
  File "", line 1, in <module>
    
NameError: name 'E' is not defined
K.<a>=NumberField(x^2-x-1) 
       
 
       
R.<A,B,C,D,E,F,G,H,I,J>=K[] 
       
Ell=EllipticCurve(K,[2+3*a,0,0,-1+a,3+8*a]) 
       
Ell.discriminant().norm() 
       
2346173225
2346173225
%cython from libc.stdint cimport int32_t, uint32_t, uint64_t cdef mpz_t A,B,C,D,E,F,G,H,I,J, b21,b22,b41,b42,b61,b62,b81,b82, temp1,temp2, sim1,sim2,sim3,sim4, disc1,disc2 mpz_init(A) mpz_init(B) mpz_init(C) mpz_init(D) mpz_init(E) mpz_init(F) mpz_init(G) mpz_init(H) mpz_init(I) mpz_init(J) mpz_init(b21) mpz_init(b22) mpz_init(b41) mpz_init(b42) mpz_init(b61) mpz_init(b62) mpz_init(b81) mpz_init(b82) mpz_init(sim1) mpz_init(sim2) mpz_init(sim3) mpz_init(sim4) mpz_init(temp1) mpz_init(temp2) mpz_init(disc1) mpz_init(disc2) cdef uint64_t quick_disc_norm_new(uint32_t a_11, uint32_t a_12, uint32_t a_21, uint32_t a_22, uint32_t a_31, uint32_t a_32, uint32_t a_41, uint32_t a_42, uint32_t a_61, uint32_t a_62): mpz_set_ui(A,a_11) mpz_set_ui(B,a_12) mpz_set_ui(C,a_21) mpz_set_ui(D,a_22) mpz_set_ui(E,a_31) mpz_set_ui(F,a_32) mpz_set_ui(G,a_41) mpz_set_ui(H,a_42) mpz_set_ui(I,a_61) mpz_set_ui(J,a_62) #b21# mpz_mul(temp1,A,A) mpz_mul(temp2,B,B) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,C,4u) mpz_add(b21,temp1,temp2) #b22# mpz_mul(temp1,A,B) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,B,B) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,D,4) mpz_add(b22,temp1,temp2) #b41# mpz_mul(temp1,A,E) mpz_mul(temp2,B,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,G,2u) mpz_add(b41,temp1,temp2) #b42# mpz_mul(temp1,A,F) mpz_mul(temp2,B,E) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,B,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,H,2u) mpz_add(b42,temp1,temp2) #b61# mpz_mul(temp1,E,E) mpz_mul(temp2,F,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,I,4u) mpz_add(b61,temp1,temp2) #b62# mpz_mul(temp1,E,F) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,F,F) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp2,J,4) mpz_add(b62,temp1,temp2) #b8 similarities# mpz_mul(temp1,A,E) mpz_mul(temp2,B,F) mpz_add(sim1,temp1,temp2) # mpz_mul(temp1,A,F) mpz_mul(temp2,B,E) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,B,F) mpz_add(sim2,temp1,temp2) # mpz_mul(temp1,E,E) mpz_mul(temp2,F,F) mpz_add(sim3,temp1,temp2) # mpz_mul(temp1,E,F) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,F,F) mpz_add(sim4,temp1,temp2) # #b81# mpz_mul(temp1,b21,I) mpz_mul(temp2,b22,J) mpz_add(temp1,temp1,temp2) #finished with b2a6 mpz_mul(temp2,sim1,G) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,sim2,H) mpz_sub(temp1,temp1,temp2) #finished with -a1a3a4 mpz_mul(temp2,sim3,C) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,sim4,D) mpz_add(temp1,temp1,temp2) #finished with a2a3^2 mpz_mul(temp2,G,G) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,H,H) mpz_sub(b81,temp1,temp2) #b81 calculated #b82# mpz_mul(temp1,b21,J) mpz_mul(temp2,b22,I) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,b22,J) mpz_add(temp1,temp1,temp2) #finished with b2a6 mpz_mul(temp2,sim1,H) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,sim2,G) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,sim2,H) mpz_sub(temp1,temp1,temp2) #finished with -a1a3a4 mpz_mul(temp2,sim3,D) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,sim4,C) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,sim4,D) mpz_add(temp1,temp1,temp2) #finished with a2a3^2 mpz_mul(temp2,G,H) mpz_mul_ui(temp2,temp2,2u) mpz_sub(temp1,temp1,temp2) mpz_mul(temp2,H,H) mpz_sub(b82,temp1,temp2) #b82 calculated #disc similarities# mpz_mul(temp1,b21,b41) mpz_mul(temp2,b22,b42) mpz_add(sim1,temp1,temp2) mpz_add(temp1,b41,b42) mpz_mul(temp1,temp1,b22) mpz_mul(temp2,b21,b42) mpz_add(sim2,temp1,temp2) mpz_mul(temp1,b21,b21) mpz_mul(temp2,b22,b22) mpz_add(sim3,temp1,temp2) mpz_mul(temp1,b21,b22) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,b22,b22) mpz_add(sim4,temp1,temp2) #disc1# mpz_mul(disc1,sim1,b61) mpz_mul(temp1,sim2,b62) mpz_add(disc1,disc1,temp1) mpz_mul_ui(disc1,disc1,9u) #this finishes +9b2b4b6 mpz_mul(temp1,sim3,b81) mpz_sub(disc1,disc1,temp1) mpz_mul(temp1,sim4,b82) mpz_sub(disc1,disc1,temp1) #this finishes -b2^2b8 mpz_mul(temp1,b41,b41) mpz_mul(temp2,b42,b42) mpz_mul_ui(temp2,temp2,3u) mpz_add(temp1,temp1,temp2) mpz_mul(temp1,temp1,b41) mpz_mul(temp2,b42,b42) mpz_mul(temp2,temp2,b42) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp1,temp1,8u) mpz_sub(disc1,disc1,temp1) #this finishes -8b4^3 mpz_mul(temp1,b61,b61) mpz_mul(temp2,b62,b62) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp1,temp1,27u) mpz_sub(disc1,disc1,temp1) #this gives disc1 #disc2# mpz_mul(disc2,sim1,b62) mpz_add(temp1,b61,b62) mpz_mul(temp1,temp1,sim2) mpz_add(disc2,disc2,temp1) mpz_mul_ui(disc2,disc2,9u) #this finishes +9b2b4b6 mpz_mul(temp1,sim3,b82) mpz_sub(disc2,disc2,temp1) mpz_add(temp1,b81,b82) mpz_mul(temp1,temp1,sim4) mpz_sub(disc2,disc2,temp1) #this finishes -b2^2b8 mpz_add(temp1,b41,b42) mpz_mul(temp1,temp1,b41) mpz_mul_ui(temp1,temp1,3u) mpz_mul(temp2,b42,b42) mpz_mul_ui(temp2,temp2,2u) mpz_add(temp1,temp1,temp2) mpz_mul(temp1,temp1,b42) mpz_mul_ui(temp1,temp1,8u) mpz_sub(disc2,disc2,temp1) #this finishes -8b4^3 mpz_mul(temp1,b61,b62) mpz_mul_ui(temp1,temp1,2u) mpz_mul(temp2,b62,b62) mpz_add(temp1,temp1,temp2) mpz_mul_ui(temp1,temp1,27u) mpz_sub(disc2,disc2,temp1) #this gives disc2 #norm# mpz_mul(temp1,disc1,disc1) mpz_mul(temp2,disc1,disc2) mpz_add(temp1,temp1,temp2) mpz_mul(temp2,disc2,disc2) mpz_sub(temp1,temp1,temp2) ret = mpz_get_si(temp1) return ret def qn_new(): return quick_disc_norm_new(1,2,3,4,9,6,7,8,9,10) 
print norm(10+53*a), 10^2+10*53-53^2 
       
Traceback (click to the left of this block for traceback)
...
NameError: name 'a' is not defined
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_6.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("cHJpbnQgbm9ybSgxMCs1MyphKSwgMTBeMisxMCo1My01M14y"),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/private/var/folders/+3/+3+bRf7bGhaMk-2vvrWEvU+++TI/-Tmp-/tmpk75WuL/___code___.py", line 3, in <module>
    exec compile(u'print norm(_sage_const_10 +_sage_const_53 *a), _sage_const_10 **_sage_const_2 +_sage_const_10 *_sage_const_53 -_sage_const_53 **_sage_const_2 
  File "", line 1, in <module>
    
NameError: name 'a' is not defined
8*10*(2*1000+1)^3 
       
640960480080
640960480080
2^64 
       
18446744073709551616
18446744073709551616
199^3 
       
7880599
7880599
qn_new() 
       
23068249636L
23068249636L
Ell=EllipticCurve(K,[1+2*a,3+4*a,9+6*a,7+8*a,9+10*a]) print Ell.b2() print Ell.b4() print Ell.b6() print Ell.b8() print Ell.discriminant()[0] print Ell.discriminant()[1] print Ell.discriminant() Ell.discriminant().norm() 
       
24*a + 17
52*a + 35
184*a + 153
1218*a + 772
-189172
-241774
-241774*a - 189172
23068249636
24*a + 17
52*a + 35
184*a + 153
1218*a + 772
-189172
-241774
-241774*a - 189172
23068249636
d1=-189172 d2=-241774 
       
d1^2+d1*d2-d2^2 
       
23068249636
23068249636
2^32 
       
4294967296
4294967296
timeit('qn_new()') 
       
625 loops, best of 3: 20.6 µs per loop
625 loops, best of 3: 20.6 µs per loop
def quick_disc(ainv): a1 = ainv[0] a2 = ainv[1] a3 = ainv[2] a4 = ainv[3] a6 = ainv[4] return -(a1^2+4*a2)^2*((a1^2+4*a2)*a6-a1*a3*a4+a2*a3^2-a4^2)-8*(a1*a3+2*a4)^3-27*(a3^2+4*a6)^2+9*(a1^2+4*a2)*(a1*a3+2*a4)*(a3^2+4*a6) 
       
timeit('quick_disc([1+2*a,3+4*a,9+6*a,7+8*a,9+10*a])') 
       
625 loops, best of 3: 281 µs per loop
625 loops, best of 3: 281 µs per loop
def f(n): Ell=EllipticCurve(K,[1+2*a,3+4*a,9+6*a,7+8*a,9+10*a]) return Ell.discriminant().norm() 
       
timeit('f(1)') 
       
125 loops, best of 3: 2.61 ms per loop
125 loops, best of 3: 2.61 ms per loop
Ell.b8() 
       
810*a + 504
810*a + 504
88**3*64+88**2*432 
       
46959616
46959616
len(46959616.bits()) 
       
26
26
%cython def boring_func(): def hi(int x): return x*x cdef int i for i in range(100): print(hi(i)) 
boring_func() 
       
0
1
4
9
16
25
36
49
64
81
100
121
144
169
196
225
256
289
324
361
400
441
484
529
576
625
676
729
784
841
900
961
1024
1089
1156
1225
1296
1369
1444
1521
1600
1681
1764
1849
1936
2025
2116
2209
2304
2401
2500
2601
2704
2809
2916
3025
3136
3249
3364
3481
3600
3721
3844
3969
4096
4225
4356
4489
4624
4761
4900
5041
5184
5329
5476
5625
5776
5929
6084
6241
6400
6561
6724
6889
7056
7225
7396
7569
7744
7921
8100
8281
8464
8649
8836
9025
9216
9409
9604
9801
0
1
4
9
16
25
36
49
64
81
100
121
144
169
196
225
256
289
324
361
400
441
484
529
576
625
676
729
784
841
900
961
1024
1089
1156
1225
1296
1369
1444
1521
1600
1681
1764
1849
1936
2025
2116
2209
2304
2401
2500
2601
2704
2809
2916
3025
3136
3249
3364
3481
3600
3721
3844
3969
4096
4225
4356
4489
4624
4761
4900
5041
5184
5329
5476
5625
5776
5929
6084
6241
6400
6561
6724
6889
7056
7225
7396
7569
7744
7921
8100
8281
8464
8649
8836
9025
9216
9409
9604
9801