login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A128050
Position of start of first occurrence of prime(n) after the decimal point in expansion of golden ratio phi.
2
19, 5, 22, 10, 34, 55, 35, 188, 131, 174, 137, 98, 90, 27, 102, 111, 166, 1, 150, 217, 479, 44, 25, 13, 81, 458, 1242, 744, 563, 96, 1602, 186, 97, 995, 259, 939, 1999, 1204, 641, 1191, 43, 833, 1682, 2833, 2708, 188, 647, 130, 62, 734, 2337, 1106, 307, 1156, 2532
OFFSET
1,1
EXAMPLE
Golden ratio phi = 1.6180339887498948482045868343656381177... (see A001622).
First occurrence of prime(1) = 2 is at the 19th digit after the decimal point, hence a(1) = 19.
First occurrence of prime(5) = 11 starts at the 34th digit after the decimal point, hence a(5) = 34.
MATHEMATICA
Module[{p = Rest[First[RealDigits[GoldenRatio, 10, 10^4]]], n = 0, a}, Reap[While[(a = SequencePosition[p, IntegerDigits[Prime[++n]], 1]) != {}, Sow[a[[1, 1]]]]][[2, 1]]] (* Paolo Xausa, Aug 01 2024 *)
PROG
(Magma) k:=3000; R := RealField(k); [ Position(IntegerToString(Round(10^k*(-1 + (Sqrt(elt<R | 5, 0>)+1) / elt<R | 2, 0>))), IntegerToString(NthPrime(n))) : n in [1..55] ]; /* Klaus Brockhaus, Feb 15 2007 */
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Gregory Allen, Feb 13 2007
EXTENSIONS
Edited, corrected and extended by Klaus Brockhaus, Feb 15 2007
STATUS
approved