OFFSET
1,1
LINKS
A. Karttunen and J. Moyer, C-program for computing the initial terms of this sequence
PROG
(Python)
from sympy import fibonacci, primerange
def a(n):
k=0
x=0
while n>0:
k=0
while fibonacci(k)<=n: k+=1
x+=10**(k - 3)
n-=fibonacci(k - 1)
return x
def ok(n): return str(a(n))[-2:]=="00"
print([n for n in primerange(1, 1001) if ok(n)]) # Indranil Ghosh, Jun 08 2017
(PARI) list(lim)=my(v=List(), w=quadgen(20), phi=(1+w)/2, p2=phi^2, x=(2*phi-2)*p2, q); lim=lim\1+1; while(x<lim, if(isprime(q=x\1), listput(v, q)); x+=p2); Vec(v) \\ Charles R Greathouse IV, Nov 10 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 01 2004
STATUS
approved