Если я правильно понял, то эти числа, которые не имеют простых делителей ровно в первой степени. Powerful иными словами. В PARI/GP есть функция ispowerful(m), которая соответствует вашей функции.
Мультипликативность для взаимнопростых ясно видна при принятии

, но может быть и проверена в диапазоне до миллиона по каждому сомножителю
{for( i=2,10^6, for( j=2,10^6,
if(gcd(i,j)==1 && ispowerful(i)*ispowerful(j)!=ispowerful(i*j),
print("wow ",i," ",j)) ))}ничего не выдала
