#include <bits/stdc++.h>
using namespace std;
long long gcd(long long a, long long b, long long& x, long long& y) {
  if (a == 0) {
    x = 0;
    y = 1;
    return b;
  }
  long long x1, y1;
  long long d = gcd(b % a, a, x1, y1);
  x = y1 - (b / a) * x1;
  y = x1;
  return d;
}
int main() {
  cin.sync_with_stdio(0);
  cin.tie(0);
  cin.exceptions(cin.failbit);
  long long n;
  cin >> n;
  long long n_ = n;
  vector<int> p;
  for (long long i = 2; i * i <= n_; i++) {
    if (n_ % i == 0) {
      p.push_back(i);
      while (n_ % i == 0) n_ /= i;
    }
  }
  if (n_ != 1) p.push_back(n_);
  if ((int)(p).size() == 1) {
    cout << "NO" << endl;
  } else {
    cout << "YES" << endl;
    cout << 2 << endl;
    long long x, y;
    gcd(p[0], p[1], x, y);
    if (x < 0) {
      swap(p[0], p[1]);
      swap(x, y);
    }
    cout << n / p[0] - x << " " << n / p[0] << endl;
    cout << -y << " " << n / p[1] << endl;
  }
  return 0;
}
