OFFSET
1,1
COMMENTS
Contains no primes. - Robert Israel, Aug 07 2018
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
It seems that a(n) is asymptotic to C*n with C=2.8...
MAPLE
filter:= proc(n) local t, p, L, f, i, c;
for f in ifactors(numtheory:-sigma(n))[2] do
p:= f[1];
L:= convert(n, base, p);
t:= 0; c:= 0;
for i from 1 to nops(L) do
if 2*L[i]+c >= p then
c:= 1; t:= t+1; if t >= f[2] then break fi;
else
c:= 0;
fi
od:
if t < f[2] then return true fi;
od;
false
end proc:
select(filter, [$1..1000]); # Robert Israel, Aug 07 2018
MATHEMATICA
fQ[n_] := Mod[ Binomial[2n, n], DivisorSigma[1, n]] > 0; Select[ Range@164, fQ] (* Robert G. Wilson v, Aug 07 2018 *)
PROG
(PARI) isok(n) = binomial(2*n, n) % sigma(n); \\ Michel Marcus, Nov 28 2013
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Aug 21 2002
STATUS
approved