初等数论程耀.pptxVIP

素数和合数算术基本定理：任意正整数n可以写成n=2a1*3a2*5a3*…,其中a1,a2,a3等为非负整数素数的判定： ①判定函数 ② 筛法求素数 ③miller_rabin素性判定素数的判定bool Is_prime( int n){ int t = (int )sqrt (n); for ( int i =2; i = t ; i ++） if(n%i ==0) return false; return true;}筛法求素数思考：如果一个数是合数，那么它的素因子都比它小？？？这样说来：如果我们的当前数是 a ,那么所有 a的倍数（当然是2倍以上啦）都不会是素数，可以这样看吧？于是，我们可以一种新的素数判定方法。筛法求素数方法：每次用一个素数，去筛掉所有它的倍数。举个例子：2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 2829 30①筛去2的倍数，剩下3 5 7 9 11 13 15 17 19 21 23 25 27 29②筛去3的倍数，剩下5 7 11 13 17 19 25 29 。。。。。。注意：开头的数一定是素数？？？筛法求素数bool prime[ maxn ];//时间复杂度O(n)???void init ( ){ memset(prime, 0 ,sizeof(prime)); for( int i=4;i maxn; i +=2) prime[i]=1; for( int i=3; i maxn; i+=2) if( ! prime[ i ]){ for ( int j = i * i ; j maxn;j+=i) prime[ j ]=1; }}miller_rabin素性判定miller_rabin是一种概率型的算法，不是确定型的。但是，只要运行次数足够多，一般出错的概率是非常小的，一般10次就好了！算法复杂度是 O((logn)^3)算法的正确性是基于费马小定理。费马小定理：对于素数p和任意整数a，有 a^(p-1)=1(mod p)。反过来，满足a^(p-1)= 1(mod p),p也几乎一定是素数。miller_rabin算法分析bool miller_rabin(LL n){ if(n2) return false; if(n==2) return true; if(!(n1)) return false; for(int i=1;i=12;i++){ LL a=random(n-2)+1; if(witness(a,n)) return false; } return true;}witness函数witness函数用来搜集n是合数的证据。bool witness(LL a,LL n){ LL m=n-1; int t=0; LL y; while(!(m1)){ m=1;t++;}//这个地方是通过一个推论来优化的，x^2=1(mod n),当那个n是合数的时候，就会出现非平凡平方根！ LL x=quick_mod(a,m,n); for(int i=1;i=t;i++){ y=muti(x,x,n); if(y==1 x!=1 x!=n-1) return true; x=y; } if(y!=1) return true; else return false; }快速幂取模题目：求出m^n ( mod p) 的值，其中 m=10000, n=100000000。分析：如果直接边乘然后边取模，直接导致 TLE！我们需要一个更优的算法。 我们可以发现：这其中包含着大量的重复的子问 题，利用分治的思想可以大大减小运算量。 举个例子：2^7 = 2^4 * 2^2 * 2^1,而2^t到2^(2t)只需要累乘就好了。 快速幂取模算法分析把指数看成一个二进制数，从低位到高位依次判断是否为1。如果是1，就要乘以2^t,否则，就不需要乘上2^t. 其中，ak等系数只能取0，或者1。如果取1的话，那么要乘以对应的a^(2^k). 算法复杂度O(log n)快速幂取模代码分析 int quick_mod ( int a , int n ){ int t =a , ret =1; while( n ){ if( n 1) ret=ret * t%n; n=1; //n右移两位，相当于除2 t = t * t %n; } return ret ; }PS：与快速幂取模类似的还有矩阵快乘，这里不展开最大公约数（gcd)普通算法：求出c=min( a , b),如果找到c|a c|b,那么c减一，然后循环继续，直到 找到 c满足条件