今天还在看数学知识这部分,已经快看完了,主要都是一些定理,性质,推论等,还是要多看看,记住。
最大公约数:在所有a和b的公约数中最大的一个,记为gcd(a,b);最小公倍数:在所有a和b的公倍数中最小的一个,记为Icm(a,b)。 定理:对任意的正整数a,b, gcd(a,b)*Icm(a,b)=a*b。
更相减损术:对任意的正整数a,b。a>=b,有gcd(a,b)=gcd(b,a-b)=gcd(a,a-b);
欧几里得算法:对任意的正整数a,b,b0,gcd(a,b)=gcd(b,amodb);
int gcd(int a,int b){
return b ? gcd(b,a%b) : a;
}
使用欧几里得算法时复杂度为O(log(a+b))。欧几里得算法是我们最常用的方法求最大公约数,但高精度取模不容易实现,需要做高精度运算时,可考虑用更相减损术。
互质即是两个正整数的最大公约数为1.
欧拉函数:1~N中与N互质的数个数为欧拉函数,记为(N).欧拉函数有很多性质,我感觉应该很少用到。我们只要掌握怎么求欧拉函数应该就行了。最快的应该就是利用线性筛法的思想在O(N)内求出1~N每个数的欧拉函数。
int v[510000],prime[510000],phi[510000],m;
void euler(int n)
{
memset(v,0,sizeof(v));
m=0;
for(int i=2;i<=n;i++){
if(v[i]==0){
v[i]=i,prime[++m]=i;
phi[i]=i-1;
}
for(int j=1;j<=m;j++){
if(prime[j]>v[i]||prime[j]>n/i) break;
v[i*prime[j]]=prime[j];
phi[i*prime[j]]=
phi[i] * (i%prime[j] ? prime[j]-1 : prime[j]);
}
}
phi[1]=1;
}
int main()
{
int n;
cin>>n;
euler(n);
for(int i=2;i<=n;i++)
cout<<phi[i]<<' ';
return 0;
}
还有一个就是通过欧拉函数性质延伸的积性函数了,如果a,b互质,有f(ab)=f(a)*f(b),即为积性函数。