Codeforces1349A Orac and LCM (数论)

题目链接:Orac and LCM

大致题意

给定长度为n的序列a, 把所有满足i < j的位置求最大公约数lcm(ai, aj), 然后放入一个集合中.

最终求集合中所有元素的最大公约数.

解题思路

首先序列长度会有1E5, 一个个暴力求取是不可行的. 不如从每个数字贡献的角度去思考.

我们现在考虑a1: 对于第一个位置的数字, 他会与后面的所有的数求lcm, 最终这些数字会参与到最后的gcd运算中.

接下来的思路以唯一分解定理所得到的质数角度去考虑, 两个数求lcm, 相当于相互补全不足的质因子(求并集), 两个数求gcd, 相当于求两个数质因子的交集.

那么对于a1这个位置, 需要与后面的所有数字进行一次并集运算, 把得到的结果, 再进行交集运算. 最终我们发现交集中, a1原有的东西一定会留下, 那么a2~an之间共有的东西也会留下.

由此我们发现, 我们可以通过先求a2~an之间共有的部分(求交集), 然后再与a1合并(求并集), 这样得到的结果是与上述相同的. 那么对于a2位置, 我们可以求a3~an之间的交集… 以此类推, 我们发现只需要维护一个后缀gcd即可.

由此我们也可以得出公式: gcd(lcm(a1,a2),lcm(a1,a3),…lcm(a1,ak))=lcm(a1,gcd(a2,a3,a4,…,ak))

AC代码

#include <bits/stdc++.h>
#define rep(i, n) for (int i = 1; i <= (n); ++i)
using namespace std;
typedef long long ll;
const int N = 1E5 + 10;
ll a[N], suf[N];
int main()
{
    int n; cin >> n;
	rep(i, n) scanf("%lld", &a[i]);

	for (int i = n; i >= 1; --i) suf[i] = gcd(a[i], suf[i + 1]);

	ll res = 0;
	rep(i, n - 1) {
		ll now = lcm(a[i], max(1ll, suf[i + 1]));
		res = gcd(res, now);
	}

	cout << res << endl;

    return 0;
}

END