OFFSET
1,1
MAPLE
N:= 1000: # to get all terms <= N
{seq(seq(p*m^2, m = 1 .. floor(sqrt(N/p))), p = select(isprime, [2*i+1 $ i = 1..floor((N-1)/2)]))};
# if using Maple 11 or previous, uncomment the next line
# sort(convert(%, list));
# Robert Israel, Oct 30 2014
MATHEMATICA
Take[Sort[Flatten[Table[Prime[n]*m^2, {n, 2, 1000}, {m, 1, 100}]]], 100]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Oct 26 2014
STATUS
approved