|
%I #19 Oct 09 2019 04:29:21
%S 0,1,0,1,0,1,1,1,1,1,3,1,2,1,3,3,3,3,1,4,4,3,2,2,4,3,5,3,2,4,5,4,5,6,
%T 1,4,2,5,4,7,4,4,3,3,6,14,3,4,10,6,3,6,4,4,4,8,7,6,8,7,10,5,11,8,5,11,
%U 4,7,7,5,8,12,5,6,9,8,11,8,5,8,9,8,10,8
%N Number of subsets of {1..n} whose sum is equal to the product of their complement.
%C Essentially the same as A178830. - _R. J. Mathar_, Jul 12 2019
%H Giovanni Resta, <a href="/A326441/b326441.txt">Table of n, a(n) for n = 0..500</a>
%e The initial terms count the following subsets:
%e 1: {1}
%e 3: {1,2}
%e 5: {3,5}
%e 6: {3,4,5}
%e 7: {2,4,5,7}
%e 8: {2,4,5,6,7}
%e 9: {2,3,5,6,7,9}
%e 10: {4,5,6,8,9,10}
%e 10: {2,3,5,6,7,8,9}
%e 10: {1,2,3,4,5,8,9,10}
%e Also the number of subsets of {1..n} whose product is equal to the sum of their complement. For example, the initial terms count the following subsets:
%e 1: {}
%e 3: {3}
%e 5: {1,2,4}
%e 6: {1,2,6}
%e 7: {1,3,6}
%e 8: {1,3,8}
%e 9: {1,4,8}
%e 10: {6,7}
%e 10: {1,4,10}
%e 10: {1,2,3,7}
%p b:= proc(n, s, p)
%p `if`(s=p, 1, `if`(n<1, 0, b(n-1, s, p)+
%p `if`(s-n<p*n, 0, b(n-1, s-n, p*n))))
%p end:
%p a:= n-> b(n, n*(n+1)/2, 1):
%p seq(a(n), n=0..100); # _Alois P. Heinz_, Jul 12 2019
%t Table[Length[Select[Subsets[Range[n]],Plus@@#==Times@@Complement[Range[n],#]&]],{n,0,10}]
%Y Cf. A028422, A053632, A059529, A063865, A178830, A301987, A325044, A325538, A326172, A326173, A326174, A326175, A326179, A326180.
%K nonn
%O 0,11
%A _Gus Wiseman_, Jul 07 2019
%E a(21)-a(83) from _Giovanni Resta_, Jul 08 2019
|
