我正在尝试使用 C++ 计算一个系列。该系列是：（
对于那些想知道的人）
我的代码如下：
#include <iostream>
#include <fstream>
#include <cmath> // exp
#include <iomanip> //setprecision, setw
#include <limits> //numeric_limits (http://en.cppreference.com/w/cpp/types/numeric_limits)
long double SminOneCenter(long double gamma)
{
using std::endl; using std::cout;
long double result=0.0l;
for (long double k = 1; k < 1000 ; k++)
{
if(isinf(pow(1.0l+pow(gamma,k),6.0l/4.0l)))
{
cout << "infinity for reached for gamma equals: " << gamma << "value of k: " << k ;
cout << "maximum allowed: " << std::numeric_limits<long double>::max()<< endl;
break;
}
// CAS PAIR: -1^n = 1
if ((int)k%2 == 0)
{
result += pow(4.0l*pow(gamma,k),3.0l/4.0l) /(pow(1+pow(gamma,k)),6.0l/4.0l);
}
// CAS IMPAIR:-1^n = -1
else if ((int)k%2!=0)
{
result -= pow(4.0l*pow(gamma,k),3.0l/4.0l) /(pow(1+pow(gamma,k)),6.0l/4.0l);
//if (!isinf(pow(k,2.0l)*zeta/2.0l))
}
// cout << result << endl;
}
return 1.0l + 2.0l*result;
}
输出将是，例如gamma = 1.7：
伽玛达到无穷大等于： 1.7值k： 892
long double由 STL 提供的 a 可以表示的最大值numeric_limits是：1.18973e+4932。
但是(1+1.7^892)= 2.19.... × 10^308，它远低于10^4932，因此不应将其视为无穷大。
如果我的代码没有错（但很可能是错的），谁能告诉我为什么所讨论的代码在不应该时会计算为无穷大？