暑假acwing算法总结17：分解质因数与筛法

• 1、分解质因数

• 用的筛法的思想

#include<iostream>
#include<vector>

using namespace std;

typedef pair<int,int>pii;

vector<pii>p;
void divid(int n)
{
    p.clear();
    for(int i=2;i<=n/i;i++)
    {
        if(n%i==0)//保证合数一定不会被整除，存储的都是质数
        {         //用到了筛法的思想比如开始能被2除尽除完之后不是2的倍数下次一定不会被4除尽
            int s=0;
            while(n%i==0)
            {
                n/=i;
                s++;
            }
            p.push_back({i,s});
        }
    }
    if(n>1)p.push_back({n,1});
}
int main()
{
    int t;cin>>t;
    while(t--)
    {
        int x;cin>>x;
        divid(x);
        for(int i=0;i<p.size();i++)
        cout<<p[i].first<<' '<<p[i].second<<endl;
        cout<<endl;
    }
    return 0;
}
//2021-07-13/16:11

• 2、筛法

• 1、暴力筛法(朴素筛)
• 从2开始循环每次是其倍数标上删除，依次循环到n-1余下的数就全是质数

void pss(int n)
{
    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;
    }
}

• 埃氏筛：把质数的倍数全部删除就是将朴素筛的

for(int j=i+i;j<=n;j++)st[j]=true;

• 移到if判断内

#include<iostream>
#include<cmath>

using namespace std;

const int N=1e6+10;
int prime[N],cnt;
bool st[N];
void ass(int n)
{
    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;
        }
    }
}
int main()
{
    int n;cin>>n;
    bls(n);
    cout<<cnt<<endl;
    return 0;
}
//2021-07-13/16:22

• 线性筛：每一个数n只会被最小质因子筛掉，在1e6的数据范围领先较小在1e7的数据范围领先一倍
• 因为只有一个最小质因子，所以每个数只筛一次所以是线性的

void xxs(int n)
{
    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;//这句话发生prime[j]一定是i的最小质因子
        }
    }
}

• 2021-07-13写于南宁