第三章-数论算法.pptVIP

III数论算法 基本概念 最大公因子及示解 孙子定理 素数判定 五道口生活网 /bbs 五道论坛 §3.1基本概念 素数：因子为1和自身 整除 ? 整除性质： （1）??? 若a|b, c为任何整数，则a|bc （2）??? 若a|b并且b|c，则a|c （3）??? 若a|b并且a|c，则a|b+c及a|b-c 线性？ 五道口生活网 /bbs 五道论坛 整除性质（续） 算术基本定理（Fundamental Theorem of Arithmetic）：任一自然数都可唯一表示成素数的积（不计因子次序） ? (4)若p为素数，p|ab，则p|a或p|b (5)若m|a并且n|a，m和n没有异于1的公因子，则mn|a. 五道口生活网 /bbs 五道论坛 最大公因子(Greatest Common Divisor) Euclid算法：gcd(a,b),a=b 即为公因子 五道口生活网 /bbs 五道论坛 一例 例 设a，b是整数，且9?a2 ? ab ? b2， 则3?(a, b)。 解： 9?(a ? b)2 ? 3ab ? 3?(a ? b)2 ? 3ab ? 3? (a ? b)2 ? 3?a ? b ? 9?(a ? b)2 9?3ab ? 3?ab 于是，3?a或3?b 结合 3?a ? b， 3?(a, b) 五道口生活网 /bbs 五道论坛 最大公因数例 例： gcd(1547，560） 1547=2*560+427 560=1*427+133 427=3*133+28 133=4*28+21 28=1*21+7 21=3*7+0 ? ∴gcd(1547，560)=7 五道口生活网 /bbs 五道论坛 欧氏算法 unsigned Euclid(unsigned int a,unsigned int b) { while(b0){ d=(unsigned int)a/b; tmp=b b=a-b*d a=tmp } return a; } 五道口生活网 /bbs 五道论坛 扩展欧氏算法 7=（1547，560） =-1*21+28 =-1*（133-4*28）+28=5*28-1*133 =5*（427-3*133）-1*133=-16*133+5*427 =-16*（560-1*427）+5*427=21*427-16*560 =21*（1547-2*560）-16*560=21*1547-58*560 7=-58*560+21*1547 x=21,y=-58 五道口生活网 /bbs 五道论坛 欧氏算法复杂度 定理：对n比特整数运行扩展的欧氏算法的时间复杂度为 定理：将扩展的欧氏算法应用于任何域上次多项式需 次域操作。 ?注意： 收敛和很快 性质： 五道口生活网 /bbs 五道论坛 性质证明 ? 证明：若 ，则 若 ，则 从而 五道口生活网 /bbs 五道论坛 孙子定理 今有物不知其數，三三數之剩二；五五數之剩三；七七數之剩二。問物幾何？ 答曰：二十三。 《孫子算經》下卷，第26題 白话： ????? 现有一堆东西不知个数，其个数被3除馀2；被5除馀3；被7除馀2；问这些东西有多少个？ 五道口生活网 /bbs 五道论坛 解法 术曰：三三数之，剩二，置一百四十；五五数之，剩三，置六十三；七七数之，剩二，置三十。并之，得二百三十三，以二百一十减之，即得。凡三三数之，剩一，则置七十五；五五数之，剩一，则置二十一；七七数之，剩一，则置十五。一百六以上，以一百五减之，即得。 五道口生活网 /bbs 五道论坛 同馀式 则下面方程组模N有唯一解： 令 ，并且当 时 。 五道口生活网 /bbs 五道论坛 应用 例：令N=15=3*5，则MOD15的数可表为： 0~（0，0），5~（2，0）， 10~（1，0） 1~（1，1），6~（0，1）， 11~（2，1） 2~（2，2），7~（1，2）， 12~（0，2） 3~（0，3），8~（2，3）， 13~（1，3） 4~（1，4），9~（0，4）， 14~（2，4） ?计算：对各分量进行： 即 五道口生活网 /bbs 五道论坛 用法 则解可表为： 实际中的用法： 求解： 及 五道口生活网 /bbs 五道论坛 例 ? 例： 扩展的欧氏算法： 五道口生活网 /bbs 五道论坛 素数的分布 ?????????????????????? ???? x ???????????????????? pi(x)----------------------------?? -------------------------?