素数的一个特殊性质及其用于伪随机数生成的方法.pdfVIP

第30卷 第 3期 北 京 化 工 大 学 学 报 Vo1．30。No．3 2003正 JOURNALOFBELLING UNIVERSITY OFCHEM ICALTECHNOLOGY 2003 素数的一个特殊性质及其用于伪随机数生成的方法 李世刚 刘 辉 陈标华 (北京化工大学可控化学反应科学与技术基础教育部重点实验室，北京 100029) 摘 要：提出素数的一个特殊性质。定义了一类超素数。证明了相关的定理。基于上述理论分析，提出一种伪随机 数生成的新方法——超素数法，统计结果表明本文方法具有 良好的统计特性，由此得到的伪随机数序列可用作伪 随机数发生器。文中给出了计算方法和数值示例。 关键词：数论 ；素数；伪随机数；MonteCarlo模拟 中囝分类号：TQ021．1 定义 1 记 M 是素数，记 zf∈N，且 0Zf 引 言 M ，如果存在既约真分数 可表为纯循环小数 0． 在计算化学和计算物理研究中，MonteCarlo模 口1口2…口 1口2… 口f…，且其循环节数为 t= (M 一 拟是一种非常重要的计算方法u．2J。MonteCarlo模 拟是以概率论和数理统计理论为基础 的计算方法。 1)，则称具有前述性质的素数 M 为超素数。 不可避免要用到随机数 u(0，1)。然而，生成真正 下面的两个例子是这类素数的例证和说明。 的 u(0，1)非常困难。一般都是按照一定的规则在 例 1，素数 7，当以1除以7，得循环小数为 r= 计算机上生成伪随机数序列 u(0，1)。在实际应用 0．142857142857…。可见，本例的循环节数为 (7— 中，随机数发生器的统计性质非常重要。Ferrenberg 1)：6。 · 等人 3【J曾报道因使用 Trausworthe序列而得到完全 例 2，并非所有素数都满足定义 1中的生成纯 错误的物理结果的算例。因此，研究性能优 良的随 循环小数的循环节数 t=(M 一1)的性质，如素数 2， 机数发生器具有重要的理论和实用价值。 3，5和 11等： 为获取伪随机数序列，提高模拟结果的有效性， ：0．5 1 ：0．333… 人们已提出多种伪随机数发生器或生成方法，这方 面近期的进展 已有文献报道 【．5]。与文献报道 的方 =0．2 =0．66090… 法不同，本文利用素数的一个特殊性质，提出一种伪 关于纯循环小数的存在条件，数论中有如下定 随机数生成的新方法——超素数法，文中给出了相 关的定理和证 明，并对生成的伪随机数序列进行了 理 统计检验，结果表明本文方法具有 良好的统计特性， 定理l… 有理数詈，0口6，(口，6)=1，则 可用作 MonteCarlo模拟的伪随机数发生器。 有理数 能表成纯循环小数的充分必要条件是 1 素数的一个特殊性质 (b，10)=1。 素数(或称质数)为只能被其本身和 1整除的整 按照本文给出的定义 1，显然 (M ，10)=1且 数[引。研究发现，在素数 中。存在一类特殊的素数 。 (M，Z)=1，可见超素数生成纯循环小数的结论是 它具有如下特殊的性质： 显然的