数论

质数

试除法判断质数

#include<iostream>
using namespace std;
int ai[105];

int main()
{
	int n;
	bool flag=true;
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>ai[i];
	for(int i=1;i<=n;i++)
	{
		flag=true;
		if(ai[i]==1)
		{
			cout<<"No"<<endl;
			continue;
		}
		for(int k=2;k<=ai[i]/k;k++)
		{
			if(ai[i]%k==0)
			{
				cout<<"No"<<endl;
				flag=false;
				break;
			}
		}
		if(flag==true)
			cout<<"Yes"<<endl;
	}
}

分解质因数

#include<cstdio>
using namespace std;
int ai[105];

int main()
{
	int n;
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	scanf("%d",&ai[i]);
	for(int i=1;i<=n;i++)
	{
		for(int k=2;k<=ai[i]/k;k++)
		{
			if(ai[i]%k==0)
			{
				int s=0;
				while(ai[i]%k==0)
				{
					ai[i]/=k;
					s++;
				}
				printf("%d %d\n",k,s);

			}
			
		}
		if(ai[i]>1)
			 printf("%d %d\n",ai[i],1);
		printf("\n");
	}
}

筛质数

埃氏筛法O(NlogN)

#include<iostream>
using namespace std;
bool prme[1000005];

int main()
{
	for(int i=2;i<=1000005;i++)
	{
		if(!prime[i])
		{
			for(int k=2*i;k<=1000005;k+=i)
				prime[k]=1;
		}
	}
}

欧拉筛法（线性筛法）O(N)

#include<iostream>
using namespace std;
int st[10000005],prime[10000005];

int main()
{
	int n,tot=0;
	cin>>n;
	st[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!st[i])
			prime[++tot]=i;
		for(int k=1;prime[k]<=n/i;k++)
		{
			st[prime[k]*i]=1;
			if(i%prime[k]==0)
				break;
		}
	}
	cout<<tot;
}

约数

试除法求约数

#include<iostream>
#include<algorithm>
using namespace std;
int ai[1005],num[2005],tot;

int main()
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>ai[i];
	for(int i=1;i<=n;i++)
	{
		tot=0;
		for(int k=1;k<=ai[i]/k;k++)
		{
			if(ai[i]%k==0)
			{
				num[++tot]=k;
				if(ai[i]/k!=k)
					num[++tot]=ai[i]/k;
			}
		}
		sort(num+1,num+1+tot);
		for(int i=1;i<=tot;i++)
			cout<<num[i]<<" ";
		cout<<endl;
	}
}

约数个数

分解质因数：n=p1^a1*p2a2…*pk^ak

公式： num=(a1+1)(a2+1)…(ak+1)

#include<iostream>
#include<unordered_map>
using namespace std;
const int MOD=1e9+7;
unordered_map<int,int> num;
int ai[105],tot;

int main()
{
	int n;
	long long res=1;
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>ai[i];
	for(int i=1;i<=n;i++)
	{
		tot=0;
		for(int k=2;k<=ai[i]/k;k++)
		{
			if(ai[i]%k==0)
			{
				while(ai[i]%k==0)
				{
					ai[i]/=k;
					num[k]++;	
				}	
			}
		}
		if(ai[i]>1)
			num[ai[i]]++;
	}
	for(auto i:num)
		res=(res*(i.second+1))%MOD;
	cout<<res;
}

约数之和

分解质因数：n=p1^a1*p2a2…*pk^ak

公式： num=(p1^0+p11+…+p1^k)(p20+p2^1+…+p2k)
…(pk^0+pk1+…+pk^k)

#include<iostream>
#include<unordered_map>
using namespace std;
const int MOD=1e9+7;
unordered_map<int,int> num;
long long ai[105],res=1;

int main()
{
	int n;
	long long res=1;
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>ai[i];
	for(int i=1;i<=n;i++)
	{
		for(int k=2;k<=ai[i]/k;k++)
		{
			if(ai[i]%k==0)
			{
				while(ai[i]%k==0)
				{
					ai[i]/=k;
					num[k]++;	
				}	
			}
		}
		if(ai[i]>1)
			num[ai[i]]++;
	}
	for(auto i:num)
	{
		int a=i.first,b=i.second;
		long long t=1;
		for(int i=1;i<=b;i++)
			t=(t*a+1)%MOD;
		res=(res*t)%MOD;
	}
	cout<<res;
}

最大公约数

欧几里得定理

#include<iostream>
using namespace std;

int gcd(int a,int b)
{
	return b?gcd(b,a%b):a;
}

int main()
{
	int n;
	cin>>n;
	while(n--)
	{
		int a,b;
		cin>>a>>b;
		cout<<gcd(a,b)<<endl;
	}
}

To be continued