#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + 5;
int exgcd(int a, int b, int &x, int &y) {
  if (b == 0) {
    x = 1;
    y = 0;
    return a;
  }
  int r = exgcd(b, a % b, x, y);
  int t = x;
  x = y;
  y = t - a / b * y;
  return r;
}
vector<int> vc;
int prime[N];
bool isprime[N];
int tot;
void work_sieve_prime() {
  tot = 0;
  for (int i = 2; i < N; i++) {
    if (!isprime[i]) prime[++tot] = i;
    for (int j = 1; j <= tot && prime[j] * i < N; j++) {
      isprime[prime[j] * i] = 1;
      if (i % prime[j] == 0) break;
    }
  }
}
int main() {
  work_sieve_prime();
  int k;
  scanf("%d", &k);
  int tp = k;
  for (int j = 1; j <= tot; j++) {
    if (tp % prime[j] == 0) {
      while (tp % prime[j] == 0) {
        tp /= prime[j];
      }
      vc.push_back(prime[j]);
    }
    if (vc.size() == 2) break;
  }
  if (tp != 1) {
    vc.push_back(tp);
  }
  if (vc.size() < 2) {
    printf("NO\n");
    return 0;
  }
  int x, y;
  int a = vc[0];
  int b = vc[1];
  exgcd(a, b, x, y);
  printf("YES\n2\n");
  if (x < 0)
    printf("%d %d\n%d %d\n", -x, k / a, k / b - y, k / b);
  else
    printf("%d %d\n%d %d\n", k / a - x, k / a, -y, k / b);
  return 0;
}
