这题看起来简单,做的时候才发现不是那么简单的。。。
M^K 在K为合数的时候都可以表示为X^Y (Y为质数),但X^Y这个数可能出现多次,例如8^2= 4^3
这是因为该数为2^6, 6等于质数2乘质数3 所以我们在计数时 ans= 指数为一个质数 - 指数为两个质数 + 指数为三个质数;
因为N最大为10^18 < 2^60 ,所以指数为 三个质数的数为:2^(2*3*5)= 2^30 , 3^(2*3*5)= 3^30, 2^(2*3*7)= 2^42
代码:
#include <stdio.h>
#include <string.h>
#include <math.h>
#define LL __int64
int main()
{
int a[]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};// 一个质因子的指数
int b[]={6,10,14,15,21,22,26,33,34,35,38,39,46,51,55,57,58};//两个质因子的指数
LL n;
while(scanf("%I64d",&n)!=EOF)
{
int ans= 0;
for(int i= 0; i<= 16; i++) //指数为一个质因子
{
int num= pow(1.0*n, 1.0/a[i]);
while(pow(num+1.0, 1.0*a[i])<= n)
num++;
num--;
if(num> 0)
ans+= num;
}
for(int i= 0; i<= 16; i++) //指数为两个质因子
{
int num= pow(1.0*n, 1.0/b[i]);
while(pow(num+1.0, 1.0*b[i])<= n)
num++;
num--;
if(num> 0)
ans-= num;
}
ans++;
if(n>= pow(2.0, 30.0)) ans++;//指数为三个质因子的只有三个数 2^30,3^30,2^42
if(n>= pow(3.0, 30.0)) ans++;
if(n>= pow(2.0, 42.0)) ans++;
printf("%d\n",ans);
}
return 0;
}