大意:
对于一个长度为n的数列,定义它的GCD Table G是一张n×n的表,其中G[i][j]=gcd(a[i],a[j]),现在乱序给出G中所有数,求原数列a。
思路:
gcd(a,b)<=min(a,b),也就是说,两个数的gcd一定是不大于其中任意一个的,换句话说,在我们得到的表中,最大的数一定是原数组的一个元素。那么我们当然可以把它加入答案之中。接下来我们就得把它对其余数的影响去掉,也就是每次这样的操作过后,我们都得遍历答案数组,把答案中已有的数与该数的gcd,让其对应的cnt-2,也就是去掉了该数的影响,那么我们之后每次都取当前最大就是合理的了。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
ll n,a;
map<ll,ll> mp;
ll sd[251000];
ll as[251000];
bool cmp(ll a,ll b)
{
return a>b;
}
ll cnt=0;//出现过的数
int main()
{
cin>>n;
for(int i=1;i<=n*n;++i)
{
cin>>a;
if(!mp[a]) sd[++cnt]=a;
mp[a]++;
}
sort(sd+1,sd+1+cnt,cmp);
// for(int i=1;i<=cnt;++i) cout<<sd[i]<<' ';
// cout<<endl;
ll ans=0;
for(int i=1;i<=cnt&&ans<n;)
{
//cout<<sd[i]<<" "<<mp[sd[i]]<<endl;
while(!mp[sd[i]])
{
// cout<<sd[i]<<' '<<mp[sd[i]]<<endl;
i++;
}
as[++ans]=sd[i];
mp[sd[i]]--;//减自己,也就是gcd(sd[i],sd[i])
for(int j=1;j<ans;++j) mp[__gcd(sd[i],as[j])]-=2;
}
for(int i=1;i<=ans;++i) cout<<as[i]<<' ';
return 0;
}