OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..2148 (terms below 2*10^10)
EXAMPLE
14 is in the sequence since bsigma(14) = bsigma(15) = 24.
MATHEMATICA
f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bsigma[m_] :=
DivisorSum[m, # &, Last@Intersection[f@#, f[m/#]] == 1 &]; a = {}; b1 = 0; Do[b2 = bsigma[k]; If[b1 == b2, a = AppendTo[a, k - 1]]; b1 = b2, {k, 1, 10^6}]; a (* after Michael De Vlieger at A188999 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 01 2017
STATUS
approved