8.12.28 ACM-ICPC数学数论 Stern–Brocot 树与 Farey 序列

本文链接：https://blog.csdn.net/tang7mj/article/details/139245870

8.12.28 ACM-ICPC数学数论 Stern–Brocot 树与 Farey 序列

引言

在数论中，Stern–Brocot 树和 Farey 序列是两个重要的概念，它们在分数的表示和数值逼近等方面具有广泛的应用。本节将详细介绍这两个概念及其在 ACM-ICPC 中的应用，并提供相关的 C++ 实现代码。

Stern–Brocot 树

Stern–Brocot 树是一种二叉树，用于生成所有正有理数的序列。树的每个节点表示一个分数，通过中位数运算生成新的分数节点。

Stern–Brocot 树的构建

Stern–Brocot 树的构建从两个初始节点开始，分别是 $0/1$ 和 $1/0$。新节点的生成遵循以下步骤：

计算两个父节点的中位数：
新分数=父节点1的分子+父节点2的分子父节点1的分母+父节点2的分母新分数=父节点1的分母+父节点2的分母父节点1的分子+父节点2的分子
将新分数插入树中，继续生成新的分数。

示例

假设初始节点为 $0/1$ 和 $1/1$，生成过程如下：

计算 $0/1$ 和 $1/1$ 的中位数：$\frac{0+1}{1+1} = \frac{1}{2}$
插入 $\frac{1}{2}$，得到新的节点序列：$0/1, 1/2, 1/1$
重复上述步骤，继续生成新的分数。

Stern–Brocot 树图示

性质

Stern–Brocot 树具有以下性质：

树中的所有节点分数都是最简分数。
树中的分数按大小有序排列。

C++ 实现

#include <iostream>
#include <vector>

struct Fraction {
    int num, den;
    Fraction(int n, int d) : num(n), den(d) {}
};

std::vector<Fraction> generateSternBrocotTree(int depth) {
    std::vector<Fraction> tree;
    tree.push_back(Fraction(0, 1)); // Initial left boundary
    tree.push_back(Fraction(1, 1)); // Initial right boundary

    for (int d = 0; d < depth; ++d) {
        std::vector<Fraction> nextTree;
        for (size_t i = 0; i < tree.size() - 1; ++i) {
            Fraction left = tree[i];
            Fraction right = tree[i + 1];
            nextTree.push_back(left);
            nextTree.push_back(Fraction(left.num + right.num, left.den + right.den));
        }
        nextTree.push_back(tree.back());
        tree = nextTree;
    }

    return tree;
}

int main() {
    int depth = 3; // Depth of the tree
    std::vector<Fraction> tree = generateSternBrocotTree(depth);

    for (const auto& frac : tree) {
        std::cout << frac.num << "/" << frac.den << " ";
    }

    return 0;
}

Farey 序列

Farey 序列是将所有在某个范围内的最简分数按大小排序后形成的序列。

Farey 序列的定义

给定一个正整数 $n$，Farey 序列 $F_n$ 包含所有分母不超过 $n$ 的最简分数，并按升序排列。

Farey 序列的构建

Farey 序列可以通过以下步骤构建：

从 $0/1$ 和 $1/1$ 开始。
按顺序插入所有符合条件的最简分数。
确保每个插入的分数与序列中的其他分数互素。

示例

对于 $n=3$ 的 Farey 序列 $F_3$，其构建过程如下：

初始分数：$0/1, 1/1$
插入 $1/2$，得到 $0/1, 1/2, 1/1$
插入 $1/3$ 和 $2/3$，最终得到 $0/1, 1/3, 1/2, 2/3, 1/1$

Farey 序列图示

性质

Farey 序列具有以下性质：

每个分数都是最简分数。
相邻分数的差分满足一定的条件：若 $a/b$ 和 $c/d$ 是相邻分数，则 $|bc - ad| = 1$。

C++ 实现

#include <iostream>
#include <vector>
#include <algorithm>

struct Fraction {
    int num, den;
    Fraction(int n, int d) : num(n), den(d) {}
    bool operator<(const Fraction& other) const {
        return num * other.den < other.num * den;
    }
};

std::vector<Fraction> generateFareySequence(int n) {
    std::vector<Fraction> sequence;
    for (int den = 1; den <= n; ++den) {
        for (int num = 0; num <= den; ++num) {
            if (std::gcd(num, den) == 1) {
                sequence.push_back(Fraction(num, den));
            }
        }
    }
    std::sort(sequence.begin(), sequence.end());
    return sequence;
}

int main() {
    int n = 3; // Maximum denominator
    std::vector<Fraction> sequence = generateFareySequence(n);

    for (const auto& frac : sequence) {
        std::cout << frac.num << "/" << frac.den << " ";
    }

    return 0;
}