|
%I #29 Dec 09 2016 10:14:42
%S 1,3,7,14,25,45,72,117,180,275,403,596,846,1206,1681,2335,3183,4342,
%T 5820,7799,10321,13622,17798,23221,30009,38706,49567,63316,80366,
%U 101805,128211,161134,201537,251495,312508,387535,478674,590072,724920,888795,1086324
%N Total sum of the odd-indexed parts of all partitions of n.
%C For more information see A206563.
%H Alois P. Heinz, <a href="/A207381/b207381.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = A066186(n) - A207382(n) = A066897(n) + A207382(n).
%e For n = 5, write the partitions of 5 and below write the sums of their odd-indexed parts:
%e . 5
%e . 3+2
%e . 4+1
%e . 2+2+1
%e . 3+1+1
%e . 2+1+1+1
%e . 1+1+1+1+1
%e . ------------
%e . 20 + 4 + 1 = 25
%e The total sum of the odd-indexed parts is 25 so a(5) = 25.
%p b:= proc(n, i) option remember; local g, h;
%p if n=0 then [1, 0$2]
%p elif i<1 then [0$3]
%p else g:= b(n, i-1); h:= `if`(i>n, [0$3], b(n-i, i));
%p [g[1]+h[1], g[2]+h[3], g[3]+h[2]+i*h[1]]
%p fi
%p end:
%p a:= n-> b(n,n)[3]:
%p seq(a(n), n=1..50); # _Alois P. Heinz_, Mar 12 2012
%t b[n_, i_] := b[n, i] = Module[{g, h}, If[n == 0 , {1, 0, 0}, If[i < 1, {0, 0, 0}, g = b[n, i - 1]; h = If[i > n, {0, 0, 0}, b[n - i, i]]; {g[[1]] + h[[1]], g[[2]] + h[[3]], g[[3]] + h[[2]] + i*h[[1]]}]]]; a[n_] := b[n, n][[3]]; Table [a[n], {n, 1, 50}] (* _Jean-François Alcover_, Dec 09 2016 after _Alois P. Heinz_ *)
%Y Cf. A066186, A066897, A066898, A181187, A194714, A206283, A206563, A207031, A207032, A207382.
%K nonn
%O 1,2
%A _Omar E. Pol_, Feb 17 2012
%E More terms from _Alois P. Heinz_, Mar 12 2012
|
