hdu2204 Eddy's 爱好【容斥原理】

最新推荐文章于 2020-02-20 22:20:18 发布

lhfl911

最新推荐文章于 2020-02-20 22:20:18 发布

阅读量401

点赞数

本文链接：https://blog.csdn.net/lhfl911/article/details/52129769

版权

数论同时被 2 个专栏收录

9 篇文章 0 订阅

订阅专栏

容斥原理

3 篇文章 0 订阅

订阅专栏

题目链接： hdu2204

题意：
给一正整数n $(1<= n <= 10^{18})$ ，求[1,n]中能表示成 $M^K(K>1)$ 的数的个数。

第一次好好看容斥，感觉智商已经不够了，看了好多题解才明白。。。

若 $M^K <= n$ 满足，则任意m <= M 都满足 $m^K<= n$ ；
满足 $M^K$ 的数中，若K不是素数，则 $M^K = {(M^a)}^p$ , 其中a*p = K且p是素数，即满足题意的数都可以表示成 $x^p$ ，p是素数，答案仅需找出满足 $x^p$ 的个数；
$x^p$ 中有许多重复的，若 $x = y^q$ (q是素数且 $p\ne q$ )， $(x^p)^q = (x^q)^p$ , 例如 $(6^2)^3 = (6^3)^2$ ，
规律满足容斥原理， $ans = f(p_1) + ... - f(p_1*p_2) -... +...$ ;
因为 $2^{60} > 10^{18}$ , 2*3*5*7 > 60, 所以所需最大素数为59，最高组合数不会过3个；
注意精度问题，虽说数据很水，不写也会过。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;


typedef long long LL;

const double eps = 1e-9;

int p[17] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59};
LL n, ans;

void dfs(int cur, int ex, int pos)
{
    if(pos > 3)
        return;
    for(int i = cur; i < 17; i++)
    {
        LL num = pow(n, 1.0/(p[i]*ex)) + eps;
        num--;
        if(num > 0)
            pos&1? ans += num: ans -= num;
        dfs(i+1, ex*p[i], pos+1);
    }
}
int main(int argc, char const *argv[])
{
    while(~scanf("%lld", &n))
    {
        ans = 0;
        dfs(0, 1, 1);
        printf("%lld\n", ans+1);
    }
    return 0;
}