bi 是n的因数，且ai<bi;
ai/bi+…+=n-1/n;
首先我们先变形一下
a1/b1+a2/b2+…+ai/bi=(n-1)/n
同乘n，（bi是n的因数，n/bi为整数）令ci=n/bi;
变为: a1c1+a2c2+…+aici=n-1;
因为ci也是n的因子，这里用到一个定理： ci 是n 的因子，则ci一定能用一个n的质因子的倍数表示出来。 因为任意一个数都能用多个质数相乘表示（除了1以外）。假设di 是n的质因子，且能表示出ci=eidi。
变为a1*(d1e1)+a2(d2e2)+…+ai(diei)=n-1;
此时上式中的每一项都是n的质因子的倍数，一个质数可以由两个小的质数倍数相加表示出来。取两个n的质因子a,b。且ab=n;
最后可以表示成ax+by=n-1
由exgcd 得知，ax+by=n-1有解的话必须是g=gcd（a,b）整倍数;
用exgcd 求出x，y ,此时求出的解是ax+by=g;的
变型：ax*(n-1)/g+by*(n-1)/g=n-1;
两边再除以n，得出
(a/n）*(n-1)x/g+(b/n）(n-1)*y/g=n-1/n
a/n=1/b; b/n=1/a;
xx=(n-1)*x/g,yy=(n-1)*y/g;
所以最后结果！！！！
xx/b+yy/a=n-1/n;
#include<queue>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<math.h>
#include<stdlib.h>
#include<stack>
#include<stdio.h>
#include<vector>
#include<map>
#include<string.h>
#define ll long long
#define db double
#define F(n) for(int i=1;i<=n;i++)
using namespace std;
const int mx=1e5+10, mod = 998244353;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;
}
ll exgcd(ll a,ll b,ll &x,ll &y){if(b==0){x=1,y=0;return a;}ll g=exgcd(b,a%b,x,y);ll tp=x;x=y;y=tp-a/b*y;return g;
}
int main()
{ll n,flag=0;scanf("%lld",&n);ll nn=sqrt(n);for(int i=2;i*i<=n;i++){if(n%i==0){ll a=i,b=n/i;ll g=gcd(a,b);if((n-1)%g==0&&a!=b){ll tx,ty;exgcd(a,b,tx,ty);tx=tx*(n-1)/g;ty=ty*(n-1)/g;while(tx<0){tx+=b;ty-=a;} while(ty<0){tx-=b;ty+=a;}flag=1;printf("YES\n2\n%lld %lld\n%lld %lld\n",ty,a,tx,b);break;}}}if(flag==0) printf("NO\n");return 0;
}