OFFSET
1,1
COMMENTS
Since m=n is allowed, perfect numbers (A000396) are terms of this sequence. - Michel Marcus, Jan 02 2015
m: 6, 24, 28, 40, 120, 216, 234, 270, 360, 496, 588, 672, 2016, ..., . - Robert G. Wilson v, Feb 28 2015
Odd members are 98553, 378225, 700245, ..., . - Robert G. Wilson v, Feb 28 2015
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..157
Robert G. Wilson v, Table of n, a(n) and m for n = 1..157
EXAMPLE
For n = 36, m = 24, 36 is the sum of the proper divisors of 24, and rad(36) = rad(24) = 6.
MATHEMATICA
rad[n_] := Times @@ (First@ # & /@ FactorInteger@ n); f[n_] := Block[{sd = DivisorSigma[1, n] - n}, If[ rad[n] == rad[sd], sd, 0]]; k = 1; lst = {}; While[k < 1000001, a = f@ k; If[a > 0, AppendTo[lst, a]]; k++]; Sort@ lst (* Robert G. Wilson v, Feb 28 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Naohiro Nomoto, Dec 15 2014
STATUS
approved
