点击上方“C语言与CPP编程”，选择“关注/置顶/星标公众号”
干货福利，第一时间送达！
大家好，我是飞宇，周末快乐，分享个有趣的编程问题。
这就是一个求阶乘的问题，大家刚刚开始学编程的时候应该都写过这样的程序。
一个求阶乘的问题，还能玩出什么样的花儿来？
我在回答区看到了一个非常有趣的回答，把各种版本的代码都举例了出来，一个比一个还离谱，后面的代码，甚至让怀疑是我学过的代码吗？
（简单利索，深藏功与名）
#include
#include
int main()
{
std::cout << std::tgamma(20 + 1) << std::endl;
}
（语言学家，你懂得，恨不得把所有语法特性都派上用场）
#include
#include
template<std::size_t...I> constexpr auto foo(std::index_sequence) { return ((I+1) * ...); }
int main()
{
std::cout << foo(std::make_index_sequence<20>()) << std::endl;
}
（void main() 有没有嗅到浓厚的历史气息？）
#include
void main(void) {
int i;
long long j;
for(i = 1, j = 1;i <= 20; j *= i++);
printf("%lld", j);
}
（可以说是非常敏捷了）
#include
int main() {
//printf("%d", 1*2*3*4*5*6*7*8*9*10);
printf("%lld", (long long)1*2*3*4*5*6*7*8*9*10*11*12*13*14*15*16*17*18*19*20);
}
（好家伙，一个简单的问题，愣是祭出了接口、继承、虚函数、虚继承、智能指针等一大堆东西出来，这很面向对象）
#include
#include
#include
struct IBaseInterface {
virtual ~IBaseInterface() = 0;
};
inline IBaseInterface::~IBaseInterface() = default;
struct IDataProvider : virtual public IBaseInterface {
virtual int first() = 0;
virtual int last() = 0;
virtual int next(int v) = 0;
};
struct ICalculator : virtual public IBaseInterface {
virtual long long calc(IDataProvider *) = 0;
};
struct IPrinter : virtual public IBaseInterface {
virtual void print(const std::string &) = 0;
};
struct ISerializer : virtual public IBaseInterface {
virtual std::string serialize(long long value) = 0;
};
struct IRunnable : virtual public IBaseInterface {
virtual void run() = 0;
};
class Foo : virtual public IRunnable {
std::shared_ptr m_dp;
std::shared_ptr m_c;
std::shared_ptr m_s;
std::shared_ptr m_p;
public:
Foo(std::shared_ptr dp, std::shared_ptr c, std::shared_ptr s, std::shared_ptr p) : m_dp(std::move(dp)), m_c(std::move(c)), m_s(std::move(s)),m_p(std::move(p)) {}
void run() override { return m_p->print(m_s->serialize(m_c->calc(m_dp.get()))); }
};
class DefaultSerializer : virtual public ISerializer {
public:
std::string serialize(long long value) override { return std::to_string(value); }
};
class StreamPrinter : virtual public IPrinter {
std::ostream &m_os;
public:
explicit StreamPrinter (std::ostream &os) : m_os(os) {}
void print(const std::string &s) override { m_os << s << std::endl; }
};
class MultiplyAccumulateCalculator : virtual public ICalculator {
public:
long long calc(IDataProvider *dp) override {
int i = dp->first();
long long j = i;
do
j *= (i = dp->next(i));
while(i != dp->last());
return j;
}
};
int main() {
struct MyDataProvider : virtual public IDataProvider {
int first() override { return 1; }
int last() override { return 20; }
int next(int v) override { return v+1; }
};
Foo foo(std::make_shared(), std::make_shared(), std::make_shared(), std::make_shared(std::cout));
foo.run();
}
（一看就是精通底层技术的大佬，把CPU拿捏得死死的）
#include
#include
double foo(int x) {
__m128 a = {1.0f, 2.0f, 3.0f, 4.0f};
__m128 b = {4.0f, 4.0f, 4.0f, 4.0f};
__m128 c = {1.0f, 1.0f, 1.0f, 1.0f};
for(int i = 0; i < x / 4; ++i, a = _mm_add_ps(a, b))
c = _mm_mul_ps(c, a);
for(int i = x % 4; i < 4; ++i)
a[i] = 1.0f;
c = _mm_mul_ps(c, a);
return (double)c[0] * (double)c[1] * (double)c[2] * (double)c[3];
}
int main() {
std::cout << foo(20) << std::endl;
}
（能看懂这段代码的，都不是普通人！）
#include
#include
#include
#include
int main() {
std::vector<int> v(std::atoi(std::end(__DATE__) - (__LINE__) / 2) - 1); // 2021年，第六行
std::iota(v.begin(), v.end(), 1);
std::cout << std::accumulate(v.begin(), v.end(), 1ull, std::multiplies<>()) << std::endl;
}
（当年看各种C++框架中，排山倒海一样的宏定义，简直令人发指）
#include
// 由于boost.preprocessor仅提供255以下的整数运算
// 所以使用sequence来 (十位个位)(千位百位)(十万位万位) 的方式来表示大整数。
// 不进位加法：(77)(66)(55) + (44)(33)(22) = (121)(99)(77)
#define PP_ADD_N_N_CARRY_OP(R, DATA, I, ELEM) (BOOST_PP_ADD(BOOST_PP_SEQ_ELEM(I, DATA), ELEM))
#define PP_ADD_N_N_CARRY(SEQ_A, SEQ_B) BOOST_PP_SEQ_FOR_EACH_I(PP_ADD_N_N_CARRY_OP, SEQ_A, SEQ_B)
// 进位加法：(121)(99)(77) = （21)(0)(78)
// 注意SEQ_A的长度要比SEQ_B长
#define PP_ADD_N_N_OP(S, STATE, ELEM_CARRY) \
BOOST_PP_SEQ_PUSH_FRONT( \
BOOST_PP_SEQ_REPLACE(STATE, 0, BOOST_PP_MOD(BOOST_PP_ADD(BOOST_PP_SEQ_HEAD(STATE), ELEM_CARRY), 100)), \
BOOST_PP_DIV(BOOST_PP_ADD(BOOST_PP_SEQ_HEAD(STATE), ELEM_CARRY), 100) \
)
#define PP_ADD_N_N(SEQ_A, SEQ_B) BOOST_PP_SEQ_REVERSE(BOOST_PP_SEQ_FOLD_LEFT(PP_ADD_N_N_OP, BOOST_PP_SEQ_NIL(0), PP_ADD_N_N_CARRY(SEQ_A, SEQ_B)))
// 没什么好说的，X*N = X+X+X+X+X+...+X
#define PP_MUL_N_1_EXP_OP(Z, I, DATA) (DATA)
#define PP_MUL_N_1_EXP(SEQ_N, N) BOOST_PP_REPEAT(N, PP_MUL_N_1_EXP_OP, SEQ_N)
#define PP_MUL_N_1_MYOP(S, STATE, ITEM) PP_ADD_N_N(STATE, ITEM)
#define PP_MUL_N_1_FWD(EXP) BOOST_PP_SEQ_FOLD_LEFT(PP_MUL_N_1_MYOP, BOOST_PP_SEQ_HEAD(EXP), BOOST_PP_SEQ_TAIL(EXP))
#define PP_MUL_N_1(SEQ_N, N) PP_MUL_N_1_FWD(PP_MUL_N_1_EXP(SEQ_N, N))
#define FACT5 PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1((1), 2), 3), 4), 5)
#define FACT10 PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(FACT5, 6), 7), 8), 9), 10)
#define FACT15 PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(FACT10, 11), 12), 13), 14), 15)
#define FACT20 PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(FACT15, 16), 17), 18), 19), 20)
#define FACT25 PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(PP_MUL_N_1(FACT20, 21), 22), 23), 24), 25)
static_assert(false, BOOST_PP_STRINGIZE(FACT10));
（泛型编程，码不惊人死不休）
···
