📄 ran_c.mht
字号:
From: <由 Microsoft Internet Explorer 5 保存>
Subject:
Date: Wed, 27 Sep 2006 22:48:32 +0800
MIME-Version: 1.0
Content-Type: text/html;
charset="gb2312"
Content-Transfer-Encoding: quoted-printable
Content-Location: http://www.oxbad.com/DSP/maindoc/arithmetic_source/c/RAN.C
X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2900.2962
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=3DContent-Type content=3D"text/html; charset=3Dgb2312">
<META content=3D"MSHTML 6.00.2900.2963" name=3DGENERATOR></HEAD>
<BODY><PRE>/* ran.c - uniform random number generator in [0, 1) */
#define a 16807 /* that is, \(a =3D =
7\sp{5}\) */
#define m 2147483647 /* that is, \(m =3D =
2\sp{31}-1\) */
#define q 127773 /* note, \(q =3D m/a\) =
=3D quotient */
#define r 2836 /* note, \(r =3D m\%a\) =
=3D remainder */
double ran(iseed) /* usage: u =3D =
ran(&iseed); */
long *iseed; /* iseed passed by =
address */
{
*iseed =3D a * (*iseed % q) - r * (*iseed / q); /* update =
seed */
if (*iseed < 0) /* wrap to positive =
values */
*iseed +=3D m;
return (double) *iseed / (double) m;
}
</PRE></BODY></HTML>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -