OFFSET
1,1
COMMENTS
The nonsquarefree unitary weird numbers that are not weird numbers are listed in A328562.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..15716 (terms below 10^10)
MATHEMATICA
weirdQ[n_, d_, s1_, m1_] := weirdQ[n, d, s1, m1] = Module[{s = s1, m = m1}, If[m == 0, False, While[d[[m]] > n, s -= d[[m]]; m--]; d[[m]] < n && If[s > n, weirdQ[n - d[[m]], d, s - d[[m]], m - 1] && weirdQ[n, d, s - d[[m]], m - 1], s < n && m < Length[d] - 1]]];
aQ[n_] := ! SquareFreeQ[n] && Module[{d = Divisors[n]}, s = Total@d - n; m = Length[d] - 1; Total@Select[d, GCD[#, n/#] == 1 &] > 2 n && weirdQ[n, d, s, m]]; Select[Range[10^7], aQ]
(* after M. F. Hasler's pari code at A006037 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 19 2019
STATUS
approved