OFFSET
2,2
COMMENTS
By Wolstenholme's theorem, if p > 3 is prime, then a(p) = p^2.
Conjecture: for n > 3, if a(n) = n^2, then n is a prime.
Note: the weak pseudoprimes n such that a(n) = n are not known.
Composite numbers m <> p^2 for which a(m) > 1 are the same as in A309391: 88, 1290, 9339, ...
LINKS
R. Meštrović, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862--2012), arXiv:1111.3057 [math.NT], 2011.
Eric Weisstein's World of Mathematics, Wolstenholme's Theorem.
Wikipedia, Wolstenholme's theorem.
FORMULA
a(n) = A309391(n) for composite n.
a(p) = p^2 for every prime p > 3.
a(p^2) = p iff p > 3 is a prime.
EXAMPLE
a(11) = gcd(11^2, A001008(11-1)) = gcd(121, 7381) = 121.
MATHEMATICA
a[n_] := GCD[n^2, Numerator[HarmonicNumber[n-1]]]; Array[a, 72, 2]
PROG
(Magma) [Gcd(k^2, Numerator(HarmonicNumber(k-1))):k in [2..80]]; // Marius A. Burtea, Jul 28 2019
(Python)
from sympy import gcd, harmonic
def A309387(n):
return gcd(n**2, harmonic(n-1).p) # Chai Wah Wu, Jul 31 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar and Thomas Ordowski, Jul 28 2019
STATUS
approved