实数的最优逼近_笔记

连分数

$x$ 为一个实数，它的连分数写成：
$$x=a_0+\frac{1}{a_1+\frac{1}{a_2+\frac{1}{a_3+\frac{1}{\ddots +\frac{1}{a_n}}}}}$$
我们简记成
$$x=[a_0,a_1,a_2,...,a_n]$$
我们约定对于 $k\ge 1,a_k\ne 0$，它满足一些显然的性质比如
$$[a_0,a_1,a_2,...,a_n]=a_0+\frac{1}{[a_1,a_2,a_3,...,a_n]}$$
记 $r_k=[a_k,a_{k+1},...,a_n]$ ，后面会用到（

对于一个有理数 $x$ ，它的连分数就直接暴力求了（也叫辗转相除法。

那么对于一个无理数呢，比如：
$$\sqrt{2}=[1,2,2,2,2,...],\sqrt{3}=[1,1,2,1,2,1,2,...]$$
有些无理数的连分数很简单的，也比较好推。像 $\pi$ 这种就只能暴力算了。

渐近分数

我们把连分数保留到 $a_k$ 的近似记作 $s_k$ ，称作一个连分数的第k个近似或者第k个截断，把 $s_k=[a0,\dots ,ak]$ 写成分数 $p_k/q_k$ ,就称为连分数的第 $k$ 个渐近分数或者 $k$ 阶渐近分数。

渐近分数可以递推来求。
$$\{\begin{array}{l}{p}_k=a_k{p}_{k-1}+p_{k-2}\\ q_k=a_k{q}_{k-1}+q_{k-2}\end{array}$$

中间分数

称
$$F_{k,n}=\frac{n{p}_k+p_k-1}{n{q}_k+q_k-1},其中(0\le n\le a_{k+1})$$
为中间分数。简单来说就是 $F_{k,n}=[a_0,...,a_k,n]$ 。

最优逼近

我们定义第一类最优逼近如下：
对实数 $x$ 与有理数 $p/q$ ，若对 $\forall 0<b\le q,a,b\in Z$ ,且$a/b\ne p/q$均有
$$|x-\frac{a}{b}|>|x-\frac{p}{q}|$$
则称 $p/q$ 是 $x$ 的一个第一类最优逼近。
简单来说，若所有分母不超过 $q$ 的有理数都比 $p/q$ 更远离 $x$ ，那么 $p/q$ 就是最优逼近。

定义第二类最优逼近 :
对实数 $x$ 与有理数 $p/q$ ，若对 $\forall 0<b\le q,a,b\in Z$，且$a/b\ne p/q$ 均有
$$|x-\frac{a}{b}|>\frac{q}{b}|x-\frac{p}{q}|$$
这说明第二类最优逼近是一个强于第一类的限制，$p/q$ 不仅要比所有分母不超过 $q$ 的有理数更接近 $x$，还要在乘上一个不小于 $1$ 的倍数后仍然最接近 $x$。

定理 1：第二类最优逼近一定是第一类最优逼近。

定理 2: 第二类最优逼近一定是渐近分数，且渐近分数一定是第二类最优逼近。
证明不难，但是很长，略了（

定理 3：第一类最优逼近一定是中间分数，且满足 $F_{k,n}\iff r_{k+2}>q_k/q_{k-1}$。
意思就是 $[a_{k+2},a_{k+3},...,a_n]>q_k/q_k-1$ 。
证明不会。

---

康康题目

题目描述

一个 $n$ 项的正实数列 $a_{1,2,...,n}$，其中 $a_k=k\sqrt{2}(k\le n,k\in \mathbb{N}\ast )$ ，将数轴上所有整数和所有 $a_k(k\le n,k\in \mathbb{N}\ast )$ 对应的点都描出。求线段中最短的一条所在的位置。

分析

这个就找距离 $\sqrt{2}$ 最近的分数，要求分母$\le n$ 。

对于任意正实数 $\alpha$，它可能的第一类最优逼近都是中间分数，而 $\sqrt{2}$ 的中间分数都不符合最优逼近的要求（或者你把中间分数的表打出来也行，反正数量级也是 $log$ 的），那么我们只考虑渐近分数，分母 $\le n$ 的渐进分数只有 $O(\log n)$ 个，直接打表即可。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll f[600],n,T;
void work(){
	cin>>n;
	for(int i=50;i>=1;i--)
		if(f[i]<=n){ cout<<f[i]<<'\n'; break; }
}
int main(){
	f[1]=1,f[2]=2;
	for(int i=3;i<=50;i++) f[i]=2ll*f[i-1]+f[i-2];
	cin>>T;
	while(T--) work();
}