cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A168033 Primes p such that floor(phi^p) is prime.

This page as a plain text file.
%I A168033 #22 Sep 08 2022 08:45:48
%S A168033 2,5,7,11,13,17,19,31,37,41,47,53,61,71,79,113,313,353,503,613,617,
%T A168033 863,1097,1361,4787,4793,5851,7741,8467,10691,12251,13963,14449,19469,
%U A168033 35449,36779,44507,51169,56003,81671,89849,94823,140057,148091,159521,183089,193201,202667,344293,387433,443609,532277,574219,616787,631181,637751,651821,692147,901657,1051849
%N A168033 Primes p such that floor(phi^p) is prime.
%C A168033 Primes in A059791. - _Charles R Greathouse IV_, Jul 29 2011
%C A168033 Also primes in A001606. - _Michel Marcus_, Oct 21 2016
%t A168033 $MaxExtraPrecision=6!; Select[Prime[Range[5! ]],PrimeQ[Floor[GoldenRatio^# ]]&]
%o A168033 (PARI) phi=(1+sqrt(5))/2;forprime(p=2,1e3,if(isprime(floor(phi^p)),print1(p", "))) \\ _Charles R Greathouse IV_, Jul 29 2011
%o A168033 (Magma) [p: p in PrimesUpTo(2000)| IsPrime(Lucas(p))]; // _Vincenzo Librandi_, Jul 11 2019
%Y A168033 Cf. A001606, A059791.
%K A168033 nonn
%O A168033 1,1
%A A168033 _Vladimir Joseph Stephan Orlovsky_, Nov 17 2009
%E A168033 a(22)-a(32) from _Charles R Greathouse IV_, Jul 29 2011
%E A168033 More terms (using A001606) from _Joerg Arndt_, Jul 11 2019