OFFSET
1,1
COMMENTS
From Michael S. Branicky, Feb 17 2026: (Start)
Also, integers k, not divisible by any of their digits, such that 1 = Sum_{j = digits of k} 1/j.
Terms cannot have any 0's or 1's and must have less than d digits d for d in 2..9. All terms must be < 999999999 since the reciprocals of the digits of any higher number will add to more than 1.
If the divisibility constraint is removed, then we recover A037268. (End)
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..744
FORMULA
EXAMPLE
a(1) = 263 -> 263/2 + 263/6 + 263/3 = 263.
a(9) = 6443 -> 6443/6 + 6443/4 + 6443/4 + 6443/3 = 6443.
MATHEMATICA
Select[Range[10^5], Function[{k, d}, If[MemberQ[d, 0], False, And[NoneTrue[Union[d], Divisible[k, #] &], k == Total[k/d] ] ] ] @@ {#, IntegerDigits[#]} &] (* Michael De Vlieger, Feb 16 2026 *)
PROG
(Python)
from collections import Counter
from fractions import Fraction
def ok(k):
digs = Counter(map(int, str(k)))
if 0 in digs or 1 in digs or any(k%d == 0 for d in digs): return False
return 1 == sum(digs[d]*Fraction(1, d) for d in digs)
print([k for k in range(10**5) if ok(k)]) # Michael S. Branicky, Feb 16 2026
(PARI) isok(k) = my(d=digits(k)); if (vecmin(d), my(v=vector(#d, i, k/d[i])); (#select(x->(type(x)!="t_INT"), v) == #d) && (vecsum(v) == k)); \\ Michel Marcus, Feb 17 2026
CROSSREFS
KEYWORD
nonn,base,fini,full
AUTHOR
Paolo P. Lava, Feb 15 2026
STATUS
approved
