OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Antti Karttunen and J. Moyer, C-program for computing the initial terms of this sequence
MAPLE
R:= NULL: count:= 0:
for n from 1 while count < 100 do
p:= floor(n*phi);
if isprime(p) then R:= R, p; count:= count+1 fi
od:
R; # Robert Israel, Jan 17 2023
MATHEMATICA
(See A184792.)
PROG
(Python)
from math import isqrt
from itertools import count, islice
from sympy import isprime
def A095280_gen(): # generator of terms
return filter(isprime, ((n+isqrt(5*n**2)>>1) for n in count(1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 04 2004
STATUS
approved