sdut3257——Cube Number(打表)

Cube Number

Time Limit: 2000MS Memory Limit: 65536KB

Submit Statistic Discuss

Problem Description

In mathematics, a cube number is an integer that is the cube of an integer. In other words, it is the product of some integer with itself twice. For example, 27 is a cube number, since it can be written as 3 * 3 * 3. 

Given an array of distinct integers (a1, a2, ..., an), you need to find the number of pairs (ai, aj) that satisfy (ai * aj) is a cube number.

Input

The first line of the input contains an integer T (1 ≤ T ≤ 20) which means the number of test cases. 

Then T lines follow, each line starts with a number N (1 ≤ N ≤ 100000), then N integers followed (all the integers are between 1 and 1000000).

Output

For each test case, you should output the answer of each case.

Example Input

1   5   1 2 3 4 9

Example Output

2

题意：给定n个数,求两两组合成立方数的情况数。

平方数升级版，分解质因数然后配对成立方数，求配对数就可以了

#include <iostream>#include <cstdio>#include <cstring>#include <queue>#include <cmath>#include <algorithm>#include <vector>#include <map>#include <string>#include <stack>using namespace std;typedef long long lsdutsdutl;#define PI 3.1415926535897932#define E 2.718281828459045#define INF 0xffffffff//0x3f3f3f3f#define mod 1000000007const int M=1005;int n,m;int cnt;int sx,sy,sz;int mp[1000][1000];int pa[M*10],rankk[M];int head[M*6],vis[M*100];int dis[M*100];int prime[M*1000];bool isprime[M*1000];int lowcost[M],closet[M];char st1[5050],st2[5050];int len[M*6];typedef pair<int ,int> ac;int dp[55][55][55][55];int has[1050000];int month[13]= {0,31,59,90,120,151,181,212,243,273,304,334,0};int dir[8][2]= {{0,1},{0,-1},{-1,0},{1,0},{1,1},{1,-1},{-1,1},{-1,-1}};void getpri(){    int i;    int j;    cnt=0;    memset(isprime,false,sizeof(isprime));    for(i=2; i<1005; i++)    {        if(!isprime[i])prime[cnt++]=i;        for(j=0; j<cnt&&prime[j]*i<10000; j++)        {            isprime[i*prime[j]]=1;            if(i%prime[j]==0)break;        }    }}struct node{    int v,w;    node(int vv,int ww)    {        v=vv;        w=ww;    }};vector<int> g[100005];char str[10005];int bit[50];int main(){    int i,j,k,t;    int u,v;    int len;    int a;    bool flag;    getpri();    scanf("%d",&t);    while(t--)    {        int ans=0;        memset(has,0,sizeof(has));        scanf("%d",&n);        for(i=0; i<n; i++)        {            scanf("%d",&a);            flag=true;   //判断能否配对成立方数            ll self=1;  //数字a多余的因数            ll wait=1;  //数字a需要配对成立方数的因数            for(j=0; j<cnt; j++)            {                if(a%prime[j]==0)                {                    k=0;                    while(a%prime[j]==0)                    {                        a/=prime[j];                        k++;                    }                    if(k%3==1)                    {                        self*=prime[j];                        wait*=prime[j]*prime[j];                    }                    else if(k%3==2)                    {                        self*=prime[j]*prime[j];                        wait*=prime[j];                    }                    if(wait>1000000||self>1000000) //如果多余因数或者需要配对的因数过大，则不可能组成立方数                    {                        flag=false;                        break;                    }                    if(a==1)                    {                        flag=true;    //被整除                        break;                    }                }            }            if(a!=1)flag=false;//防止出现整除不了质因数的特例            if(flag)//&&a==1)//可以配对(a必须等于1,防止出现整除不了质因数的特例)            {                ans+=has[self];//统计值为余数(非立方数)的数字的个数                has[wait]++;//相当于给上一句

                //这里好难解释，比如3 9 两个数,3的self是3，wait是9，has[9]++；9的self是9,正好是3的wait,这不就配对了么，ans+=has[9]=1，而9 wait是3，has[3]++，等它的盖世英雄驾着七彩祥云来找他配对~

            }        }        printf("%d\n",ans);    }    return 0;}