OFFSET
1,1
COMMENTS
Records in A002618. - Artur Jasinski, Jan 23 2008
Also records in A174857. - Vladimir Shevelev, Mar 31 2010
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = prime(n) * (prime(n) - 1).
a(n) = prime(n) * phi(prime(n)). - Artur Jasinski, Jan 23 2008
From Reinhard Zumkeller, Sep 17 2011: (Start)
a(n) = A009262(prime(n)). - Enrique Pérez Herrero, May 12 2012
a(n) = prime(n)! mod (prime(n)^2). - J. M. Bergot, Apr 10 2014
a(n) = 2*A008837(n). - Antti Karttunen, May 01 2015
Sum_{n>=1} 1/a(n) = A136141. - Amiram Eldar, Nov 09 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(2)*zeta(3)/zeta(6) (A082695).
Product_{n>=1} (1 - 1/a(n)) = A005596. (End)
EXAMPLE
2*1, 3*2, 5*4, 7*6, 11*10, 13*12, 17*16, ...
MAPLE
A036689 := proc(n) local p ; p := ithprime(n) ; p*(p-1) ; end proc: # R. J. Mathar, Apr 11 2011
MATHEMATICA
Table[Prime[n] EulerPhi[Prime[n]], {n, 100}] (* Artur Jasinski, Jan 23 2008 *)
Table[Prime[n] (Prime[n] - 1), {n, 1, 50}] (* Bruno Berselli, Apr 22 2014 *)
#(#-1)&/@Prime[Range[50]] (* Harvey P. Dale, Sep 08 2019 *)
PROG
(Magma) [ n*(n-1): n in PrimesUpTo(220) ]; // Bruno Berselli, Apr 11 2011
(PARI) forprime(p=2, 1e3, print1(p^2-p", ")) \\ Charles R Greathouse IV, Jun 10 2011
(PARI) A036689(n) = ((p->(p-1)*p)(prime(n))); \\ Antti Karttunen, Dec 14 2024
(Haskell)
a036689 n = a036689_list !! (n-1)
a036689_list = zipWith (*) a000040_list $ map pred a000040_list
-- Reinhard Zumkeller, Sep 17 2011
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Deleted two incorrect comments. - N. J. A. Sloane, May 07 2020
STATUS
approved