题意:
给你一个数列,数列中互相之间满足 l c m ( x , y ) / g c d ( x , y ) lcm(x,y)/gcd(x,y) lcm(x,y)/gcd(x,y)开方后是整数的可以构成一个多重集合,每过一秒数列中的所有数都会变为他们的乘积,问第w秒的时候,数列中最大的多重集合中的元素个数。
思路:
通过观察式子 l c m ( x , y ) / g c d ( x , y ) lcm(x,y)/gcd(x,y) lcm(x,y)/gcd(x,y)可以知道同属一个集合内的元素需要满足何种性质。对每个数分解质因数,则 x = p 1 k 1 ∗ p 2 k 2 ∗ ∗ ∗ p m K m x=p_1^{k_1}*p_2^{k_2}***p_m^{K_m} x=p1k1∗p2k2∗∗∗pmKm,y同理也可以表示为这种形式,那么对于x和y的各个质因数来说,lcm就等于拥有最大次数的各个质因数的最大幂次的乘积,gcd同理取最小次数,不难得出,只有当lcm/gcd的各个质因数的次方数为偶数时,才能在开方后是整数,于此我们可以对每一个数分解质因数,并统计每个数的奇数次质因子的乘积,乘积相同的才能被放在同一个集合;于是,在第0秒时,数列被分成了:只含有偶数次质因子的数和奇数次质因子乘积相同的数,统计他们的最大值就是第0秒的答案;而第一秒及之后,偶数个集合的奇数次质因子乘积相同的数会变为偶数次,然后融入到之前的偶数次质因子集合中,而奇数次质因子的则会保持原样,这样,答案就出来了。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
#define pb push_back
ll a[300010];
int pri[1000010],cnt,vis[1000010];
void oula()
{
for(int i = 2; i <= 1000000; i++){
if(!vis[i])pri[++cnt] = i;
for(int j = 1; j <= cnt && pri[j]*i <= 1000000; j++){
vis[pri[j]*i] = 1;
if(i%pri[j] == 0)break;
}
}
}
void solve()
{
int n;
cin>>n;
vector<int>v(1,0);
int sum = 0,sum2 = 0;
for(int i = 1; i <= n; i++){
cin>>a[i];
ll k = 1,f = 0;
for(int j = 1; j <= cnt && pri[j] <= sqrt(a[i]); j++){
int num = 0;
while(a[i]%pri[j] == 0){
num++;
a[i]/=pri[j];
}
if(num%2)k*=pri[j],f = 1;
}
if(a[i] > 1)k*=a[i],f = 1;
if(k>1)v.pb(k);
if(!f)sum++;
}
v[0] = 0;
sort(v.begin(),v.end());
int c = 1,maxx = 0;
v.pb(-1);
for(int i = 1; i < v.size(); i++){
if(v[i]==v[i-1])c++;
else{
if(c%2 == 0)sum2+=c;
c = 1;
}
maxx = max(maxx,c);
}
int q;
cin>>q;
while(q--){
int w;
scanf("%d",&w);
if(w)cout<<max(sum+sum2,maxx)<<endl;
else cout<<max(sum,maxx)<<endl;
}
}
int main(){
oula();
int t;
cin>>t;
while(t--){
solve();
}
}