OFFSET
1,2
PROG
(PARI) lista(nn) = {vec = vector(nn, i, i); powpr = select((i->((omega(i)==1) || (i==1))), vec); for (i = 1, #powpr, if (issquare(powpr[i]), print1(i, ", ")); ); } \\ Michel Marcus, Oct 02 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved