login
a(n) = binomial(bigOmega(n),omega(n)), where omega = A001221 and bigOmega = A001222.
15

%I #19 Jul 10 2022 09:43:46

%S 1,1,1,2,1,1,1,3,2,1,1,3,1,1,1,4,1,3,1,3,1,1,1,6,2,1,3,3,1,1,1,5,1,1,

%T 1,6,1,1,1,6,1,1,1,3,3,1,1,10,2,3,1,3,1,6,1,6,1,1,1,4,1,1,3,6,1,1,1,3,

%U 1,1,1,10,1,1,3,3,1,1,1,10,4,1,1,4,1,1,1,6,1,4,1,3,1,1,1,15,1,3,3,6

%N a(n) = binomial(bigOmega(n),omega(n)), where omega = A001221 and bigOmega = A001222.

%H Reinhard Zumkeller, <a href="/A181591/b181591.txt">Table of n, a(n) for n = 1..10000</a>

%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>

%F a(n) = A007318(A001222(n),A001221(n)).

%F a(n) = A088529(n) = A136565(n) for 2 <= n < 24.

%F a(n) = A136277(A156552(n)). - _Antti Karttunen_, May 29 2017

%t a[n_] := Binomial[PrimeOmega[n], PrimeNu[n]]; Table[a[n], {n, 1, 100}] (* _Jean-François Alcover_, Jun 29 2013 *)

%o (PARI) a(n) = binomial(bigomega(n), omega(n)); \\ _Michel Marcus_, Jul 10 2022

%Y Cf. A001221, A001222, A007318, A136277.

%Y Cf. A088529, A136565 (equal for n = 2..23).

%K nonn

%O 1,4

%A _Reinhard Zumkeller_, Nov 01 2010