A249951 Numbers n such that A113630(n) is prime.
3, 5, 8, 23, 42, 62, 63, 75, 90, 98, 177, 192, 207, 213, 222, 228, 233, 263, 288, 297, 317, 320, 402, 453, 455, 528, 570, 602, 620, 650, 672, 752, 780, 797, 810, 863, 867, 870, 875, 912, 932, 935, 978, 1010, 1043, 1065, 1067, 1070, 1110, 1125, 1133, 1142, 1190
Offset: 1
Keywords
Examples
See example section of A113630.
Links
- Chai Wah Wu, Table of n, a(n) for n = 1..2434
Programs
-
Haskell
a249951 n = a249951_list !! (n-1) a249951_list = filter ((== 1) . a010051'' . a113630) [1..] -- Reinhard Zumkeller, Nov 22 2014
-
Magma
[n: n in [0..1500] | IsPrime(1+2*n+3*n^2+4*n^3+ 5*n^4+6*n^5+7*n^6+8*n^7+9*n^8)]; // Vincenzo Librandi, Nov 09 2014
-
PARI
for(n=1,5000,s=sum(i=1,9,i*n^(i-1));if(ispseudoprime(s),print1(n,", "))) \\ Derek Orr, Nov 09 2014
-
Python
from sympy import isprime A249951_list, m = [], [362880, -1229760, 1607760, -1011480, 309816, -40752, 1584, -4, 1] for n in range(1,10**5+1): for i in range(8): m[i+1]+= m[i] if isprime(m[-1]): A249951_list.append(n)
Comments