数学常用模板

数学常用模板

Euler筛（线性筛）

#include<bits/stdc++.h>
using namespace std;
int p_table[10010];
bool isPrime[10010];
int n;
void Euler(int n){
	memset(isPrime,true,sizeof(isPrime));
	memset(p_table,0,sizeof(p_table));
	p_table[0]=0;//p_table[0]代表素数的个数
	isPrime[0]=isPrime[1]=false;//0和1不是素数特判
	for (int i=2;i<=n;i++){
		if (isPrime[i]) p_table[++p_table[0]]=i;
		for (int j=1;j<=p_table[0] && i*p_table[j]<=n; j++){//i*p_table[j]即要筛的数
			isPrime[i*p_table[j]]=false;//筛掉i*p_table[j]
			if (i%p_table[j]==0) break;//注释[1]
		}
	}
}
int main(){
	scanf("%d",&n);
	Euler(n);
	for (int i=1;i<=n;i++){
		if (isPrime[i]) printf("%d ",i);
	}
}

注释

[1]
对比与埃氏筛，欧拉筛最大的优势在于每一个合数只被筛了一次，如何做到只被筛一次呢？
考虑对于一个合数n (n>=2) ，n=factor_max * p 此处factor_max 代表n的最大因子，那么易知p是一个素数，而且p比factor_max的任何一个因子都要小，即p是n的最小素因子。
综上，我们关于p可以总结出两个性质：

1. p是一个素数
2. p比factor_max中的任何一个因子都要小（或等于）

当我们在上面的程序实现的时候，如果当i%p_table[j]==0时（i 对应factor_max，p_table[j]对应p），由于我们是从小到大枚举，那么此时的p_table[j]就是性质2的零界状态，即取等，再往后p_table[j]只会越来越大，肯定比factor_max中的最小因子要大，所以此时就可以break了。（可以举例试试看）

---

快速幂

#include<bits/stdc++.h>
using namespace std;
long long b,p,k,base,ans=1,ap;
int main()
{
    cin>>b>>p>>k;
    base=b;
    ap=p;//ap为输出需要，存初始的p
    while(p>0)
    {
        if (p&1)
        {
            ans=ans*base;
            ans=ans%k;
        }
        base=base*base;
        base=base%k;
        p>>=1;
    }
    cout<<b<<"^"<<ap<<" mod "<<k<<"="<<ans%k;
}

---

最大公约数

#include<bits/stdc++.h>
using namespace std;
int x,y;
int gcd(int a,int b){
	return b?gcd(b,a%b):a;
}
int main()
{
	scanf("%d%d",&x,&y);
	printf("%d",gcd(x,y));
	return 0;
}

---

扩展欧几里得算法

求关于x的同余方程 ax≡1 (mod b) 的最小正整数解。