オイラー機能---------最大公約数

シーク、N整数所与1 <= X、Y <= N およびGCD（x、y）は素数の数（x、y）のペアの数です。
X、Y最大公約数を見つけることであるGCD（X、Y）。
入力フォーマット
入力N整数
出力フォーマット
番号で表される条件を満足する出力整数。
データ範囲
1≤N≤1071≤N≤107
入力サンプル：
4

サンプル出力：
4

アイデア：
（X、Y）= P
（X / P、Y / P）= 1。
列挙子すべての素数、検索の予選のどのように多くの数

#include <cstring>
#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1e7 + 10;
typedef long long LL;
int primes[N], cnt;
bool st[N];
int phi[N];
LL s[N];
void init(int n){
 for (int i = 2; i <= n; i ++){
  if (!st[i])   primes[cnt ++] = i, phi[i] = i - 1;
  for (int j = 0; i * primes[j] <= n; j ++){
   st[i * primes[j]] = 0;
   if (i % primes[j] == 0){
    phi[i * primes[j]] = primes[j] * phi[i];
    break; 
   }
   phi[i * primes[j]] = phi[i] * (primes[j] - 1);
   }
 }
 for (int i = 1; i <= n; i ++)   s[i] = s[i - 1] + phi[i];
} 
int main(){
 int n;
 cin >> n;
 init(n);
  LL res = 0;
 for (int i = 0; i < cnt; i ++){
  int p = primes[i];
  res += s[n / p] * 2 + 1;
 }
  cout << res << endl;
  return 0;
}