#include #include #define rep(i,n) for(int i=0;i vi; typedef vector vl; typedef vector> vvi; typedef vector> vvl; typedef long double ld; typedef pair P; template ostream& operator<<(ostream& os, const static_modint& a) {os << a.val(); return os;} template ostream& operator<<(ostream& os, const dynamic_modint& a) {os << a.val(); return os;} template istream& operator>>(istream& is, static_modint& a) {long long x; is >> x; a = x; return is;} template istream& operator>>(istream& is, dynamic_modint& a) {long long x; is >> x; a = x; return is;} template istream& operator>>(istream& is, vector& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;} template ostream& operator<<(ostream& os, const pair& p){os << p.first << ' ' << p.second; return os;} template ostream& operator<<(ostream& os, const vector& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); return os;} template ostream& operator<<(ostream& os, const vector>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : ""); return os;} template ostream& operator<<(ostream& os, const set& se){for(T x : se) os << x << " "; os << "\n"; return os;} template ostream& operator<<(ostream& os, const unordered_set& se){for(T x : se) os << x << " "; os << "\n"; return os;} template ostream& operator<<(ostream& os, const atcoder::segtree& seg){int n = seg.max_right(0, [](S){return true;}); rep(i, n) os << seg.get(i) << (i == n - 1 ? "\n" : " "); return os;} template ostream& operator<<(ostream& os, const atcoder::lazy_segtree& seg){int n = seg.max_right(0, [](S){return true;}); rep(i, n) os << seg.get(i) << (i == n - 1 ? "\n" : " "); return os;} template void chmin(T& a, T b){a = min(a, b);} template void chmax(T& a, T b){a = max(a, b);} // Sieve of Eratosthenes // https://youtu.be/UTVg7wzMWQc?t=2774 struct Sieve { int n; vector f, primes; Sieve(int n=1):n(n), f(n+1) { f[0] = f[1] = -1; for (long long i = 2; i <= n; ++i) { if (f[i]) continue; primes.push_back(i); f[i] = i; for (long long j = i*i; j <= n; j += i) { if (!f[j]) f[j] = i; } } } bool isPrime(int x) { return f[x] == x;} vector factorList(int x) { vector res; while (x != 1) { res.push_back(f[x]); x /= f[x]; } return res; } vector> factor(ll x) { vector> res; for (int p : primes) { int y = 0; while (x%p == 0) x /= p, ++y; if (y != 0) res.emplace_back(p,y); } if (x != 1) res.emplace_back(x,1); return res; } int num_of_divisors(long long x){ auto f = factor(x); long long res = 1; for(auto [p, val] : f){ res *= val + 1; } return res; } set divisors(long long x){ set se = {1}; auto f = factor(x); for(auto [key, value] : f){ set old = se; for(long long y : old){ rep(i, value + 1){ se.insert(y); y *= key; } } } return se; } }; const int M = 200000; int main(){ long long x; cin >> x; if(x == 1){ cout << "2\n"; cout << "1 2\n"; cout << "b g\n"; return 0; } Sieve sieve(M); vector vec; auto f = sieve.factor(x); long long tot = 0; for(auto [prime, cnt] : f){ if(prime == 2 and cnt >= 2){ vec.emplace_back(4); tot += 5; cnt -= 2; } rep(_, cnt){ vec.emplace_back(prime); tot += prime + 1; } } if(tot > M){ cout << "-1\n"; return 0; } cout << tot << "\n"; vector> G; int idx1 = 0; int idx2 = 0; for(long long val : vec){ idx2 += val; } vector color(tot); for(long long val : vec){ rep(_, val){ G.emplace_back(idx1, idx2); color[idx1] = 'g'; idx1++; } color[idx2] = 'b'; idx2++; } while(idx1 < idx2 - 1){ G.emplace_back(idx1, idx1 + 1); idx1++; } for(auto [u, v] : G){ cout << u + 1 << ' ' << v + 1 << "\n"; } cout << color; return 0; }