国际大学生程序设计竞赛--数论与算法.pptVIP

数论基本知识 信息学中应用的关于素数和整除的知识并不多，关键在于灵活应用。 素数和整除问题 如果a和b是整数，a≠0，若有整数c使b＝ac，就说a整除b。在a整除b时，记a是b的一个因子，b是a的倍数。用符号a∣b表示a整除b，a不能整除b记为a ⊥b。 整除基本性质有： （1）若a∣b， a∣c，则a∣（b＋c） （2）若a∣b，则对所有整数c， a∣bc （3）若a∣b， b∣c，则a∣c （传递性） 自反，反对称，传递性，偏序关系。 ∣,Z 是一个格。 素数（prime）和合数（compound），如果一个整数p只有1和p两个因子，则p为素数，不为素数的其它数为合数。如果n为合数，则n必有一个小于或等于n的平方根的数因子。 算术基本定理：每个正整数都可以唯一地表示成素数的乘积。其中素数因子从小到大依次出现。 数论基本知识 除法的定义和同余 a＝dq＋r。 同余：a≡b(mod c) 最大公约数gcd(a,b) 最小公倍数lcm(a,b) ab＝gcd(a,b)×lcm(a,b) 如果gcd(a,b)=1,则a与b互素。 数论基本知识 1。如何求出1～n中的所有素数？ Eraosthenes氏筛法：每次求出一个新的素数，就把n以内的它的所有倍数都筛去。 2。给出一个数n，如何判断它是不是素数？ 1）朴素的判别法 从2开始试除小于n的所有自然数，时间复杂度为O(n). 2） 如果a是n的因子，那么n/a也是n的因子，所以如果n有一个大于1的真因子，则它必有一个不大于n1/2的因子，时间复杂度O(n1/2)。 等等。这些方法太慢。 伪素数：如果n是一个正整数，并且存在和n互素的正整数a满足an-1 ≡ 1(mod n), 我们说n是基于a的伪素数。如果一个数是伪素数，它几乎肯定是素数。另一方面，如果一个数不是伪素数，它一定不是素数。 数论基本知识 Miller-Rabbin测试 function Miller-Rabbin(n:longint):boolean; begin for i:=1 to s do Begin a:=random(n-2)+2; If modular_exp(a,n-1,n)1 then return false; end; End; return true; End; 这是一个概率型算法，属于Monte-Carlo算法系列。 时间复杂度为O（s×log3n）。 最大公约数的求法 欧几里德辗转相除法 function gcd(a,b:longint):longint; Begin if b=0 then gcd:=a; else gcd:=gcd(b, a mod b); End; 扩展的欧几里德算法 如果（a，b）＝d，那么一定存在x，y满足ax＋by＝d。 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 extended_gcd:=extended_gcd(b, a mod b); t:=x; x:=y; y:=t-(a div b) * y; end; End; 佳佳的困惑 题目：给出一个数N，含数字1，2，3，4，把N的所有数字重新排列一下组成一个新数，使它是7的倍数。 分析：把数字1，2，3，4从中抽出，然后把其他数字按照原顺序排列组成的自然数w，w×10000整除7取余有7种可能，即是0，1，2，3，4，5，6。如果能把1，2，3，4排列出7个数，使它们整除7取余的值分别为0，1，2，3，4，5，6，把这4位数接在w后面即为问题的解。 除法表达式 题目：除法表达式有如下的形式：X1/X2/X3/…/Xk.其中Xi是正整数且X≤1000 000 000（1≤i≤k，k≤10 000）。除法表达式应当按照从左到右的顺序计算。可以在表达式中嵌入顺序。现在给一个除法表达式E要求告诉是否可以通过增加括号使表达式为E‘，E’是整数。 分析： 显然X1为分子，X2比为分母，其它数可以为分子也可以为分母。 方法一：将分母X2分解质因数，由于X2≤1000000000，所以质因数个数不超过29个。逐一扫描X1，X3，X4…，Xk，