OFFSET
1,1
COMMENTS
Even numbers n such that n divides binomial(n, [n/2]) and A010551(n) does not divide j!*(n-j)! exactly 7 times for j = 0..n. - Peter Luschny, Aug 04 2017
LINKS
David A. Corneth, Table of n, a(n) for n = 1..274 (Terms <= 60000)
MAPLE
isa := proc(n) local bn, bm;
if n mod 2 = 0 then bn := binomial(n, iquo(n, 2)):
if modp(bn, n) = 0 then
bm := (n, j) -> `if`(modp(bn, binomial(n, j)) = 0, 1, 0):
return 1 <> add(bm(n, j), j=2..iquo(n, 2)-1)
fi fi; false end:
select(isa, [$1..5000]); # Peter Luschny, Aug 04 2017
MATHEMATICA
Do[s=Count[Table[IntegerQ[Binomial[n, Floor[n/2]]/ Binomial[n, j]], {j, 0, n}], True]; s1=IntegerQ[Binomial[n, n/2]/n]; If[ !Equal[s, 7] && Equal[s1, True], Print[n]], {n, 1, 10000}]
(* Second program: *)
Select[Range@ 5000, Function[n, And[Divisible[Binomial[n, n/2], n], Count[Table[Divisible[Binomial[n, Floor[n/2]], Binomial[n, j]], {j, 0, n}], True] != 7]]] (* Michael De Vlieger, Jul 30 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Mar 17 2003
EXTENSIONS
More terms from Michael De Vlieger, Jul 30 2017
STATUS
approved