HDU 2204 Eddy's爱好（容斥原理）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2204

解题报告：输入一个n让你求出[1，n]范围内有多少个数可以表示成形如m^k的样子。

不详细说了，自己一开始也忽略了三个素数的乘积的乘方的情况。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long INT;
INT prim[70];

int dabiao()
{
    int f = 0;
    for(int i = 2;i < 60;++i)
    {
        int flag = 1;

        for(int j = 2;j < i;++j)
        if(i % j == 0)
        {
            flag = 0;
            break;
        }
        if(flag) prim[f++] = i;
    }
    return f;
}

int calc(INT a,INT b,INT n)
{
    INT s = 1;
    while(b--)
    {
        if((INT)(n / a) <= s) return 0;
        s *= a;
    }
    return 1;
}
int main()
{
    INT n;
    int num = dabiao();
    while(scanf("%I64d",&n)!=EOF)
    {
        INT tot = 1;
        for(int i = 0;i < num;++i)
        {
            INT temp = (INT)pow((double)n,1.0 / prim[i]);
            if(temp == 1 ) break;   //开方之后不足1的话，则退出
            tot += temp - 1;   //减1是因为每次都会统计到1
        }
        for(int i = 0;i < num;++i)
        for(int j = i+1;j < num;++j)
        {
            INT temp = pow((double)n,1.0/(prim[i] * prim[j]));
            if(temp == 1) break;
            tot -= temp - 1;
        }
        for(int i = 0;i < num;++i)
        for(int j = i + 1;j < num;++j)
        for(int k = j + 1;k < num;++k)
        {
            INT temp = pow((double)n,1.0/(prim[i] * prim[j] * prim[k]));
            if(temp == 1) break;
            tot += temp - 1;
        }
        printf("%I64d
",tot);
    }
    return 0;
}