AtCoder 2004 Anticube

题面

题意

给出n个数，求其中最多选择几个数使选择的数中任意选两个数，它们的乘积都不是立方数。

做法

首先思路很好想，将每个数根据质因数的个数模3分组，也就是除以它因数中的立方数使它的因数中不含有立方数，然后再顺便算出与它相乘后得到立方数的最小数，然后存入map，对于每个数（1要特判）在它和与它相乘得立方数的最小数之间选择出现次数较多的即可。
难点在于分解质因数，如果暴枚质数，暴力分解，可以发现100000内有9000多个质数，复杂度为O（n*9000），稳T。因此要优化一下，仅暴枚到约2300以内的质数即可，因为剩下的数必定是1或者是一个大于2300的质数或是它的平方，然后判断一下这个数是不是完全平方数就可以快速分解了。

代码

#include<iostream>
#include<cstdio>
#include<map>
#include<cmath>
#define ll long long
#define N 100100
#define MZ 333
using namespace std;

ll n,a[N],b[N],zs[N],zz,ans;
bool P[N];
map<ll,ll>cnt;
map<ll,bool>vis;

inline void get()
{
    ll i,j;
    for(i=2; i<=100000; i++)
    {
        if(!P[i]) zs[++zz]=i;
        for(j=1; j<=zz&&zs[j]*i<=100000; j++) P[zs[j]*i]=1;
    }
}

int main()
{
    ll i,j,k,p,q,t,tmp;
    get();
    cin>>n;
    for(i=1; i<=n; i++)
    {
        scanf("%lld",&a[i]);
    }
    for(i=1; i<=n; i++)
    {
        t=a[i];
        b[i]=q=1;
        for(j=1; j<=MZ&&t!=1; j++)
        {
            p=0;
            for(; t%zs[j]==0; t/=zs[j]) p++;
            p%=3;
            if(!p) continue;
            for(k=1; k<=p; k++) q*=zs[j];
            for(k=p+1; k<=3; k++) b[i]*=zs[j];
        }
        if(t)
        {
            tmp=sqrt(t);
            if(tmp*tmp==t)
            {
                q*=t;
                b[i]*=tmp;
            }
            else
            {
                q*=t;
                b[i]*=t*t;
            }
        }
        a[i]=q;
        cnt[a[i]]++;
    }
    for(i=1; i<=n; i++)
    {
        if(vis[a[i]] || vis[b[i]] || a[i]==1) continue;
        vis[a[i]]=1;
        ans+=max(cnt[a[i]],cnt[b[i]]);
    }
    cout<<ans+(!!cnt[1]);
}