OFFSET
1,2
COMMENTS
Each term is a product of a finite subsequence of A005478.
PROG
(Python)
import itertools, math, sympy
def fibprimegen(limit): # Generate Fibonacci primes <= limit
a, b = 1, 2
while b <= limit:
if sympy.isprime(b):
yield b
a, b = b, a+b
LIMIT=1000000
fibprimes=list(fibprimegen(LIMIT))
fibprimeseqs=itertools.chain.from_iterable(
itertools.combinations(fibprimes, n) for n in range(len(fibprimes)+1))
print(sorted(a for a in map(math.prod, fibprimeseqs) if a <= LIMIT))
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jack Brennen, Oct 04 2024
STATUS
approved