NOIP2014复习
数论
(一)
上述代码可以进行优化,就在于i的范围,因为当i*i>n时,其余的比较全是浪费- -
上述代码时间复杂度O(logn)
-----------------------------------------------------未完待续.....----------------------------------------------
------------------------------------------------------------------------------------------------by 大神MPS
快要NOIP了,决定在家里补习一些知识点,数论是首当其冲。
§1.1质数
判断一个数是否是素数
idea 1:
void prime(int n){
if(n==2)return true;
if(n<2)return false;
int i;
for(i=1;i<=n;i++)
if(n%i==0)return false;
return true;
}
上述代码可以进行优化,就在于i的范围,因为当i*i>n时,其余的比较全是浪费- -
idea 2:
void prime(int n){
if(n==2)return true;
if(n<2)return false;
int i;
for(i=1;i<=sqrt(n);i++)
if(n%i==0)return false;
return true;
}
上述代码时间复杂度O(logn)
这里再来点优化,因为2的倍数比较多,所以- -,哪种优化你懂的。
当如果我们对n个数判断是否为素数时怎么办?
《1》直接枚举n个数,每个数用一次prime函数,时间复杂度O(nlogn)
《2》方案1显然不能使用在大数据上面,所以我们用筛法,时间复杂度O(logn)
memset(f,1,sizeof(f));
for(i=1;i*i<=n;i++)
if(f[i])
for(j=i*i;j<=n;j+=i)f[j]=0;
//0为合数,1为质数
§1.2最大公约数和最小公倍数
根据欧几里得的算法,最小公约数算法如下:
int gcd(int a,int b){
if(b==0)return a;
else return gcd(b,a%b);
}
根据数论定理,最小公倍数=a*b/gcd(a,b)