山东省第六届省赛 - Cube number （立方数）

题目

题意：

给你一组数，问你两个数乘起来是立方数的对数。

POINT：

先降立方操作。除掉立方数。

之后得到的数肯定为.  x=a1^1*a2^1*a3^2*a4^2……

a1,a2,a3都是它的质因数。

即 这个数是一些质数的平方或者本身乘起来。

所以要找到对应的数和他乘起来刚好为a1^3*a2^3*a3^3*a4^3.而且要快。

假设to就是我们要找的数。

遍历for(int i=1;pri[i]*pri[i]<=x;i++)可以找到 每次x%(pri*pri)==0就是要乘一次的（to*=pri）. 当前to=a3*a4

怎么找乘两次的。temp=x/to/to,就是a1^1*a2^1。把to=to*temp*temp，就是a1^2*a2^2*a3^1*a4^1。

就找到了。

枚举 for(int i=1;pri[i]<=x;i++)会超时哦。

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<cmath>
using namespace std;
const int maxn = 100000+55;
#define LL long long

LL num[1000100];
LL a[100010];
LL is[1000100];
LL pri[maxn];
LL cnt=0;


void init()
{
	for(LL i=2;i<=1000000;i++){
		if(is[i]==0){
			pri[++cnt]=i;
			for(LL j=2*i;j<=1000000;j+=i)
				is[j]=1;
		}
	}
}
LL ans=0;

LL f(LL x)
{
	LL to = 1;
	for(LL i=1;pri[i]*pri[i]<=x;i++){
		while(x%(pri[i]*pri[i]*pri[i])==0)
			x/=pri[i]*pri[i]*pri[i];
		if(x%(pri[i]*pri[i])==0)
			to=to*pri[i];
	}
	LL temp=x/to/to;
	to=to*temp*temp;
	if(to!=1&&to<=1000000)
		ans+=num[to];
	return x;
}


int main()
{
	init();
	LL T;
	scanf("%lld",&T);
	while(T--){
		LL n;
		scanf("%lld",&n);
		memset(num,0,sizeof num);
		ans=0;
		for(LL i=1;i<=n;i++){
			scanf("%lld",&a[i]);
			a[i]=f(a[i]);
			num[a[i]]++;
		}
		ans+=num[1]*(num[1]-1)/2;
		printf("%lld\n",ans);
	}
	return 0;
}