OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
56347 is a term since 56347 = 29^2 * 67, 56347 + 1 = 56348 = 2^2 * 14087, 56347 + 2 = 56349 = 3^3 * 2087 and 56347 + 3 = 56350 = 2 * 5^2 * 7^2 * 23 all have the same number of unitary and nonunitary prime divisors.
MATHEMATICA
q[n_] := n == 1 || Count[(e = FactorInteger[n][[;; , 2]]), 1] == Length[e]/2; v = q /@ Range[4]; seq = {}; Do[v = Append[Drop[v, 1], q[k]]; If[And @@ v, AppendTo[seq, k - 3]], {k, 5, 1.3*10^6}]; seq
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Sep 30 2021
STATUS
approved