竞赛辅导——数论问题.pptVIP

数论问题 Lecture by Hao Wu 常见的数论问题 1.数的幂运算 2.欧拉定理的应用 3.素数测试 4.Pell方程 5.其它 案例1：模的运算 描述 人与人是不同的,有些人喜欢阅读满是女孩子的杂志，有些人喜欢在地下室引爆炸弹，而还有一些却喜欢一些麻烦的数字游戏。比如ESSE论坛的一次活动： 每个人选择两个数字Ai和Bi写在纸上，其他人不能看见。过了一段时间后，每个人说出自己纸上的数字，然后每个人的目标是求出所有的AiBi的和模M的值，最先算出结果的，就是胜利者。 作为一个程序员，你当然有办法编一个程序，以最快的速度算出结果，赢得比赛。 相当于求 A1B1+A2B2+...+AHBH) mod M. 解题思路（1） 取模运算可以与加、减、乘、乘方交换次序； (a+b)%m=a%m+b%m a*b%m=(a%m)*(b%m) ab%m=(a%m)b 解题思路（2） 有一类常见的运算，就是求形如的值，一般情况下，我们可以用幂的定义来做：通过乘法的倍乘来实现，需要做b次乘法和模运算。如果b比较大时，这样的方法就不行了，本节介绍一种只需要O(log2b)次乘法和模运算的方法。 解题思路（3） 首先将指数变成二的幂次和的形式。如5=1+4=20+22，11=1+2+8=20+21+23。如果b可以表示为的形式，那么ab可以表示为 解题思路（4） 表面上看这样的表达式更复杂了，但注意到上面的表达式的乘积项不超过(log2b)，同时，形如（i=0,1….）有如下的递推式： 解题思路（5） 后项可以通过前项的平方来得到。所以可以先的值求出来，不超过(log2b)次乘法。 核心代码及解释 （1） voidcal(unsigned*r,longunsigneda,unsignedm){ //计算a的2j次方对m取模的结果,用数组r保存 r[0]=a%m; for(i=1;i32;++i) r[i]=(long)r[i-1]*r[i-1]%m;//ri+1=ri*ri } 核心代码及解释 （2） unsignedcalc(unsignedm,longunsigneda,longunsignedb){ //函数计算ab%m的值 intbl[32],r[32];//bl是b的二进制结果，r[i]=a的2j次方对m取模的结果 cal(r,a,m);//计算r[i] while(b){//将b用二进制形式表示,结果保存在bl中 bl[i++]=b%2; b=1;} for(j=0,k=1;ji;++j)//计算乘方，时间复杂度为i=O(log2b) if(bl[j])//如果b对应二进制位bj为1,则对应的rbj在幂的乘积表达式中，否则不在 k=k*r[j]%m; Return k;} 核心代码及解释 （3） Int main(){ scanf(%d,z); while(z--){ scanf(%u%d,m,i); for(j=s=0;ji;++j){ scanf(%lu%lu,a,b); s=(s+calc(m,a,b))%m;}//计算AjBj的值，并将其和求出来，结果用变量s保存 printf(%u\n,s);}} 欧拉定理1 如果p是素数，则对任意的a，有 欧拉定理2 如果p不是素数，则对任意的a，有 案例2：快乐2004 Time Limit: 1000ms Special Time Limit:2500ms Memory Limit:32768KB 描述 考虑到一个正整数X，S是所有2004X的因子的和，你的目标是求出S除以29的余数。 看看X=1的例子，20041的所有因子为1,2,3,4,6,12,167,334,501,668,1002和2004，因此S=4704，而4704 mod 29 =6。 输入 输入有多个测试序列，每个测试序列一行一个正整数X，(1=X。 X=0表示输入结束，并且不需要处理。 输出 对每个测试序列，输出一行就是S除以29的余数。 输入样例 1 10000 0 输出样例 6 10 解题思路（1） 题目要求的是（2004X）的因子和对29取模的结果 求某个数A的所有因子和，首先将该数分解： 其中p1,p2,...,pn都是素数，则A的任意因子可以表示为： 解题思路（2） 如：2004=22*3*167，则 解题思路（3） 解题思路（4） X还很大( )，根据Euler定理 解题思路（5） 解题思路（6） 当X=10000时，计算200410000=220000*310000所有因子和对29取模的结果 核心代码及解释 （1） int mod29(int x) { for