OFFSET
1,2
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local F, i, q, r;
F:= sort(ifactors(n)[2], (a, b) -> a[1]<b[1] );
if F[1][1] > 3 then return false fi;
for i from 1 to nops(F)-1 do
q:= nextprime(F[i][1]);
if F[i+1][1] = q then if F[i][2]=F[i+1][2] then return false fi
elif F[i+1][1] > nextprime(q) then return false
fi;
od;
true
end proc:
filter(1):= true:
select(filter, [$1..500]); # Robert Israel, Apr 01 2026
PROG
(PARI) is(n) = { my (v=-1); forprime (p=2, oo, if (n==1, return (1), v==v=valuation(n, p), return (0), n\=p^v)) }
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Jul 23 2022
STATUS
approved
