Fast lattice reduction for -linear pseudorandom number generators

Authors:
Shin Harase, Makoto Matsumoto and Mutsuo Saito

Journal:
Math. Comp. **80** (2011), 395-407

MSC (2010):
Primary 11K45; Secondary 65C10

Published electronically:
June 18, 2010

MathSciNet review:
2728986

Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: Sequences generated by an -linear recursion have wide applications, in particular, pseudorandom number generation. The dimension of equidistribution with -bit accuracy is a most important criterion for the uniformity of the generated sequence. The fastest known method for computing these dimensions is proposed by Couture and L'Ecuyer, based on Lenstra's lattice basis reduction and the dual lattice to the lattice of vector-valued generating functions (with components in the formal power series ) associated to the output -vector sequence. In this paper we propose a similar but faster algorithm, where (1) the state space is used to represent vectors with components in the formal power series, (2) the dual lattice is not necessary, and (3) Lenstra reduction is replaced with a simpler basis reduction. The computational complexity of our method is smaller than for the Couture-L'Ecuyer method. Experiments show that our method improves the speed by a factor of 10 for Mersenne Twister MT19937 and for WELL generators with state sizes of 19937 bits and 44497 bits.

**1.**Raymond Couture and Pierre L’Ecuyer,*Lattice computations for random numbers*, Math. Comp.**69**(2000), no. 230, 757–765. MR**1651748**, 10.1090/S0025-5718-99-01112-6**2.**Raymond Couture, Pierre L’Ecuyer, and Shu Tezuka,*On the distribution of 𝑘-dimensional vectors for simple and combined Tausworthe sequences*, Math. Comp.**60**(1993), no. 202, 749–761, S11–S16. MR**1176708**, 10.1090/S0025-5718-1993-1176708-4**3.**M. Fushimi and S. Tezuka,*The k-distribution of generalized feedback shift register pseudorandom numbers*, Commun. ACM**26**(1983), no. 7, 516-523.**4.**Shin Harase,*Maximally equidistributed pseudorandom number generators via linear output transformations*, Math. Comput. Simulation**79**(2009), no. 5, 1512–1519. MR**2488100**, 10.1016/j.matcom.2008.06.006**5.**Pierre L’Ecuyer and Raymond Couture,*An implementation of the lattice and spectral tests for multiple recursive linear random number generators*, INFORMS J. Comput.**9**(1997), no. 2, 206–217. MR**1477315**, 10.1287/ijoc.9.2.206**6.**P. L'Ecuyer and F. Panneton,*-linear random number generators*, Advancing the Frontiers of Simulation: A Festschrift in Honor of George Samuel Fishman (C. Alexopoulos, D. Goldsman, and J. R. Wilson, eds.), Springer-Verlag, 2009, pp. 169-193.**7.**A. K. Lenstra,*Factoring multivariate polynomials over finite fields*, J. Comput. System Sci.**30**(1985), no. 2, 235–248. MR**801825**, 10.1016/0022-0000(85)90016-9**8.**Kurt Mahler,*An analogue to Minkowski’s geometry of numbers in a field of series*, Ann. of Math. (2)**42**(1941), 488–522. MR**0004272****9.**Kurt Mahler,*On a theorem in the geometry of numbers in a space of Laurent series*, J. Number Theory**17**(1983), no. 3, 403–416. MR**724538**, 10.1016/0022-314X(83)90057-4**10.**James L. Massey,*Shift-register synthesis and 𝐵𝐶𝐻 decoding*, IEEE Trans. Information Theory**IT-15**(1969), 122–127. MR**0242556****11.**M. Matsumoto and T. Nishimura,*Mersenne twister: a -dimensionally equidistributed uniform pseudo-random number generator*, ACM Trans. Model. Comput. Simul.**8**(1998), no. 1, 3-30.**12.**François Panneton, Pierre L’ecuyer, and Makoto Matsumoto,*Improved long-period generators based on linear recurrences modulo 2*, ACM Trans. Math. Software**32**(2006), no. 1, 1–16. MR**2272349**, 10.1145/1132973.1132974**13.**Sachar Paulus,*Lattice basis reduction in function fields*, Algorithmic number theory (Portland, OR, 1998) Lecture Notes in Comput. Sci., vol. 1423, Springer, Berlin, 1998, pp. 567–575. MR**1726102**, 10.1007/BFb0054893**14.**Wolfgang M. Schmidt,*Construction and estimation of bases in function fields*, J. Number Theory**39**(1991), no. 2, 181–224. MR**1129568**, 10.1016/0022-314X(91)90044-C**15.**Shu Tezuka,*The 𝑘-dimensional distribution of combined GFSR sequences*, Math. Comp.**62**(1994), no. 206, 809–817. MR**1223233**, 10.1090/S0025-5718-1994-1223233-9**16.**Li-Ping Wang and Harald Niederreiter,*Successive minima profile, lattice profile, and joint linear complexity profile of pseudorandom multisequences*, J. Complexity**24**(2008), no. 2, 144–153. MR**2400313**, 10.1016/j.jco.2007.07.001**17.**Liping Wang and Yuefei Zhu,*𝐹[𝑥]-lattice basis reduction algorithm and multisequence synthesis*, Sci. China Ser. F**44**(2001), no. 5, 321–328. MR**1895107****18.**Li-Ping Wang, Yue-Fei Zhu, and Ding-Yi Pei,*On the lattice basis reduction multisequence synthesis algorithm*, IEEE Trans. Inform. Theory**50**(2004), no. 11, 2905–2910. MR**2097012**, 10.1109/TIT.2004.836670

Retrieve articles in *Mathematics of Computation*
with MSC (2010):
11K45,
65C10

Retrieve articles in all journals with MSC (2010): 11K45, 65C10

Additional Information

**Shin Harase**

Affiliation:
Graduate School of Mathematical Sciences, The University of Tokyo, Tokyo, Japan

Email:
sharase@orange.ocn.ne.jp

**Makoto Matsumoto**

Affiliation:
Graduate School of Mathematical Sciences, The University of Tokyo, Tokyo, Japan

Email:
matumoto@ms.u-tokyo.ac.jp

**Mutsuo Saito**

Affiliation:
Department of Mathematics, Hiroshima University, Hiroshima, Japan

Email:
saito@math.sci.hiroshima-u.ac.jp

DOI:
https://doi.org/10.1090/S0025-5718-2010-02391-9

Received by editor(s):
July 15, 2009

Received by editor(s) in revised form:
October 18, 2009

Published electronically:
June 18, 2010

Additional Notes:
The first author was partially supported by Grant-in-Aid for JSPS Fellows 21$·$4427

The second author was partially supported by JSPS Grant-in-Aid for Challenging Exploratory Research 21654017, Scientific Research (A) 19204002 and JSPS Core-to-Core Program 18005.

The third author was partially supported by JSPS Grant-in-Aid for Challenging Exploratory Research 21654004

Article copyright:
© Copyright 2010
American Mathematical Society