NOIP中的数学---第3课约数

信奥教练Andy

已于 2022-06-29 14:21:22 修改

阅读量789

点赞数 3

分类专栏： NOIP中的数学文章标签：算法线性代数 c++

于 2021-01-18 21:52:07 首次发布

本文链接：https://blog.csdn.net/qq_32431299/article/details/112791704

版权

NOIP中的数学专栏收录该内容

13 篇文章 19 订阅

订阅专栏

约数定义

思考三秒钟，注意约数是成对出现的
如果数a能被数b整除， a就叫做 b 的倍数， b就叫做a的约数。约数和倍数都表示一个整数与另一个整数的关系，不能单独存在。
最大公约数是指两个或多个整数共有约数中最大的一个。例如：24 ，18 共同的约数为 1 2 3 6 ，其中6是最大的。所以GCD(24,18)=6 。

一级题目

试除法求所有约数

给定n个正整数ai，对于每个整数ai,请你按照从小到大的顺序输出它的所有约数。
输入格式
第一行包含整数n。
接下来n行，每行包含一个整数ai。
输出格式
输出共n行，其中第 i 行输出第 i 个整数ai的所有约数。
数据范围
1≤n≤100,
2≤ai≤2∗109
输入样例：
2
6
8
输出样例：
1 2 3 6
1 2 4 8

解法
1、大暴力

#include<iostream>
using namespace std;
int main()
{
    int n,x;
    cin >> n;
    while(n --)
    {
        cin >> x;
        for(int i = 1;i <= x;i ++)
        {
            if(x % i == 0)cout << i<<" "; 
           
        }
        cout <<endl;
    }
}

Time Limit Exceeded
当n==100的时候就已经…

2、代码优化
考虑约数是成对出现的，所以我们只要求出其中那个小的i,而和它一对的约数n/i也可以得到了，但要注意两点：
（1）题目要求约数从小到大排序，所以要建立数组暂存一下约数，然后排序输出
（2）注意输出n/i的时候要考虑i的平方的情况，也就是i==n/i的时候就不用重复记录约数了
上“酸菜”

#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
int a[10000];
int main()
{
    int n,x;
    cin >> n;
   
    while(n --)
    {
        cin >> x;
        int j = 0;//约数数组的下标，也是约数的总个数
        memset(a,0,sizeof(a));
        for(int i = 1;i <= x / i;i ++)
        {
            if(x % i == 0)
            {
                a[j++] = i;
                if(i != x / i) a[j++] = x / i;//防止i^2的出现
            }
        }
        sort(a,a+j);
        for(int i = 0;i <j;i++)cout<<a[i]<<" ";
        cout <<endl;
    }
}

二级题目：约数个数

公式：对于任意一个整数N：
N的分解质因数
N=（p1^a1）*（p2^a2)*……*(pn^an)其p1 …pn为质数，括号纯为了好看
则约数的个数是(a1+1)(a2+1)……*(an+1)
例如72=2³ * 3²,则72的约数个数为（3+1）×（2+1）=12个。
看例题：约数个数

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

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

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

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

数据范围 1≤n≤100, 1≤ai≤2∗10⁹
输入样例：
3
2
6
8
输出样例：
12

题解按上面的公式，对每个乘数进行分解质因数，并利用map<int ,int>记录对应的质数的指数，然后遍历map进行指数求和即可。

#include<iostream>
#include<map>
using namespace std;
const int M = 1e9+7;
long long ans = 1;
int main()
{
    int n;
    cin >> n;
    map<int,int> mp;
    while(n--)
    {
        int x;
        cin >> x;
        //下面进行分解X的质因数形式
        for(int i = 2;i <= x/i;i++)
        {
            while(x % i ==0)
            {
                x /= i;
                mp[i] ++;
            }
        }
         if(x>1) mp[x]++;
        //接下来遍历已经处理好的map，进行累加求和即可
        
    }
    for(auto i:mp) ans = ans * (i.second +1) % M;
        cout << ans  << endl;//注意cout<<ans %M；最后取模会炸掉，思考模的性质
    
}

关于for(auto i:mp) ans = ans * (i.second +1) % M;做一个解释，这是C++11的新特性。

//迭代器版本遍历map
auto i = mp.begin();
while(i !=mp.end())
{
	cout<<i -> first<<":"<<i -> second<<endl;
	i++;
	}
另一种写法
for(auto i = mp.begin();i != mp.end();i++)
	cout<<i -> first<<":"<<i -> second<<endl;
//C++11新特性
for(auto i:mp)
	cout<<i -> first<<":"<<i -> second<<endl;

最大公约数

