OFFSET
1,1
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..1000
FORMULA
Primes p(j) such that (p(j)+p(j+1))/2 = q(m)^w, where q(m) is a prime.
EXAMPLE
n = prime(9750374) = 174689077, next prime = 174689101, mean = 174689089 = 13217^2, a prime power. The arithmetic mean of two consecutive primes is never prime, while between two consecutive primes, prime powers occur. These prime powers are in the middle of gap: p+d/2 = q^w. The prime power is most often square and very rarely occurs more than once (see A053706).
MATHEMATICA
fi[x_] := FactorInteger[x] ff[x_] := Length[FactorInteger[x]] Do[s=(Prime[n]+Prime[n+1])/2; s1=ff[s]; If[Equal[s1, 1], Print[{n, p=Prime[n], s, fi[s], s-p, s1}]], {n, 1, 10000000}]
Select[Partition[Prime[Range[25*10^5]], 2, 1], PrimePowerQ[Mean[#]]&][[;; , 1]] (* Harvey P. Dale, Oct 15 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, May 26 2003
STATUS
approved