开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Internet/网络编程 > 查看源码
XRAND.C
资源名称：mercwsrc.zip　[点击查看]
上传用户：kesirui
上传日期：2007-01-07
资源大小：263k
文件大小：2k
源码类别：

Internet/网络编程

开发平台：

WINDOWS

XRAND.C：源码内容
  1. #include <math.h>
  2. #include <stdio.h>
  3. #include <stdlib.h>
  5. typedef unsigned char bool;
  6. #define TRUE 1
  8. #define TRIAL_MAX 1000000
  9. #define SQUARE(x) ( (x) * (x) )
  11. int number_range ( int from, int to );
  12. int number_mm ( void );
  16. int main( argc, argv )
  17.     int argc;
  18.     char *argv[];
  19. {
  20.     int nBin;
  22.     for ( nBin = 2; nBin < 1000; nBin++ )
  23.     {
  24. int * prgiBin;
  25. int iRand;
  26. int iBin;
  27. int nMean;
  28. int nSigma2;
  29. double dSigma;
  31. prgiBin = calloc( nBin, sizeof(int) );
  32. for ( iRand = 0; iRand < TRIAL_MAX; iRand++ )
  33.     prgiBin[number_range(0, nBin-1)] += 1;
  35. nMean = TRIAL_MAX / nBin;
  36. nSigma2 = 0;
  37. for ( iBin = 0; iBin < nBin; iBin++ )
  38.     nSigma2 += SQUARE( prgiBin[iBin] - nMean );
  39. dSigma   = sqrt( (double) (nSigma2 / (nBin - 1)) );
  40. printf( "%6d %6.4fn", nBin, dSigma / nMean );
  41. free( prgiBin );
  42.     }
  44.     return 0;
  45. }
  49. /*
  50.  * Generate a random number.
  51.  */
  52. int number_range( int from, int to )
  53. {
  54.     int power;
  55.     int number;
  57.     if ( ( to = to - from + 1 ) <= 1 )
  58. return from;
  60.     for ( power = 2; power < to; power <<= 1 )
  61. ;
  63.     while ( ( number = number_mm( ) & (power - 1) ) >= to )
  64. ;
  66.     return from + number;
  67. }
  71. int number_mm( void )
  72. {
  73.     static int rgiState[3+55];
  74.     int *piState;
  75.     int iState1;
  76.     int iState2;
  77.     int iRand;
  79.     piState = &rgiState[0];
  81.     if ( piState[0] == 0 )
  82.     {
  83. int iState;
  85. piState[3+0] = 1;
  86. piState[3+1] = 1;
  87. for ( iState = 3+2; iState < 3+55; iState++ )
  88. {
  89.     piState[iState] = (piState[iState-1] + piState[iState-2])
  90.     & ((1 << 30) - 1);
  91. }
  93. piState[0] = 1;
  94. piState[1] = 3+55 - 55;
  95. piState[2] = 3+55 - 24;
  96.     }
  98.     iState1   = piState[1];
  99.     iState2   = piState[2];
  100.     iRand   = (piState[iState1] + piState[iState2])
  101. & ((1 << 30) - 1);
  102.     piState[iState1] = iRand;
  103.     if ( ++iState1 == 3+55 )
  104. iState1 = 3+0;
  105.     if ( ++iState2 == 3+55 )
  106. iState2 = 3+0;
  107.     piState[1] = iState1;
  108.     piState[2] = iState2;
  109.     return iRand >> 6;
  110. }