|
%I #13 Apr 07 2023 04:28:36
%S 1554070,1596070,1725430,1859830,1952230,2095030,2242870,2293270,
%T 2553670,2607430,2716630,2772070,3116470,3481030,3607030,3670870,
%U 3800230,3998470,4065670,4410070,4623430,4841830,5065270,5140870,5371030,5527270,5606230,6009430,6597430
%N Nonsquarefree unitary weird numbers that are also weird numbers.
%C The nonsquarefree unitary weird numbers that are not weird numbers are listed in A328562.
%H Amiram Eldar, <a href="/A328563/b328563.txt">Table of n, a(n) for n = 1..15716</a> (terms below 10^10)
%t 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]]];
%t 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]
%t (* after _M. F. Hasler_'s pari code at A006037 *)
%Y Subsequence of A064114.
%Y Intersection of A006037 and A292705.
%Y Cf. A328562.
%K nonn
%O 1,1
%A _Amiram Eldar_, Oct 19 2019
|
