算法学习——数论

最新推荐文章于 2024-05-04 21:37:16 发布

lies@

最新推荐文章于 2024-05-04 21:37:16 发布

阅读量159

点赞数

分类专栏：算法学习文章标签：算法

本文链接：https://blog.csdn.net/ironman321/article/details/108171038

版权

本文介绍了数论在算法学习中的应用，包括等差数列的求和问题、因子链的计算、寻找特定数字的因子序列、等比数列的比例推算、C语言循环语句分析、正则表达式的处理以及糖果分配问题的解决策略。通过实例分析和代码展示，阐述了数论在各种算法问题中的核心作用。

摘要由CSDN通过智能技术生成

1. 等差数列

题目

数学老师给小明出了一道等差数列求和的题目。

但是粗心的小明忘记了一部分的数列，只记得其中 N 个整数。

现在给出这 N 个整数，小明想知道包含这 N 个整数的最短的等差数列有几项？

输入格式

输入的第一行包含一个整数 N。

第二行包含 N 个整数 A1,A2,⋅⋅⋅,AN。(注意 A1∼AN 并不一定是按等差数
列中的顺序给出)

输出格式

输出一个整数表示答案。

数据范围

2≤N≤100000,
0≤Ai≤109

输入样例：

5
2 6 4 10 20

输出样例：

样例解释

包含 2、6、4、10、20 的最短的等差数列是 2、4、6、8、10、12、14、16、18、20。

分析

求最小公倍数的代码：

int gcd(int a,int b)
{
   
    return b ? gcd(b,a % b) : a;
}

代码

#include<iostream>
#include<algorithm>
using namespace std;

const int N = 100010;

int n;
int a[N];

int gcd(int a,int b)
{
   
    return b ? gcd(b,a % b) : a;
}

int main()
{
   
    cin >> n;
    for(int i = 0;i  < n;i++) cin >> a[i];
    
    sort(a,a + n);
    
    int d = 0;
    for(int i = 0;i < n;i++) d = gcd(d,a[i] - a[0]);
    
    if(!d) //表示数字全一样
        cout << n << endl;
    else
        cout << (a[n - 1] - a[0]) / d + 1 << endl;
    return 0;
}

2.X的因子链

题目

输入正整数 X，求 X 的大于 1 的因子组成的满足任意前一项都能整除后一项的严格递增序列的最大长度，以及满足最大长度的序列的个数。

输入格式

输入包含多组数据，每组数据占一行，包含一个正整数表示 X。

输出格式

对于每组数据，输出序列的最大长度以及满足最大长度的序列的个数。

每个结果占一行。

数据范围

1≤X≤220

输入样例：

2
3
4
10
100

输出样例：

1 1
1 1
2 1
2 2
4 6

分析

在这里插入图片描述

算数基本定理：
在这里插入图片描述

最大长度

数列的个数

多重集合的排列数问题。

分解质因数的方法：筛法求素数（线性筛法）

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;

typedef long long LL;

const int N = (1 << 20 ) + 10;

int primes[N];  //所有的质因子
bool st[N];
int cnt;
int minp[N];  //最小质因子

void get_primes(int n)
{
   
    for(int i = 2;i <= n;i++)
    {
   
        if(!st[i])
        {
   
            minp[i] = i;
            primes[cnt++] = i;
        }
        
        for(int j = 0;primes[j] * i <= n;j++)
        {
   
            st[primes[j] * i] = true;
            minp[primes[j] * i] = primes[j];
            if(i % primes[j] == 0) break;
        }
    }
}

int main()
{
   
    get_primes(N - 1);
    int sum[N];  //sum存放的是一串 a 
    int x;
    while(scanf("%d",&x) != -1)
    {
   
        int k = 0,tot = 0;
        while(x > 1)
        {
   
            int p = minp[x];  //p是x的最小质因子
            sum[k] = 0;
            while(x % p == 0)
            {
   
                x /= p;
                sum[k]++;
                tot++;
            }
            k++;
        }
    
    LL res = 1;
    for(int i = 1;i <=tot;i++) res *