Deje q=1−p.
La aproximación de Stirling afirma n!∼(ne)n√2πn. Esto nos da:
(npn)∼(n/e)n(np/e)np(nq/e)nq√2πn√2πnp√2πnq
El (n/e)n cancela con (n/e)np(n/e)nq, por lo que esto se simplifica a:
(npn)∼1pnpqnq1√2πnpq
Al tomar el logaritmo en base 2, se obtiene:
log2(npn)∼−n(plog2p+qlog2q)−log2(2πnpq)/2
Dividiendo por n, se transforma a:
1nlog2(npn)∼−(plog2p+qlog2q)−log2(2πnpq)/(2n)=H(p)+O(lnn/n)
Y al n tiende a infinito, se obtiene la deseada límite.
EDIT: no Se necesita toda la fuerza de la aproximación de Stirling. n!=(n/e)ng(n) donde nα≤g≤nβ es suficiente.