4 votos

(Promedio) Número de ciclos de longitud m en permutaciones de N con k ciclos

Supongamos que tenemos permutaciones en $[1,2,...,n]$ que tienen exactamente $k$ ciclos (que no se $|s(n,k)|$ donde $s(n,k)$ es el número de Stirling de primera especie).

¿Cuál es el número promedio de ciclos de longitud $m$ para una de estas permutaciones? Vamos a llamar a este número $C(m,n,k)$.

Así, por $m=1$, este es el número total de puntos fijos en todas aquellas permutaciones dividido por |s(n,k)|. Incluso para este caso estoy teniendo problemas para la formulación de una relación de recurrencia.

Sabemos por definición que $\sum_m C(m,n,k)=k$$\sum_m m.C(m,n,k)=n$.

Por ejemplo, para $n=4,k=2$, $|s(4,2)|=11$ y tenemos, $$C(1,4,2) = 8/11, \qquad C(2,4,2) = 6/11, \qquad C(3,4,2) = 8/11$$ y para $n=9, k=3$, $|s(9,3)|=118124$ y tenemos, $$\;\;\;C(1,9,3)=117612/|s(9,3)|=0.9957..., \\\; C(2,9,3)= 63504/|s(9,3)|=0.5376...,\\\; C(3,9,3)=46032/|s(9,3)|=0.3897..., \\\; C(4,9,3)=37800/|s(9,3)|=0.3200...,\\\; C(5,9,3)=33264/|s(9,3)|=0.2816...,\\\; C(6,9,3)=30240/|s(9,3)|=0.2560...,\\\;C(7,9,3)=25920/|s(9,3)|=0.2194...$$

Por el momento, estoy interesada sobre todo en el caso de que $n=k^2$. Desde el anterior experimento parece posiblemente \lim_{k\to\infty} C(1,k^2,k)=1. edit: De Marko la respuesta de que esto no es cierto.

En última instancia, estoy interesada en la distribución acumulativa, decir, $$\lim_{k\to\infty} \sum_{m=0}^{kx} C(m,k^2,k)$$

Cualquier referencia o ayuda sería muy apreciada.

6voto

Marko Riedel Puntos 19255

Para la primera parte de la pregunta con $Q(n, k, m)$ el número de ciclos de tamaño $m$ entre las permutaciones de $[n]$ tener $k$ ciclos obtener la especie

$$\mathfrak{P}_{=k}(\mathfrak{C_{\lt m}}(\mathcal{Z}) + \mathcal{U}\mathfrak{C_{= m}}(\mathcal{Z}) + \mathfrak{C_{\gt m}}(\mathcal{Z})).$$

Esto produce que la bivariante de generación de función

$$G(z, u) = \frac{1}{k!}\a la izquierda(\log\frac{1}{1-z}-\frac{z^m}{m} + u\frac{z^m}{m}\right)^k.$$

La estadística deseada ha de generación de función

$$\left.\frac{\partial}{\partial u} G(z, u)\right|_{u=1} \\ = \left. \frac{1}{(k-1)!}\a la izquierda(\log\frac{1}{1-z}-\frac{z^m}{m} + u\frac{z^m}{m}\right)^{k-1} \frac{z^m}{m} \right|_{u=1} \\ = \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \frac{z^m}{m}.$$

La extracción de los coeficientes obtenemos

$$n! [z^n] \frac{z^m}{m} \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \frac{n!}{m} [z^{n-m}] \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \frac{n!}{m (n-m)!} (n-m)! [z^{n-m}] \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1}.$$

Dividir por ${n\brack k}$ por la expectativa

$$\bbox[5px,border:2px solid #00A000]{ \frac{n!}{m (n-m)!} {n-m\brack k-1} {n\brack k}^{-1}.}$$

Para la prueba de que estos se suma a $k$ que debe sostener por primera principios observamos que necesitamos $n-m\ge k-1$ o $n+1-k\ge m$ y obtener la reclamación

$${n\brack k}^{-1} \sum_{m=1}^{n+1-k} \frac{n!}{m (n-m)!} {n-m\brack k-1} = k.$$

El FEAG de la suma es

$$\sum_{n\ge k} \frac{w^n}{n!} \sum_{m=1}^{n+1-k} \frac{n!}{m (n-m)!} (n-m)! [z^{n-m}] \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \sum_{n\ge k} w^n \sum_{m=1}^{n+1-k} \frac{1}{m} [z^{n}] z^m \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1}.$$

