数论-欧几里得算法-第十届蓝桥杯省赛C++B组-等差数列

数论-欧几里得算法-第十届蓝桥杯省赛C++B组-等差数列

题目：

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

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

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

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

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

输出格式
输出一个整数表示答案。

数据范围
2≤N≤100000,
0≤Ai≤10⁹
输入样例：
5
2 6 4 10 20
输出样例：
10
样例解释
包含 2、6、4、10、20 的最短的等差数列是 2、4、6、8、10、12、14、16、18、20。

题解：

$给出一个等差数列的一部分序列，要求该完整序列的最小长度，也就是要公差尽量的大。$

$就是要找所有数差之间的最大公约数，即该等差数列的最大公差。$

$首先每一项都是大于等于零的，所以公差是大于等于零的，因此将序列先按升序排序，再遍历求公差的最大公约数。$

$最终答案项数，就是最大项和最小项之差除公差再加一，即\frac{A_n-A_1}{d}+1，显然要特判d=0的情况。$

欧几里得算法说明：

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

$$\begin{gathered} 不用考虑a与b的大小关系。 \\ 因为，如果a<b，在第一次递归后将返回gcd(b,a\%b)，而因为a<b，所以a\%b=a，即gcd(b,a)。 \end{gathered}$$

---

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#define P pair<int,int>
#define x first
#define y second
#define ll long long
#define inf 0x7fffffff
using namespace std;
const int N=1e5+10;
int n,a[N];
int gcd(int a,int b)
{
    return b ? gcd(b,a%b) : a;
}

int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);

    sort(a+1,a+n+1);
    int d=0;
    for(int i=2;i<=n;i++)
        d=gcd(d,a[i]-a[i-1]);

    if(!d) cout<<n<<endl;
    else cout<<(a[n]-a[1])/d+1<<endl;

    return 0;
}