记录编号 477715 评测结果 AAAAAAAAAA
题目名称 [vijos1889]天真的因数分解 最终得分 100
用户昵称 GravatarHallmeow 是否通过 通过
代码语言 C++ 运行时间 0.458 s
提交时间 2017-12-06 15:53:45 内存使用 6.04 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
using namespace std;
#define pos(i,a,b) for(int i=(a);i<=(b);i++)
#define N 500010
#define LL long long
LL k;
int mu[N],notprime[N],prime[N];
void getmu(){
	mu[1]=1;notprime[1]=1;
	pos(i,2,N-10){
		if(!notprime[i]){
			mu[i]=-1;
			prime[++prime[0]]=i;
		}
		for(int j=1;j<=prime[0]&&prime[j]*i<=N-10;j++){
			notprime[i*prime[j]]=1;
			if(i%prime[j]==0){
				mu[i*prime[j]]=0;break;
			}
			mu[i*prime[j]]=-mu[i];
		}
	}
}
LL l,r;
LL tot(LL x){
	LL top=(LL)sqrt(x);
	LL sum(0);
	for(LL i=1;i<=top;i++){
		sum+=x/(i*i)*mu[i];
	}
	return sum;
}
int main(){
	freopen("naive.in","r",stdin);
	freopen("naive.out","w",stdout);
	scanf("%lld",&k);
	getmu();
	l=k;r=(LL)1e11;
	while(l<=r){
		LL mid=(l+r)>>1;
		if(mid-tot(mid)<k) l=mid+1;
		else r=mid-1;
	}
	cout<<l;
	return 0;
}