A033195 Erroneous version of A133675.
3, 4, 7, 8, 11, 16, 19, 27, 28, 43, 67, 163
Offset: 1
Keywords
References
- David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989, p. 260.
Crossrefs
Cf. A133675.
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
ok(n)={(-n)%4<2 && quadclassunit(-n).no > 1 && !#select(k->k<>2, quadclassunit(-n).cyc)} \\ Andrew Howroyd, Jul 20 2018
for(n=1, 500, if((-n)%4<=1&&quadclassunit(-n)[1]==2, print1(n, ", ")))
isA328825(d) = (d>0) && ((d%4==0)||(d%4==3)) && quadclassunit(-d)[2]==[3] \\ Corrected by Jianing Song, May 17 2021
a(15)=2 because the forms x^2 + xy + 4y^2 and 2x^2 + xy + 2y^2 have discriminant -15.
dLim=150; cnt=Table[0, {dLim}]; nn=Ceiling[dLim/4]; Do[d=b^2-4a*c; If[GCD[a, b, c]==1 && 0<-d<=dLim, cnt[[ -d]]++ ], {b, 0, nn}, {a, b, nn}, {c, a, nn}]; cnt
{a(n)=local(m); if(n<3, 0, forvec(v=vector(3,k,[0,(n+1)\4]), if( (gcd(v)==1)&(-v[1]^2+4*v[2]*v[3]==n), m++ ), 1); m)} /* Michael Somos, May 31 2005 */
isA329182(d) = (d>0) && ((d%4==0)||(d%4==3)) && quadclassunit(-d)[2]==[2,2] \\ Jianing Song, May 17 2021
isA330219(d) = (d>0) && ((d%4==0)||(d%4==3)) && quadclassunit(-d)[2]==[4] \\ Jianing Song, May 17 2021
n = 27 is not a member of this sequence because N = 27 = (0/2)^2 + 27*(2/2)^2 is of the form (x/2)^2+27*(y/2)^2, but for the prime factor 3 which appears to an odd power in N = 27, 3 is not of the form (x/2)^2+27*(y/2)^2 because 4 * 3 = 12 is not of the form x^2+27*y^2 (also 3 itself is not of the form x^2+27*y^2).
Comments