cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A279785 Number of ways to choose a strict partition of each part of a strict partition of n.

This page as a plain text file.
%I A279785 #23 Sep 18 2019 17:45:17
%S A279785 1,1,1,3,4,7,11,18,28,47,71,108,166,252,382,587,869,1282,1938,2832,
%T A279785 4153,6148,8962,12965,18913,27301,39380,56747,81226,115907,166358,
%U A279785 236000,334647,475517,671806,947552,1335679,1875175,2630584,3687589,5150585,7183548
%N A279785 Number of ways to choose a strict partition of each part of a strict partition of n.
%C A279785 This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = -1, g(n) = -A000009(n). - _Seiichi Manyama_, Nov 14 2018
%H A279785 Seiichi Manyama, <a href="/A279785/b279785.txt">Table of n, a(n) for n = 0..10000</a> (terms 0..5000 from Alois P. Heinz)
%H A279785 Gus Wiseman, <a href="/A063834/a063834.txt">Sequences enumerating triangles of integer partitions</a>
%F A279785 G.f.: Product_{k>0} (1 + A000009(k)*x^k). - _Seiichi Manyama_, Nov 14 2018
%e A279785 The a(6)=11 twice-partitions are:
%e A279785 ((6)), ((5)(1)), ((51)), ((4)(2)), ((42)), ((41)(1)),
%e A279785 ((3)(2)(1)), ((31)(2)), ((32)(1)), ((321)), ((21)(2)(1)).
%p A279785 with(numtheory):
%p A279785 g:= proc(n) option remember; `if`(n=0, 1, add(add(
%p A279785       `if`(d::odd, d, 0), d=divisors(j))*g(n-j), j=1..n)/n)
%p A279785     end:
%p A279785 b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
%p A279785       `if`(n=0, 1, b(n, i-1)+`if`(i>n, 0, g(i)*b(n-i, i-1))))
%p A279785     end:
%p A279785 a:= n-> b(n$2):
%p A279785 seq(a(n), n=0..70);  # _Alois P. Heinz_, Dec 20 2016
%t A279785 nn=20;CoefficientList[Series[Product[(1+PartitionsQ[k]x^k),{k,nn}],{x,0,nn}],x]
%t A279785 (* Second program: *)
%t A279785 g[n_] := g[n] = If[n==0, 1, Sum[Sum[If[OddQ[d], d, 0], {d, Divisors[j]}]* g[n - j], {j, 1, n}]/n]; b[n_, i_] := b[n, i] = If[n > i*(i + 1)/2, 0, If[n==0, 1, b[n, i-1] + If[i>n, 0, g[i]*b[n-i, i-1]]]]; a[n_] := b[n, n]; Table[a[n], {n, 0, 70}] (* _Jean-François Alcover_, Feb 07 2017, after _Alois P. Heinz_ *)
%Y A279785 Cf. A000009, A063834, A270995, A279375, A327553, A327605.
%K A279785 nonn
%O A279785 0,4
%A A279785 _Gus Wiseman_, Dec 18 2016