OFFSET
0,1
COMMENTS
If a(n) > -1 then a(n)*3^n is the least amicable number k such that A007949(k) = n.
LINKS
EXAMPLE
a(1) = 4865, because 4865 is an odd number and 4865 % 3 > 0 and 4865 * 3 = 14595 is an amicable number, and no lesser number has this property.
PROG
(PARI)
sigmap(k)=if(k, sigma(k)-k, 0)
cycle(k, TT=2)=my(x=k, T=1); while(x>0&&T<=TT, x=sigmap(x); if(x==k, return(T)); T++)
a(n, TT=2)=my(p3n=3^n); forstep(m=1, +oo, 2, if(m%3&&cycle(p3n*m, TT)==2, return(m)))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jean-Marc Rebert, Nov 21 2022
STATUS
approved