12 和 18 的最大公约数是多少？

最大公约数：最大公约数，也称最大公因数、最大公因子，指两个或多个整数共有约数中最大的一个。例如：18 与 12 的最大公约数为 6 。

短除法

短除法是求最大公因数的一种方法：先把每个数的因数找出来，然后再找出公因数，最后在公因数中找出最大公因数。
在这里插入图片描述

因式分解法
在这里插入图片描述
在初中数学题中，基本上我们就是采取因式分解或者短除法的形式来求最大公约数。

但是它们存在的问题是：当公共素因子较小时，通过观察可以很快找出；但是当公因子较大时，仅仅通过观察已经很难找出甚至在一定时间内找不出。

比如求 22008 和 655 的最大公约数时，很难直接找到其公因子。
那么有没有更好的方法来求解最大公约数呢？答案是有的，就是接下来要介绍的欧几里得算法。

欧几里得算法

欧几里得算法（英语：Euclidean algorithm），又称辗转相除法，是求最大公约数的算法。

辗转相除法基于如下原理：两个整数的最大公约数等于其中较小的数和两数的差的最大公约数。

辗转相除法首次出现于欧几里得的《几何原本》（第VII卷，命题i和ii）中，而在中国则可以追溯至东汉出现的《九章算术》。

辗转相除法基于如下原理：两个整数的最大公约数等于其中较小的数和两数的差的最大公约数。例如，252 和 105 的最大公约数是 21；因为 252 − 105 = 21 × (12 − 5) = 147 ，所以 147 和 105 的最大公约数也是 21。

在这个过程中，较大的数缩小了，所以继续进行同样的计算可以不断缩小这两个数直至其中一个变成零。

这时，所剩下的还没有变成零的数就是两数的最大公约数。

将上面的较大的数缩小的过程中往往使用的是 MOD 操作。

MOD，是一个数学运算符号-----求余运算符。

例如 a mod b = c，表明 a 除以 b 余数为 c 。在整数的除法中，只有能整除与不能整除两种情况，所以当不能整除时，就会产生余数。

我们借助于 MOD 使用辗转相除法的概念来求数字 1112和数字 695的最大公约数。

当余数变为 0 的时候，最后一个操作的除数是最大公约数，即 139是数字 1112和数字 695的最大公约数。
在这里插入图片描述

设计来源于算法动画讲解

一般算法流程如下：

在这里插入图片描述

动画理解

通过动画来理解一下为什么使用辗转相除法可以找到最大公约数。

将最大的公约数设置为 n，当然虽然一开始对于每个整数是不知道可以分段成多少个 n 的，但是，可以知道 1112 和 695 都是最大公约数 n 的倍数。

在这里插入图片描述

通过 mod 操作，不停的找余数。

在这里插入图片描述

最后两个数是倍数关系，可以整除，余数为 0 ，结束了操作。

此时剩下在这里插入图片描述
的一条线段的长度就是 1112 和 695 的最大公因数。
**

C++代码实现

递归实现

#include <iostream>

using namespace std;

int gcd(int a, int b){
    if (a%b == 0) {
        return b;
    }
    return gcd(b, a%b);
}

int main() {
    int a = 14, b = 18;
    printf("%d\n",gcd(a,b));
    return 0;
}

非递归实现

int gcd(int a, int b){
    int temp = a;
    while(a%b != 0){
        a = b;
        b = temp%b;
        temp = a;
    }
    return b;
}

小赠品-最小公倍数

最小公倍数，就是a b的乘积除以它们两个的最大公约数，就是它们的最小公倍数。代码如下：

int MinMultiple( int a, int b){
    return (a * b)/gcd(a, b);
}

这样子就可以了。
练习题目
1、最大公约数
2、最小公倍数
3、小b与灯泡
4、分数约分
5、最大公约数之和V2
6、最大公约数之和V3
7、最小公倍数之和V2
8、最小公倍数之和V3
9、最小公倍数计数

信奥教练Andy

关注

3
点赞
踩
3

收藏

觉得还不错? 一键收藏
打赏
0
评论
NOIP中的数学---第3课约数

约数定义思考三秒钟，注意约数是成对出现的一级题目试除法求所有约数给定n个正整数ai，对于每个整数ai,请你按照从小到大的顺序输出它的所有约数。输入格式第一行包含整数n。接下来n行，每行包含一个整数ai。输出格式输出共n行，其中第 i 行输出第 i 个整数ai的所有约数。数据范围1≤n≤100,2≤ai≤2∗109输入样例：268输出样例：1 2 3 61 2 4 8解法1、大暴力#include<iostream>using namespace
复制链接

扫一扫