bzoj 2986: Non-Squarefree Numbers

Description

一个正整数K被称为squarefree,如果它没有一个D^2(D>1)这样的约数。

Input

读入一个正整数N

Output

找出第N个不是squarefree的数。1<=N<=10^10

Sample Input

10

Sample Output

27

Hint
前10个非squarefree的数
4 8 9 12 16 18 20 24 25 27

说好双倍经验系列，数据读入大了10倍结果要longlong了是个什么鬼(╯°口°)╯(┴—┴

改了以后还T了(╯°口°)╯(┴—┴

-----------------------------------------------------------------------------------------------------------------------------------

这题可以直接用莫比乌斯函数做。和2440的思路相反

#include<cstdio>#include<string>#include<cstring>using namespace std;int mu[10000001],prime[10000001];bool check[10000001]; int tot;inline void findmu(){     memset(check,false,sizeof(check));     mu[1]=1;     int i,j;     for(i=2;i<=10000000;i++)     {          if(!check[i])          {               tot++;               prime[tot]=i;               mu[i]=-1;          }          for(j=1;j<=tot;j++)          {               if(i*prime[j]>10000000)                    break;               check[i*prime[j]]=true;               if(i%prime[j]==0)               {                    mu[i*prime[j]]=0;                    break;               }               else                    mu[i*prime[j]]=-mu[i];          }     }}inline long long gets(long long x){     long long i;     long long t=0;     for(i=1;i*i<=x;i++)          t+=mu[i]*(x/(i*i));     return t;}int main(){     findmu();     long long n;     scanf("%lld",&n);     long long l=1,r=(long long)2100000000*(long long)100;     while(l<=r)     {          long long mid=(l+r)/(long long)2;          long long t=gets(mid);          if(mid-t>=n)               r=mid-1;          else               l=mid+1;     }     printf("%lld\n",l);     return 0;}