【HDU】2204 Eddy's爱好容斥原理

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

poursoul

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

阅读量653

点赞数

分类专栏：数学文章标签： HDU

本文链接：https://blog.csdn.net/u013368721/article/details/39102321

版权

数学专栏收录该内容

25 篇文章 0 订阅

订阅专栏

传送门：【HDU】2204 Eddy's爱好

题目分析：首先，对于所有形如M^K的数我们都可以转化成M^(p1^k1 * p2^k2 * p3^k3 * ... )的形式，其中p1,p2,p3..为素数。则所有的M^K都可以转化成M'^p，其中p为素数。我们意识到2^60>10^18，所以只要求出前60以内的所有素数即可。然后，由于2*3*5*7>60，其中一个K的质因数最多只有三种。

令K'为每种质因数只有一个的数。

我们可以通过x = pow ( n , 1 / K' )求得指数为K'的数在[ 1 , n ]中出现的次数。

注意到可能存在x^p1 = y^p2（p1，p2为素数），那么我们就可以用容斥原理来去除重复。

则ans = sum { 指数为素数的个数x - 1 } - sum { 指数为两个质因数之积的个数x - 1 } + sum {指数为三个质因数之积的个数x - 1 } + 1。（x-1是因为不能包括这个数本身，否则K = 1不合题意；最后+1是因为还有个1）

这题精度要尤其注意。

代码如下：

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

#define REP( i , a , b ) for ( int i = ( a ) ; i <  ( b ) ; ++ i )
#define FOR( i , a , b ) for ( int i = ( a ) ; i <= ( b ) ; ++ i )
#define REV( i , a , b ) for ( int i = ( a ) ; i >= ( b ) ; -- i )
#define CLR( a , x ) memset ( a , x , sizeof a )
#define CPY( a , x ) memcpy ( a , x , sizeof a )

typedef long long LL ;

const double eps = 1e-8 ;

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

LL n ;

int dfs ( LL n , int cur , int num ) {
	int ans = 0 ;
	REP ( i , cur , 17 ) {
		int tmp = pow ( n , 1.0 / ( num * prime[i] ) ) + eps ;
		if ( pow ( 1.0 * tmp , 1.0 * num * prime[i] ) > n + eps ) -- tmp ;
		-- tmp ;
		if ( tmp <= 0 ) break ;
		ans += tmp - dfs ( n , i + 1 , num * prime[i] ) ;
	}
	return ans ;
}

int main () {
	while ( ~scanf ( "%I64d" , &n ) ) printf ( "%d\n" , dfs ( n , 0 , 1 ) + 1 ) ;
	return 0 ;
}