信息学竞赛中数论常见问题.docVIP

Prime 概述: 1的数,除了1和本身没有其他因子. 1既不是素数也不是合数,0和所有的负整数同样如此. 100以内的素数 {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97} 1.打表(1) #define n 1000003 int prime[n]; //要判断的范围是多大,就要设置多大的数组,对应于某个数,看它下标对应的数是1，则是素数,否则,不是素数 void getPrime(){ int i,j; int bound=sqrt((double)n); for(i=2;in;i++)prime[i]=1;//将所有的数置1，表示这些数都是素数. for(i=2;ibound;i++){//注意从2开始 for(j=i+i;jn;j+=i)//将2的倍数置成0，表示不是素数 prime[j]=0; } } 打表(2) int prim[50],pn=0; bool pp[101]; void pre() { int i,j; for (i=2;i=100;i++) { if (!pp[i]) prim[pn++]=i; for (j=0;jpni*prim[j]=100;j++) { pp[prim[j]*i]=1; if (i%prim[j]==0) break; } } } Pn是1至100范围内素数的个数,prim[i]一次存放这pn个素数.这种方法的好处是存储素数时节省空间,输出素数时节省时间 直接判断 bool isPrime2(int x){ if(x=1)return false; if(x==2)return true; if(x%2==0)return false; double bound=sqrt((double)x); for(int i=3;i=bound;i+=2){ if(x%i==0){ return false; } } return true; } GCD(最大公约数)---欧几里得算法 a,b为正整数，设集合A = {x*a+y*b|x,y是整数}，则A中最小正元素是gcd(a,b) long kgcd(long a,long b) { if (a==0)return b; if (b==0)return a; if (!(a1)!(b1))return kgcd(a1,b1)1; else if (!(b1))return kgcd(a,b1); else if (!(a1))return kgcd(a1,b); else return kgcd(abs(a-b),min(a,b)); } LCM (最小公倍数) LCM ( a, b ) = a * b / GCD ( a, b ) 实际上最好写成a/GCD(a,b)*b long lcm(long a,long b) { long c,d,sw; c=(a=b)?a:b; d=(a=b)?a:b; while (c%d!=0) { sw=c%d; c=d; d=sw; } return (a/d)*b; } 求多个数的lcm,需要将res初始化为1 相关题目: 1./showproblem.php?pid=1019 Least Common Multiple 2./showproblem.php?pid=1222 wolf and rabbit 扩展欧几里得(可以求解模线性方程) a*x+b*y=gcd(a,b)一系列解是x=x+b,y=y-a int extgcd(int a,int b,int x,int y){ if(b==0){ x=1; y=0; return a; } int d=extgcd(b,a%b,x,y); int t=x; x=y; y=t-a/b*y; return d; } a*x+n*y=b==ax=b(mod n)模线性方程 定理：方程ax=b(mod n)对于未知量x有解，当且仅当gcd(a,n)|b. 令d=gcd(a,n). 设为x ,y为所求 满足ax+ny=gcd(a,n). 则原方程有一解x0=x*(b/d)