std::accumulate使用的一个小细节

今天使用std::accumulate模板函数的时候出现了一个错误，特此记录一下。

#include <iostream>

#include <numeric>

int main()

{

	int LevelColRow[][3] =

	{

		1,  2,	   1,

		2,  4,	   2,

		3,  8,	   4,

		4,  16,	   8,

		5,  32,	   16,

		6,  64,	   32,

		7,  128,   64,

		8,  256,   128,

		9,  512,   256,

		10, 1024,  512,

		11, 2048,  1024,

		12, 4096,  2048,

		13, 8192,  4096,

		14, 16384, 8192,

		15, 32768, 16384,

		16, 65536, 32768 /*,

		17, 131072, 65536/*,

		18, 262144, 131072,

		19, 524288, 262144*/

	};

	double xx = std::accumulate( std::begin( LevelColRow ), std::end( LevelColRow ),

				     0, [] (double val, int lcr[3]) { return(val + double(lcr[1]) * double(lcr[2]) ); } );

	double xx2 = 0;

	for ( auto x : LevelColRow )

	{

		xx2 += double(x[1]) * double(x[2]);

	}

	std::cout << (int64_t) xx << std::endl;

	std::cout << (int64_t) xx2 << std::endl;

	std::cout << "Hello, world!\n";

}

这个代码是用于求这个行列组总共有多少个格子的，但是算出的结果总是为负数。

经过排查，当格子数大于2^31个时候，就出现问题了，这就应该是计算的结果实际是int32的。

查看accumulate(http://en.cppreference.com/w/cpp/algorithm/accumulate)的声明如下：

template< class InputIt, class T >

T accumulate( InputIt first, InputIt last, T init );

template< class InputIt, class T, class BinaryOperation >

T accumulate( InputIt first, InputIt last, T init,

              BinaryOperation op );

那么原因就很简单了，因为我在使用的时候，第三个参数写的是0（默认为int），改为0.0就可以了。

巴特西

std::accumulate使用的一个小细节

最新文章

热门文章