OFFSET
1,1
LINKS
K. Brockhaus, Table of n, a(n) for n = 1..108.
Florian Luca and Igor E. Shparlinski, On the largest prime factor of n! + 2^n - 1, Journal de Theorie des Nombres de Bordeaux 17 (2005), 859-870.
MAPLE
# this code will produce FAIL for terms that can't be computed within a
# 30 second time limit
spf:= proc(n) local F;
F:= select(f -> f[1]::integer, ifactors(n, easy)[2]);
if F = [] then F:= traperror(timelimit(30, ifactors(n)[2])) fi;
if F = lasterror then FAIL else min(seq(f[1], f=F)) fi
end proc:
seq(spf(n!+2^n-1), n=1..100); # Robert Israel, Sep 18 2016
MATHEMATICA
a = {}; Do[AppendTo[a, n! + 2^n - 1], {n, 1, 40}]; b = {}; Do[c = FactorInteger[a[[n]]]; d = c[[1]]; AppendTo[b, d[[1]]], {n, 1, Length[a]}]; b
Table[FactorInteger[n! + 2^n - 1] [[1, 1]], {n, 1, 40}] (* Vincenzo Librandi, Sep 18 2016 *)
PROG
(Magma) trialdiv:=function(n, P) val:=0; for p in P do if n mod p eq 0 then val:=p; break; end if; end for; return val; end function; P:=PrimesUpTo(300000000); [ trialdiv(a, P) where a is Factorial(n)+2^n-1: n in [1..70] ]; //a(61) requires a separate computation. // Klaus Brockhaus, Nov 06 2009
(PARI) a(n)=factor(n!+2^n-1)[1, 1] \\ Charles R Greathouse IV, Feb 01 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Apr 06 2008, corrected Apr 22 2008
EXTENSIONS
a(1) - a(40) verified and a(41) - a(70) added by Klaus Brockhaus, Nov 06 2009
STATUS
approved