%I #12 Sep 08 2022 08:46:03
%S 1,3,5,9,65,71,225,269,651,843,1971,4851,13667,13961,33309,168111
%N Numbers k such that 3^k + k^3 - 1 is prime.
%C a(17) > 2*10^5. - _Robert Price_, Apr 15 2014
%t Select[Range[0, 3000], PrimeQ[3^# + #^3 - 1] &]
%o (Magma) [n: n in [0..800] | IsPrime(3^n + n^3 - 1)]
%o (PARI) is(n)=ispseudoprime(3^n+n^3-1) \\ _Charles R Greathouse IV_, Jun 13 2017
%Y Cf. A215435.
%K nonn
%O 1,2
%A _Vincenzo Librandi_, Sep 03 2012
%E a(12)-a(16) from _Robert Price_, Apr 15 2014