OFFSET
1,1
COMMENTS
Equivalently, let D_{p}(k) = {1 < d < L : d | k} where p is the smallest prime divisor of k and L = k/p, then k is a term if there exists a subset S of D_{p}(k) such that L = Sum(S).
We will refer to these numbers as 'p-composable'. Each p-composable number belongs to a unique class determined by its smallest prime factor p. For example, 135 is 3-composable, while 2275 is 5-composable. The smallest p-composable numbers appear in A392492. An integer k is a term of this sequence if and only if k is p-composable for some p. Composite numbers that are not p-composable are listed in A392500.
By definition, the index p is the smallest prime divisor of k. Consequently, all 2-composable numbers are even, while all p-composable numbers where p > 2 are odd. No prime number can be p-composable because if n is prime, then L = n/n = 1, making D_{p}(n) empty, and the empty sum cannot equal 1. Furthermore, if the predecessor of an odd prime is composable, it must be a 2-composable number.
LINKS
Peter Luschny, Composable, rigid, and sparse numbers. A Python notebook.
EXAMPLE
The first few p-classes of composable numbers start:
p = 2 | 12, 18, 24, 30, 36, 40, 42, 48, ...
p = 3 | 135, 225, 315, 345, 405, 495, ...
p = 5 | 1645, 1925, 2275, 3185, 4025, ...
p = 7 | 14651, 17017, 19019, 20111, ...
p = 11 | 8041, 23881, 46189, 62491, ...
PROG
(Python) # See also links.
# The function 'subset_sum_exists' is defined in A392652.
from sympy import divisors
def is_composable(n):
if n < 4: return 0
divs = divisors(n)[1:-1]
if not divs: return 0
if len(divs) == 2: return 0
L = divs.pop()
if sum(divs) < L: return 0
return divs[0] if subset_sum_exists(divs, L) else 0
print([n for n in range(253) if is_composable(n)])
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Jan 14 2026
STATUS
approved
