[go: up one dir, main page]

login
A156165
Numbers k such that k![7]+1 is prime (n![7] = A114799(n) = septuple factorial).
11
0, 1, 2, 4, 6, 9, 10, 12, 13, 24, 25, 26, 29, 31, 35, 36, 47, 49, 57, 58, 64, 71, 73, 75, 78, 80, 97, 123, 125, 129, 131, 135, 147, 150, 159, 183, 201, 250, 251, 255, 298, 336, 337, 458, 467, 556, 570, 657, 743, 801, 908, 925, 1003, 1209, 1473, 1524, 1716, 1881, 1926
OFFSET
1,3
COMMENTS
a(103) > 50000. - Robert Price, Sep 03 2012
LINKS
MATHEMATICA
mf[n_, k_] := Product[n - i k, {i, 0, Quotient[n - 2, k]}];
Reap[For[k = 0, k <= 2000, k++, If[PrimeQ[mf[k, 7] + 1], Sow[k]]]][[2, 1]] (* Jean-François Alcover, Feb 26 2019 *)
Select[Range[0, 2000], PrimeQ[Times@@Range[#, 1, -7]+1]&] (* Harvey P. Dale, Aug 21 2021 *)
PROG
(PARI) mf(n, k=7)=prod(i=0, (n-2)\k, n-i*k)
for( n=0, 9999, ispseudoprime(mf(n)+1) & print1(n", "))
KEYWORD
nonn
AUTHOR
M. F. Hasler, Feb 10 2009
STATUS
approved