A037031 Number of combinations of n objects taken pi(n) at a time.

1, 1, 2, 3, 6, 10, 20, 35, 70, 126, 210, 462, 792, 1716, 3003, 5005, 8008, 19448, 31824, 75582, 125970, 203490, 319770, 817190, 1307504, 2042975, 3124550, 4686825, 6906900, 20030010, 30045015, 84672315, 129024480, 193536720, 286097760, 417225900, 600805296

Offset

0,3

Comments

Number of subsets of [n] in which exactly half of the elements are primes: a(6) = 20: {}, {1,2}, {1,3}, {1,5}, {2,4}, {2,6}, {3,4}, {3,6}, {4,5}, {5,6}, {1,2,3,4}, {1,2,3,6}, {1,2,4,5}, {1,2,5,6}, {1,3,4,5}, {1,3,5,6}, {2,3,4,6}, {2,4,5,6}, {3,4,5,6}, {1,2,3,4,5,6}. - Alois P. Heinz, Oct 11 2022

Links

Alois P. Heinz, Table of n, a(n) for n = 0..2000 (terms n = 1..1000 from Vincenzo Librandi)

Formula

a(n) = C(n,pi(n)).

Example

C(16,6) = 8008.

Maple

with(numtheory); A037031:=n->binomial(n, pi(n)); seq(A037031(k), k=0..100); # Wesley Ivan Hurt, Oct 19 2013

Mathematica

Table[Binomial[n, PrimePi[n]], {n, 40}] (* Harvey P. Dale, Mar 22 2012 *)

Prog

(Python)
from math import comb
from sympy import primepi
def A037031(n): return comb(n, primepi(n)) # Chai Wah Wu, Oct 14 2022

Crossrefs

Cf. A000040, A000720.

Sequence in context: A061551 A026034 A178381 * A336675 A336678 A056202

Adjacent sequences: A037028 A037029 A037030 * A037032 A037033 A037034

Keyword

nonn

Author

G. L. Honaker, Jr.

Extensions

a(0)=1 prepended by Alois P. Heinz, Oct 11 2022

Status

approved