# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a365073 Showing 1-1 of 1 %I A365073 #20 Dec 13 2024 09:42:16 %S A365073 1,1,3,6,14,26,60,112,244,480,992,1944,4048,7936,16176,32320,65088, %T A365073 129504,261248,520448,1046208,2090240,4186624,8365696,16766464, %U A365073 33503744,67064064,134113280,268347392,536546816,1073575936,2146703360,4294425600,8588476416,17178349568 %N A365073 Number of subsets of {1..n} that can be linearly combined using nonnegative coefficients to obtain n. %H A365073 Andrew Howroyd, Table of n, a(n) for n = 0..100 %H A365073 Steven R. Finch, Monoids of natural numbers, March 17, 2009. %e A365073 The subset {2,3,6} has 7 = 2*2 + 1*3 + 0*6 so is counted under a(7). %e A365073 The a(1) = 1 through a(4) = 14 subsets: %e A365073 {1} {1} {1} {1} %e A365073 {2} {3} {2} %e A365073 {1,2} {1,2} {4} %e A365073 {1,3} {1,2} %e A365073 {2,3} {1,3} %e A365073 {1,2,3} {1,4} %e A365073 {2,3} %e A365073 {2,4} %e A365073 {3,4} %e A365073 {1,2,3} %e A365073 {1,2,4} %e A365073 {1,3,4} %e A365073 {2,3,4} %e A365073 {1,2,3,4} %t A365073 combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]]; %t A365073 Table[Length[Select[Subsets[Range[n]],combs[n,#]!={}&]],{n,0,5}] %o A365073 (PARI) %o A365073 a(n)={ %o A365073 my(comb(k,b)=while(b>>k, b=bitor(b, b>>k); k*=2); b); %o A365073 my(recurse(k,b)= %o A365073 if(bittest(b,0), 2^(n+1-k), %o A365073 if(2*k>n, 2^(n+1-k) - 2^sum(j=k, n, !bittest(b,j)), %o A365073 self()(k+1, b) + self()(k+1, comb(k,b)) ))); %o A365073 recurse(1, 1<