3.- Generación de números aleatorios
3.- Generación de números aleatorios
3.- Generación de números aleatorios
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>Generación</strong> <strong>de</strong> <strong>números</strong> <strong>aleatorios</strong><br />
Las posibilida<strong>de</strong>s <strong>de</strong> elección <strong>de</strong> x0 para que se cumpla la primera condición <strong>de</strong>l teorema son<br />
Teorema <strong>3.</strong><br />
n<br />
( ) ∏ ( 1)<br />
− = ϕ m Pi<br />
P<br />
i=<br />
1<br />
Sea m=p e , a es raíz primitiva módulo m sii<br />
1) si p=2 ⇒ si:<br />
• e=1⇒ a <strong>de</strong>be ser impar.<br />
• e=2⇒ a mod 4 <strong>de</strong>be ser <strong>3.</strong><br />
• e=3⇒ a mod 8 <strong>de</strong>be ser 3, 5 ó 7.<br />
• e≥4⇒ a mod 8 <strong>de</strong>be ser 3 ó 5.<br />
2) si p es impar ⇒ si:<br />
•<br />
•<br />
⎧a<br />
≠ 0 mod p<br />
⇒#<br />
⎨ p 1<br />
⎩a<br />
≠ 1mod<br />
p<br />
e = 1 − / q<br />
<strong>3.</strong>4.- Otros métodos<br />
⎧ #<br />
⎨ 1<br />
⎩a<br />
≠ 1mod<br />
p<br />
e >⇒ p − 2<br />
con<br />
ei<br />
−1<br />
i<br />
<strong>3.</strong>4.1.- Generadores Cogruenciales Cuadráticos<br />
q divisor primo <strong>de</strong><br />
p −1<br />
2<br />
Son <strong>de</strong>l tipo = ( dx + ax + c)<br />
mod m . La longitud <strong>de</strong> periodo máxima que se alcanza<br />
xn+ 1 n n<br />
con ellos es m, igual que en el caso <strong>de</strong>l G.G.L mixto y sin embargo ha <strong>de</strong> realizar más operaciones que<br />
éste.<br />
Teorema 4.<br />
Para obtener longitud máxima <strong>de</strong> periodo (m) en un generador congruencial cuadrático, se ha<br />
<strong>de</strong> cumplir:<br />
1) c y m <strong>de</strong>ben ser primos relativos.<br />
2) d y a-1 han <strong>de</strong> ser múltiplos <strong>de</strong> todos los factores primos impares <strong>de</strong> m. (Si m=2 e ⇒a-<br />
1=d=1).<br />
3) Si:<br />
•<br />
⎧d ≡ ( a −1)<br />
mod 4<br />
m = k4<br />
⇒ ⎨<br />
⎩d<br />
<strong>de</strong>be ser par<br />
59