OFFSET
1,2
COMMENTS
Numbers k such that A020733(2*k) = 1.
LINKS
Robert Israel, Table of n, a(n) for n = 1..250
EXAMPLE
a(4)=6 is in the sequence because binomial(12,6) = 924 = 2^2*3*7*11 has 4 distinct prime factors while binomial(12,0) to binomial(12,5) all have at most 3.
7 is not in the sequence because binomial(14,7) = 3432 = 2^3*3*11*13 and binomial(14,6) = 3003 = 3*7*11*13 both have 4 distinct prime factors.
MAPLE
filter:= proc(n) local t, v, i, m;
m:= 0: t:= 1:
for i from 1 to n-1 do
t:= t * ifactor(2*n-i+1)/ifactor(i);
if type(t, `*`) then v:= nops(t) else v:= 1 fi;
if v > m then m:= v fi;
od;
t:= t*ifactor(n+1)/ifactor(n);
type(t, `*`) and nops(t) > m
end proc:
filter(1):= true:
select(filter, [$1..2500]); # Robert Israel, May 26 2020
MATHEMATICA
Select[Range@ 1001, Max@ Most@ # < Last@ # &@ PrimeNu@ Binomial[2 #, Range[0, #]] &] (* Michael De Vlieger, May 26 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, May 26 2020
STATUS
approved