OFFSET
1,1
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000
EXAMPLE
22625 is a term since 22625, 22626 and 22627 each have an equal number of unitary and nonunitary divisors. 22625 has 4 unitary divisors (1, 125, 181 and 22625) and 4 nonunitary divisors (5, 25, 905 and 4525), 22626 has 8 unitary divisors and 8 nonunitary divisors, and 22627 has 4 unitary divisors and 4 nonunitary divisors.
MATHEMATICA
q[n_] := DivisorSigma[0, n] == 2^(PrimeNu[n] + 1); v = q /@ Range[3]; seq = {}; Do[v = Append[Drop[v, 1], q[k]]; If[And @@ v, AppendTo[seq, k - 2]], {k, 4, 2 * 10^5}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov and Amiram Eldar, Jun 06 2020
STATUS
approved