sunnyso

¦^´_ 10# sjgau

thanks in advance

sjgau

¦^´_ 11# sunnyso


    integer ªº½d³ò¡A¥u¦³¨ì +32767, ¤£¦æ¨Ï¥Î

long ªº½d³ò¡AÁöµM¦³¨ì +2147483647
¦ý¬O¡A¤U¤@­Ó ¶Ã¼Æ¡Aµ¥©ó «e¤@­Ó¶Ã¼Æ¡A­¼¤W a, a= 7^5= 16807
´NÃz±¼¤F

©Ò¥H¡A¥²¶·¨Ï¥Î double,

single ªººë½T«×¡A¥u¦³ 2^24, ¤j¬ù¬O 10^7, ¤£°÷¨Ï¥Î
©Ò¥H¡A¥²¶·¨Ï¥Î double

·íµM¡A¥²¶· ¤ÏÂдú¸Õ ¨C¤@­Ó ²Ó¸`

stillfish00

¦^´_ 8# sjgau
¦Û¤v¼g¯uªº¥i¥Hµo²{«Ü¦h°ÝÃD¡A¾Ç¨ì«Ü¦h¡A
§Æ±æ¯à©ß¿j¤Þ¥É¡A§Ú¬O±Ä¥ÎMicrosoft Visual Basic ªº¨º²Õ°Ñ¼Æ¡A¥é³yrnd¨ç¼Æ¤@¼Ë¶Ç¦^single¡C

1. Function myRnd(Optional Number) As Single
   
2.   Dim tmp As Double
   
3.   'LCG Algorithm¡G
   
4.   '  X(n+1) = (a * X(n) + c) mod m
   
5.   Const a As Long = 1140671485
   
6.   Const c As Long = 12820163
   
7.   Const m As Long = 16777216  '2^24
   
8.   Const x0 As Long = 1     '®Ú¾Ú§A¨Ï¥Î­þ²Õ°Ñ¼Æ¡Ax0¦³®É¦³»Ý­n²Å¦Xªº±ø¥ó¡C
   
9.   
   
10.   Static seed As Variant
    
11.   '§Q¥ÎIsEmpty§PÂ_¡A¨Óªì©l¤Æ¤@ staticÅܼƬ°«D¹s­È¡C
    
12.   If IsEmpty(seed) Then seed = x0
    
13.   
    
14.   '§Q¥ÎIsMissing§PÂ_ optional °Ñ¼Æ¬O§_¶Ç¤J¡A­Y¨ç¼Æ¦³¶Ç¤J°Ñ¼Æ®É¡A³]¬°ºØ¤l¡C
    
15.   If Not IsMissing(Number) Then seed = Number
    
16.    
    
17.   'myRnd = (a * seed + c) Mod m
    
18.   '(a * seed + c)¹L¤j®É¡Amod ·|µo¥Í·¸¦ì¡A¥u¯à¥Î¨ä¥L¤èªk¨D¡C
    
19.   
    
20.   tmp = a * seed + c  'double
    
21.   seed = tmp - m * Fix(tmp / m)   '§Y­pºâ (a * seed + c) Mod m¡A¦¹¬°²£¥Íªº¶Ã¼Æ­È¡A¨Ã³]¬°¤U¦¸ºØ¤l
    
22.   myRnd = seed / m  '¨Ï¶Ç¦^¶Ã¼Æ­È¤¶©ó0¨ì1¤§¶¡
    
23. End Function

½Æ»s¥N½X

stillfish00

¦^´_ 13# stillfish00
¹ê´ú«á¤´¤£º¡·N¡A©I¥sm¦¸«á¡A¤£·|¦^¨ì²Ä¤@²Õ²£¥ÍªºÀH¾÷¼Æ¡A
¥D¦]¬O tmp = a * seed + c §Y«K¬O double ¼Æ¤´¦s¦³ºë½T«×»~®t...   

LongLong¤S¥u¦³¦b64¦ì¤¸¥­¥x¤W¥i¥H¨Ï¥Î

c_c_lai

¦^´_ 14# stillfish00
³o¬O§Ú´ú¸Õªºµ²ªG¡G
A.  rng = myRnd(sngVal)

¤U¸ü (18.04 KB)

2013-11-1 09:09

B.  rng = myRnd(sngVal, 2 ^ 32)

¤U¸ü (17.97 KB)

2013-11-1 09:09

²¤¥[¼W­×³¡¤À¡G

1. Option Explicit
   
2. 
   
3. Sub Ex()
   
4.     Dim sngVal As Single, rng As Range, xi As Long
   
5.    
   
6.     With Sheets(1)
   
7.         .Cells.Clear
   
8.         
   
9.         For xi = 1 To 5
   
10.             Set rng = .Range(Chr(64 + xi) & "1")
    
11.             For sngVal = 1 To 20
    
12.                 '  rng = myRnd(sngVal)
    
13.                 rng = myRnd(sngVal, 2 ^ 32)
    
14.                 Set rng = rng.Offset(1)
    
15.             Next sngVal
    
16.         Next xi
    
17.         .Columns.AutoFit
    
18.     End With
    
19. End Sub
    
20. 
    
21. Function myRnd(Optional Number, Optional m As Double = 2 ^ 24) As Single
    
22.     Dim tmp As Double
    
23.     '  LCG Algorithm¡G
    
24.     '      X(n+1) = (a * X(n) + c) mod m
    
25.     Const a As Long = 1140671485
    
26.     Const c As Long = 12820163
    
27.     '  Const m As Long = 16777216    '  2^24
    
28.     Const x0 As Long = 1         '  ®Ú¾Ú§A¨Ï¥Î­þ²Õ°Ñ¼Æ¡Ax0¦³®É¦³»Ý­n²Å¦Xªº±ø¥ó¡C

½Æ»s¥N½X

stillfish00

¦^´_ 15# c_c_lai
¨ä¹ê§Úªº·N«ä¬O³o¼Ë¡G

1. Sub EX()
   
2.   Dim a As Single, b As Single
   
3.   Dim i
   
4.   
   
5.   a = Rnd
   
6.   For i = 2 To 16777216
   
7.     Rnd
   
8.   Next
   
9.   b = Rnd
   
10.   Debug.Print a, b
    
11.   
    
12.   a = myRnd
    
13.   For i = 2 To 16777216
    
14.     myRnd
    
15.   Next
    
16.   b = myRnd
    
17.   Debug.Print a, b
    
18.   
    
19. End Sub

½Æ»s¥N½X

¹³rnd¤@¼Ë¡A³om¦¸¨úªº¶Ã¼Æ³£·|¤£¤@¼Ë¡A
¦ýrnd¨ú§¹m¦¸«á¡A¦A¨ú¤U¤@¦¸®É¡A·|¦^¨ì²Ä¤@¦¸¨úªº­È¡A¦A¤@¦¸´`Àô¡C
´N¹³LCG¸Ì­±´£¨ìªº¥Í©R¶g´Á

