OFFSET
1,1
COMMENTS
LINKS
PROG
(PARI)
A156552(n) = {my(f = factor(n), p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ From A156552
isA332228(n) = ((n%2)&&!isprimepower(n)&&2==(sigma(n)%8));
isA332229(n) = isA332228(A156552(n));
(PARI)
v156552sigs = readvec("a156552.txt"); \\ Factorization file for A156552 prepared by Hans Havermann, available at https://oeis.org/A156552/a156552.txt
isA156552not_a_primepower(n) = if(n<=2, 0, my(prsig=v156552sigs[n]); length(prsig[1])>1);
A323243(n) = if(n<=2, n-1, my(prsig=v156552sigs[n], ps=prsig[1], es=prsig[2]); prod(i=1, #ps, ((ps[i]^(1+es[i]))-1)/(ps[i]-1)));
isA332229(n) = (!(n%2)&&isA156552not_a_primepower(n)&&(2==(A323243(n)%8)));
k=0; for(n=1, 10000, if(isA332229(n), k++; print1(n, ", ")));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 13 2020
STATUS
approved