Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #35 Nov 09 2023 20:13:20
%S 0,1,3,4,7,8,9,12,13,15,27,31,49,57,60,75,139,147,283,327,488,604,700,
%T 825,908,1051,1064,1215,5319,9669,10136,16675,25656,28933,35864,47671,
%U 68028,73380,186223,194965,221649,233059,240644,513007,543128,551491,648872,989124,994536
%N Numbers k such that 3^k + 16 is prime.
%C Indices of primes in A205646.
%C a(50) > 10^6. - _Robert Price_, Oct 28 2020
%F {n: 3^n + 16 is in A000040} = {n: 3^n + 16 is prime} = {n: A000244(n) is prime} = {n: A205646(n) is prime}.
%e 57 is in the sequence because (3^57) + 16 = 1570042899082081611640534579 is prime.
%t Select[Range[0, 2000], PrimeQ[3^# + 16] &] (* _T. D. Noe_, Jan 30 2012 *)
%o (PARI) is(n)=isprime(3^n+16) \\ _Charles R Greathouse IV_, Feb 17 2017
%Y Cf. Sequences of numbers k such that 3^k + m is prime:
%Y (m = 2) A051783, (m = -2) A014224, (m = 4) A058958, (m = -4) A058959,
%Y (m = 8) A217136, (m = -8) A217135, (m = 10) A217137, (m = -10) A217347,
%Y (m = 14) A219035, (m = -14) A219038, (m = 16) A205647, (m = -16) A219039,
%Y (m = 20) A219040, (m = -20) A219041, (m = 22) A219042, (m = -22) A219043,
%Y (m = 26) A219044, (m = -26) A219045, (m = 28) A219046, (m = -28) A219047,
%Y (m = 32) A219048, (m = -32) A219049, (m = 34) A219050, (m = -34) A219051. Note that if m is a multiple of 3, 3^k + m is also a multiple of 3 (for k greater than 0), and as such isn't prime.
%K nonn,hard,more
%O 1,3
%A _Jonathan Vos Post_, Jan 30 2012
%E 5319 from _Nicolas M. Perrault_, Nov 10 2012
%E a(30)-a(40) from _Robert Price_, Oct 23 2013
%E a(41) discovered by Lelio R Paula, November 2016
%E a(42)-a(43) from _Robert Price_, Apr 02 2020
%E a(44)-a(46) from _Robert Price_, May 14 2020
%E a(47)-a(49) from _Robert Price_, Oct 28 2020