数学问题（数论）试除法做质数的判断、分解质因数，筛质数

质数

在大于1的整数中，只有1和本身两个约数，则称为素数（质数）。

1.质数的判断（试除法）

int prime(long long n)

{

    if(n<2) return 0;//要求大于1 的正整数

    else

    {

        for(int i=2;i<n;i++)//遍历一遍，时间复杂度妥妥的O(n)超时

        {

            if(n%i==0) return 0;

        }

    }

    return 1;

}

优化一下：

由于n/d=x,同时n/x=d也成立，所以可以得出约数是成对存在的。

由此我们只需要枚举一下较小的的那个，即d<=n/d，把循环条件改写成i<=n/i，这样时间复杂度就为O(sqrt(n))了。

题目链接：866. 试除法判定质数 - AcWing题库

 题面：

2.分解质因数（试除法）

质因数是指能整除给定整除的质数因子。

【若两个数只有1一个共同的质因子，称为互质】

【1和任何一个正整数都是互为质数】

【任何一个数都可以看成多个质因子相乘，n=p1^k1*p2^k2*…*pn^kn】

链接：AcWing 867. 分解质因数 - AcWing

题面：

 代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    ll t;
    cin>>t;
    while(t--)
    {
        ll n;
        cin>>n;
        for(int i=2;i<=n;i++)//时间复杂度是O(n),超时咯
        {
            if(n%i==0)//看似是找出来的是n的因数，而非质因数。
/*但其实这里没问题，由于下面一旦找到了能被整除的i，n会进入循环一直除i知道除不了为止。拿4来举例子，在i=2的时候就n/=2 ——>n=2，又n/=2——>n=0.没有i=4什么事了，所以把因数为合数的情况全部除尽了，保证都是质因数。*/
            {
                int s=0;
                while(n%i==0)//找出指数为多少
                {
                    n/=i;
                    s++;
                }
                cout<<i<<" "<<s<<endl;
            }
        }
        cout<<endl;//输出有每个数据要换行的要求
    }
    return 0;
}

优化一下：

由于n中最多只有一个大于sqrt(n)的质因子

因此循环条件换成i<=n/i，同时加上一个特判if(n>1) cout<<n<<" "<<1<<endl;

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    ll t;
    cin>>t;
    while(t--)
    {
        ll n;
        cin>>n;
        for(int i=2;i<=n/i;i++)//时间复杂度为O(sqrt(n))，最多为sqrt(n)
        {
            if(n%i==0)
            {
                int s=0;
                while(n%i==0)
                {
                    n/=i;
                    s++;
                }
                cout<<i<<" "<<s<<endl;
            }
        }
        if(n>1) cout<<n<<" "<<1<<endl;
        cout<<endl;
    }
    return 0;
}

3.筛质数

题目链接：868. 筛质数 - AcWing题库

题面：

代码：

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

const int N=1e6+10;

int primes[N],cnt,n;

bool st[N];

int main()

{

    cin>>n;

    for(int i=2;i<=n;i++)

    {

        if(!st[i])//依次从小到大质数遍历，最小的为2

            primes[cnt++]=i;

        for(int j=i+i;j<=n;j+=i)//把i的倍数都删掉

            st[j]=true;

    }

    cout<<cnt<<endl;

    return 0;

}

时间复杂度计算：

n/2+n/3+n/4+…+n/n

=n(1/2+1/3+1/4+…+1/n)

当n趋于正无穷时，后面为调和级数，值为ln n+c （c为0.57）

又因为ln n<log2 n

所以时间复杂度为：O（nlog2 n）

优化一下：把第二个循环放到if里面去，就是埃氏筛法

具体的埃氏筛法传送门：

埃拉托色尼筛选法巧解质数问题（埃氏筛法求解素数问题）_吴雨4的博客-CSDN博客

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int primes[N],cnt,n;
bool st[N];
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])//依次从小到大质数遍历，最小的为2
        {
            primes[cnt++]=i;
            for(int j=i+i;j<=n;j+=i)//把i的倍数都删掉
                st[j]=true;
        }
    }
    cout<<cnt<<endl;
    return 0;
}

 两次都能AC，但是第二个筛法比第一个快了五倍。

 还有一种线性筛法：（数据范围为1e7时比埃氏筛法又快了一倍）

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int primes[N],cnt,n;
bool st[N];
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!st[i])//依次从小到大质数遍历，最小的为2
            primes[cnt++]=i;
            for(int j=0;primes[j]<=n/i;j++)
            {
                st[primes[j]*i]=true;
                if(i%primes[j]==0) break;
            }
    }
    cout<<cnt<<endl;
    return 0;
}