Agc003_D AntiCube
题目大意
给定$N$个数,求一个最大的子集,使得任意两两的乘积不是一个完全立方数。
$n\leq 10^5 A_i\leq 10^{10}$
题解
考虑两两乘积为$x^3$,由于$x^3\leq 10^{20}$,那么$x\leq 10^{\frac{20}{3}}$,那么$x$最多出现一个超过$10^{\frac{10}{3}}$的质因子,且这一质因子不会超过$10^{\frac{20}{3}}$。又由于$A_i\leq 10^{10}$,所以至多出现一个$> 10^5$的质因子,由于$x^3$由两个$A_i$相乘,那么$x$一定不会出现超过$>10^5$的质因子。
所以$x$由$[2,10^5]$内的质数幂构成,且$[10^{\frac{10}{3}},10^5]$内质数至多出现$1$个,且指数不超过$1$。
那么考虑筛出$10^5$以内的质数。
对于每一个$A_i$我们直接暴力枚举$<10^{\frac{10}{3}}$的质数(大约$300$个左右),直接筛掉,过滤掉质数是$3$的倍数的质因子,记录它由哪几个剩下的质因子构成,并将它每一个质因子质数模$3$的余数($1$或$2$)压进一个二进制状态中。对于剩下的数,如果它不为$1$,我们判断它是不是一个$<10^5$的质数或者一个$<10^5$的质数的平方,如果不是,那么它一定无法和别的任何一个数组成立方数,因为它肯定是一个$>10^5$的质数。否则,我们直接把这个稍大的质数当做之前过滤的数一样压入状态中。
对于两个过滤完的质因子的集合完全相同的数,它们能组成立方数当且仅当它们的二进制状态是恰好相反的,这样就能凑出每一个质因子的质数都是$3$的倍数。
然后就可以直接按照质因子集合放在一起,记录每一种状态和它的补集包含的数的数量,选集合大小更大的那一堆即可。
复杂度$O(n\log n+300n)$。
#include<algorithm> #include<iostream> #include<cstring> #include<cstdio> #include<cmath> #define LL long long #define M 300020 #define MAXN 100000 using namespace std; int n,m,sz[M],K[M],res,mp[200002]; int od[M],ans; LL p[M],pri[M],t[M]; bool isp[M]; bool cmp(int x,int y){return t[x]<t[y];} LL read(){ LL nm=0; char cw=getchar(); for(;!isdigit(cw);cw=getchar()); for(;isdigit(cw);cw=getchar()) nm=nm*10+(cw-'0'); return nm; } int main(){ memset(isp,true,sizeof(isp)); for(int i=2;i<MAXN;i++){ if(!isp[i]) continue; pri[++m]=i; for(int j=(i<<1);j<MAXN;j+=i) isp[j]=false; } n=read(); for(int j,i=1;i<=n;i++){ for(p[i]=read(),K[i]=sz[i]=0,j=t[i]=1;j<=308;j++){ if(p[i]%pri[j]) continue; int cnt=0; while(!(p[i]%pri[j])) p[i]/=pri[j],++cnt; cnt%=3; if(cnt) K[i]|=((cnt&1)<<sz[i]),sz[i]++,t[i]*=pri[j]; } if(p[i]>1ll){ if(p[i]<MAXN) K[i]|=(1<<sz[i]),sz[i]++,t[i]*=p[i]; else { LL qq=sqrt(p[i]); if(qq*qq!=p[i]) ans++,i--,n--; else sz[i]++,t[i]*=qq; } } } for(int i=1;i<=n;i++) od[i]=i; sort(od+1,od+n+1,cmp); for(int l=1,r;l<=n;l=r+1){ for(r=l;r<n&&t[od[r]]==t[od[r+1]];r++); if(t[od[l]]==1){ans++;continue;} int maxn=(1<<sz[od[l]])-1; for(int i=l;i<=r;i++) mp[K[od[i]]]++; for(int i=l;i<=r;i++){ int x=K[od[i]],ot=(maxn^K[od[i]]); if(mp[x]>=mp[ot]) ans+=mp[x],mp[x]=mp[ot]=0; } } printf("%d\n",ans); return 0; }