Question

Given k positive integer numbers a₁ < a₂ < a₃ < … < a_k, and all integers greater than a_k, we want to check whether the set A = {a_i : i ∈ [1,k]} ∪ {n : n > a_k, n ∈ ℕ } = {a₁, a₂, a₃, … , a_k, a_k+1, a_k+2, …} is closed under addition. This means that:

∑_{1 ≤ i ≤ k} a_i*b_i ∈ A, for any non-negative integers b_i.

For example, {2,4,6,7,8,….} is closed under addition.

Is there any easy way to do this? Any functions that we could use in Mathematica or Matlab?

Answer 1

If the discontinuous portion less than k of the sets are not large I believe you can approach it directly like this:

a = {2, 4, 6};
Tr /@ Subsets[a, {2}];
TakeWhile[%, # < Last@a &];
Complement[%, a] === {}