OFFSET
1,1
COMMENTS
If n is a term of the sequence, then the parallelogram formed by the vectors {n,phi(n)},{n+1,phi(n+1)} has the same area as that of an integral square.
EXAMPLE
Det[{{26,phi(26)},{27,phi(27)}}] = Det[{{26,12},{27,18}}] = 12^2, so 26 is a term of the sequence.
MATHEMATICA
f[n_] := Det[{{n, EulerPhi[n]}, {n + 1, EulerPhi[n + 1]}}]; Do[If[IntegerQ @ Sqrt @ f[n], Print[n]], {n, 1, 10^5}]
PROG
(PARI) isok(n) = issquare(matdet([n, eulerphi(n); n+1, eulerphi(n+1)])); \\ Michel Marcus, Sep 26 2019
(Magma) [k:k in [1..5000000]|IsSquare(k*EulerPhi(k+1)-(k+1)*EulerPhi(k))]; // Marius A. Burtea, Sep 26 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Jan 30 2002
EXTENSIONS
a(18)-a(37) from Amiram Eldar, Sep 26 2019
STATUS
approved