OFFSET
1,2
COMMENTS
Indices of primes occurring in A107900.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Continued Fraction.
EXAMPLE
17 is in this sequence because the 17th prime is 59, 59/17 = 3 + 1/(2 + 1/8), the terms of the continued fraction are [3; 2, 8], and the sum of the terms is the prime number 13.
MAPLE
filter:= n -> isprime(convert(Term(NumberTheory:-ContinuedFraction(ithprime(n)/n), all), `+`)):
select(filter, [$1..300]); # Robert Israel, Jan 08 2024
MATHEMATICA
s={}; Do[If[PrimeQ[Plus@@ContinuedFraction[Prime[n]/n]], s=Append[s, n]], {n, 300}]; A107899=s
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, May 26 2005
STATUS
approved