1. もし入力値が0だったら（具体的には、7FFFFFFFとのANDをとって0だったら）5までとぶ
2. 入力値を浮動小数として解釈し、0と比較することで符号を見る。（上の0かどうか生のビットパターンを見るのに比べ、こちらは実際の小数の値で比較する）
3. もし0未満、符号が負であれば以下の処理を行う。（小数でなく、ビットパターンのままで行う）
   『 現在の値 = ((入力値 >> 24) + 入力値 ) And &H00FFFFFF) 』
4. 次の式で「現在の値」を更新する。
   『 現在の値 = (現在の値 * &HFFFD43FD - &H003C613D) And &H00FFFFFF 』
5. 現在の値を小数に変換、ビットパターン33800000から生成される小数2^-24を掛けて、そのまま関数の返り値として返す。

• 入力値が0なら、上の1の比較によって5の処理までとぶ。
  「現在の値」は前回Rndを呼んだ時と変わらないので、前回と同じように2^-24を掛けて返す。
  そのため、前回と同じ値が返る。
  
  
• 入力値が正なら、「現在の値」に対して4の計算を行った後、2^-24を掛けて返す。
  
  
• 入力値が負なら、「現在の値」を3の式で生成した後、4の計算を行い、2^-24を掛けて返す。

/* ------------------------------------------- */
/* RndTest.c */
#include <windows.h>
static long x=327680;

__declspec(dllexport) float CALLBACK Rnd2(float f)
{
    unsigned u=*(unsigned*)&f;
    if (u&2146828287) {
        if (f<0) x=(u>>24)+u;
        x=x*16598013+12820163&16777215;
    }
    return x*(1.0/16777216.0);
}

__declspec(dllexport) void CALLBACK SRnd1(double d)
{
    unsigned short*s=(unsigned short*)&d;
    x=((s[2]^s[3])<<8)|(x&255);
}

__declspec(dllexport) void CALLBACK SRnd2(float f)
{
    unsigned char*c=(unsigned char*)&f;
    x=((c[0]^c[2])<<8)|((c[1]^c[3])<<16)|(x&255);
}

int WINAPI DllMain(HINSTANCE a,DWORD b,PVOID c) { return TRUE; }
/* ------------------------------------------- */

