HDU 2204 Eddy's爱好(数论，容斥原理)

最新推荐文章于 2019-05-25 11:07:07 发布

acm_Gavin

最新推荐文章于 2019-05-25 11:07:07 发布

阅读量727

点赞数

分类专栏：数学文章标签：算法容斥原理数学

本文链接：https://blog.csdn.net/perAres/article/details/51044319

版权

数学专栏收录该内容

1 篇文章 0 订阅

订阅专栏

Eddy's爱好

Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2162 Accepted Submission(s): 939

Problem Description

Ignatius 喜欢收集蝴蝶标本和邮票，但是Eddy的爱好很特别，他对数字比较感兴趣，他曾经一度沉迷于素数，而现在他对于一些新的特殊数比较有兴趣。
这些特殊数是这样的：这些数都能表示成M^K，M和K是正整数且K>1。
正当他再度沉迷的时候，他发现不知道什么时候才能知道这样的数字的数量，因此他又求助于你这位聪明的程序员，请你帮他用程序解决这个问题。
为了简化，问题是这样的：给你一个正整数N，确定在1到N之间有多少个可以表示成M^K（K>1)的数。

Input

本题有多组测试数据，每组包含一个整数N，1<=N<=1000000000000000000(10^18).

Output

对于每组输入，请输出在在1到N之间形式如M^K的数的总数。
每组输出占一行。

Sample Input

  
  
   
   10
36
1000000000000000000

Sample Output

Author

Eddy

Recommend

lcy

容斥原理求解，x^p<=n,开根号求大于1的x的个数，对于指数使用容斥原理，注意精度问题。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
#include<vector>
#include<cmath>
#include<string>
using namespace std;
typedef long long LL;
const int MAXN=1000008;
int prime[MAXN];
bool vis[MAXN];
LL ans,n;
int tot;
int getPrime(int n){
	tot=0;
	memset(vis,0,sizeof(vis));
	for(int i=2;i<=n;i++){
		if(!vis[i]){
            prime[tot++]=i;
		}
		for(int j=0;j<tot&&i*prime[j]<=n;j++){
			vis[i*prime[j]]=true;
			if(i%prime[j]==0){
                break;
			}
		}
	}
	return tot;
}
void dfs(int dep,int mul,int num){
    for(int i=num;i<tot&&prime[i]*mul<=(int)(log(n)/log(2.0));i++){
        int tp=mul*prime[i];
        if(dep&1){
            ans-=(LL)(pow(n,1/(double)tp)+0.000001)-1;//1不算，注意精度
        }else{
            ans+=(LL)(pow(n,1/(double)tp)+0.000001)-1;
        }
        dfs(dep+1,tp,i+1);
    }
}
int main(){
    getPrime(1000000);
    while(~scanf("%I64d",&n)){
        ans=1;//1满足条件
        dfs(0,1,0);
        printf("%I64d\n",ans);
    }

    return 0;
}