OFFSET
1,1
COMMENTS
Factorizations: 13*19, 2^3*3^2*5*29, 2^4*3^3*7*29, 2^5*3^3*5^2*7, 2^2*3^3*29*61, 2^5*3^3*5*61, 2^7*3*5*11*19*29*1021
FORMULA
PROG
(PARI) A049060(n)={ local(i, resul, rmax, p) ; if(n==1, return(1) ) ; i=factor(n) ; rmax=matsize(i)[1] ; resul=1 ; for(r=1, rmax, p=0 ; for(j=1, i[r, 2], p += i[r, 1]^j ; ) ; resul *= p-1 ; ) ; return(resul) ; } isA051153(r)={ local(s, t) ; s=A049060(r) ; t=A049060(s) ; if( 2*r == t, return(1), return(0) ) ; } { for(n=1, 30000000, if( isA051153(n), print(n, ", ") ) ; ) ; } \\ R. J. Mathar, Oct 12 2006
CROSSREFS
KEYWORD
nonn,more
AUTHOR
EXTENSIONS
Corrected and extended by R. J. Mathar, Oct 12 2006
STATUS
approved