数学小记:1988IMOP6_随笔
JueFan 一只绝帆

数学小记:1988IMOP6

改编版 OI 题:给定 nn,对满足 (xy+1)(x2+y2),x,y[1,n](xy+1)\mid(x^2+y^2),x,y\in[1,n](x,y)(x,y) 计数,n1018n\le 10^{18}

对其分析,首先注意到 xy,x2,y2xy,x^2,y^2 都是二次项,令 xyx\le y 则可以推出 x2xyy2x^2\le xy\le y^2,则 x2+y2xy+1yx\frac{x^2+y^2}{xy+1}\approx \frac{y}{x}

这是初步的估计,接下来我们分析一下:

  • (xy+1)(yx1)=y2xy+yx1<y2+x2(xy+1)(\frac y x-1)=y^2-xy+\frac yx -1<y^2+x^2
  • (xy+1)(yx+1)=y2+xy+yx+1>y2+x2(xy+1)(\frac y x+1)=y^2+xy+\frac y x+1> y^2+x^2

所以答案要么是 yx\lfloor\frac y x\rfloor 要么是 yx\lceil\frac y x\rceil

y=qx+ry=qx+r,若 qqx2+y2xy+1\frac{x^2+y^2}{xy+1}r(x,x)r\in(-x,x),我们先讨论 r=0r=0 的特殊情况。

x2+(qx)2x(qx)+1=qx2+q2x2=q2x2+qx2=q,y=x3\begin{aligned}&\frac{x^2+(qx)^2}{x(qx)+1}=q \\&x^2+q^2x^2=q^2x^2+q \\&x^2=q,y=x^3 \end{aligned}

代入验证发现 y=x3y=x^3 确实是一组特解。

讨论 r>0r>0

x2+(qx+r)2x(qx+r)+1=qx2+q2x2+2qxr+r2=q2x2+qxr+qx2+r2+qxr=q\begin{aligned}&\frac{x^2+(qx+r)^2}{x(qx+r)+1}=q \\&x^2+q^2x^2+2qxr+r^2=q^2x^2+qxr+q \\&x^2+r^2+qxr=q\end{aligned}

显然 qxrq,x2+r2+qxr>qqxr\ge q,x^2+r^2+qxr>q,该情况不成立。

讨论 r<0r<0,另设 y=qxr,r(0,x)y=qx-r,r\in(0,x)

x2+(qxr)2x(qxr)+1=qx2+q2x22qxr+r2=q2x2qxr+qx2+r2=qxr+q\begin{aligned}&\frac{x^2+(qx-r)^2}{x(qx-r)+1}=q \\&x^2+q^2x^2-2qxr+r^2=q^2x^2-qxr+q \\&x^2+r^2=qxr+q \end{aligned}

看起来我们一筹莫展了,但是观察形式:x2+r2=qxr+q=q(xr+1)x^2+r^2=qxr+q=q(xr+1),我们由 (x,y)(x,y) 递归到了 (xymodx,x)(x-y\bmod x,x)

而递归终止条件是我们推过的特解,y=x3y=x^3

注意这个递归的形式十分苛刻,不仅要保证两参数的形式,还要保证 yx\lceil\frac{y}{x}\rceil 是固定的值 qq,它等于递归终点的 x2x^2

因此我们会计数了,枚举递归终点,反向生成即可,log\log 层之内结束。

 评论
评论插件加载失败
正在加载评论插件
由 Hexo 驱动 & 主题 Keep
总字数 231.7k 访客数 访问量