OFFSET
0,3
LINKS
Antti Karttunen, Table of n, a(n) for n = 0..4096
EXAMPLE
PROG
(PARI)
A249151(n) = { my(uplim, padicvals, b); uplim = (n+3); padicvals = vector(uplim); for(k=0, n, b = binomial(n, k); for(i=1, uplim, padicvals[i] += valuation(b, prime(i)))); k = 1; while(k>0, for(i=1, uplim, if((padicvals[i] -= valuation(k, prime(i))) < 0, return(k-1))); k++); };
\\ Alternative implementation:
A001142(n) = prod(k=1, n, k^((k+k)-1-n));
A055881(n) = { my(i); i=2; while((0 == (n%i)), n = n/i; i++); return(i-1); }
for(n=0, 4096, write("b249151.txt", n, " ", A249151(n)));
CROSSREFS
One more than A249150.
Cf. A249423 (numbers k such that a(k) = k+1).
Cf. A249429 (numbers k such that a(k) > k).
Cf. A249433 (numbers k such that a(k) < k).
Cf. A249434 (numbers k such that a(k) >= k).
Cf. A249424 (numbers k such that a(k) = (k-1)/2).
Cf. A249425 (record positions).
Cf. A249427 (record values).
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 25 2014
STATUS
approved