数论(算法基础)

一、判断质数        ( sqrt ( N ) )

bool isPrime(int n)
{
	if(n<2) return false;
	for(int i=2;i<=n/i;i++)
		if(n%i==0) return false;
	return true;
}

二、分解质因子        ( N )和 ( Log(N) ~ sqrt ( N ) ) 

void divide (int n)
{
	for(int i=2;i<=n;i++)
		if(n%i==0)
		{
			int s=0;
			while(n%i==0){
				s++;
				n/=i;
			}
			printf("%d %d\n",i,s);
		}
}
void divide (int n)
{
	for(int i=2;i<=n/i;i++)
		if(n%i==0)
		{
			int s=0;
			while(n%i==0){
				s++;
				n/=i;
			}
			printf("%d %d\n",i,s);
		}
	
	if(n>1) printf("%d 1\n",n);
}	

三、朴素筛质数        ( n Log ( n ) )

const int N = 10e4+10;
bool st[N];
int prime[N];

int getPrime(int n)
{
	int cnt=0;
	for(int i=2;i<=n;i++){
		if(!st[i]){
			prime[cnt++]=i;
		}
		for(int j=i+i;j<=n;j+=i) st[j]=true;
	}
	return cnt;                        //返回 2~n 中质数个数
}

四、埃氏筛法                (n loglog n)

const int N = 10e4+10;
bool st[N];
int prime[N];

int getPrime(int n)
{
	int cnt=0;
	for(int i=2;i<=n;i++){
		if(!st[i]){
			prime[cnt++]=i;
            for(int j=i+i;j<=n;j+=i) st[j]=true;
		}
		
	}
	return cnt;
}

五、线性筛法                ( n ) 

const int N = 10e8;
bool st[N];
int prime[N];	
	
int getPrime(int n)
{
	int cnt=0;
	for(int i=2;i<=n;i++)
	{
		if(!st[i]) prime[cnt++]=i;
		for(int j=0;prime[j]<=n/i;j++){
			st[prime[j]*i]=true;
			if(i%prime[j]==0) break;
		}
	}
	return cnt;
}

六、试除法求约数

vector <int> getDivisors(int n)
{
	vector <int> res;
	for(int i=1;i<=n/i;i++)
	{
		if(n%i==0){
			res.push_back(i);
			if(i!=n/i) res.push_back(n/i);
		}
	}
	sort(res.begin(),res.end());
	return res;
}

七、约数个数  && 约数之和

//约数个数
#include <bits/stdc++.h>
using namespace std;

const long long MOD = 10e9+7;
map <int,int> mp;

int main()
{
	int n;
	long long res=1;
	cin>>n;
	for(int i=2;i<=n/i;i++)
	{
		while(n%i==0){
			n/=i;
			mp[i]++;
		}
	}
	if(n>1) mp[n]++;
	for(auto p:mp){
		res=res*(p.second+1)%MOD;
	}
	cout<<res;
	return 0;
}
//约数求和 
#include <bits/stdc++.h>
using namespace std;

const long long MOD = 10e9+7;
map <int,int> mp;

int main()
{
	int n;
	long long res=1;
	cin>>n;
	for(int i=2;i<=n/i;i++)
	{
		while(n%i==0){
			n/=i;
			mp[i]++;
		}
	}
	if(n>1) mp[n]++;		//分解质因数
	 
	for(auto t:mp){
		int p=t.first,a=t.second;
		long long k=1;
		while(a--) k=(k*p+1)%MOD;
		res=res*k%MOD;
	} 
	cout<<res;
	return 0;
}

八、辗转相除法

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

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值