A328563 Nonsquarefree unitary weird numbers that are also weird numbers.
1554070, 1596070, 1725430, 1859830, 1952230, 2095030, 2242870, 2293270, 2553670, 2607430, 2716630, 2772070, 3116470, 3481030, 3607030, 3670870, 3800230, 3998470, 4065670, 4410070, 4623430, 4841830, 5065270, 5140870, 5371030, 5527270, 5606230, 6009430, 6597430
Offset: 1
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 1..15716 (terms below 10^10)
Programs
-
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 *)
Comments