std::vector<std::vector<int>> combination(int n) {
        std::vector<std::vector<int>> dp(n + 1, std::vector<int>(n + 1, 0));

        for (int i = 0; i <= n; ++i) {
            dp[i][0] = 1; // C(n, 0) = 1
            dp[i][i] = 1; // C(n, n) = 1
        }
        for (int i = 1; i <= n; ++i) {
            for (int j = 1; j < i; ++j) {
                dp[i][j] = dp[i - 1][j - 1] + dp[i - 1][j];
            }
        }
        return dp;
    }

    glm::vec2 CalculateBezierPoint(
        std::span<glm::vec2> points,
        float const          t) {
        // your code here:
        glm::vec2 result = glm::vec2 {0, 0};

        int n = static_cast<int>(points.size()) - 1;
        auto comb = combination(n);

        for (int i = 0; i <= n; ++i) {
            float coeff = static_cast<float>(comb[n][i]) * std::pow(t, i) * std::pow(1 - t, n - i);
            result += coeff * points[i];
        }
        return result;
    }