randlc2s.c
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:3k
源码类别:

数学计算

开发平台:

Unix_Linux

  1. /* gmp_randinit_lc_2exp_size -- initialize a random state with a linear
  2.    congruential generator of a requested size.
  3. Copyright 1999, 2000, 2001 Free Software Foundation, Inc.
  4. This file is part of the GNU MP Library.
  5. The GNU MP Library is free software; you can redistribute it and/or modify
  6. it under the terms of the GNU Lesser General Public License as published by
  7. the Free Software Foundation; either version 3 of the License, or (at your
  8. option) any later version.
  9. The GNU MP Library is distributed in the hope that it will be useful, but
  10. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  11. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  12. License for more details.
  13. You should have received a copy of the GNU Lesser General Public License
  14. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  15. #include <stdio.h> /* for NULL */
  16. #include "gmp.h"
  17. #include "gmp-impl.h"
  18. /* Array of LC-schemes, ordered in increasing order of the first
  19.    member (the 'm2exp' value).  The end of the array is indicated with
  20.    an entry containing all zeros.  */
  21. /* All multipliers are in the range 0.01*m and 0.99*m, and are
  22. congruent to 5 (mod 8).
  23. They all pass the spectral test with Vt >= 2^(30/t) and merit >= 1.
  24. (Up to and including 196 bits, merit is >= 3.)  */
  25. struct __gmp_rand_lc_scheme_struct
  26. {
  27.   unsigned long int m2exp; /* Modulus is 2 ^ m2exp. */
  28.   const char *astr; /* Multiplier in string form. */
  29.   unsigned long int c; /* Addend. */
  30. };
  31. static const struct __gmp_rand_lc_scheme_struct __gmp_rand_lc_scheme[] =
  32. {
  33.   {32, "29CF535",      1},
  34.   {33, "51F666D",      1},
  35.   {34, "A3D73AD",      1},
  36.   {35, "147E5B85",      1},
  37.   {36, "28F725C5",      1},
  38.   {37, "51EE3105",      1},
  39.   {38, "A3DD5CDD",      1},
  40.   {39, "147AF833D",      1},
  41.   {40, "28F5DA175",      1},
  42.   {56, "AA7D735234C0DD",  1},
  43.   {64, "BAECD515DAF0B49D", 1},
  44.   {100, "292787EBD3329AD7E7575E2FD", 1},
  45.   {128, "48A74F367FA7B5C8ACBB36901308FA85", 1},
  46.   {156, "78A7FDDDC43611B527C3F1D760F36E5D7FC7C45", 1},
  47.   {196, "41BA2E104EE34C66B3520CE706A56498DE6D44721E5E24F5", 1},
  48.   {200, "4E5A24C38B981EAFE84CD9D0BEC48E83911362C114F30072C5", 1},
  49.   {256, "AF66BA932AAF58A071FD8F0742A99A0C76982D648509973DB802303128A14CB5", 1},
  50.   {0, NULL, 0} /* End of array. */
  51. };
  52. int
  53. gmp_randinit_lc_2exp_size (gmp_randstate_t rstate, mp_bitcnt_t size)
  54. {
  55.   const struct __gmp_rand_lc_scheme_struct *sp;
  56.   mpz_t a;
  57.   /* Pick a scheme.  */
  58.   for (sp = __gmp_rand_lc_scheme; sp->m2exp != 0; sp++)
  59.     if (sp->m2exp / 2 >= size)
  60.       goto found;
  61.   return 0;
  62.  found:
  63.   /* Install scheme.  */
  64.   mpz_init_set_str (a, sp->astr, 16);
  65.   gmp_randinit_lc_2exp (rstate, a, sp->c, sp->m2exp);
  66.   mpz_clear (a);
  67.   return 1;
  68. }