IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> C++知识库 -> 数学问题(数论)试除法做质数的判断、分解质因数,筛质数 -> 正文阅读

[C++知识库]数学问题(数论)试除法做质数的判断、分解质因数,筛质数

质数

在大于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;
}

  C++知识库 最新文章
【C++】友元、嵌套类、异常、RTTI、类型转换
通讯录的思路与实现(C语言)
C++PrimerPlus 第七章 函数-C++的编程模块(
Problem C: 算法9-9~9-12:平衡二叉树的基本
MSVC C++ UTF-8编程
C++进阶 多态原理
简单string类c++实现
我的年度总结
【C语言】以深厚地基筑伟岸高楼-基础篇(六
c语言常见错误合集
上一篇文章      下一篇文章      查看所有文章
加:2022-07-03 10:32:27  更:2022-07-03 10:34:58 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年5日历 -2024/5/11 17:15:15-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码