#include <bits/stdc++.h>
using namespace std;
int k, pcnt, tot;
int pri[100500], p[20];
bool vis[100500];
int exgcd(int a, int b, int &x, int &y) {
  if (!b) return x = 1, y = 0, a;
  int ret = exgcd(b, a % b, x, y);
  int buf = x;
  x = y, y = buf - a / b * y;
  return ret;
}
inline void get() {
  int bd = std::sqrt(k), tar = k;
  for (register int i = 2; i <= bd; ++i) {
    if (!vis[i]) {
      pri[++pcnt] = i;
      if (!(tar % i)) {
        tar /= i, p[++tot] = i;
        while (!(tar % i)) tar /= i;
      }
    }
    for (register int j = 1; j <= pcnt; ++j) {
      if (pri[j] * i > bd) break;
      vis[pri[j] * i] = true;
      if (!(i % pri[j])) break;
    }
  }
  if (tar != 1 && tar > bd) p[++tot] = tar;
}
void solve() {
  int x, y;
  if (tot == 1) return puts("NO"), void();
  puts("YES");
  int gcd = exgcd(p[1], p[2], x, y);
  int fac1 = k / p[1], fac2 = k / p[2];
  if (x < 0)
    printf("2\n%d %d\n%d %d", -x, fac1, fac2 - y, fac2);
  else
    printf("2\n%d %d\n%d %d", fac1 - x, fac1, -y, fac2);
}
int main(void) {
  scanf("%d", &k);
  get();
  solve();
}
