




OFFSET

1,1


COMMENTS

Members of a pair (m,k) such that m = sum of aliquot parts of k minus the sum of the prime distinct divisors of k and k = sum of aliquot parts of m minus the sum of the prime distinct divisors of m.
We introduce in a single sequence the “pseudoamicable numbers” and the “pseudoperfect numbers” in the case m = k. In this sequence 42 and 1316 are pseudoperfect numbers: for example, the sum of the aliquot parts of 42 is 1 + 2 + 3 + 6 + 7 + 14 + 21 = 54 and 54  (2 + 3 + 7) = 54  12 = 42.
The pseudoamicable numbers are (36, 50), (12540240, 29559057). (see the example).
Is the sequence finite?
a(8) > 5*10^9.  Donovan Johnson, Apr 25 2013
a(8) <= 37778715690312487141376. It is easy to verify that if p = (2^k1) and q = 4^k2*2^k1 are two primes, then n = 2^(k1)*p*q is in the sequence, because A206773(n) = n. This happens for k = 2, 3, 7, and 19, which give the terms 42, 1316, 131080256, and 3777871569031248714137, respectively. Up to 3*10^12 there are no other n for which A206773(n) = n.  Giovanni Resta, May 03 2016
a(8) <= 72872313094554244192.  Giovanni Resta, Jan 28 2020


LINKS

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


EXAMPLE

The divisors of m = 50 are {1, 2, 5, 10, 25, 50} and (1 + 2 + 5 + 10 + 25)  (2 + 5) = 43  7 = 36;
The divisors of k = 36 are {1, 2, 3, 4, 6, 9, 12, 18, 36} and (1+2+3+4+6+9+12+18)  (2 + 3) = 55  5 = 50.


MAPLE

with(numtheory):for n from 1 to 10^7 do:x:=factorset(n):n1:=nops(x):s:=sum('x[i] ', 'i'=1..n1):s1:=sigma(n)ns: y:=factorset(s1):n2:=nops(y): ss:=sum('y[i] ', 'i'=1..n2):s2:=sigma(s1)s1ss:if s2=n then printf(`%d, `, n):else fi:od:


CROSSREFS

Cf. A023890, A001065, A000203, A000396, A063990, A206773.
Sequence in context: A044862 A162526 A078299 * A239802 A266242 A274241
Adjacent sequences: A225025 A225026 A225027 * A225029 A225030 A225031


KEYWORD

nonn,hard


AUTHOR

Michel Lagneau, Apr 24 2013


EXTENSIONS

a(7) from Donovan Johnson, Apr 25 2013


STATUS

approved



