%I #21 Sep 08 2022 08:46:11
%S 331,82171,3751197451
%N Primes of the form 166^k - 165^k.
%C a(4) has 1641 digits and a(5) has 1717 digits. - _Vincenzo Librandi_, Feb 22 2015
%C The values of k are 2, 3, 5, 739, 773, ... - _Jinyuan Wang_, May 30 2020
%H Michael P. May, <a href="https://arxiv.org/abs/1510.01028">On the Existence and Frequency Distribution of the Shell Primes</a>, arXiv preprint arXiv:1510.01028 [math.CO], 2015.
%t Select[Table[166^n - 165^n, {n, 1000}], PrimeQ] (* _Vincenzo Librandi_, Feb 22 2015 *)
%o (Magma) [a: n in [0..600] | IsPrime(a) where a is 166^n-165^n]; // _Vincenzo Librandi_, Feb 22 2015
%o (PARI) select(m->ispseudoprime(m), vector(10, k, 166^k-165^k)) \\ _Jinyuan Wang_, May 30 2020
%Y Cf. A254298.
%K nonn,bref
%O 1,1
%A _Michael P. May_, Feb 21 2015