OFFSET
1,2
LINKS
Michael S. Branicky, Table of n, a(n) for n = 1..101
MAPLE
f:= proc(n) option remember; `if`(n<1, 0, f(n-1)+(n+1)/n) end:
q:= n-> isprime(numer(f(n))):
select(q, [$1..2000])[]; # Alois P. Heinz, Aug 27 2024
PROG
(PARI) ss(n)=sum(i=1, n, (i+1)/i); for (i=1, 1000, if(isprime(numerator(ss(i))), print1(i, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 07 2004
EXTENSIONS
a(30) and beyond from Michael S. Branicky, Aug 26 2024
STATUS
approved