# 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<