【Foreign】Uria [欧拉函数]

Uria

Time Limit: 20 Sec  Memory Limit: 512 MB

Description

　　从前有个正整数 n。

　　对于一个正整数对 (a,b)，如果满足 a + b ≤ n 且 a + b 是 a * b 的因子，则成为神奇的数对。

　　求神奇的数对的个数。

Input

　　一行一个正整数 n。

Output

　　一行一个整数表示答案，保证不会超过 64 位有符号整数类型的范围。

Sample Input

　　21

Sample Output

　　11

HINT

　　n ≤ 1e14

Solution

　　

Code

 #include<iostream>
 #include<string>
 #include<algorithm>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<cmath>
 using namespace std;
 typedef long long s64;
 typedef unsigned int u32;

 const int ONE = 1e7 + ;
 const u32 MOD = ;

 s64 n, Q;
 s64 Ans;
 bool isp[ONE];
 s64 phi[ONE], prime[ONE], p_num;

 int get()
 {
         int res=,Q=;char c;
         while( (c=getchar())< || c> )
         if(c=='-')Q=-;
         res=c-;
         while( (c=getchar())>= && c<= )
         res=res*+c-;
         return res*Q;
 }

 void Get_phi(int MaxN)
 {
         phi[] = ;
         for(int i = ; i <= MaxN; i++)
         {
             if(!isp[i])
                 prime[++p_num] = i, phi[i] = i - ;
             for(int j = ; j <= p_num && i * prime[j] <= MaxN; j++)
             {
                 isp[i * prime[j]] = ;
                 if(i % prime[j] == )
                 {
                     phi[i * prime[j]] = phi[i] * prime[j];
                     break;
                 }
                 phi[i * prime[j]] = phi[i] * phi[prime[j]];
             }
         }
 }

 int main()
 {
         cin>>n;    Q = sqrt(n);
         Get_phi(Q);
         for(int i = ; i <= Q; i++)
             Ans += (n / i / i) * phi[i];
         printf("%lld", Ans);
 }