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.
%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