求使 y=sqrt(x+a)+sqrt(x+b) 成立的正整数对 (x,y) 的数量这一类的题如何解？ 第1页

       #约数列表 dlist<-function(n) {     ans=c()     for(j in 1:n)     {         if((n%%j)==0){ans=c(ans,j)}     }     ans }  #求因子p的指数 index<-function(n,p) {     r=n;m=log(n,p)     if(m%%1==0)return(m)     for(j in 0:m)     {         if(r%%p!=0)return(j)         r=r/p     }     return(j) }  #非负取整函数 pf<-function(x) {     if(x>0){y=floor(x)}else{y=0}     y }  #方程y=sqrt(x+a)+sqrt(x+b)整数解个数（理论计算） Key<-function(a,b) {     K=abs(a-b)      alpha=index(K,2)      L=K/2^alpha     dL=length(dlist(L))      if(alpha%%2==0)     {         ans=pf(alpha/2-1)*dL+floor(dL/2)     }else{         ans=floor(alpha/2)*floor(dL)          }     ans }  fx<-function(a,b,x) {     y=sqrt(x+a)+sqrt(x+b)     y }  #方程y=sqrt(x+a)+sqrt(x+b)整数解个数（数值验证） key<-function(a,b) {     K=abs(a-b);m=min(a,b)     for(j in -m:(K^2-m))     {         if(fx(a,b,j)%%1==0)print(j)     } }