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.

A259145 Numbers k such that k^2 - phi(k) is prime, where phi() is A000010.

This page as a plain text file.
%I A259145 #30 Sep 08 2022 08:46:13
%S A259145 2,3,7,13,33,35,65,67,77,79,91,133,139,151,163,193,221,247,249,287,
%T A259145 299,321,337,341,349,377,379,437,457,481,533,541,551,561,581,591,595,
%U A259145 611,613,643,721,727,763,769,779,789,803,817,843,851,869,917,919,991
%N A259145 Numbers k such that k^2 - phi(k) is prime, where phi() is A000010.
%C A259145 Conjecture: a(n) is a cyclic number (see A003277) for all n.
%C A259145 A065508 is the subsequence of prime terms. - _Michel Marcus_, Jun 19 2015
%H A259145 Amiram Eldar, <a href="/A259145/b259145.txt">Table of n, a(n) for n = 1..10000</a>
%e A259145 a(1) = 2, since phi(2) = 1, thus 2^2 - 1 = 3 (prime).
%e A259145 a(3) = 7, since phi(7) = 6, thus 7^2 - 6 = 43 (prime).
%e A259145 a(5) = 33, since phi(33) = 20, thus 33^2 - 20 = 1069 (prime).
%t A259145 Select[Range[2000], PrimeQ[#^2 - EulerPhi[#]] &]
%o A259145 (Magma) [n: n in [1..1000] | IsPrime(n^2 - EulerPhi(n))]; // _Vincenzo Librandi_, Jun 21 2015
%o A259145 (PARI) main(size)={ v=vector(size); i=0; m=1; while(i<size, if(isprime(m^2-eulerphi(m)),v[i++]=m);m++;); return(v);} /* _Anders Hellström_, Jul 08 2015 */
%Y A259145 Cf. A000010, A003277, A065508, A258435.
%K A259145 nonn,easy
%O A259145 1,1
%A A259145 _Carlos Eduardo Olivieri_, Jun 19 2015