A241004 Numbers n such that anti-phi(sigma*(n)) = n, where anti-phi is A066452 and sigma* is the sum of anti-divisors of n (A066417).

90, 137, 162, 581, 714, 773, 3735, 4557, 71028

Offset

1,1

Comments

Like A001229 but using anti-phi(n) (A066452) and sigma*(n) (A066417).

Links

Table of n, a(n) for n=1..9.

Example

90 is in the sequence: Anti-divisors of 90 are 4, 12, 20, 36, 60 and their sum is 132. Anti-phi of 132 is 90.

Maple

isA241004 := proc(n)
    simplify( n = A066452(A066417(n))) ;
end proc:
for n from 1 do
    if isA241004(n) then
        printf("%d\n", n) ;
    end if;
end do: # R. J. Mathar, Aug 07 2014

Prog

(PARI) sad(n) = my(k); if(n>1, k=valuation(n, 2); sigma(2*n+1)+sigma(2*n-1)+sigma(n/2^k)*2^(k+1)-6*n-2, 0);
antidiv(n) = {my(v = []); for (k=2, n-1, if (abs((n % k) - k/2) < 1, v = concat(v, k)); ); v; }
antiphi(n) = {my(vad = antidiv(n)); my(nbad = 0); for (j=1, n-1, isad = 1; for (k=1, #vad, if ((j % vad[k]) == 0, isad = 0; break); ); nbad += isad; ); nbad; }
isok(n) = n == antiphi(sad(n)); \\ Michel Marcus, Feb 25 2016

Crossrefs

Cf. A001229, A066417, A066452.

Sequence in context: A352231 A249553 A125010 * A294310 A048705 A119895

Adjacent sequences: A241001 A241002 A241003 * A241005 A241006 A241007

Keyword

nonn,more

Author

Paolo P. Lava, Aug 07 2014

Extensions

a(9) from Michel Marcus, Feb 25 2016

Status

approved