sjgau

¤ñ¸û¦nªº ¶Ã¼Æ²£¥Í¾¹ªº ±ø¥ó
¤@ ²£¥Íªº­pºâ³t«× ¶V§Ö¶V¦n

¤G ¶Ã¼Æªº½d³ò©M ­Ó¼Æ¡A¶V¤j¶V¦n


¤T ºÉ¥i¯àªº ±µªñ§¡¤Ã¤À°t

¥| ¦b²£¥Í¤U¤@­Ó ¬Û¦Pªº¶Ã¼Æ¤§«e¡A­n§â©Ò¦³ªº¶Ã¼Æ³£¶]¹L¤@¹M


¤­ ´`Àôªº¶g´Á­n ©T©w

¤» ´`Àôªº ¶g´Á­n ¶V¤j ¶V¦n

sunnyso

¦^´_ 10# sjgau

¤£ª¾¤j¤jªº¤é¥»¤§®È¦p¦ó¡H

¦³ªÅ ³Â·Ð«ü¾É³o­Ó¶Ã¼Æ

ÁÂÁÂ

sjgau

¦^´_ 18# sunnyso


    11/8 ¦^°ê
³o¨â¤Ñ¡A¦b¥x¤j¦³±Â½Òªº ¦æµ{
»Ý­n ©ú¤Ñ ¤~¯à ¶}©l

sjgau

¦^´_ 19# sjgau


    http://www.csie.ntu.edu.tw/train/teacher_display.php?num=18
§Ú¦b ¥x¤j£¾§c ½Òµ{¤¶²Ðºô­¶