Codeforces Round #825 (Div. 2)B. Playing with GCD

最新推荐文章于 2022-10-20 16:38:06 发布

2,4(1H,3H)-PD are mine

最新推荐文章于 2022-10-20 16:38:06 发布

阅读量550

点赞数

分类专栏： codeforces刷题文章标签： codeforces 数学最大公因数

本文链接：https://blog.csdn.net/qq_60256199/article/details/127269849

版权

codeforces刷题专栏收录该内容

18 篇文章 1 订阅

订阅专栏

题目链接：Problem - B - Codeforces

题目大意：

t(1≤t≤10e5)组测试，每一组给定长度为n(1≤n≤10e5)的数组a，问能否构造出一个数组b，使得ai=gcd( bi,b(i+1) )(1≤i≤n)

gcd为最大公因数。

每组输出YES或NO即可。

思路：

从条件出发，考虑相邻的两个ai和a(i+1)：

$a_{i}=(b_{i},b_{i+1})\rightarrow a_{i}|b_{i}, a_{i}|b_{i+1}$

$a_{i+1}=(b_{i+1},b_{i+2})\rightarrow a_{i+1}|b_{i+1}, a_{i+1}|b_{i+2}$

由上式可以看出b(i+1)同时被ai和a(i+1)整除。

故对于bi需满足下式：

$a_{i}|b_{i}; a_{i-1}|b_{i}$

考虑构造bi=[ai,a(i-1)](最小公倍数)，构造出bi后逐一检查ai=gcd( bi,b(i+1) )是否成立。

如果都成立，说明可以找到b数组直接输出YES。

如果ai不等于gcd( bi,b(i+1) )，由于：

$b_{i}=[a_{i-1},a_{i}]$

$b_{i+1}=[a_{i},a_{i+1}]$

所以一定有gcd( bi,b(i+1) )=k*ai (k为整数)

考虑重新构造一个c数组，使得gcd(ci,c(i+1))=ai。在之前的构造中，我们考虑的是最小的bi满足bi=[ai,a(i-1)]，现在ci仍需要满足：

$a_{i}|c_{i}; a_{i-1}|c_{i}$

故有 $c_{i}>=b_{i}\rightarrow (c_{i},c_{i+1})>=(b_{i},b_{i+1})$

( bi,b(i+1) )已经比ai大了，这样的ci是不可能构造出来的。

所以如果有一个不成立，直接输出NO。

边界情况a[0]=a[n+1]=1即可。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cstdlib>
#include <cmath>
#include <stack>
#include <queue>
#include <vector>
#include <map>
#include <set>
#include <algorithm>
#define pii pair<int,int>
#define pll pair<LL,LL>
#define pil pair<int,LL>
#define pli pair<LL,int>
#define se second 
#define fi first
#define endl '\n'
#define rep(i,a,b) for (register int i=a;i<b;++i)
#define per(i,a,b) for (register int i=a;i>b;--i)
typedef long long LL;
using namespace std;
const LL MOD=998244353;
const int N=1e5+10;
int a[N],b[N];
int lcm(int x,int y)
{
	return x/__gcd(x,y)*y;
}
void solve()
{
	int n;
	cin>>n;
	rep(i,1,n+1) cin>>a[i];
	a[0]=a[n+1]=1;
	rep(i,1,n+2) b[i]=lcm(a[i-1],a[i]);
	rep(i,1,n+1)
	if(a[i]!=__gcd(b[i],b[i+1])){
		cout<<"NO"<<endl;
		return;
	}
	cout<<"YES"<<endl;
}
int main()
{
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	int _=1;
	cin>>_;
	while(_--){
		solve();
	}
	return 0;
}