46

背景

考虑以下:

template <unsigned N>
struct Fibonacci
{
    enum
    {
        value = Fibonacci<N-1>::value + Fibonacci<N-2>::value
    };
};

template <>
struct Fibonacci<1>
{
    enum
    {
        value = 1
    };
};

template <>
struct Fibonacci<0>
{
    enum
    {
        value = 0
    };
};

这是一个常见的例子,我们可以将斐波那契数的值作为编译时常量:

int main(void)
{
    std::cout << "Fibonacci(15) = ";
    std::cout << Fibonacci<15>::value;
    std::cout << std::endl;
}

但是您显然无法在运行时获得该值:

int main(void)
{
    std::srand(static_cast<unsigned>(std::time(0)));

    // ensure the table exists up to a certain size
    // (even though the rest of the code won't work)
    static const unsigned fibbMax = 20;
    Fibonacci<fibbMax>::value;

    // get index into sequence
    unsigned fibb = std::rand() % fibbMax;

    std::cout << "Fibonacci(" << fibb << ") = ";
    std::cout << Fibonacci<fibb>::value;
    std::cout << std::endl;
}

因为fibb不是编译时常量。

问题

所以我的问题是:

在运行时查看此表的最佳方法是什么?最明显的解决方案(“解决方案”应该轻描淡写)是有一个大的 switch 语句:

unsigned fibonacci(unsigned index)
{
    switch (index)
    {
    case 0:
        return Fibonacci<0>::value;
    case 1:
        return Fibonacci<1>::value;
    case 2:
        return Fibonacci<2>::value;
    .
    .
    .
    case 20:
        return Fibonacci<20>::value;
    default:
        return fibonacci(index - 1) + fibonacci(index - 2);
    }
}

int main(void)
{
    std::srand(static_cast<unsigned>(std::time(0)));

    static const unsigned fibbMax = 20;    

    // get index into sequence
    unsigned fibb = std::rand() % fibbMax;

    std::cout << "Fibonacci(" << fibb << ") = ";
    std::cout << fibonacci(fibb);
    std::cout << std::endl;
}

但是现在表格的大小是非常硬编码的,将其扩展为 40 并不容易。

我想出的唯一一个具有类似查询方法的方法是:

template <int TableSize = 40>
class FibonacciTable
{
public:
    enum
    {
        max = TableSize
    };

    static unsigned get(unsigned index)
    {
        if (index == TableSize)
        {
            return Fibonacci<TableSize>::value;
        }
        else
        {
            // too far, pass downwards
            return FibonacciTable<TableSize - 1>::get(index);
        }
    }
};

template <>
class FibonacciTable<0>
{
public:
    enum
    {
        max = 0
    };

    static unsigned get(unsigned)
    {
        // doesn't matter, no where else to go.
        // must be 0, or the original value was
        // not in table
        return 0;
    }
};

int main(void)
{
    std::srand(static_cast<unsigned>(std::time(0)));

    // get index into sequence
    unsigned fibb = std::rand() % FibonacciTable<>::max;

    std::cout << "Fibonacci(" << fibb << ") = ";
    std::cout << FibonacciTable<>::get(fibb);
    std::cout << std::endl;
}

这似乎工作得很好。我看到的唯一两个问题是:

  • 调用堆栈可能很大,因为计算 Fibonacci<2> 需要我们通过 TableMax 一直到 2,并且:

  • 如果该值在表之外,则返回零而不是计算它。

那么我有什么遗漏吗?似乎应该有更好的方法在运行时挑选出这些值。

一个 switch 语句的模板元编程版本可能会生成一个特定数量的 switch 语句?

提前致谢。

4

7 回答 7

30
template <unsigned long N>
struct Fibonacci
{
    enum
    {
        value = Fibonacci<N-1>::value + Fibonacci<N-2>::value
    };
    static void add_values(vector<unsigned long>& v)
    {
        Fibonacci<N-1>::add_values(v);
        v.push_back(value);
    }
};

template <>
struct Fibonacci<0>
{
    enum
    {
        value = 0
    };
    static void add_values(vector<unsigned long>& v)
    {
        v.push_back(value);
    }

};

template <>
struct Fibonacci<1>
{
    enum
    {
        value = 1
    };
    static void add_values(vector<unsigned long>& v)
    {
        Fibonacci<0>::add_values(v);
        v.push_back(value);
    }
};



int main()
{
    vector<unsigned long> fibonacci_seq;
    Fibonacci<45>::add_values(fibonacci_seq);
    for (int i = 0; i <= 45; ++i)
        cout << "F" << i << " is " << fibonacci_seq[i] << '\n';
}

