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 A113501 #29 Feb 16 2025 08:32:59
%S A113501 1,2,3,9,14,23,29,81,128,210,468,473,746,950,3344,4043,4839,14376,
%T A113501 39521,64563,72984,82899,84338,85206,86121,139160
%N A113501 Indices of prime NSW numbers A088165.
%C A113501 Very closely related to indices of prime Pell-Lucas numbers (A099088).
%C A113501 a(27) > 221400. - _Robert Price_, Mar 29 2019
%H A113501 H. Li, T. MacHenry, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL16/MacHenry/machenry7.html">Permanents and Determinants, Weighted Isobaric Polynomials, and Integer Sequences</a>, J. Int. Seq. 16 (2013) #13.3.5, example 47.
%H A113501 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/NSWNumber.html">NSW Number</a>
%H A113501 Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/IntegerSequencePrimes.html">Integer Sequence Primes</a>
%e A113501 NSW(1) = 7, NSW(2) = 41, NSW(3) = 239, NSW(9) = 9369319, ...
%t A113501 nsw = LinearRecurrence[{6, -1}, {7, 41}, 1000]; Position[nsw, _?(PrimeQ[#] &)] // Flatten (* _Amiram Eldar_, Dec 07 2018 *)
%o A113501 (PARI) isok(n) = my(w=3+quadgen(32)); isprime(imag((1+w)*w^n)); \\ _Michel Marcus_, Dec 07 2018
%Y A113501 Cf. A002315, A088165, A099088.
%K A113501 nonn,more
%O A113501 1,2
%A A113501 _Eric W. Weisstein_, Jan 09 2006
%E A113501 a(19)-a(20) from _Eric W. Weisstein_, May 22 2006
%E A113501 a(21) from _Eric W. Weisstein_, Aug 29 2006
%E A113501 a(22) from _Eric W. Weisstein_, Nov 11 2006
%E A113501 a(23) from _Eric W. Weisstein_, Nov 26 2006
%E A113501 a(24) from _Eric W. Weisstein_, Dec 10 2006
%E A113501 a(25) from _Eric W. Weisstein_, Jan 25 2007
%E A113501 a(26) from _Robert Price_, Dec 07 2018