Ver el código de «Problema de Flavio Josefo»

[[Archivo:Josephusbust.jpg|miniaturadeimagen|Busto de Flavio Josefo quién vivió en la región de Alejandría, en el primer siglo de la era común. A dicho personaje se le atribuye la historia que dio origen al problema que lleva su nombre.]]

[[File:JosephusProblemDrawing.png|thumb|right|Gráfico de la evolución del problema de Flavio Josefo para 500 personas y un valor de salto de 6. El eje horizontal es el número de persona. El eje vertical (de arriba hacia abajo) es el tiempo (el número de ciclo). Una persona viva es representada como un punto verde, una muerta como un punto negro.<ref name="formulae">{{cite web |url=https://formulae.org/?example=Josephus_problem |title=Josephus problem in Fōrmulæ programming language |first=Laurence |last = R.Ugalde |website=Fōrmulæ |access-date=2021-07-26}}</ref>]]

En [[matemáticas]] y en las [[ciencias de la computación]], el '''problema de Flavio Josefo''' (o '''permutación de Josefo''') es un problema teórico relacionado con un cierto [[Fórmulas de echar a suertes|problema de echar suertes.]]

Hay gente de pie en un [[círculo]] a la espera de ser ejecutada. La cuenta comienza en un punto y dirección específica del círculo. Después de que se haya salteado a un número determinado de personas, la siguiente persona es ejecutada. El procedimiento se repite con las personas restantes, a partir de la siguiente persona, que va en la misma dirección y omitiendo el mismo número de personas, hasta que solo una persona permanece y se libra.

El problema (dado el número de personas, punto de partida, dirección, y el número de personas a saltar) es elegir la posición en el círculo inicial para evitar la ejecución.

== Historia ==
El problema recibe el nombre de [[Flavio Josefo]], historiador judío que vivió en el {{siglo|I||s}}. Según el relato del [[asedio de Yodfat]] de Josefo, él y sus 40 soldados quedaron atrapados en una cueva por los [[soldados romanos]]. Eligieron el suicidio durante la captura, y establecieron un método para suicidarse en serie por sorteo. Josefo afirma que, por suerte o, posiblemente de la mano de Dios, él y otro hombre se mantuvieron hasta el final y se rindieron a los romanos en lugar de matarse a sí mismos. Esta es la historia dada en el libro 3, capítulo 8, parte 7 de [[La guerra de los judíos|La Guerra de los Judíos]] de Josefo (escritura de sí mismo en tercera persona):

{{quotation|Sin embargo, en esta angustia extrema, él no estaba desprovisto de su sagacidad habitual; pero confiando a sí mismo a la providencia de Dios, puso su vida en peligro: "Y ahora," dijo él, "ya que se resuelve entre nosotros que vamos a morir, vamos, encomendemos nuestras mutuas muertes a la determinación por suertes. Aquel a quien el azar caiga primero, que sea degollado por quien reciba el segundo azar, y por lo tanto la fortuna hará su progreso a través de todos nosotros; y ninguno de nosotros será asesinado por su propias manos, porque no sería justo si, cuando el resto haya muerto, alguien se arrepintiese y se salvase". Esta propuesta les pareció muy justa; y cuando él se había impuesto con ellos para determinar este asunto por azar, participó en la rifa también. El que obtuvo la primera suerte puesta en su cuello desnudo por quien obtuvo la siguiente, como supuso que el general iba a morir entre ellos, de inmediato pensó que la muerte, si Josefo moría entre ellos, sería más dulce que la vida; sin embargo, fue con el compañero de su izquierda hasta el final. Digamos que tal cosa ocurrió por casualidad, o por la providencia de Dios. Y como él estaba muy deseoso de no ser condenado por la suerte y de no de empapar su mano derecha con la sangre de sus compatriotas, lo convenció a confiar en él, y en vivir para él como a sí mismo.<ref name="Josephus18">[[Flavio Josefo]], [http://www.gutenberg.org/files/2850/2850-h/2850-h.htm#link32HCH0008 Guerra de los judíos III. 8. 7.] (Traducción por William Whiston).</ref>}}

Los detalles del mecanismo utilizado en esta hazaña son bastante vagos. Según Dowdy y Mays,<ref>{{harvnb|Dowdy|Mays|1989|page=125}}</ref> en 1612 Bachet sugirió que el mecanismo específico fue disponer a los hombres en un círculo y contar de tres en tres para determinar el orden de eliminación.<ref>{{citation|first=C. G.|last=Bachet|year=1612|title=Problemes Plaisants ed Delectables qui se font par les Nombres|page=174}}</ref> Esta historia se ha repetido con frecuencia y los detalles específicos varían considerablemente de una fuente a otra. Por ejemplo, Herstein y Kaplansky (1974) tienen a Josefo y a sus 39 compañeros colocados en círculo con cada séptimo hombre eliminado.<ref>{{citation|first1=I.N.|last1=Herstein|first2=I.|last2=Kaplansky|year=1974|title=Matters Mathematical|publisher=Harper and Row|pages=121–126}}</ref> Una historia del problema se puede encontrar en la carta de S. L. Zabell al editor de la Fibonacci Quarterly.<ref>{{citation|first=S. L.|last=Zabell|year=1976|title=Letter to the editor|journal=Fibonacci Quarterly|volume=14|pages=48 & 51}}</ref>

Josefo tenía un cómplice; el problema era entonces encontrar los lugares de los dos últimos supervivientes (de cuya conspiración aseguraría sus supervivencias). Alega que se colocó junto al otro hombre en las posiciones 31 y 16 respectivamente.<ref name=Ball1896>{{cita libro|last1=Rouse Ball|first1=W. W.|título=Mathematical Recreations and Essays|fecha=1896|editorial=Macmillan|edición=2nd|url=https://www.gutenberg.org/ebooks/26839}}</ref>

== Variantes y generalizaciones ==
Una versión medieval del problema de Josefo involucra a 15 cristianos y a 15 turcos a bordo de un barco en una tormenta que se hundirá a menos que la mitad de los pasajeros sean arrojados por la borda. Los 30 se disponen en un círculo y cada novena persona ha de ser arrojada al mar. ¿Dónde deben los cristianos colocarse para garantizar que solo los turcos se lanzan?<ref>{{citation|first=J.R.|last=Newman|title=The World of Mathematics|volume=4|year=1988|publisher=Tempus|pages=2403–2405}}</ref> En otras versiones los roles de los turcos y los cristianos están intercambiados.

En ''Matemáticas concretas: Una Fundación de Ciencias Informáticas'', Graham, Knuth y Patashnik describen y estudian una variante "estándar":<ref>{{citation|first1=R.L.|last1=Graham|first2=D.E.|last2=Knuth|first3=O.|last3=Patashnik|title=Concrete Mathematics: A Foundation for Computer Science|publisher=Addison Wesley|year=1989|page=8|isbn=978-0-201-14236-5}}</ref> determinan dónde se disponen los últimos supervivientes si hay {{mvar|n}} personas inicialmente y cada segunda persona ({{mvar|k}} menor o igual a 2) es eliminada.

Una generalización del problema es la siguiente. Suponemos que cada {{mvar|m}}-ésima persona será ejecutada de un grupo de tamaño {{mvar|n}}, en el cual la {{mvar|n}}-ésima persona es la superviviente. Si hay una adición de {{mvar|x}} personas al círculo, entonces el superviviente está en la ({{math|''p'' + ''mx''}})-ésima posición si es menor o igual a {{math|(''p'' + ''mx'') − (''n'' + ''x'')}}<ref name=Robinson1960>{{cita publicación|last1=Robinson|first1=W. J.|título=The Josephus Problem|url=https://archive.org/details/sim_mathematical-gazette_1960-02_44_347/page/47|publicación=The Mathematical Gazette|fecha=1960|volumen=44|número=347|páginas=47–52|doi=10.2307/3608532}}<!--|accessdate=1 August 2015--></ref>

== Solución ==
Se resuelve el problema de manera explícita cuando cada segunda persona muere, es decir, <math>k=2</math>. (Para un caso más general donde <math>k\neq 2</math>, puede verse la solución más abajo.) Expresamos la solución recursiva. Sea <math>f(n)</math> la posición de la superviviente cuando hay inicialmente <math>n</math> personas (y <math>k=2</math>).En la primera vuelta alrededor del círculo, todas las personas en posición par mueren. En la segunda vuelta al círculo, la nueva segunda persona muere, después la nueva 4.ª persona, etc .; es como si no existiera primera vez alrededor del círculo.

Si el número inicial de la gente es par, entonces la persona en la posición <math>x</math> durante la segunda vuelta alrededor del círculo estaba originalmente en la posición <math>2x - 1</math> (para cualquier valor de <math>x</math>). Sea <math>n=2j</math> la persona en <math>f(j)</math> que ahora sobrevivirá, estaba originalmente en la posición <math>2f(j) - 1</math>. Esto nos da la recurrencia

:<math>f(2j)=2f(j)-1\;.</math>

Si el número inicial de la gente es impar, entonces pensamos que la persona en la posición primera morirá al final de la primera vuelta alrededor del círculo. Una vez más, durante la segunda vuelta alrededor del círculo, la nueva segunda persona muere, después la nueva 4.ª persona, etc. En este caso, la persona en la posición <math>x</math> estaba originalmente en la posición <math>2x+1</math>. Esto nos da la recurrencia

:<math>f(2j+1)=2f(j)+1\;.</math>

Cuando contabilizamos los valores de <math>n</math> and <math>f(n)</math>, vemos el patrón:

{| class="wikitable"
|<math>n</math> || 1 || 2 || 3 || 4 || 5 || 6 || 7 || 8 || 9 || 10 || 11 || 12 || 13 || 14 || 15 || 16
|-
|<math>f(n)</math> || 1 || 1 || 3 || 1 || 3 || 5 || 7 || 1 || 3 || 5 || 7 || 9 || 11 || 13 || 15 || 1
|}

Esto sugiere que <math>f(n)</math> es una secuencia creciente impar que retorna con <math>f(n)=1</math>, siempre que el índice ''n'' sea una potencia de 2. Por lo tanto, si elegimos ''m'' y ''l'' de manera que <math>n=2^m+l</math> y <math>0\leq l<2^m</math>, entonces <math>f(n)=2 \cdot l+1</math>.
Está claro que los valores de la tabla satisfacen esta ecuación. O podemos pensar que después de que <math>l</math> personas hayan muerto, solo hay <math>2^m</math> personas y nos vamos a la <math>2l+1</math>-ésima persona. Él debe ser el sobreviviente. De tal modo que <math>f(n)=2l+1</math>. A continuación, damos una demostración por inducción.

'''Teorema:''' Si <math>n=2^m+l</math> y <math>0\leq l<2^m</math>, entonces <math>f(n) = 2l+1</math>.

'''Prueba:''' Usamos la inducción sobre <math>n</math>. El caso base <math>n=1</math> es correcto. Consideramos separadamente los casos cuando <math>n</math> es par y cuando <math>n</math> es impar.

Si <math>n</math> es par, entonces se elige <math>l_1</math> y <math>m_1</math> tal que <math>n/2 = 2^{m_1}+l_1</math> y <math>0\leq l_1 < 2^{m_1}</math>. Notar que <math>l_1 = l/2</math>.
Tenemos <math>f(n) = 2f(n/2)-1=2((2l_1)+1) - 1=2l+1</math>, donde la segunda igualdad se deduce de la hipótesis de inducción.

Si <math>n</math> es impar, entonces se elige <math>l_1</math> y <math>m_1</math> tal que <math>(n-1)/2 = 2^{m_1}+l_1</math> y <math>0\leq l_1 < 2^{m_1}</math>. Notar que <math>l_1 = (l-1)/2</math>.
Tenemos <math>f(n) = 2f((n-1)/2)+1=2((2l_1)+1) + 1=2l+1</math>, donde la segunda igualdad se deduce de la hipótesis de inducción. Esto completa la demostración.

La podemos resolver para <math>l</math> para obtener una expresión específica para <math>f(n)</math>:

:<math>f(n) = 2(n-2^{\lfloor \log_2(n) \rfloor})+1</math>

La expresión más elegante recurre a la representación binaria de tamaño <math>n</math>: <math>f(n)</math> se pueden obtener por un desplazamiento cíclico de un bit a la izquierda de <math>n</math> por sí solo. Si representamos <math>n</math> en binario como <math>n=1 b_1 b_2 b_3\dots b_m</math>, entonces la solución está dada por<math>f(n)=b_1 b_2 b_3 \dots b_m 1</math>. La prueba de esto se desprende de la representación de <math>n</math> como <math>2^m+l</math> o de la expresión de arriba para <math>f(n)</math>.

'''Implementación:''' Si n indica el número de personas, la posición de seguridad está dada por la función <math>f(n) = 2l+1</math> ,donde  <math>n=2^m+l</math> y <math>0\leq l<2^m</math>.

Ahora si representamos el número en formato binario, el primer bit indica  <math>2^m</math> y los bits restantes se denotarán por <math>l</math>. Por ejemplo, cuando n = 41, su representación binaria es

'''n    = 1 0   1 0   0 1'''

'''2<sup>m</sup>  = 1 0   0 0   0 0'''

'''l     =       0 1   0 0   1'''<syntaxhighlight lang="java">
    /**
	 * 
	 * el parámetro n el número de personas en el círculo
	 * devuelve la posición de aquellos que sobrevivirán la ejecución
	 *   f(N) = 2L + 1 donde N =2^M + L y 0 <= L < 2^M
	 */
	public int getSafePosition(int n) {
		// find value of L for the equation
		int valueOfL = n - Integer.highestOneBit(n);
		int safePosition = 2 * valueOfL  + 1;
		
		return safePosition;
	}
</syntaxhighlight>

=== El caso general ===
La forma más sencilla de resolver este problema en el caso general usando la [[programación dinámica]] mediante la realización de la primera etapa y luego utilizando la solución del problema restante. Cuando el índice comienza desde 1, la persona en la posición <math>s</math> se desplaza desde donde está la primera persona en la posición <math>((s-1)\bmod n)+1</math>, donde n es el número total de personas. Sea <math>f(n,k)</math> la posición del superviviente, después de que la <math>k</math>-ésima persona haya muerto, nos quedamos con un círculo de <math>n-1</math> personas, y se inicia el siguiente conteo con la persona cuyo número en la situación original fue <math>(k\bmod n)+1</math>. La posición del superviviente en el círculo restante sería <math>f(n-1,k)</math> si empezamos a contar en <math>1</math>; cambiando al hecho de que estamos a partir de <math>(k\bmod n)+1</math>, se obtiene la recurrencia.

: <math>f(n,k)=((f(n-1,k)+k-1) \bmod n)+1,\text{ con }f(1,k)=1\,,</math>

que toma una forma más simple

: <math>g(n,k)=(g(n-1,k)+k) \bmod n,\text{ con }g(1,k)=0</math>

si numeramos las posiciones de <math>0</math> a <math>n-1</math> en su lugar.

Esta aproximación tiene un [[Cota superior asintótica|tiempo de ejecución]] <math>O(n)</math>, pero para una <math>k</math> pequeña y una <math>n</math> grande, hay otra aproximación. El segundo enfoque también utiliza la programación dinámica, pero tiene tiempo de ejecución <math>O(k\log n)</math>. Está basado en que se considera matar a la ''k''-ésima, 2''k''-ésima,..., <math>(\lfloor n/k \rfloor k)</math>-ésima persona una a una, e ir cambiando la numeración.

Este enfoque mejorado toma la forma

:<math>g(n,k) = \begin{cases}
0 & \text{si } n=1\\
(g(n-1,k)+k) \bmod n & \text{si } 1 < n < k\\
\left \lfloor \frac{k((g(n',k)-n \bmod k) \bmod n')}{k-1} \right \rfloor \text{donde } n'= n- \left \lfloor \frac{n}{k} \right \rfloor & \text{si } k \le n\\
\end{cases} </math>

== Referencias ==
{{listaref|2}}
== Enlaces externos ==
* {{cita libro |autor=[[Thomas H. Cormen]], [[Charles E. Leiserson]], [[Ronald L. Rivest]], y [[Clifford Stein]] |título=[[Introducción a los algoritmos]] |edición=Segunda |editorial=[[MIT Press]] and McGraw-Hill |año=2001 |isbn=0-262-03293-7 |capítulo=Capítulo 14: Aumentar estructuras de datos |página=318}}
* {{IMMP |flavius|El juego de Flavio Josefo|clase=recurrence}} (inglés)
* [[Armin Shams-Baragh]] [http://www.cs.man.ac.uk/~shamsbaa/Josephus.pdf Formulando el problema extendido de Josefo] (inglés)
* [http://mathworld.wolfram.com/JosephusProblem.html El problema de Josefo en la enciclopedia MathWorld] (inglés)
* {{cita web|last1=[[Daniel Erman]]|título=El problema de Josefo - Numberphile|url=https://www.youtube.com/watch?v=uCsD3ZGzMgE|editorial=[[Brady Haran]]|fechaacceso=29 de octubre de 2016|ubicación=YouTube|formato=video|fecha=28 de octubre de 2016}} (inglés)
* {{Traducido ref|en|Josephus problem|inglés|oldid=775434834|trad=total}}
{{Traducido ref|en|Josephus problem|oldid=775434834|trad=parcial|fecha= 14 de abril de 2017}}
{{Control de autoridades}}
[[Categoría:Problemas computacionales]]
[[Categoría:Permutaciones]]
[[Categoría:Combinatoria]]
[[Categoría:Problemas matemáticos]]
[[Categoría:Flavio Josefo]]