#include <cmath>
#include <deque>
#include <algorithm>
#include <iterator>
#include <list>
#include <tuple>
#include <map>
#include <unordered_map>
#include <queue>
#include <set>
#include <unordered_set>
#include <stack>
#include <string>
#include <vector>
#include <fstream>
#include <iostream>
#include <functional>
#include <numeric>
#include <iomanip> 
#include <stdio.h>
#include <assert.h>
//eolibraries
#define lnf 3999999999999999999
#define inf 999999999
#define fi first
#define se second
#define pb push_back
#define ll long long
#define ld long double
#define all(c) (c).begin(),(c).end()
#define sz(c) (int)(c).size()
#define make_unique(a) sort(all(a)),a.erase(unique(all(a)),a.end())
#define pii pair <int,int>
#define ftpii pair<pii,pii>
#define rep(i,n) for(int i = 0 ; i < n ; i++) 
#define drep(i,n) for(int i = n-1 ; i >= 0 ; i--)
#define crep(i,x,n) for(int i = x ; i < n ; i++)
#define vi vector <int> 
#define vec(...) vector<__VA_ARGS__>
#define fcin ios_base::sync_with_stdio(false),cin.tie(0),cout.tie(0);
//eodefine
using namespace std;

const int _maxn = 5e7;

vi prms;

void fdprm(int nupto) {

	//glblhere
	vi prm;
	//glblhere

	prm.resize(nupto+3,0);

	crep(i,2,sqrt(nupto)+1){
		if(prm[i] == 0) {
			for(int j = i*i ; j <= nupto ; j += i) prm[j] = i;
		}
	}

	crep(i,2,nupto+1) if(prm[i]==0) prms.pb(i); 

}

int d[_maxn+100];

int main(){
fcin;
	fdprm(_maxn);

	// vi d(_maxn+10),cnt(_maxn+10);
	rep(i,_maxn) d[i]=i;

	// cout<<sz(prms)<<"\n";
	for(auto p : prms){
		if(p*p>_maxn) break;
		for(ll j=p*p;j<=_maxn;j+=p*p){
			// d[j]=1;
		}
	}
	
	// cout<<sz(prms)<<"\n";

	// rep(i,_maxn){
	// 	cnt[d[i]]++;
	// }

	// ll ans=0;
	// crep(i,1,_maxn){
	// 	ans+=(ll)cnt[i]*cnt[i];
	// }

	// cout<<ans<<"\n";
/*
*/
	return 0;
}