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 A153440 #17 Sep 08 2022 08:45:39 %S A153440 1,2,11,44,45,56,62,63,110,170,219,234,245,261,263,333,395,398,402, %T A153440 413,428,434,437,498,557,558,578,633,692,695,723,731,750,761,774,794, %U A153440 797,804,806,846,854,855,863,906,923,926,977,1046,1085,1086 %N A153440 Numbers k such that k^9*(k^9+1)+1 is prime. %C A153440 It seems numbers of the form k^n*(k^n+1)+1 with n > 0, k > 1 may be primes only if n has the form 3^j. When n is even, k^(4*n)+k^(2*n)+1=(k^(2*n)+1)^2-(k^n)^2=(k^(2*n)+k^n+1)*(k^(2*n)-k^n+1) so composite. But why if n odd > 3 and not a power of 3, k^n*(k^n+1)+1 is always composite? %H A153440 Pierre CAMI, <a href="/A153440/b153440.txt">Table of n, a(n) for n=1,...,38019</a> %t A153440 k9pQ[n_]:=Module[{c=n^9},PrimeQ[c(c+1)+1]]; Select[Range[1200],k9pQ] (* _Harvey P. Dale_, Oct 29 2014 *) %t A153440 Select[Range[1100], PrimeQ[(#^9 (#^9 + 1)) + 1] &] (* _Vincenzo Librandi_, Jan 17 2015 *) %o A153440 (Magma) [n: n in [0..1100] | IsPrime(n^9*(n^9+1)+1)]; // _Vincenzo Librandi_, Jan 17 2015 %Y A153440 Cf. A153438. %K A153440 nonn %O A153440 1,2 %A A153440 _Pierre CAMI_, Dec 26 2008