antiprimes := proc(n) local ap, i, max_divisors, num_divisors; max_divisors := 0; ap := []; for i from 1 while numelems(ap) < n do num_divisors := numelems(NumberTheory:-Divisors(i)); if num_divisors > max_divisors then ap := [op(ap), i]; max_divisors := num_divisors; end if; end do; return ap; end proc: antiprimes(20);