cp's OEIS Frontend

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.

Showing 1-10 of 38 results. Next

A055026 Number of Gaussian primes of successive norms (indexed by A055025).

Original entry on oeis.org

4, 8, 4, 8, 8, 8, 8, 8, 4, 8, 8, 8, 8, 8, 8, 8, 8, 4, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 4, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 4, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 4, 8, 8
Offset: 1

Views

Author

N. J. A. Sloane, Jun 09 2000

Keywords

Comments

These are the primes in the ring of integers a+bi, a and b rational integers, i = sqrt(-1).

Examples

			There are 8 Gaussian primes of norm 5, +-1+-2i and +-2+-i, but only two inequivalent ones (2+-i).
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, A16.
  • L. W. Reid, The Elements of the Theory of Algebraic Numbers, MacMillan, NY, 1910, see Chap. V.

Crossrefs

Programs

  • Mathematica
    m = 32; Length /@ Split[Sort[Select[Flatten[Table[{a^2 + b^2, a + b*I}, {a, -m, m}, {b, -m, m}], 1], PrimeQ[#[[2]], GaussianIntegers -> True] & ]], #1[[1]] == #2[[1]] & ][[1 ;; 87]] (* Jean-François Alcover, Apr 08 2011 *)

Extensions

More terms from Reiner Martin, Jul 20 2001

A055027 Number of inequivalent Gaussian primes of successive norms (indexed by A055025).

Original entry on oeis.org

1, 2, 1, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2
Offset: 1

Views

Author

N. J. A. Sloane, Jun 09 2000

Keywords

Comments

These are the primes in the ring of integers a+bi, a and b rational integers, i = sqrt(-1).
Two primes are considered equivalent if they differ by multiplication by a unit (+-1, +-i).

Examples

			There are 8 Gaussian primes of norm 5, +-1+-2i and +-2+-i, but only two inequivalent ones (2+-i).
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, A16.
  • L. W. Reid, The Elements of the Theory of Algebraic Numbers, MacMillan, NY, 1910, see Chap. V.

Crossrefs

Programs

  • Mathematica
    norms = Union[ #*Conjugate[#]& [ Select[ Flatten[ Table[a + b*I, {a, 0, 31}, {b, 0, 31}]], PrimeQ[#, GaussianIntegers -> True] &]]]; f[norm_] := (Clear[a, b]; primes = {a + b*I} /. {ToRules[ Reduce[a^2 + b^2 == norm, {a, b}, Integers]]}; primes //. {p1___, p2_, p3___, p4_, p5___} /; MatchQ[p2, (-p4 | I*p4 | -I*p4)] :> {p1, p2, p3, p5} // Length); A055027 = f /@ norms (* Jean-François Alcover, Nov 30 2012 *)

Extensions

More terms from Reiner Martin, Jul 20 2001

A300588 Imaginary part y of the n-th Gaussian prime x + i*y, x >= y >= 0, ordered by norm x^2 + y^2 = A055025(n)^2.

Original entry on oeis.org

1, 1, 0, 2, 1, 2, 1, 4, 0, 2, 5, 3, 5, 4, 1, 3, 7, 0, 4, 7, 6, 2, 9, 7, 1, 2, 8, 4, 1, 10, 9, 5, 2, 12, 11, 9, 5, 8, 0, 7, 10, 6, 1, 3, 14, 12, 7, 4, 10, 5, 11, 0, 10, 14, 13, 1, 8, 5, 17, 16, 4, 13, 6, 12, 1, 5, 15, 2, 9, 19, 12, 17, 11, 5, 14, 10, 18, 4, 6, 16, 20, 19, 10, 13, 0, 4, 6
Offset: 1

Views

Author

M. F. Hasler, Mar 09 2018

Keywords

Comments

With the restriction Re(z) >= Im(z) >= 0 used here and in A239621, there is exactly one Gaussian prime z for each possible norm |z|^2 in A055025. Sequence A239397 lists both, (x, y) and (y, x), for each of these having x > y (i.e., except for x = y = 1).
According to the graph, the values seem rather uniformly distributed between 0 and the upper bound sqrt(A055025(n)/2) ~ sqrt(n log n), in contrast to the values of the real parts A300587(n).

Examples

			The smallest Gaussian primes with Re(z) >= Im(z) >= 0, ordered by norm, are 1+i, 2+i, 3, 3+i, ...
Their imaginary parts, listed here, are a(1) = 1, a(2) = 1, a(3) = 0, a(4) = 1,
		

Crossrefs

Even bisection of A239621. See A300587 for real parts, A055025 for the norms.

Programs

  • PARI
    c=1; for(n=1,oo, matsize(f=factor(n*I))[1]<=2 && vecsum(f[,2])==2+(f[1, 1]==I) && !write("/tmp/b300588.txt",c" "min(real(f=f[3-f[1,2],1]),imag(f))) && c++>1e4 && break) \\ Replace write("/tmp/b300588.txt",c" by print1(" to print the values.

Formula

a(n) = A239621(2n) = A239397(4n-3) (= A239397(4n-4) for n > 1).
a(n) = sqrt(A055025(n) - A300587(n)^2).

A218217 a(n) = (x(n+1) - x(n))^2 + (y(n+1) - y(n))^2, where x(n)^2 + y(n)^2 = A055025(n) is norm of Gaussian prime and x(n) >= y(n) >= 0.

Original entry on oeis.org

1, 2, 4, 2, 2, 2, 10, 20, 4, 10, 8, 4, 2, 10, 4, 20, 58, 16, 10, 2, 20, 58, 8, 40, 2, 40, 20, 10, 90, 2, 20, 10, 116, 2, 8, 20, 10, 68, 50, 10, 20, 26, 4, 146, 8, 34, 10, 40, 34, 40, 130, 104, 20, 2, 160, 50, 10, 180, 2, 180, 90, 58, 40, 130, 16, 116, 194, 50
Offset: 1

Views

Author

Thomas Ordowski, Oct 23 2012

Keywords

Comments

We walk through the Gaussian primes in the first octant of the Gaussian plane along increasing norm: 1+i, 2+i, 3, 3+2i, 4+i, 5+2i, 6+i, 5+4i, 7, 7+2i etc. The sequence lists the squared distance between consecutive Gaussian primes along this walk.

Examples

			The squared distance between 5+4i and 6+i is (6-5)^2+(4-1)^2 =10 = a(7).
		

Crossrefs

Cf. A055025.

Programs

  • Mathematica
    nMx = 1000; modLst = {2}; Do[p = Prime[n]; If[Mod[p, 4] == 1, AppendTo[modLst, p], If[p^2 < nMx, AppendTo[modLst, p^2]]], {n, 2, PrimePi[nMx]}]; modLst = Union[modLst]; last = {1, 1}; Table[pr = PowersRepresentations[n, 2, 2][[1]]; dist = SquaredEuclideanDistance[last, pr]; last = pr; dist, {n, Rest[modLst]}] (* T. D. Noe, Oct 29 2012 *)

A001481 Numbers that are the sum of 2 squares.

Original entry on oeis.org

0, 1, 2, 4, 5, 8, 9, 10, 13, 16, 17, 18, 20, 25, 26, 29, 32, 34, 36, 37, 40, 41, 45, 49, 50, 52, 53, 58, 61, 64, 65, 68, 72, 73, 74, 80, 81, 82, 85, 89, 90, 97, 98, 100, 101, 104, 106, 109, 113, 116, 117, 121, 122, 125, 128, 130, 136, 137, 144, 145, 146, 148, 149, 153, 157, 160
Offset: 1

Views

Author

Keywords

Comments

Numbers n such that n = x^2 + y^2 has a solution in nonnegative integers x, y.
Closed under multiplication. - David W. Wilson, Dec 20 2004
Also, numbers whose cubes are the sum of 2 squares. - Artur Jasinski, Nov 21 2006 (Cf. A125110.)
Terms are the squares of smallest radii of circles covering (on a square grid) a number of points equal to the terms of A057961. - Philippe Lallouet (philip.lallouet(AT)wanadoo.fr), Apr 16 2007. [Comment corrected by T. D. Noe, Mar 28 2008]
Numbers with more 4k+1 divisors than 4k+3 divisors. If a(n) is a member of this sequence, then so too is any power of a(n). - Ant King, Oct 05 2010
A000161(a(n)) > 0; A070176(a(n)) = 0. - Reinhard Zumkeller, Feb 04 2012, Aug 16 2011
Numbers that are the norms of Gaussian integers. This sequence has unique factorization; the primitive elements are A055025. - Franklin T. Adams-Watters, Nov 25 2011
These are numbers n such that all of n's odd prime factors congruent to 3 modulo 4 occur to an even exponent (Fermat's two-squares theorem). - Jean-Christophe Hervé, May 01 2013
Let's say that an integer n divides a lattice if there exists a sublattice of index n. Example: 2, 4, 5 divide the square lattice. The present sequence without 0 is the sequence of divisors of the square lattice. Say that n is a "prime divisor" if the index-n sublattice is not contained in any other sublattice except the original lattice itself. Then A055025 (norms of Gaussian primes) gives the "prime divisors" of the square lattice. - Jean-Christophe Hervé, May 01 2013
For any i,j > 0 a(i)*a(j) is a member of this sequence, since (a^2 + b^2)*(c^2 + d^2) = (a*c + b*d)^2 + (a*d - b*c)^2. - Boris Putievskiy, May 05 2013
The sequence is closed under multiplication. Primitive elements are in A055025. The sequence can be split into 3 multiplicatively closed subsequences: {0}, A004431 and A125853. - Jean-Christophe Hervé, Nov 17 2013
Generalizing Jasinski's comment, same as numbers whose odd powers are the sum of 2 squares, by Fermat's two-squares theorem. - Jonathan Sondow, Jan 24 2014
By the 4 squares theorem, every nonnegative integer can be expressed as the sum of two elements of this sequence. - Franklin T. Adams-Watters, Mar 28 2015
There are never more than 3 consecutive terms. Runs of 3 terms start at 0, 8, 16, 72, ... (A082982). - Ivan Neretin, Nov 09 2015
Conjecture: barring the 0+2, 0+4, 0+8, 0+16, ... sequence, the sum of 2 distinct terms in this sequence is never a power of 2. - J. Lowell, Jan 14 2022
All the areas of squares whose vertices have integer coordinates. - Neeme Vaino, Jun 14 2023
Numbers represented by the definite binary quadratic forms x^2 + 2nxy + (n^2+1)y^2 for any integer n. This sequence contains the even powers of any integer. An odd power of a number appears only if the number itself belongs to the sequence. The equation given in the comment by Boris Putievskiy 2013 is Brahmagupta's identity with n = 1. It proves that any set of numbers of the form a^2 + nb^2 is closed under multiplication. - Klaus Purath, Sep 06 2023

References

  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 106.
  • David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989.
  • L. Euler, (E388) Vollständige Anleitung zur Algebra, Zweiter Theil, reprinted in: Opera Omnia. Teubner, Leipzig, 1911, Series (1), Vol. 1, p. 417.
  • S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 98-104.
  • G. H. Hardy, Ramanujan, pp. 60-63.
  • P. Moree and J. Cazaran, On a claim of Ramanujan in his first letter to Hardy, Expos. Math. 17 (1999), pp. 289-312.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Disjoint union of A000290 and A000415.
Complement of A022544.
A000404 gives another version. Subsequence of A091072, supersequence of A046711.
Column k=2 of A336820.

Programs

  • Haskell
    a001481 n = a001481_list !! (n-1)
    a001481_list = [x | x <- [0..], a000161 x > 0]
    -- Reinhard Zumkeller, Feb 14 2012, Aug 16 2011
    
  • Magma
    [n: n in [0..160] | NormEquation(1, n) eq true]; // Arkadiusz Wesolowski, May 11 2016
    
  • Maple
    readlib(issqr): for n from 0 to 160 do for k from 0 to floor(sqrt(n)) do if issqr(n-k^2) then printf(`%d,`,n); break fi: od: od:
  • Mathematica
    upTo = 160; With[{max = Ceiling[Sqrt[upTo]]}, Select[Union[Total /@ (Tuples[Range[0, max], {2}]^2)], # <= upTo &]]  (* Harvey P. Dale, Apr 22 2011 *)
    Select[Range[0, 160], SquaresR[2, #] != 0 &] (* Jean-François Alcover, Jan 04 2013 *)
  • PARI
    isA001481(n)=local(x,r);x=0;r=0;while(x<=sqrt(n) && r==0,if(issquare(n-x^2),r=1);x++);r \\ Michael B. Porter, Oct 31 2009
    
  • PARI
    is(n)=my(f=factor(n));for(i=1,#f[,1],if(f[i,2]%2 && f[i,1]%4==3, return(0))); 1 \\ Charles R Greathouse IV, Aug 24 2012
    
  • PARI
    B=bnfinit('z^2+1,1);
    is(n)=#bnfisintnorm(B,n) \\ Ralf Stephan, Oct 18 2013, edited by M. F. Hasler, Nov 21 2017
    
  • PARI
    list(lim)=my(v=List(),t); for(m=0,sqrtint(lim\=1), t=m^2; for(n=0, min(sqrtint(lim-t),m), listput(v,t+n^2))); Set(v) \\ Charles R Greathouse IV, Jan 05 2016
    
  • PARI
    is_A001481(n)=!for(i=2-bittest(n,0),#n=factor(n)~, bittest(n[1,i],1)&&bittest(n[2,i],0)&&return) \\ M. F. Hasler, Nov 20 2017
    
  • Python
    from itertools import count, islice
    from sympy import factorint
    def A001481_gen(): # generator of terms
        return filter(lambda n:(lambda m:all(d & 3 != 3 or m[d] & 1 == 0 for d in m))(factorint(n)),count(0))
    A001481_list = list(islice(A001481_gen(),30)) # Chai Wah Wu, Jun 27 2022

Formula

n = square * 2^{0 or 1} * {product of distinct primes == 1 (mod 4)}.
The number of integers less than N that are sums of two squares is asymptotic to constant*N/sqrt(log(N)), hence lim_{n->infinity} a(n)/n = infinity.
Nonzero terms in expansion of Dirichlet series Product_p (1 - (Kronecker(m, p) + 1)*p^(-s) + Kronecker(m, p)*p^(-2s))^(-1) for m = -1.
a(n) ~ k*n*sqrt(log n), where k = 1.3085... = 1/A064533. - Charles R Greathouse IV, Apr 16 2012
There are B(x) = x/sqrt(log x) * (K + B2/log x + O(1/log^2 x)) terms of this sequence up to x, where K = A064533 and B2 = A227158. - Charles R Greathouse IV, Nov 18 2022

Extensions

Deleted an incorrect comment. - N. J. A. Sloane, Oct 03 2023

A055664 Norms of Eisenstein-Jacobi primes.

Original entry on oeis.org

3, 4, 7, 13, 19, 25, 31, 37, 43, 61, 67, 73, 79, 97, 103, 109, 121, 127, 139, 151, 157, 163, 181, 193, 199, 211, 223, 229, 241, 271, 277, 283, 289, 307, 313, 331, 337, 349, 367, 373, 379, 397, 409, 421, 433, 439, 457, 463, 487, 499, 523, 529, 541, 547, 571
Offset: 1

Views

Author

N. J. A. Sloane, Jun 09 2000

Keywords

Comments

These are the norms of the primes in the ring of integers a+b*omega, a and b rational integers, omega = (1+sqrt(-3))/2.
Let us say that an integer n divides a lattice if there exists a sublattice of index n. Example: 3 divides the hexagonal lattice. Then A003136 (Loeschian numbers) is the sequence of divisors of the hexagonal lattice. Say that n is a "prime divisor" if the index-n sublattice is not contained in any other sublattice except the original lattice itself. The present sequence gives the prime divisors of the hexagonal lattice. Similarly, A055025 (Norms of Gaussian primes) is the sequence of "prime divisors" of the square lattice. - Jean-Christophe Hervé, Dec 04 2006

Examples

			There are 6 Eisenstein-Jacobi primes of norm 3, omega-omega^2 times one of the 6 units [ +-1, +-omega, +-omega^2 ] but only one up to equivalence.
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, A16.
  • L. W. Reid, The Elements of the Theory of Algebraic Numbers, MacMillan, NY, 1910, see Chap. VI.

Crossrefs

Cf. A055665-A055668, A055025-A055029, A135461, A135462. See A004016 and A035019 for theta series of Eisenstein (or hexagonal) lattice.
The Z[sqrt(-5)] analogs are in A020669, A091727, A091728, A091729, A091730 and A091731.

Programs

  • Mathematica
    Join[{3}, Select[Range[600], (PrimeQ[#] && Mod[#, 6] == 1) || (PrimeQ[Sqrt[#]] && Mod[Sqrt[#], 3] == 2) & ]] (* Jean-François Alcover, Oct 09 2012, from formula *)
  • PARI
    is(n)=(isprime(n) && n%3<2) || (issquare(n,&n) && isprime(n) && n%3==2) \\ Charles R Greathouse IV, Apr 30 2013

Formula

Consists of 3; rational primes == 1 (mod 3) [A002476]; and squares of rational primes == -1 (mod 3) [A003627^2].

Extensions

More terms from David Wasserman, Mar 21 2002

A103431 Subsequence of the Gaussian primes, where only Gaussian primes a+bi with a>0, b>=0 are listed. Ordered by the norm N(a+bi)=a^2+b^2 and the size of the real part, when the norms are equal. a(n) is the real part of the Gaussian prime. Sequence A103432 gives the imaginary parts.

Original entry on oeis.org

1, 1, 2, 3, 2, 3, 1, 4, 2, 5, 1, 6, 4, 5, 7, 2, 7, 5, 6, 3, 8, 5, 8, 4, 9, 1, 10, 3, 10, 7, 8, 11, 4, 11, 7, 10, 6, 11, 2, 13, 9, 10, 7, 12, 1, 14, 2, 15, 8, 13, 4, 15, 1, 16, 10, 13, 9, 14, 5, 16, 2, 17, 12, 13, 11, 14, 9, 16, 5, 18, 8, 17, 19, 7, 18, 10, 17, 6, 19, 1, 20, 3, 20, 14, 15, 12, 17
Offset: 1

Views

Author

Sven Simon, Feb 05 2005; corrected Feb 20 2005 and again on Aug 06 2006

Keywords

Comments

Definition of Gaussian primes (Pieper, Die komplexen Zahlen, p. 122): 1) i+i, norm N(i+i) = 2. 2) Natural primes p with p = 3 mod 4, norm N(p) = p^2. 3) primes a+bi, a>0, b>0 with a^2 + b^2 = p = 1 mod 4, p natural prime. Norm N(a+bi) = p. b+ai is a different Gaussian prime number, b+ai cannot be factored into a+bi and a unit. 4) All complex numbers from 1) to 3) multiplied by the units -1,i,-i, these are the associated numbers. The sequence contains all the Gaussian primes mentioned in 1) - 3).
Every complex number can be factored completely into the Gaussian prime numbers defined by the sequence, an additional unit as factor can be necessary. This factorization can be used to calculate the complex sigma, as defined by Spira. The elements a(n) are ordered by the size of their norm. If the two different primes a+bi and b+ai have the same norm, they are ordered by the size of the real part of the complex prime number. So a+bi follows b+ai in the sequence, if a > b.
Of course this is not the only possible definition. As primes p = 1 mod 4 can be factored in p = (-i)(a+bi)(b+ai) and the norm N(a+bi) = N(b+ai) = p, these primes a+bi occur much earlier in the sequence than p does in the sequence of natural primes. 4+5i with norm 41 occurs before prime 7 with norm 49.

References

  • H. Pieper, "Die komplexen Zahlen", Verlag Harri Deutsch, p. 122

Crossrefs

Programs

  • Maple
    N:= 100: # to get all terms with norm <= N
    p1:= select(isprime,[seq(i,i=3..N,4)]):
    p2:= select(isprime,[seq(i,i=1..N^2,4)]):
    p2:= map(t -> GaussInt:-GIfactors(t)[2][1][1],p2):
    p3:= sort( [1+I, op(p1),op(p2)],(a,b) -> Re(a)^2 + Im(a)^2  < Re(b)^2 + Im(b)^2):
    g:= proc(z)
        local a,b;
        a:= Re(z); b:= Im(z);
        if b = 0 then z
        else
          a:= abs(a);
          b:= abs(b);
          if a = b then a
          elif a < b then a,b
          else b,a
          fi
        fi
    end proc:
    map(g, p3); # Robert Israel, Feb 23 2016
  • Mathematica
    maxNorm = 500;
    norm[z_] := Re[z]^2 + Im[z]^2;
    m = Sqrt[maxNorm] // Ceiling;
    gp = Select[Table[a + b I, {a, 1, m}, {b, 0, m}] // Flatten, norm[#] <= maxNorm && PrimeQ[#, GaussianIntegers -> True]&];
    SortBy[gp, norm[#] maxNorm + Abs[Re[#]]&] // Re (* Jean-François Alcover, Mar 04 2019 *)

Extensions

Edited (mostly to correct meaning of norm) by Franklin T. Adams-Watters, Mar 04 2011
a(48) corrected by Robert Israel, Feb 23 2016

A055029 Number of inequivalent Gaussian primes of norm n.

Original entry on oeis.org

0, 0, 1, 0, 0, 2, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0
Offset: 0

Views

Author

N. J. A. Sloane, Jun 09 2000

Keywords

Comments

These are the primes in the ring of integers a+bi, a and b rational integers, i = sqrt(-1).
Two primes are considered equivalent if they differ by multiplication by a unit (+-1, +-i).

Examples

			There are 8 Gaussian primes of norm 5, +-1+-2i and +-2+-i, but only two inequivalent ones (2+-i).
		

References

  • R. K. Guy, Unsolved Problems in Number Theory, A16.
  • L. W. Reid, The Elements of the Theory of Algebraic Numbers, MacMillan, NY, 1910, see Chap. V.

Crossrefs

Programs

Formula

a(n) = A055028(n)/4.
a(n) = 2 if n is a prime = 1 (mod 4); a(n) = 1 if n is 2, or p^2 where p is a prime = 3 (mod 4); a(n) = 0 otherwise. - Franklin T. Adams-Watters, May 05 2006
a(n) = if n = 2 then 1 else 2*A079260(n) + A079261(A037213(n)). - Reinhard Zumkeller, Nov 11 2012

Extensions

More terms from Reiner Martin, Jul 20 2001

A091727 Norms of prime ideals of Z[sqrt(-5)].

Original entry on oeis.org

2, 3, 5, 7, 23, 29, 41, 43, 47, 61, 67, 83, 89, 101, 103, 107, 109, 121, 127, 149, 163, 167, 169, 181, 223, 227, 229, 241, 263, 269, 281, 283, 289, 307, 347, 349, 361, 367, 383, 389, 401, 409, 421, 443, 449, 461, 463, 467, 487
Offset: 1

Views

Author

Paul Boddington, Feb 02 2004

Keywords

Comments

Consists of primes congruent to 1, 2, 3, 5, 7, 9 (mod 20) together with the squares of all other primes.
From Jianing Song, Feb 20 2021: (Start)
The norm of a nonzero ideal I in a ring R is defined as the size of the quotient ring R/I.
Note that Z[sqrt(-5)] has class number 2.
For primes p == 1, 9 (mod 20), there are two distinct ideals with norm p in Z[sqrt(-5)], namely (x + y*sqrt(-5)) and (x - y*sqrt(-5)), where (x,y) is a solution to x^2 + 5*y^2 = p.
For p == 3, 7 (mod 20), there are also two distinct ideals with norm p, namely (p, x+y*sqrt(-5)) and (p, x-y*sqrt(-5)), where (x,y) is a solution to x^2 + 5*y^2 = p^2 with y != 0; (2, 1+sqrt(-5)) and (sqrt(-5)) are respectively the unique ideal with norm 2 and 5.
For p == 11, 13, 17, 19 (mod 20), (p) is the only ideal with norm p^2. (End)

Examples

			From _Jianing Song_, Feb 20 2021: (Start)
Let |I| be the norm of an ideal I, then:
|(2, 1+sqrt(-5))| = 2;
|(3, 2+sqrt(-5))| = |(3, 2-sqrt(-5))| = 3;
|(sqrt(-5))| = 5;
|(7, 1+3*sqrt(-5))| = |(7, 1-3*sqrt(-5))| = 7;
|(23, 22+3*sqrt(-5))| = |(23, 22-3*sqrt(-5))| = 23;
|(3 + 2*sqrt(-5))| = |(3 - 2*sqrt(-5))| = 29;
|(6 + sqrt(-5))| = |(6 - sqrt(-5))| = 41. (End)
		

References

  • David A. Cox, Primes of the form x^2+ny^2, Wiley, 1989.
  • A. Frohlich and M. J. Taylor, Algebraic number theory, Cambridge university press, 1991.

Crossrefs

Cf. A091728.
The number of distinct ideals with norm n is given by A035170.
Norms of prime ideals in O_K, where K is the quadratic field with discriminant D and O_K be the ring of integers of K: A055673 (D=8), A341783 (D=5), A055664 (D=-3), A055025 (D=-4), A090348 (D=-7), A341784 (D=-8), A341785 (D=-11), A341786 (D=-15*), A341787 (D=-19), this sequence (D=-20*), A341788 (D=-43), A341789 (D=-67), A341790 (D=-163). Here a "*" indicates the cases where O_K is not a unique factorization domain.

Programs

  • PARI
    isA091727(n) = { my(ms = [1, 2, 3, 5, 7, 9], p, e=isprimepower(n,&p)); if(!e || e>2, 0, bitxor(e-1,!!vecsearch(ms,p%20))); }; \\ Antti Karttunen, Feb 24 2020

Extensions

Offset corrected by Jianing Song, Feb 20 2021

A341784 Norms of prime elements in Z[sqrt(-2)], the ring of integers of Q(sqrt(-2)).

Original entry on oeis.org

2, 3, 11, 17, 19, 25, 41, 43, 49, 59, 67, 73, 83, 89, 97, 107, 113, 131, 137, 139, 163, 169, 179, 193, 211, 227, 233, 241, 251, 257, 281, 283, 307, 313, 331, 337, 347, 353, 379, 401, 409, 419, 433, 443, 449, 457, 467, 491, 499, 521, 523, 529, 547, 563
Offset: 1

Views

Author

Jianing Song, Feb 19 2021

Keywords

Comments

Also norms of prime ideals in Z[sqrt(-2)], which is a unique factorization domain. The norm of a nonzero ideal I in a ring R is defined as the size of the quotient ring R/I.
Consists of the primes congruent to 1, 2, 3 modulo 8 and the squares of primes congruent to 5, 7 modulo 8.
For primes p == 1, 3 (mod 8), there are two distinct ideals with norm p in Z[sqrt(2)], namely (x + y*sqrt(-2)) and (x - y*sqrt(-2)), where (x,y) is a solution to x^2 + 2*y^2 = p; for p = 2, (sqrt(-2)) is the unique ideal with norm p; for p == 5, 7 (mod 8), (p) is the only ideal with norm p^2.

Examples

			norm(1 + sqrt(-2)) = norm(1 + sqrt(-2)) = 3;
norm(3 + sqrt(-2)) = norm(3 + sqrt(-2)) = 11;
norm(3 + 2*sqrt(-2)) = norm(3 + 2*sqrt(-2)) = 17;
norm(1 + 3*sqrt(-2)) = norm(1 + 3*sqrt(-2)) = 19.
		

Crossrefs

The number of nonassociative elements with norm n (also the number of distinct ideals with norm n) is given by A002325.
The total number of elements with norm n is given by A033715.
Norms of prime ideals in O_K, where K is the quadratic field with discriminant D and O_K be the ring of integers of K: A055673 (D=8), A341783 (D=5), A055664 (D=-3), A055025 (D=-4), A090348 (D=-7), this sequence (D=-8), A341785 (D=-11), A341786 (D=-15*), A341787 (D=-19), A091727 (D=-20*), A341788 (D=-43), A341789 (D=-67), A341790 (D=-163). Here a "*" indicates the cases where O_K is not a unique factorization domain.

Programs

  • PARI
    isA341784(n) = my(disc=-8); (isprime(n) && kronecker(disc,n)>=0) || (issquare(n, &n) && isprime(n) && kronecker(disc,n)==-1)
Showing 1-10 of 38 results. Next