[go: up one dir, main page]

login
A136683
Numbers k such that A136675(k) is prime.
8
2, 3, 4, 5, 6, 9, 20, 21, 29, 119, 132, 151, 351, 434, 457, 462, 572, 611, 930, 1107, 1157, 1452, 1515, 2838, 3997, 5346, 6463, 6725, 7664, 10234, 14168, 14299
OFFSET
1,1
COMMENTS
A136675(k) = Numerator of Sum_{j=1..k} (-1)^(j+1)/j^3.
LINKS
Eric Weisstein's World of Mathematics, Harmonic Number
MATHEMATICA
Do[ f=Numerator[ Sum[ (-1)^(k+1)*1/k^3, {k, 1, n} ] ]; If[ PrimeQ[f], Print[ {n, f} ] ], {n, 1, 151} ]
Flatten[Position[Numerator[Accumulate[Table[(-1)^(k+1) 1/k^3, {k, 3000}]]], _?PrimeQ] ] (* Harvey P. Dale, Feb 12 2013 *)
PROG
(PARI) isok(n) = ispseudoprime(numerator(sum(k=1, n, (-1)^(k+1) / k^3))); \\ Daniel Suteu, Mar 15 2019
KEYWORD
nonn,more
AUTHOR
Alexander Adamchuk, Jan 16 2008
EXTENSIONS
More terms from Harvey P. Dale, Feb 12 2013
a(25)-a(28) from Amiram Eldar, Mar 15 2019
a(29)-a(32) from Robert Price, Apr 22 2019
STATUS
approved