c++rand函数源代码

glibc的rand源码如下所示：

int

__random_r (buf, result)

            struct random_data *buf

            int32_t *result

{

  int32_t *state

  if (buf == NULL || result == NULL)

    goto fail

  state = buf->state

  if (buf->rand_type == TYPE_0)

    {

      int32_t val = state[0]

      val = ((state[0] * 1103515245) + 12345) & 0x7fffffff

      state[0] = val

      *result = val

    }

  else

    {

      int32_t *fptr = buf->fptr

      int32_t *rptr = buf->rptr

      int32_t *end_ptr = buf->end_ptr

      int32_t val

      val = *fptr += *rptr

      /* Chucking least random bit.  */

      *result = (val >> 1) & 0x7fffffff

      ++fptr

      if (fptr >= end_ptr)

        {

          fptr = state

          ++rptr

        }

      else

        {

          ++rptr

          if (rptr >= end_ptr)

            rptr = state

        }

      buf->fptr = fptr

      buf->rptr = rptr

    }

  return 0

 fail:

  __set_errno (EINVAL)

  return -1

}

参考： https://github.com/lattera/glibc/blob/master/stdlib/random_r.c

return一句的功能只是把得数的绝对值控制在0~32767之间，关键是next = next * 1103515245 + 12345一句。next变量是静态的，变化以后的值是不消失的，下一次还可用，所以每执行一次next = next * 1103515245 + 12345就获得一个新值，这个新值被return取整除、取余除后控制在0~32767之间就形成了0~32767的随机数。比如第一次，next=1*1103515245 + 12345=1103527590，经return中的对65536取整、对32768取余后就是16838；若再来一次，next已经是1103527590了，那么next = next * 1103515245 + 12345就为next = 1103527590 * 1103515245 + 12345=1217759518843121895；但这个数已经溢出了，实际上表示成了-1770082073，经return中取整取余后返回的就是-27009；可以算出来，再下一次是10113……不过，你这个函数并不能产生真正意义上的“随机”数，因为作为基数的1103515245在每次开始时是不变的，所以只能得出同样的序列，即每次执行都产生16838、-27009、10113……这样一组无限多的数。所以实际的C随机函数rand()的形参并不是void而是一个长整型变量，通常来调用实时时间函数获取实时时间值来得到，因为时间是每时每刻都在变化的，所以充当"1103515245”角色的数就每调用一次的值都不同。这样每次执行就都能得到不重复的序列。至于12345，我想随便取个数都行！仅供参考……

---