Ahora podemos extender el interior de la suma de los infinitos términos, porque cuando $m\gt n+1-k$ tenemos $m+k-1\gt n$ y no hay ninguna contribución a la $[z^n].$ Tenemos

$$\sum_{n\ge k} w^n \sum_{m\ge 1} \frac{1}{m} [z^{n}] z^m \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \sum_{n\ge k} w^n [z^n] \sum_{m\ge 1} \frac{1}{m} z^m \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \sum_{n\ge k} w^n [z^n] \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k} \\ = \frac{1}{(k-1)!} \left(\log\frac{1}{1-w} \right)^{k}.$$

La extracción de los coeficientes encontramos

$${n\brack k}^{-1} n! [w^n] k \frac{1}{k!} \left(\log\frac{1}{1-w} \right)^{k} = {n\brack k}^{-1} k {n\brack k} = k$$

como se reivindica. Para comprobar el número dos de la reclamación es

$${n\brack k}^{-1} \sum_{m=1}^{n+1-k} \frac{n!}{(n-m)!} {n-m\brack k-1} = n.$$

Re-uso de la computación en el primero de los rendimientos

$$\sum_{n\ge k} w^n [z^n] \sum_{m\ge 1} z^m \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \sum_{n\ge k} w^n [z^n] \frac{z}{1-z} \frac{1}{(k-1)!} \left(\log\frac{1}{1-z} \right)^{k-1} \\ = \frac{w}{1-w} \frac{1}{(k-1)!} \left(\log\frac{1}{1-w} \right)^{k-1} = w \frac{d}{ps} \frac{1}{k!} \left(\log\frac{1}{1-w} \right)^{k}.$$

La extracción de los coeficientes encontramos

$${n\brack k}^{-1} n! [w^n] w \frac{d}{ps} \frac{1}{k!} \left(\log\frac{1}{1-w} \right)^{k} = {n\brack k}^{-1} n {n\brack k} = n$$

de nuevo como se reivindica. La primera comprobación de validez debe mantener porque cuando nos suma las expectativas de que el número de ciclos de todos los tamaños posibles nos debe obtener $k$ ciclos, la cual es constante en el problema. El segundo debe mantener porque si se suma las longitudes veces el número de ciclos de todos los tamaños posibles debemos cubrir todos los de $n,$ también una constante aquí.

Estas fórmulas fueron verificadas mediante el ciclo de índice $Z(S_n)$ de la grupo simétrico con la siguiente secuencia de comandos de Maple que es posible alrededor de $n=45.$ $n=20$ $k=15$ (permutaciones de veinte elementos que tengan quince ciclos) para el recuento total de ciclos de longitudes de uno a seis los valores

$$10995785640, 2640350580, 737990400, 191280600, 39070080, 4651200$$

y por la expectativa encontramos

$$11.28998110, 2.710993931, 0.7577355487, 0.1963983683, \\ 0.04011541140, 0.004775644215$$

y estas hecho se suma a los quince años.

con(planta);


pet_cycleind_symm :=
proc(n)
opción de recordar;

 si n=0 entonces devolver 1; fi;

 ampliar(1/n*añadir(a[l]*pet_cycleind_symm(n-l), l=1..n));
end;



Q :=
proc(n, k)
opción de recordar;
local idx, conjclass, cyc, conde;

 count := tabla([seq(p=0, q=1..n+1-k)]);

 si n=1 entonces
 idx : = [[1]]
otra cosa
 idx := pet_cycleind_symm(n);
fi;

 para conjclass en idx ¿
 si grado(conjclass) = k entonces
 para cyc en indets(conjclass) ¿
 conde[op(1, cyc)] :=
 conde[op(1, cyc)]
 + grado(conjclass, cyc)
 * lcoeff(conjclass);
od;
fi;
od;

 [seq(n!*recuento de[p], p=1..n+1-k)];
end;

QX :=
(n, k), m) - > n!/m/(n-m)! * abs(stirling1(n-m, k-1));

QXS := (n, k) -> [seq(QX(n, k), m), m=1..n+1-k)];

Observación. El anteriormente utilizado la técnica de la aniquilado coeficiente de extractores (ACE), también conocida como la sustitución de la regla formal el poder de la serie. Hay varios ejemplos en este MSE enlace Yo y en este MSE enlace II y también aquí en este MSE enlace III.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X