|
|
A055666
|
|
Number of inequivalent Eisenstein-Jacobi primes of successive norms (indexed by A055664).
|
|
5
|
|
|
1, 1, 2, 2, 2, 1, 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, 1, 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, 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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,3
|
|
COMMENTS
|
These are the primes in the ring of integers a+b*omega, a and b rational integers, omega = (1+sqrt(-3))/2.
Two primes are considered equivalent if they differ by multiplication by a unit (+-1, +-omega, +-omega^2).
|
|
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.
|
|
LINKS
|
|
|
EXAMPLE
|
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.
|
|
MATHEMATICA
|
norms = Join[{3}, Select[Range[2000], (PrimeQ[#] && Mod[#, 6] == 1) || (PrimeQ[Sqrt[#]] && Mod[Sqrt[#], 3] == 2) &]]; r[n_] := Length[Reduce[n == a^2 - a*b + b^2, {a, b}, Integers]]/6; A055666 = r /@ norms (* Jean-François Alcover, Oct 24 2013 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|