OFFSET
1,1
REFERENCES
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Jonathan Webster, Table of n, a(n) for n = 1..375 (terms 1..200 from H. C. Williams)
E. F. Ecklund, Jr. et al., A new function associated with the prime factors of C(n,k), Math. Comp., 28 (1974), 647-649.
R. F. Lukes, R. Scheidler, and H. C. Williams, Further Tabulation of the Erdős-Selfridge Function, Math. Comput. 66 (1997) 1709-1717.
R. Scheidler and H. C. Williams, A method for tabulating the number-theoretic function g(k), Math. Comp. 59 (199) (1992) 251-257.
Brianna Sorenson, Jonathan P Sorenson, and Jonathan Webster, An Algorithm and Estimates for the Erdős-Selfridge Function (work in progress), arXiv:1907.08559 [math.NT], 2019.
Eric Weisstein's World of Mathematics, Erdős-Selfridge function.
MAPLE
A003458 := proc(n)
local m;
for m from n+2 do
if A020639( binomial(m, n)) > n then
return m ;
end if;
end do:
end proc:
seq(A003458(n), n=1..16) ; # R. J. Mathar, Mar 27 2024
MATHEMATICA
f[n_] := Block[{k = n + 2, p = Table[Prime[i], {i, 1, PrimePi[n]}]}, While[ First[ Sort[ Mod[ Binomial[k, n], p]]] == 0, k++ ]; k]; Table[ f[n], {n, 1, 40}]
(* Alternative: *)
esf[n_]:=Module[{m=n+2}, While[FactorInteger[Binomial[m, n]][[1, 1]]<=n, m++]; m]; Array[esf, 50] (* Harvey P. Dale, Nov 03 2013 *)
PROG
(PARI) a(n) = local(m, i, f); m=0; i=n+1; while(m<=n, i=i+1; m=factor(binomial(i, n))[1, 1]); i /* Ralf Stephan */
(PARI) primorial(x)=prod(i=1, primepi(x), prime(i))
a(n)=my(G=primorial(n), N=n+1); while(gcd(binomial(N++, n), G)!=1, ); N \\ Charles R Greathouse IV, Apr 28 2026
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
EXTENSIONS
Extended by Robert G. Wilson v, Dec 01 2002
STATUS
approved
