Description

众所周知,lala是一名非常喜欢看天线宝宝拉拉吃面包的高中生,在经过无数次的重播那堪称经典的画面之后,lala终于感到看不下去了,决定自己也去吃回面包。
于是lala来到了一家面包店,但由于这家面包店的面包实在太多了,lala难以选择该买哪些面包,但这自然难不倒机智的lala,他先给每种面包定了一个美味度,并只购买美味度能表示成x^p(x,p均为大于等于2的正整数)的面包,但即便如此,lla还是难以快速地知道该买哪些面包,会买多少面包,你能帮帮他吗?

Solution

方法

输入一个k,题目要求判断满足一个xp=k
转换一下logxk=p,那么就是x开p次方等于k=x1p,我们发现k其实并不大,可以直接枚举。
然后判断一下就好了。

小优化

因为次数p是可以因式分解的,所有只用枚举质数为次数就好了。

注意

打c++如果用pow可能会有精度问题,我们可以在pow外面打一个round四舍五入。
注意如果x=1那么就不用管它了。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fod(i,a,b) for(i=a;i>=b;i--)
typedef long long ll;
using namespace std;
const int maxn=1000007;
const ll da=100000000000000;
ll i,j,k,l,t,n,m,ans,tot,o,p;
int zhi[100];
bool bz,az[maxn];
ll a[50000];
int main(){
freopen("bread.in","r",stdin);
freopen("bread.out","w",stdout);
fo(i,2,50){
if(!az[i]){
zhi[++tot]=i;
fo(j,2,50/i){
az[i*j]=1;
}
}
}
scanf("%lld",&n);
fo(i,1,n){
scanf("%lld",&k);
if(k==1)continue;
fo(j,1,tot){
double u=round(pow(k,1.0/zhi[j]));
if(round(pow(1.0*u,1.0*zhi[j]))==k){
++ans;
break;
}
}
}
printf("%lld\n",ans);
}