cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A335145 Numbers that are both unitary and nonunitary Zumkeller numbers.

This page as a plain text file.
%I A335145 #6 Aug 04 2025 07:42:43
%S A335145 150,294,630,726,750,840,1014,1050,1470,1650,1734,1890,1950,2058,2166,
%T A335145 2550,2850,2940,2970,3174,3234,3450,3510,3630,3750,3822,4350,4410,
%U A335145 4650,4998,5046,5070,5082,5250,5550,5586,5670,5766,6150,6450,6762,6930,7050,7098,7260
%N A335145 Numbers that are both unitary and nonunitary Zumkeller numbers.
%e A335145 150 is a term since its unitary divisors, {1, 2, 3, 6, 25, 50, 75, 150} can be partitioned in two disjoint sets of equal sum: 1 + 2 + 3 + 25 + 50 + 75 = 6 + 150, and so are its nonunitary divisors, {5, 10, 15, 30}: 5 + 10 + 15 = 30.
%t A335145 zumQ[n_] := Module[{d = Divisors[n], ud, nd, sumUd, sumNd, x},ud = Select[d, CoprimeQ[#, n/#] &]; nd = Complement[d, ud]; sumUd = Plus @@ ud; sumNd = Plus @@ nd; sumUd >= 2*n && sumNd > 0 && EvenQ[sumUd] && EvenQ[sumNd] && CoefficientList[Product[1 + x^i, {i, ud}], x][[1 + sumUd/2]] > 0 && CoefficientList[Product[1 + x^i, {i, nd}], x][[1 + sumNd/2]] > 0]; Select[Range[10000], zumQ]
%Y A335145 Intersection of A290466 and A335142.
%K A335145 nonn
%O A335145 1,1
%A A335145 _Amiram Eldar_, May 25 2020