reviewed
approved
proposed
editing
(PARI) lista(nn)=vec = vector) = my(nn, i, i); powpr = select((i->((omega(i)==1) || (i==1))), [1))), vec); ..nn])); for (i = 1, #powpr, if (isprime(powpr[i]), print1(i, ", ")); ); \\, ", ")); ); \\ _Michel Marcus_, Jun 03 2021
a(n) = A027883(n) + 1;
a(n) = A025528(A000040(n)) + 1;
Cf. A025528, A065515, A000040.