OFFSET
1,1
LINKS
PROG
(PARI)
is_A228058(n) = if(!(n%2)||(omega(n)<2), 0, my(f=factor(n), y=0); for(i=1, #f~, if(1==(f[i, 2]%4), if((1==y)||(1!=(f[i, 1]%4)), return(0), y=1), if(f[i, 2]%2, return(0)))); (y));
is_A388275(k) = if(1==k, 1, my(f=factor(k), x=f[#f~, 1]^f[#f~, 2], s = sigma(f)); if(s%x, 0, my(f=factor(s), y=f[#f~, 1]^f[#f~, 2]); (x==y)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Sep 20 2025
STATUS
approved
