题目链接: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):
由上式可以看出b(i+1)同时被ai和a(i+1)整除。
故对于bi需满足下式:
考虑构造bi=[ai,a(i-1)](最小公倍数),构造出bi后逐一检查ai=gcd( bi,b(i+1) )是否成立。
如果都成立,说明可以找到b数组直接输出YES。
如果ai不等于gcd( bi,b(i+1) ),由于:
所以一定有gcd( bi,b(i+1) )=k*ai (k为整数)
考虑重新构造一个c数组,使得gcd(ci,c(i+1))=ai。在之前的构造中,我们考虑的是最小的bi满足bi=[ai,a(i-1)],现在ci仍需要满足:
故有
( 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;
}