' ---------------------------------------------
' 標準モジュール
Declare Function Rnd2! Lib "RndTest" (Optional ByVal a! = 1)
Declare Sub SRnd1 Lib "RndTest" (ByVal a#)
Declare Sub SRnd2 Lib "RndTest" (ByVal a!)
Sub Randomize2(Optional ByVal a)
    If IsMissing(a) Then SRnd2 (Timer()) Else SRnd1 (a)
End Sub
' ---------------------------------------------

Public Declare Sub CopyValtoVal Lib "kernel32" Alias "RtlMoveMemory" _
  (Destination As Any, Source As Any, ByVal Length As Long)
Dim num As Variant
Dim first As Boolean

Private Sub CheckFirst()

'初期値の設定
If Not first Then
  num = CDec(327680)
  first = True
End If
End Sub

Function MyRnd(ByVal v As Single) As Single
Dim n As Long, d As Variant

Call CheckFirst

If v <> 0 Then
  If v < 0 Then     '入力が負なら、その値で現在の値を更新。
    
    'singleをlongとして読む
    Call CopyValtoVal(n, v, 4)
    
    'VBにはUnsigned longがないので、Dec型に変換
    dec = CDec(n And &H7FFFFFFF) + CDec(&H10000000) * 8 '最上位ビットは別扱い
    
    num = dec + Int(dec / &H1000000)    'num = dec + (dec>>24)
  End If
  
  'num = (num * &HFD43FD + &HC39EC3) mod &h1000000 をmodナシで行う
  num = num * CDec(&HFD43FD) + &HC39EC3
  num = num - CDec(Int(num / &H1000000)) * &H1000000
End If

MyRnd = num / &H1000000

End Function

Sub MyRandomize(Optional ByVal v As Double)
Dim s As Single, i As Integer, l As Long
Dim b(3) As Byte, inte(3) As Integer

Call CheckFirst

If IsMissing(v) Then   '引数なしの場合

  'バイト配列に変換
  s = Timer
  Call CopyValtoVal(b(0), s, 4)
  
  num = ((b(0) Xor b(2)) * &H100&) Or ((b(1) Xor b(2)) * &H10000) Or (CLng(num) And &HFF)
        
Else
  '引数ありの場合、doubleをintege配列に変換
  Call CopyValtoVal(inte(0), v, 8)
  
  i = inte(2) Xor inte(3)
  ' unsigned int→unsigned long化
  l = (i And &H7FFF&) + IIf(i And &H8000, &H8000&, 0)
  
  num = (l * &H100) Or (CLng(num) And &HFF)
End If

End Sub

パーツ
商品名
店名
値段
(ケース)
TUKUMO ST-5688
TUKUMO ex
7,980円
(VGA)
玄人志向 GFX5700-A128CE
TUKUMO ex
13,650円
(M/B)
AOpen AK77-600N
TUKUMO ex
7,547円
(CPU)
Athlon 2500+ BOX
PC-Success 2号店
9,030円
(OS)
Windows XP Home OEM版
PC-Success 2号店
11,580円
(MEM)
SanMax=Hynix PC3200 CL3 512MB（相性保証付）
OVERTOP
10,390円
(HDD)
Maxtor M6Y160P0 160GB
クレバリー
9,628円
計
69,805円

   ALL  Integer   Float  MemoryR MemoryW MemoryRW  DirectDraw
 45539    77806   94882    26359   37569    43310          60

Rectangle   Text Ellipse  BitBlt    Read   Write    Copy  Drive
    56760  44991    7554      56   51638   55231   41967  C:\100MB

   ALL  Integer   Float  MemoryR MemoryW MemoryRW  DirectDraw
 21198    40726   49479    15606   19447    24663          59

Rectangle   Text Ellipse  BitBlt    Read   Write    Copy  Drive
    37792  19980    5413     997   18927   17885   22346  C:\20MB

'OpenGLのコンテキストハンドル
Dim hRC As Long

'OpenGLの初期化を行う
Function InitOpenGL(ByVal hdc As Long) As Boolean
  Dim pfd As PIXELFORMATDESCRIPTOR
  Dim res&

  'ピクセルフォーマットの指定
  With pfd
    .nSize = Len(pfd)
    .nVersion = 1
    
    'OpenGLでダブルバッファを使用
    .dwFlags = PFD_DRAW_TO_WINDOW Or PFD_SUPPORT_OPENGL Or PFD_DOUBLEBUFFER
    
    'パレットでなくRGBAで。
    .iPixelType = PFD_TYPE_RGBA
    .cColorBits = 32
    
    'RGBの各ビットは指定しなくてもよい
    'アキュムレーションバッファやステンシルバッファも今回は使用しない
    'Zバッファだけ使用する
    .cDepthBits = 16
    .iLayerType = PFD_MAIN_PLANE
  End With
  
  'ピクセルフォーマットが利用可能か？
  res = ChoosePixelFormat(hdc, pfd)
  
  '0なら失敗
  If res = 0 Then Exit Function
    
  'ピクセルフォーマットの指定
  res = SetPixelFormat(hdc, res, pfd)
  If res = 0 Then Exit Function
  
  'OpenGLのコンテキストを生成する
  hRC = wglCreateContext(hdc)
  If hRC = 0 Then Exit Function
  
  '今作成したコンテキストをカレントにする
  Call wglMakeCurrent(hdc, hRC)
  
  '成功
  InitOpenGL = True
  
End Function


Sub TerminateOpenGL()
  'OpenGLの終了

  'コンテキストを選択していない状態にする
  Call wglMakeCurrent(0, 0)
  If hRC > 0 Then Call wglDeleteContext(hRC)

End Sub

Sub DrawFrame()
  Static cnt&
  Dim deg As Single
  
  'メインの描画処理
  cnt = cnt + 1

  'バッファのクリア
  Call glClear(clrColorBufferBit Or clrDepthBufferBit)

  'ビューの指定
  Call glViewport(0, 0, Picture1.ScaleWidth, Picture1.ScaleHeight)

  '描画

  'X軸・Y軸を示す線を描く
  Call glColor3f(1, 1, 1) '白
  Call glLineWidth(2) '線の太さ
  
  Call glBegin(bmLines) '直線群
    Call glVertex2f(-1, 0)
    Call glVertex2f(1, 0)
    Call glVertex2f(0, 1)
    Call glVertex2f(0, -1)
  Call glEnd
  
  '適当に三角形を描画
  Call glBegin(bmTriangles) '三角形
    '最初の頂点
    Call glColor3f(1, 0, 0) '赤
    deg = cnt * 0.1
    Call glVertex2f(0.8 * Cos(deg), 0.8 * Sin(deg))
    
    Call glColor3f(0, 1, 0) '緑
    deg = cnt * 0.1 + 2 / 3 * 3.1415926535 '前の点より120度=2/3π余分に回転
    Call glVertex2f(0.8 * Cos(deg), 0.8 * Sin(deg))
    
    Call glColor3f(0, 0, 1) '青
    deg = cnt * 0.1 + 4 / 3 * 3.1415926535 '前の点より120度=2/3π余分に回転
    Call glVertex2f(0.8 * Cos(deg), 0.8 * Sin(deg))
  Call glEnd
  
  '最後にバッファを入れ替えて画面に反映
  Call SwapBuffers(Picture1.hdc)

End Sub