OFFSET
1,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
Robert P. Schneider, A golden product identity for e, Mathematics Magazine, Vol. 87, No. 2 (2014), pp. 132-134.
FORMULA
a(n) = n if and only if n is a prime.
Product_{n>=1} (1 - 1/phi^n)^(-a(n)/n) = e, where phi is the golden ratio (A001622) (Schneider, 2014). - Amiram Eldar, Jun 11 2022
MATHEMATICA
a[ n_] := EulerPhi @ n - MoebiusMu @ n; Table[a[n], {n, 70}] (* Michael Somos, Jul 19 2011 *)
PROG
(PARI) {a(n) = if( n==0, 0, eulerphi( n) - moebius( n))} /* Michael Somos, Jul 19 2011 */
(Magma) [EulerPhi(n) - MoebiusMu(n): n in [1..100]]; // G. C. Greubel, Sep 03 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Mar 25 2000
STATUS
approved