OFFSET
1,1
COMMENTS
MAPLE
a := proc(n) local x, x0, y, k, bound; bound := 1000;
x := 2*n + 1; x0 := x;
y := binomial(4*n + 1, 2*n);
for k from 0 to bound while y mod x = 1 do
x := x * x0 od;
if k < bound then k else print("No k below ", bound) fi end:
seq(a(n), n = 1..100); # Peter Luschny, Apr 22 2021
PROG
(PARI) a(n) = my(x=2*n+1, b=binomial(2*x-1, x-1)); for(k=1, oo, if(Mod(b, x^k)!=1, return(k-1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Apr 18 2021
STATUS
approved