NOIP数论.pptxVIP

数学 素数和整除问题 素数筛法 O(NlogN)算法不再赘述 介绍一种O(N)算法： 将i从2到N扫描 从前i未被筛过则存入prime数组 用j从1循环，用i*prime[j]筛素数，若i mod prime[j] = 0 跳出 不难发现，每个合数被删仅一次 x = a^k1 * b^k2 ( a,b 为素数,ab) i = a^(k1-1) * b^k2 prime[j] = a 时 x 被筛去 复杂度为O(N) for (i=2;i=n;i++) { if (!hash[i]) prime[++tot]=i; for (j=1;prime[j]*i=n;j++) { hash[prime[j]*i]=1; if (i%prime[j]==0) break; } } 素数和整除问题 最大公约数与最小公倍数 高精度GCD? 模拟单精度实现过程太过麻烦 设a=b，则 Gcd(a,b)= ①Gcd(a/2,b/2) *2 ……a%2=0 b%2=0 ②Gcd(a,b/2) ……a%2=1 b%2=0 ③Gcd(a/2,b) ……a%2=0 b%2=1 ④Gcd(b,a-b) ……a%2=1 b%2=1 为什么是log级别的？ 素数和整除问题 (NOIp 2009) 已知 Gcd(x,a0)=a1 Lcm(x,b0)=b1 求满足条件的x个数 (a0,a1,b0,b1=2,000,000,000，2000组数据) x = a^kx * tx a0= a^ka0 * ta0 a1= a^ka1 * ta1 b0= a^kb0 * tb0 b1= a^kb1 * tb1 min (kx,ka0)=ka1 max(kx,kb0)=kb1 若 ka0=ka1 则 kx=ka1 否则 kx=ka1 若 kb0=kb1 则 kx=kb1 否则 kx=kb1 可由此确定每个素数因子的范围 乘法原理求总解数 sqrt(2,000,000,000)内仅有4000+素数 扩展欧几里得 扩展欧几里得： 如何求解 ax+by=gcd(a,b) 回顾欧几里得算法 int gcd(int a,int b) { return b?gcd(b,a%b):a; } 递归到底时 a=gcd(a,b)，b=0，令x=1，y=0，满足ax+by=gcd(a,b) a%b=a-a/b*b ax+by=c - bx’+(a-a/b*b)y’=c xa+yb=c - y’a+(x’-a/b*y’)b=c 回溯时令 x=y’ y=x’-a/b*y’ 扩展欧几里得 代码求ax+by=gcd(a,b) #include stdio.h int extended_gcd(int a, int b, int *x, int *y) { int ret, tmp; if (!b) { *x = 1; *y = 0; return a; } ret = extended_gcd(b, a % b,x, y); tmp = *x; *x = *y; *y = tmp - a / b * (*y); return ret; } int main() { int a,b,x, y, z; scanf(%d%d,a,b); z=extended_gcd(a,b, x, y); printf(%d %d %d\n,z,x,y); system(pause); return 0; } 扩展欧几里得 求解 ax+by=c ? 判断c是否被gcd(a,b)整除，整除则有解，反之无解 求解ax+by=gcd(a,b) 若ax+by=gcd(a,b)的解为(x’,y’) 则ax+by=c的解(x,y)=(c/gcd(a,b)*x’, c/gcd(a,b)*y’) 如何求出不定方程的所有整数解？ (x+k*b/gcd(a,b)，y-k*a/gcd(a,b))为所有合法解 其中k∈Z 简单排列组合计数 排列与组合 排列数： N个不同物体不重复地取M个做排列的方法数P(N,M)=N!/(N-M)! 组合数： N个不同物体不重复地选取M个的方法数C(N,M)=N!/((N-M)!*M!) 可重排列： K种元素，重数分别为n1,n2,…,nk，所有n=Σni个元素的全排列数为n!/(n1!n2!..nk!) K种元素，重数