最近看到一个代码，感觉十分有趣， M a r k Mark Mark 一下，来源不太清楚了。

首先说一下代码的功能：通过若干次随机使得随机数的和超过 R A N D _ M A X RAND\_MAX RAND_MAX，求这个随机次数的期望。

代码里通过 1 e 8 1e8 1e8 次操作，最后求随机次数的平均值的方法来得到这个期望。

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

const int MAXN = 1e8;

int main() {
   
    printf("rand: 0-%d:\n", RAND_MAX);
    
    long long cnt1 = 0;
    for (int i = 0; i < MAXN; i++)
        for (long long j = 0; j < RAND_MAX; j += rand())
            cnt1++;
    
    printf("%f\n", cnt1 * 1.0 / MAXN);
    
    long long cnt2;
    for (int k = 0; k < 10; k++)
    {
   
        cnt2 = 0;
        srand((unsigned)time(NULL));
        for (int i = 0; i < MAXN; i++)
            for (long long j = 0; j < RAND_MAX; j += rand())
                cnt2++;
        
        printf("%d: %f\n", k, cnt2 * 1.0 / MAXN);
    }
    
    return 0;
}

这里不管通过什么样的随机种子( s e e d seed seed)，最后得到的期望都是接近于 2.71828 … 2.71828… 2.71828… 的，也就是接近于自然常数 e e e，充分的体现了 r a n d ( ) rand() rand() 的稳定性。

不过，我的概率论学得有些差劲，并不知道如何科学的去解释期望是 e e e 就是对的。求大佬们给一个解释……谢谢。