概率算法新.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * §5.2 素数测定(数的素性测定) 当n为合数，对?a∈[1，n-1]，有 吗？ 若成立，则只要n是合数，? 否则n为素数。遗憾的是此命题亦不成立。 ?结论：我们不能通过验证an-1modn是否为1来判定n是否为素数 例如： ∵ ∴ 若是 时呢？ 设n=15(合数)， * * §5.2 素数测定(数的素性测定) 3.伪素数和素性伪证据 设2≤a≤n-2，一个满足 （即n可整除an-1-1）的合数n称为以a为底的伪素数，a称为n的素性伪证据。 实际上，符合费马小定理逆命题的数，我们称为拟素数(Probable Prime) 。 拟素數一是素數，二是伪素數。 后来数学家渐把符合一些素数性质的逆命题的合数称为伪素数 伪素数有多少？ 在前10亿个自然数中，有50,847,534个素数，而以2为底的伪素数则只有5597个。 因此在n整除2n-1 -1的情况下，出现合数的机会仅有5597/(5597 = 0.00011 。 而我们同时考虑以2和3为底的伪素数，则只有1272个。 因此n同时可以整除2n-1 -1和整除3n-1 -1的情况下，碰上合数的机会便更低了，仅有1272/(1272 = 0.000025。 若我们把该测试扩张至其他底，自然会将找到伪素数的机会降低，会降至0吗？否！ * * §5.2 素数测定(数的素性测定) 若将Fermat测试改为从2~n-2之间随机选a，则只有选到一个伪证据时，对合数的测试失败(返回true). 伪证据有多少？ 总体情况是伪证据相当少 1000之内的奇合数测试误差概率3.3%。n较大时，概率更小。 有些合数伪证据比例相当高 如561，有318个伪证据，超过证据数的一半（2~559）。极端情况：Fermat(651693055693681)返回true的概率99.9965% 一般地，对 ，存在无穷多个合数，使得Fermat测试发现他们是合数的概率小于δ 即：对任意的p0， Fermat测试都不是p-正确的。因此，以前将Fermat测试重复固定次数，并不能将误差降到任意小的ε内。 * * §5.2 素数测定(数的素性测定) Fermat测试改进 强伪素数 设n是一个大于4的奇整数，s和t是使得 的正整数，其中t为奇数，设B(n)是如下定义的整数集合： 当且仅当2≤a≤n-2且满足下述2个条件之一： 或 ? 当n为素数时， ，均有 当n为合数时，若 ，则称n为一个以a为底的强伪素数，称a为n素性的强伪证据。 n为素数，说明它对所有底均为强伪素数。 * * §5.2 素数测定(数的素性测定) Btest(a，n){//n为奇数， 返回 。即返回//真说明n是强伪素数或素数 s←0； t ←n-1； // t开始为偶数 repeat s++；t ← t÷2； until t mod 2 = 1； //n-1=2st ， t为奇数 x ←at mod n； if x=1 or x=n-1 then return true；//满足①or②， for i ←1 to s-1 do{ //验证 x ← x2 mod n； if x=n-1 then return true； //满足②， } return false； } *