【题单 - 数学专题】最大公约数

最新推荐文章于 2024-06-09 11:24:10 发布

繁凡さん

最新推荐文章于 2024-06-09 11:24:10 发布

阅读量689

点赞数 2

CC 4.0 BY-SA版权

分类专栏： # 第五章数学知识数学 - 约数

本文链接：https://blog.csdn.net/weixin_45697774/article/details/110956197

这篇博客整理了关于最大公约数的算法模板，包括基于值域预处理的快速GCD算法，以及在不同题目中应用GCD解决问题的方法，如缘分问题、最大公约数和最小公倍数问题等。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

整理的算法模板合集： ACM模板

Part 6.2.2 最大公约数

如果两个数有一个共同的约数，那么这个约数就被称为公约数。最大公约数就是指这两个数的所有公约数中，最大的一个。

求解两个数的最大公约数，可以采用欧几里得算法解决。

P5435 基于值域预处理的快速 GCD算法

给定 $n$ 个正整数 $a_1,a_2,\dots,a_n$ ，再给定 $n$ 个正整数 $b_1,b_2,\dots,b_n$ ，你需要对每对 $(i, j)$ 求出 $a_i$ 与 $b_j$ 的最大公因数。

不难发现你的输出应有 $n^2$ 个正整数。为了减少输出对程序的运行效率的影响，你只需要输出 nn 行，每行一个整数 $A_i$ 。

其中对于 $i\in[1,n]i∈[1,n]，A_i=\sum_{j=1}^{n}i^j\gcd(a_i,b_j)A$ 。由于答案可能过大，你只需要输出模 $998, 244, 353$ 后的结果即可。

还是直接上模板。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstring>

using namespace std;
typedef long long ll;

const int N = 5000, M = 1000000, T = 1000, Mod = 998244353;

int pre[T + 2][T + 2];
int a[N + 2], b[N + 2];
int fac[M + 2][3];
bool isp[M + 2];
int pri[M / 10], tot;
int n;

inline int read()
{
   
   
	int re=0,k=1;char ch=getchar();
	while(ch>'9'||ch<'0'){
   
   if(ch=='-')k=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){
   
   re=re*10+ch-48;ch=getchar();}<