Advertisement
Morass

Number Theory

Jan 20th, 2016
97
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.18 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<30)
  6. #define CL(A,I) (memset(A,I,sizeof(A)))
  7. #define DEB printf("DEB!\n");
  8. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  9. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  10. typedef long long ll;
  11. typedef pair<ll,ll> pll;
  12. typedef vector<int> vi;
  13. typedef pair<int,int> ii;
  14. #define IN(n) int n;scanf("%d",&n);
  15. #define FOR(i, m, n) for (int i(m); i < n; i++)
  16. #define REP(i, n) FOR(i, 0, n)
  17. #define F(n) REP(i, n)
  18. #define FF(n) REP(j, n)
  19. #define FT(m, n) FOR(k, m, n)
  20. #define aa first
  21. #define bb second
  22. void ga(int N,int *A){F(N)scanf("%d",A+i);}
  23. #define MX (1<<17)
  24. ll P[MX],L,cn[MX]={1,1};
  25. void gen(void){
  26.     for(int i(2);i<MX;++i)
  27.         if(!cn[i]){
  28.             P[L++]=i;
  29.             for(int j(i<<1);j<MX;j+=i)
  30.                 cn[j]=1;
  31.         }
  32. }
  33. ll sl(ll n){
  34.     ll t(1),h(1);
  35.     for(ll i(0),ct(0),d(1);P[i]*P[i]<=n;++i,ct=0,d=1){
  36.         while(!(n%P[i]))
  37.             n/=P[i],++ct,d*=P[i];
  38.         t*=ct+1,h*=(d-d/P[i]);
  39.     }
  40.     if(n!=1)t<<=1,h*=n-1;
  41.     return t+h-1;
  42. }
  43. int main(void){
  44.     ll N;
  45.     gen();
  46.     while(~scanf("%lld",&N))printf("%lld\n",N-sl(N));
  47.     return 0;
  48. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement