%I #20 Mar 19 2019 12:54:50
%S 1,1,3,6,15,28,66,120,253,465,903,1596,3003,5151,9180,15576,26796,
%T 44253,74305,120295,196878,314028,502503,788140,1241100,1917861,
%U 2968266,4531555,6913621,10421895,15705210,23409903,34857075,51445296,75774205,110759286
%N a(n) = A000217(A000041(n)).
%C a(n) is the number of partitions of 2n that are sum-symmetric. That is, a(n) is the number of partitions of 2n that can be divided into two subsequences (no central summand) that each total to n. Example: Of the 11 partitions of 6, there are 6 that are sum-symmetric (partition subsequences bracketed [] and listed in descending order for clarity:) [3][3], [3][2,1], [3][1,1,1], [2,1][2,1], [2,1][1,1,1], [1,1,1][1,1,1]. As this example suggests, a(n) = p(n)*(p(n)+1)/2. - _Gregory L. Simay_, Oct 26 2015
%H Alois P. Heinz, <a href="/A086737/b086737.txt">Table of n, a(n) for n = 0..10000</a> (terms n = 1..1000 from Robert Israel)
%p f:= proc(n) local p;
%p p:= combinat:-numbpart(n);
%p p*(p+1)/2
%p end proc:
%p map(f, [$1..100]); # _Robert Israel_, Oct 26 2015
%t pp = Array[PartitionsP, 40, 0]; pp (pp + 1)/2 (* _Jean-François Alcover_, Mar 19 2019 *)
%o (PARI) a(n) = apply(x->x*(x+1)/2, numbpart(n)); \\ _Michel Marcus_, Oct 26 2015
%Y Cf. A000041, A000217, A108796.
%K nonn
%O 0,3
%A _Jon Perry_, Jul 29 2003
%E a(0)=1 prepended by _Alois P. Heinz_, Mar 25 2017
|