OFFSET
1,1
COMMENTS
The bi-unitary version of A171641.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bdiv[n_] := Select[Divisors[n], Last@Intersection[f@#, f[n/#]] == 1 &]; fun[p_, e_] := If[OddQ[e], (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1)-p^(e/2)]; bsigma[n_] := If[n==1, 1, Times @@ (fun @@@ FactorInteger[n])]; seq={}; Do[s=bsigma[n]; If[OddQ[s] || s<=2n, Continue[]]; div = bdiv[n]; If[Coefficient[Times @@ (1 + x^div) // Expand, x, s/2] == 0, AppendTo[seq, n]], {n, 1, 10000}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jan 11 2019
STATUS
approved