ll@p,@q,b[2d5],v[2d5],n=Divisor(q*q,b),t=0;
rep(i,n){
	ll x=b[i];
	if((x+q)%p+(q*q/x+q)%p<1){
		v[t++]=(x+q)/p;
	}
}
wt(t);
rep(i,t)wt(v[i],v[t+~i]);