初等数论 (4)课件.pptVIP

初等数论 雅礼 朱全民 内容介绍 素数及其性质 欧几里得公式gcd(a,b) 计算n的最大互质数 欧几里得公式推广：d=gcd(a,b)=ax+by 计算同余方程 ax≡b(mod n)(n0) 求解同余式组a≡ai(mod ni）(1≤i≤k) 计算不定方程ax+by=c 数论的应用 初等数论 基本概念 整除和约数 素数和合数 唯一分解定理 除法的定义和同余 最大公约数(gcd)和最小公倍数(lcm) 互素、两两互素 素数测试 问题一：求1~n的所有素数 Eraosthenes的筛子 对于素数p，删除p*p, p*(p+1), p*(p+2)… 例：求[n, n+m]内的所有素数 n很大 m很小 如何筛？ 如何加速？ 素数测试 问题二：给数n，它是素数吗？ 朴素算法：枚举sqrt(n)内约数 改进算法：枚举sqrt(n)内素数 预处理出sqrt(n)内素数：筛？ Monte-Carlo算法 不断选基b，重复以下测试 测试越多，正确概率越大 欧几里得公式gcd(a,b) gcd(a,b)=gcd(b, a mod b) 程序：递归求解 时间复杂度O(logb) 传球游戏 N个人围圈玩传球游戏，开始时第一个人拿着球，每个人把球传给左手的第K个人。满足1≤K≤N/2。求K的最大值，使得第一个人重新拿到球之前，每个人都拿过球。 输入：数串n（3 ≤ N ≤ 10200） 输出：K的最大值 分析 如果n是奇数，那么最大的可能的[n/2]。由于k*2=n-1，所以gcd(k*2,n)=1，可以得到gcd(k,n)=1，所以直接输出[n/2]即可。 如果N是偶数且为4的倍数：则N/2显然不可能是符合要求的K，接着试验N/2-1，显然有gcd(N/2-1, N/2)=1，而由于N是4的倍数，所以N/2是偶数，N/2-1是奇数，所以gcd(N/2-1, N)=1。此时满足题意的K就是N/2-1。 N是偶数但不是4的倍数。此时N/2-1也是偶数，显然不可能是正确答案。那么就来考察一下N/2-2：这是一个奇数。根据刚才的介绍，有gcd(N/2-2,N/2)≤2，但是N/2是奇数，所以gcd(N/2-2,N/2)=1。最后同样可以得到gcd(N/2-2,N)=1。所以此时的答案就是N/2-2。 欧几里得推广 用一个线性组合来表示最大公约数，即d=gcd(a,b)=ax+by。要求计算出其中的整系数X和Y (X与Y可能为0或负数) 分析： 若b=0,则gcd(a,b)=a，x=1，y=0； 若b0,则首先计算出d’=gcd(b,a mod b)和满足d’=bx’+(a mod b)y’的x’、y’。在这种情况下, 有d=gcd(a,b)=d’=gcd(b,a mod b)。 为了得到满足d=ax+by的x和y，我们将d’=gcd(b,a mod b)转化成线性组合d’=bx’+(a mod b)y’=bx’+(a-[a／b]b)y’ =ay’+b(x’-[a／b]y’)。对应d=ax+by，选择x=y’;y=x’-[a／b]*y’就可以满足等式。 扩展的Euclid算法 function extended_gcd(a,b:longint; var x,y:longint):longint; begin if b=0 then begin extended_gcd:=a; x:=1; y:=0; end else begin //计算a与a mod b的最大公约数d’和满足d’=ax’+a mod b*y’的x’、y’值。 //x=y’;y=x’-[a／b]*y’ extended_gcd:=extended_gcd(b, a mod b, x, y); t:=x; x:=y; y:=t-(a div b)*y end; end; 满足ax+by=d的数对（x,y）不是惟一的 因为当x增加b，y减少a，和不变。 除法表达式 除法表达式有如下的形式： X1 / X2 / X3 / … / Xk 其中Xi是正整数且Xi≤109 (1≤i≤k，k≤10 000)。 除法表达式应当按照从左到右的顺序求和， 例如表达式1/2/1/2的值为1/4。 可以在表达式中嵌入括号以改变计算顺序， 例如表达式(1 / 2) / (1 / 2)的值为1。 现在给一个除法表达式E要求告诉是否可以通过增加括号使表达式值为E′，E′是整数。 分析 X2必须在分母 其他都可以在分子 最后结果是整数吗？ X2分解因数？ 复杂度？？？ 每次约掉X2和Xi的最大公约数！ 线性同余方程ax≡b(mod n)(n0) 已知同余方程ax≡b(mod n)(