AtCoder Grand Contest 003 D - Anticube 贪心

题意

给出n个数a[1..n]，要求选出最多的数使得其中任意两个数的乘积不是完全立方数。
n<=100000,a[i]<=10^10

分析

我们把一个数质因数分解后的每个指数都模3，不难发现每一个数都唯一对应一个和它相乘为完全立方数的数，设为其补数。然后在每个数和其补数中选一个较大的加入答案即可。
现在问题在于怎么分解把每个数化简并得到其补数。
不难注意到一个数的立方因子一定在$10^{\frac{10}{3}}$以内，那么就可以用这里的素数去将一个数化简。
求补数的话，可以用$10^{\frac{10}{3}}$以内的素数先去分解x，然后x最多只会剩下两个素因子，只要判断一下x是否是完全平方数即可。
因为时限较大所以不筛素数也可以跑过去。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<map>
#include<cmath>
using namespace std;

typedef long long LL;

const int N=100005;

int n;
LL a[N],b[N];
map<LL,int> ma;

LL sqr(LL x)
{
    return x*x;
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        LL x;scanf("%lld",&x);
        for (LL j=2;j*j*j<=x;j++)
            while (x%(j*j*j)==0) x/=j*j*j;
        ma[x]++;a[i]=x;
        LL y=1;
        for (LL j=2;j*j*j<=x;j++)
            if (x%j==0)
            {
                y*=(x%(j*j)==0)?j:j*j;
                while (x%j==0) x/=j;
            }
        if (sqr((LL)sqrt(x))==x) y*=(LL)sqrt(x);
        else y*=x*x;
        b[i]=y;
    }
    int ans=0;
    if (ma[1]) ans++,ma[1]=0;
    for (int i=1;i<=n;i++) ans+=max(ma[a[i]],ma[b[i]]),ma[a[i]]=ma[b[i]]=0;
    printf("%d",ans);
    return 0;
}