1. 程式人生 > >HDU2204 Eddy's愛好

HDU2204 Eddy's愛好

bit break pow can sca c++ 素數 eps 復雜

題意:問[1,n]有幾個是m^k的形式(m>0,k>1, n<=1e18)

題解:枚舉m,會發現很多,所以枚舉k,可以發現k<=63那麽只要每個k計算pow(n , 1/k),接下來去重,可以用容斥來做,2^61的復雜度,這裏考慮m^k = m^(k1+k2)一個合數可以分解為素數,那麽只要對素數容斥就可以了

#include <bits/stdc++.h>
#define maxn 101000
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
const double eps = 1e-6
; ll a[100]; int main(){ ll num, sum, cnt, ans, n, t; num = 0; for(ll i=2, j;i<=63;i++){ for(j=2;j*j<=i;j++){ if(i%j == 0) break; } if(j*j>i) a[num++] = i; } while(~scanf("%lld", &n)){ ans = t = cnt = 0; while
((1LL<<a[t]) <= n) t++; for(ll i=1;i<(1LL<<t);i++){ sum = 1;cnt = 0; for(ll j=0;j<t;j++){ if(i&(1LL<<j)) sum *= a[j],cnt++; } if(sum<=63&&(1LL<<sum)<=n){ if
(cnt&1) ans += (int)(eps+pow(n, 1.0/sum))-1; else ans -= (int)(eps+pow(n, 1.0/sum))-1; } if(cnt%2 == 0) ans--; else ans++; } printf("%lld\n", ans); } return 0; }

HDU2204 Eddy's愛好