A035664 Number of partitions of n into parts 7k+2 and 7k+5 with at least one part of each type.
0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 3, 1, 3, 2, 3, 4, 3, 7, 4, 8, 6, 8, 10, 9, 15, 11, 17, 15, 19, 22, 21, 29, 25, 35, 33, 39, 44, 44, 57, 52, 67, 65, 76, 84, 86, 103, 101, 122, 124, 139, 153, 158, 185, 185, 216, 222, 247, 268, 282, 317, 327, 369, 387, 422, 458, 482
Offset: 1
Keywords
Links
- Alois P. Heinz, Table of n, a(n) for n = 1..1000 (first 125 terms from Robert Price)
Programs
-
Maple
b:= proc(n, i, t, s) option remember; `if`(n=0, t*s, `if`(i<1, 0, b(n, i-1, t, s)+(h-> `if`(h in {2, 5}, add(b(n-i*j, i-1, `if`(h=2, 1, t), `if`(h=5, 1, s)), j=1..n/i), 0))(irem(i, 7)))) end: a:= n-> b(n$2, 0$2): seq(a(n), n=1..75); # Alois P. Heinz, Aug 14 2020
-
Mathematica
nmax = 69; s1 = Range[0, nmax/7]*7 + 2; s2 = Range[0, nmax/7]*7 + 5; Table[Count[IntegerPartitions[n, All, s1~Join~s2], x_ /; ContainsAny[x, s1 ] && ContainsAny[x, s2 ]], {n, 1, nmax}] (* Robert Price, Aug 14 2020 *) nmax = 69; l = Rest@CoefficientList[Series[(-1 + 1/Product[(1 - x^(7 k + 2)), {k, 0, nmax}])*(-1 + 1/Product[(1 - x^(7 k + 5)), {k, 0, nmax}]), {x, 0, nmax}], x] (* Robert Price, Aug 16 2020 *)
Formula
G.f.: (-1 + 1/Product_{k>=0} (1 - x^(7 k + 2)))*(-1 + 1/Product_{k>=0} (1 - x^(7 k + 5))). - Robert Price, Aug 16 2020