OFFSET
0,2
COMMENTS
Sets of this type may be called "positive combination-free".
Also subsets of {1..n} such that no element can be written as a (strictly) positive linear combination of the others.
LINKS
S. R. Finch, Monoids of natural numbers, March 17, 2009.
FORMULA
a(n) = 2^n - A365043(n).
EXAMPLE
The subset S = {3,5,6,8} has 6 = 2*3 + 0*5 + 0*8 and 8 = 1*3 + 1*5 + 0*6 but neither of these is strictly positive, so S is counted under a(8).
The a(0) = 1 through a(5) = 20 subsets:
{} {} {} {} {} {}
{1} {1} {1} {1} {1}
{2} {2} {2} {2}
{3} {3} {3}
{2,3} {4} {4}
{2,3} {5}
{3,4} {2,3}
{2,3,4} {2,5}
{1,2,3,4} {3,4}
{3,5}
{4,5}
{2,3,4}
{2,4,5}
{3,4,5}
{1,2,3,4}
{1,2,3,5}
{1,2,4,5}
{1,3,4,5}
{2,3,4,5}
{1,2,3,4,5}
MATHEMATICA
combp[n_, y_]:=With[{s=Table[{k, i}, {k, y}, {i, 1, Floor[n/k]}]}, Select[Tuples[s], Total[Times@@@#]==n&]];
Table[Length[Select[Subsets[Range[n]], And@@Table[combp[Last[#], Union[Most[#]]]=={}, {k, Length[#]}]&]], {n, 0, 10}]
PROG
(Python)
from itertools import combinations
from sympy.utilities.iterables import partitions
def A365044(n):
mlist = tuple({tuple(sorted(p.keys())) for p in partitions(m, k=m-1)} for m in range(1, n+1))
return n+1+sum(1 for k in range(2, n+1) for w in combinations(range(1, n+1), k) if w[:-1] not in mlist[w[-1]-1]) # Chai Wah Wu, Nov 20 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 26 2023
EXTENSIONS
a(15)-a(34) from Chai Wah Wu, Nov 20 2023
STATUS
approved