经过深思熟虑,我想出了这个解决方案。当然,您仍然必须在运行时将值添加到容器中,但(重要的是)它们不是在运行时计算的。

附带说明一下,重要的是不要在Fibonacci<1>上面定义Fibonacci<0>,否则编译器在解析对 的调用时会感到非常Fibonacci<0>::add_values困惑,因为Fibonacci<0>没有指定模板特化。

当然,TMP 有其局限性:您需要预先计算的最大值,并且在运行时获取值需要递归(因为模板是递归定义的)。

于 2009-05-25T23:16:14.090 回答
17

我知道这个问题很老,但它引起了我的兴趣,我不得不尝试在运行时不填充动态容器的情况下进行操作:

#ifndef _FIBONACCI_HPP
#define _FIBONACCI_HPP


template <unsigned long N>
struct Fibonacci
{
    static const unsigned long long value = Fibonacci<N-1>::value + Fibonacci<N-2>::value;

    static unsigned long long get_value(unsigned long n)
    {
        switch (n) {
            case N:
                return value;
            default:
                return n < N    ? Fibonacci<N-1>::get_value(n)
                                : get_value(n-2) + get_value(n-1);
        }
    }
};

template <>
struct Fibonacci<0>
{
    static const unsigned long long value = 0;

    static unsigned long long get_value(unsigned long n)
    {
        return value;
    }
};

template <>
struct Fibonacci<1>
{
    static const unsigned long long value = 1;

    static unsigned long get_value(unsigned long n)
    {
        return value;
    }
};

#endif

这似乎可行,并且当使用优化编译时(不确定您是否要允许),调用堆栈不会深入 - 当然对于值(参数)n> N,堆栈上有正常的运行时递归,其中 N 是模板实例化中使用的 TableSize。但是,一旦您低于 TableSize,生成的代码将替换在编译时计算的常量,或者最坏的情况是通过跳过跳转表“计算”的值(在 gcc 中使用 -c -g -Wa,-adhlns=main. s 并检查了清单),与我认为您的显式 switch 语句会导致的结果相同。

像这样使用时:

int main()
{
    std::cout << "F" << 39 << " is " << Fibonacci<40>::get_value(39) << '\n';
    std::cout << "F" << 45 << " is " << Fibonacci<40>::get_value(45) << '\n';
}

在第一种情况下根本没有计算调用(在编译时计算的值),在第二种情况下调用堆栈深度最差:

fibtest.exe!Fibonacci<40>::get_value(unsigned long n=41)  Line 18 + 0xe bytes    C++
fibtest.exe!Fibonacci<40>::get_value(unsigned long n=42)  Line 18 + 0x2c bytes    C++
fibtest.exe!Fibonacci<40>::get_value(unsigned long n=43)  Line 18 + 0x2c bytes    C++
fibtest.exe!Fibonacci<40>::get_value(unsigned long n=45)  Line 18 + 0xe bytes    C++
fibtest.exe!main()  Line 9 + 0x7 bytes    C++
fibtest.exe!__tmainCRTStartup()  Line 597 + 0x17 bytes    C

即它递归直到它在“表”中找到一个值。(通过在调试器中逐行执行反汇编验证,也通过用随机数 <= 45 替换测试整数来验证)

递归部分也可以用线性迭代解决方案代替:

static unsigned long long get_value(unsigned long n)
{
    switch (n) {
        case N:
            return value;    
        default:
            if (n < N) {
                return Fibonacci<N-1>::get_value(n);
            } else {
                // n > N
                unsigned long long i = Fibonacci<N-1>::value, j = value, t;
                for (unsigned long k = N; k < n; k++) {
                    t = i + j;
                    i = j;
                    j = t;
                }
                return j;
            }
    }
}
于 2011-05-18T06:11:16.620 回答
4

如果您有支持可变参数模板(C++0x 标准)的 C++ 编译器,您可以在编译时将 fibonacii 序列保存在元组中。在运行时,您可以通过索引访问该元组中的任何元素。

#include <tuple>   
#include <iostream>

template<int N>
struct Fib
{
    enum { value = Fib<N-1>::value + Fib<N-2>::value };
};

template<>
struct Fib<1>
{
    enum { value = 1 };
};

template<>
struct Fib<0>
{
    enum { value = 0 };
};

// ----------------------
template<int N, typename Tuple, typename ... Types>
struct make_fibtuple_impl;

template<int N, typename ... Types>
struct make_fibtuple_impl<N, std::tuple<Types...> >
{
    typedef typename make_fibtuple_impl<N-1, std::tuple<Fib<N>, Types... > >::type type;
};

