OFFSET
0,1
COMMENTS
a(n) > n, with a(n)/n -> 1 as n -> infinity. - Robert Israel, May 28 2018
LINKS
Robert Israel, Table of n, a(n) for n = 0..7485
FORMULA
Let S(n, m) = Sum_{k=1..m} k^n. Define a(n) = min{ m | S(n, m)<m! }.
EXAMPLE
a(3)=6 because S(3,6)=441<720=6! but S(3,5)=225>=120=5! and so for S(3,j), j=0,1,2,3,4.
MAPLE
f:= proc(n) local k, L;
L:= sum(k^n, k=1..n);
for k from n+1 do
L:= L + k^n;
if L < k! then return k fi
od
end proc:
map(f, [$0..100]); # Robert Israel, May 28 2018
PROG
(PARI) a(n) = {my(s = 0, ok = 0, m = 1); until (ok, s += m^n; if (s < m!, ok = 1, m++); ); return (m); } \\ Michel Marcus, Jul 15 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Hieronymus Fischer, Jan 19 2006
STATUS
approved