ACM程序设计12.pptVIP

ACM程序设计 第十二讲 若干数论问题 剩余类 完全剩余系 简化剩余系 欧拉函数 素数判定 大整数分解 剩余类 按照某一模m的余数将全体整数进行分类 例如 奇数和偶数 就是模2的剩余类 完全剩余系 在每个对模m同余的剩余类中选出一个数构成具有m个元素的集合，称为完全剩余系 有下面结论： 任意m个连续整数构成对模m的剩余系 0，1，2，…，m-1构成对模m的最小非负剩余系 m=2k+1时，-k,-k+1,…,k-1,k构成对模m的最小绝对值剩余系；m=2k时， -k,-k+1,…,k-1构成对模m的最小绝对值剩余系 设gcd(b,m)=1，c为任意整数，a1…am是对模m的一个剩余系，那么ba1+c,ba2+c,…,bam+c也是对模m的一个剩余系 简化剩余系 在模m的每个剩余类中取出一个和m互质的数构成的集合，称为简化剩余系 一个剩余系里的元素要么都和m互质，要么都不和m互质 欧拉函数 把对模m的简化剩余系的元素个数称为m的欧拉函数，记做Φ(m) 欧拉函数表示0,1,…,m-1中和m互质的数的个数 Φ(1)的取值是0还是1，根据具体题目要求决定 求欧拉函数的方法： 将m分解质因数，设其所有两两不同的质因数为p1p2…pn,则Φ(m)=m*∏(1-1/pi) 求一个连续区间数的欧拉函数 方法1：一个一个按上面公式求 太慢了！ 方法2：筛法 设一个数组phi[11], 赋初值phi[1]=1,phi[2]=2,…,phi[10]=10 对数组下标是2的倍数的每个元素，分别乘上（1-1/2） 对数组下标是3的倍数的每个元素，分别乘上（1-1/3） 对5对7（所有素数）都类似操作，就得到了欧拉函数值 注意：一般需要先求出此区间的所有的素数 非常精简的筛法求欧拉函数的模板 const int Max=10001; int eu[Max] void eular() { int i,j; for(i=2;iMax;i++)if(!eu[i]) for(j=i;jMax;j+=i){ if(!eu[j])eu[j]=j; eu[j]=eu[j]/i*(i-1); } } 例题 Relatives（POJ2407） 题意： 给出n，求小于等于n且与n互质的元素个数 n1,000,000,000 思路： 显然，本题是求n的欧拉函数 素数判定 试除法 米勒罗宾测试 试除法 筛法求1，000，000以内的素数表 int prime[200000]; bool bol[1000001]; int prime_table() { int I,j,flag=0; memset(bo,0,sizeof(bo)); bo[0]=bo[1]=1; for(i=2;i=1000;i++) if(!bo[i]){ j=i*i; for(;j=1000000;j+=i)bo[j]=1; } for(i=0;i=1000000;i++)if(!bo[i])prime[flag++]=i; return flag; } 米勒罗宾测试 人们发现伪素数非常少，因此对费马小定理中a随机取若干个值，如果每次答案都是素数，那么它不是素数的概率就大大减少了，这就是米勒罗宾算法的原理，也是随机算法的一个著名的应用例子。 米勒罗宾算法实现 bool Miller_Rabin(long long n) { if(n2)return false; if(n==2)return true; if(!(n1))return false; long long k=0,i,j,m,a; m=n-1; while(!(m1))m=1,k++; for(i=0;i10;i++){ if(pri[i]=n)return true; a=multimod(pri[i],m,n); if(a==1)continue; for(j=0;jk;j++){ if(a==n-1)break; a=multi(a,a,n); } if(j==k)return false; } return true; } Java语言中的Miller-Rabin函数 Java的BigInteger类提供了一个利用Miller-Rabin算法判断素数的函数 boolean isProbablePrime(int certainty) 如果此大整数可能为素数，返回true 如果一定是合数，