template<typename ... Types>
struct make_fibtuple_impl<0, std::tuple<Types...> >
{
    typedef std::tuple<Fib<0>, Types... > type;
};

template<int N>
struct make_fibtuple : make_fibtuple_impl<N, std::tuple<> >
{};

int main()
{
   auto tup = typename make_fibtuple<25>::type();
   std::cout << std::get<20>(tup).value;  
   std::cout << std::endl; 

   return 0;
}
于 2011-06-18T21:04:21.040 回答
4

使用 C++11:您可以创建std::array一个简单的 getter:https ://ideone.com/F0b4D3

namespace detail
{

template <std::size_t N>
struct Fibo :
    std::integral_constant<size_t, Fibo<N - 1>::value + Fibo<N - 2>::value>
{
    static_assert(Fibo<N - 1>::value + Fibo<N - 2>::value >= Fibo<N - 1>::value,
                  "overflow");
};

template <> struct Fibo<0u> : std::integral_constant<size_t, 0u> {};
template <> struct Fibo<1u> : std::integral_constant<size_t, 1u> {};

template <std::size_t ... Is>
constexpr std::size_t fibo(std::size_t n, index_sequence<Is...>)
{
    return const_cast<const std::array<std::size_t, sizeof...(Is)>&&>(
        std::array<std::size_t, sizeof...(Is)>{{Fibo<Is>::value...}})[n];
}

template <std::size_t N>
constexpr std::size_t fibo(std::size_t n)
{
    return n < N ?
        fibo(n, make_index_sequence<N>()) :
        throw std::runtime_error("out of bound");
}
} // namespace detail

constexpr std::size_t fibo(std::size_t n)
{
    // 48u is the highest
    return detail::fibo<48u>(n);
}

在 C++14 中,您可以简化一些功能:

template <std::size_t ... Is>
constexpr std::size_t fibo(std::size_t n, index_sequence<Is...>)
{
    constexpr std::array<std::size_t, sizeof...(Is)> fibos{{Fibo<Is>::value...}};
    return fibos[n];
}
于 2014-03-20T14:35:29.093 回答
1

我的想法是递归地将斐波那契序列保存在可变参数模板中,然后将其转换为数组。所有这些都是在编译时完成的。例如,当 n = 5 时,我们有:

F<5>::array
= F<4, 0>::array
= F<3, 0, 1>::array
= F<2, 0, 1, 1>::array
= F<1, 0, 1, 1, 2>::array
= F<0, 0, 1, 1, 2, 3>::array
= { 0, 1, 1, 2, 3 }

然后我们可以在运行时索引数组。

我的 C++14 实现:

#include <cstdint>
#include <array>
#include <iostream>


template<uint64_t n>
struct Helper { static constexpr uint64_t value = Helper<n - 1>::value + Helper<n - 2>::value; };

template<>
struct Helper<0> { static constexpr uint64_t value = 0; };

template<>
struct Helper<1> { static constexpr uint64_t value = 1; };


template<u_int64_t x>
class Fib {
private:
    template<u_int64_t n, u_int64_t...rest>
    struct Get {
        static constexpr std::array<u_int64_t, n + sizeof...(rest)> value = Get<n - 1, rest..., Helper<sizeof...(rest)>::value>::value;
    };

    template<u_int64_t...rest>
    struct Get<0, rest...> {
        static constexpr std::array<u_int64_t, sizeof...(rest)> value{rest...};
    };
public:
    static constexpr std::array<u_int64_t, x> sequence = Get<x>::value;
};

template<u_int64_t x>
constexpr std::array<u_int64_t, x> Fib<x>::sequence;


int main() {
    for (int i = 0; i < 45; i++) std::cout << "F" << i << " = " << Fib<45>::sequence[i] << std::endl;
}
于 2021-05-01T07:48:25.570 回答
0

C(以及大多数情况下是 C++)的基本特性之一是您不需要为不需要的东西付费。

查找表的自动生成并不是编译器需要为您做的事情。即使您需要该功能,也不是其他所有人都需要。

如果您想要一个查找表,请编写一个程序来制作一个。然后在您的程序中使用该数据。

如果您希望在运行时计算值,请不要使用模板元程序,只需使用常规程序来计算值。

于 2009-05-25T23:29:47.243 回答
0

您可以使用预处理器元编程技术生成开关或静态数组。如果复杂性不超过该方法的限制,并且您不希望通过生成代码或数据的额外步骤来扩展您的工具链,那么这是一个不错的决定。

于 2013-05-05T02:39:42.510 回答