#include using namespace std; void _main(); int main() { _main(); return 0; } using ll = long long; using ull = unsigned long long; using ld = long double; using vi = std::vector; using vvi = std::vector; using vl = std::vector; using vii = std::vector >; using vvl = std::vector; using vll = std::vector >; using vd = std::vector; using vvd = std::vector; using vs = std::vector; using vvs = std::vector; using vb = std::vector; using vvb = std::vector; using vc = std::vector; using vvc = std::vector; using pii = std::pair; using pll = std::pair; using piil = std::pair, ll>; using mii = std::map; using mll = std::map; using pql = std::priority_queue; using pqi = std::priority_queue; using pqiil = std::priority_queue, ll> >; using pqii = std::priority_queue >; #define pb push_back #define ps push #define eb emplace_back #define is insert #define er erase #define f first #define s second #define lb lower_bound #define ub upper_bound #define sf(i) sizeof(i) #define endl "\n" #define sz(v) ((int)(v).size()) #define all(v) (v).begin(), (v).end() #define rep(i, L, R) for(ll i = L;i<=R;i++) #define pcis precision template struct infinity { static constexpr T max=std::numeric_limits::max(); static constexpr T min=std::numeric_limits::min(); static constexpr T value=std::numeric_limits::max()/2; static constexpr T mvalue=std::numeric_limits::min()/2; }; templateconstexpr T INF=infinity::value; constexpr ll lINF=INF; constexpr int iINF = INF; constexpr ld PI = 3.1415926535897932384626; double is[10000001]; ll phi(ll N) { if (is[N] != 0)return is[N]; ll res = N; for (long i = 2; i <= sqrt(N); i++) { if (N % i == 0) { res = res - res / i; while (N % i == 0) N /= i; } } if (N > 1) res = res - res / N; return is[N] = res; } ll gcd(ll a, ll b) { if (b==0)return a; else return gcd(b, a%b); } struct ra { ll up = 1; ll down = 1; }; ra radd(ra a, ra b) { ra c; c.up = a.up*b.down+b.up*a.down; c.down = a.down*b.down; return c; } ra remp(ra a, ll em) { ra c; c.up = a.up*em; c.down = a.down; return c; } ra rdemp(ra a, ll em) { ra c; c.up = a.up; c.down = a.down*em; return c; } ra rmul(ra a, ra b) { ra c; c.up = a.up*b.up; c.down = a.down*b.down; return c; } ra rdiv(ra a, ra b) { ra c; c.up = a.up*b.down; c.down = a.down*b.up; return c; } ra rsub(ra a, ra b) { return radd(a, remp(b, -1)); } bool rsame(ra a, ra b) { return (a.up*b.down==a.down*b.up); } ra mr(ll a, ll b) { ra c; c.up = a; c.down = b; return c; } ra DP[10000001]; // N일때 기댓값 반환 ra rei = {0,1}; ra GET(ll N) { if (N==1)return {0, 1}; if (!rsame(DP[N], rei))return DP[N]; ra kitai = rei; for (ll i = 1; i*i<=N; i++) { // i가 되는 경우 if (N%i == 0) { kitai = radd(kitai, remp(GET(i) , phi(N/i))); if (i != 1 && i*i != N) { kitai = radd(kitai, remp(GET(N/i) , phi(i))); } } } kitai = radd(kitai, mr(N, 1)); kitai = rmul(kitai, mr(1, N-1)); return DP[N] = kitai; } void _main() { ll N; cin >> N; for (ll i = 1; i<=N; i++)DP[i]=rei; cout.pcis(10); cout << GET(N).up/GET(N).down; // cout << GET(N).up << "/" << GET(N).down << endl; }