irbit1.cpp

来自「C++数值算法原书代码」· C++ 代码 · 共 14 行

#include "nr.h"

int NR::irbit1(unsigned long &iseed)
{
	unsigned long newbit;

	newbit =  ((iseed >> 17) & 1)
		^ ((iseed >> 4) & 1)
		^ ((iseed >> 1) & 1)
		^ (iseed & 1);
	iseed=(iseed << 1) | newbit;
	return int(newbit);
}