关于质数和约数问题

质数和约数问题

试除法

试除法用于判断一个数是否为质数

#include<iostream>
using namespace std;
bool f(int n);
bool f(int m)
{
    for(int i=2;i<m;i++)
    {
        if(m%i==0)
            return false;
    }
    return true;
}
int main()
{
    int n;
    cin>>n;
    if(f(n))
        cout<<"Yes"<<endl;
    else
        cout<<"No"<<endl;
    return 0;

}

埃氏筛法

筛取n内的全部质数

#include<iostream>
using namespace std;
const int N=1e6+10;//随便定义了范围，看题目来
int prime[N],cnt;
bool st[N];//判断是否为质数，false为质数，true为合数；
int main()
{
    int n;
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])//当bool数组定义为全局变量默认定义为false
        {
            prime[cnt++]=i;
            for(int j=i+i;j<=n;j+=i)
                st[j]=true;//如果i为质数，则它的倍数一定是合数；
        }
    }
    cout<<cnt<<endl;
    return 0;
}

但是回看这个代码，6就会被2和3重复筛；
当数值越大，越多的数被重复筛，筛的次数也会越多，因此可以继续优化

线性筛法（可将时间复杂度降到O(n))

由于上述筛法会将部分数字重复筛取；
采用线性筛法通过质因数达到每个数只会筛取一次

#include<iostream>
using namespace std;
const int N=1e6+10;//随便定义了范围，看题目来
int prime[N],cnt;
bool st[N];//判断是否为质数，false为质数，true为合数；
int main()
{
    int n;
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])//当bool数组定义为全局变量默认定义为false
            prime[cnt++]=i;
        for(int j=0;prime[j]<=n/i;j++)
        {
            st[i*prime[j]]=true;
            if(i%prime[j]==0)
                break;
        }
    }
    cout<<cnt<<endl;
    return 0;
}

可以自己模拟一下，理解一下这个方法

分解质因数

可将n分解成这样
n=p1^a1+p2^a2+p3^a3+.....+pm^am;
例如：

给定n个正整数ai，将每个数分解质因数， 并按照质因数从小到大的顺序输出每个质因数的底数和指数。

#include<iostream>
using namespace std;
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n;
        cin>>n;
        for(int i=2;i<=n/i;i++)
        {
            int cnt=0;
            if(n%i==0)
            {
                while(n%i==0)
                {
                    cnt++;
                    n/=i;
                }
                cout<<i<<" "<<cnt<<endl;
            }
        }
        if(n!=1) cout<<n<<" "<<1<<endl;
        cout<<endl;
    }
    return 0;
}

试除法求约数

这个思路很简单，就是一个数i能被n整除
，那么n/i也是n的约数

#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
int main()
{
    int n;
    cin>>n;
    vector<int> s;
    //当n过大时，可利用vector函数防止数组空间浪费
    for(int i=1;i<=n/i;i++)
    {
        if(n%i==0)
        {
            s.push_back(i);
            if(i!=n/i)
            s.push_back(n/i);
        }
    }
    sort(s.begin(),s.end());
    for(int i=0;i<s.size();i++)
    cout<<s[i]<<" ";
    cout<<endl;
    return 0;
}

求约数个数

给定 n 个正整数 ai ，请你输出这些数的乘积的约数个数，答案对 109+7 取模。

输入格式 第一行包含整数 n 。

接下来 n 行，每行包含一个整数 ai 。

输出格式 输出一个整数，表示所给正整数的乘积的约数个数，答案需对 109+7 取模。

数据范围
1≤n≤100 ,
1≤ai≤2×109

	n=p1^a1+p2^a2+p3^a3+.....+pm^am;
	(p为质因数，a为指数）
由于有一个公式n的正约数个数=（a1+1)*(a2+1)*.....*(am+1)

#include<iostream>
#include<unordered_map>
using namespace std;
unordered_map<int,int> h;//第一个数存底数，第二个数存指数
const int mod=1e9+7;
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n;
        cin>>n;
        for(int i=2;i<=n/i;i++)
        {
            if(n%i==0)
            {
                while(n%i==0)
                {
                    h[i]++;
                    n/=i;
                }
            }
        }
        if(n!=1)
        h[n]++;
    }
    unordered_map<int,int>::iterator it=h.begin();
   // 利用迭代器来遍历h
    long long ans=1;
    while(it!=h.end())
    {
        int a=it->first,b=it->second;
        ans=ans*(b+1)%mod;
        it++;
    }
    cout<<ans;
    return 0;
}

开始不知道这个公式觉得挺难弄；
蓝桥杯十一届国赛c++组出过一个关于约数个数的题

约数和

约数和的公式是
f(n)=(p1^0+p1^1+....+p1^a1)*(p2^0+p2^1+....+p2^a2)*.....(pn^0+...pn^an)
因此简化代码可写成
int ans=1;

int b=(p1,p2…pn),a=(a1,a2…an);
//就是前面分解质因数的模板，提供思路
for(int i=0;i<a;i++)
{
ans=ans*(b+1);
b++;
}

最大公约数

最小公倍数

欧拉函数