Codeforces1036F Relatively Prime Powers 【容斥原理】

题目分析:

这种题目标题写莫比乌斯反演会不会显得太恐怖了，那就容斥算了。

gcd不为1的肯定可以开根。所以把根式结果算出来就行了。

辣鸡题目卡我精度。

代码：

 #include<bits/stdc++.h>
 using namespace std;

 const long long LMAX = ;

 long long n;
 int mu[];

 void init(){
     for(int i=;i<=;i++){
     int p = i;
     mu[i] = -;
     for(int j=;j*j<=p;j++){
         int cnt = ; while(p%j == ) p/=j,cnt++;
         if(cnt !=  && cnt != ) mu[i] = ;
         else if(cnt == ) mu[i]*=-;
     }
     if(p != ) mu[i]*=-;
     }
 }

 long long fast_pow(int now,int pw){
     long long ans = ,dt = now;
     int bit = ;
     while(bit <= pw){
     if(bit & pw){
         if(ans < LMAX/dt) ans *= dt;
         else ans = LMAX;
     }
     if(dt < LMAX/dt) dt *= dt;
     else dt = LMAX;
     bit<<=;
     }
     return ans;
 }

 void work(){
     long long ans = ;
     for(int i=;i<=;i++){
     if(mu[i] == ) continue;
     int z = pow((long double)n,1.0/(long double)i);
     if(z == ) break;
     if(fast_pow(z,i) > n) z--;
     if(fast_pow(z+,i) <= n) z++;
     z--; ans += z*mu[i];
     }
     n -= ans;n--;
     printf("%I64d\n",n);
 }

 int main(){
     int Tmp; scanf("%d",&Tmp);
     init();
     while(Tmp--){
     scanf("%I64d",&n);
     work();
     }
     return ;
 }