extern void sprand(seed, rngstate) long seed;
long rngstate[57];
{
  int i, ii;
  long last, next, v[55];
  rngstate[0] = last = seed;
  next = 1;
  for (i = 1; i < 55; i++) {
    ii = (21 * i) % 55;
    rngstate[ii] = next;
    next = last - next;
    if (next < 0)
      next += 1000000000;
    last = rngstate[ii];
  }
  (rngstate[57 - 2]) = 0;
  (rngstate[57 - 1]) = 24;
  for (i = 0; i < 165; i++) {
    lprand(rngstate);
  }
}


x86 assembly (clang, -O0):
<sprand>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1e0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x8(%rbp),%rcx
mov    %rcx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
movq   $0x1,-0x28(%rbp)
movl   $0x1,-0x14(%rbp)
cmpl   $0x37,-0x14(%rbp)
jge    11ad <sprand+0x9d>
imul   $0x15,-0x14(%rbp),%eax
mov    $0x37,%ecx
cltd
idiv   %ecx
mov    %edx,-0x18(%rbp)
mov    -0x28(%rbp),%rdx
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
sub    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
cmpq   $0x0,-0x28(%rbp)
jge    118f <sprand+0x7f>
mov    -0x28(%rbp),%rax
add    $0x3b9aca00,%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1141 <sprand+0x31>
mov    -0x10(%rbp),%rax
movq   $0x0,0x1b8(%rax)
mov    -0x10(%rbp),%rax
movq   $0x18,0x1c0(%rax)
movl   $0x0,-0x14(%rbp)
cmpl   $0xa5,-0x14(%rbp)
jge    11f8 <sprand+0xe8>
mov    -0x10(%rbp),%rdi
mov    $0x0,%al
call   1030 <lprand@plt>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11d2 <sprand+0xc2>
add    $0x1e0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sprand>:
push   %r15
push   %r14
push   %rbx
mov    %rsi,%r14
mov    %rdi,(%rsi)
lea    0x150(%rsi),%r8
lea    0xa8(%rsi),%r9
mov    $0x1,%ebx
mov    $0x2a,%ecx
mov    $0x36,%r10d
mov    $0x15,%esi
movabs $0x29e4129e4129e413,%r11
12 e4 29
nopl   0x0(%rax,%rax,1)
mov    %rcx,%rax
mul    %r11
mov    %rcx,%rax
sub    %rdx,%rax
shr    %rax
add    %rdx,%rax
shr    $0x5,%rax
imul   $0xfffffffffffffe48,%rax,%r15
mov    %rsi,%rax
mul    %r11
mov    %rsi,%rax
sub    %rdx,%rax
shr    %rax
add    %rdx,%rax
shr    $0x5,%rax
imul   $0xfffffffffffffe48,%rax,%rax
sub    %rbx,%rdi
lea    0x3b9aca00(%rdi),%rdx
test   %rdi,%rdi
cmovs  %rdx,%rdi
mov    %rbx,(%r9,%rax,1)
mov    %rdi,(%r8,%r15,1)
sub    %rdi,%rbx
lea    0x3b9aca00(%rbx),%rax
test   %rbx,%rbx
cmovs  %rax,%rbx
add    $0x150,%r8
add    $0x2a,%rcx
add    $0x150,%r9
add    $0x2a,%rsi
add    $0xfffffffe,%r10d
jne    1150 <sprand+0x40>
movaps 0xe29(%rip),%xmm0
movups %xmm0,0x1b8(%r14)
mov    $0xa5,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r14,%rdi
xor    %eax,%eax
call   1030 <lprand@plt>
add    $0xffffffff,%ebx
jne    11f0 <sprand+0xe0>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @sprand(i64 %0, i64* %1) {
%3 = alloca i64
%4 = alloca i64*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i64
%8 = alloca i64
%9 = alloca [55 x i64]
store i64 %0, i64* %3
store i64* %1, i64** %4
%10 = load i64, i64* %3
store i64 %10, i64* %7
%11 = load i64*, i64** %4
%12 = getelementptr inbounds i64, i64* %11, i64 0
store i64 %10, i64* %12
store i64 1, i64* %8
store i32 1, i32* %5
br label %13
13: ; preds = %39, %2
%14 = load i32, i32* %5
%15 = icmp slt i32 %14, 55
br i1 %15, label %16, label %42
16: ; preds = %13
%17 = load i32, i32* %5
%18 = mul nsw i32 21, %17
%19 = srem i32 %18, 55
store i32 %19, i32* %6
%20 = load i64, i64* %8
%21 = load i64*, i64** %4
%22 = load i32, i32* %6
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i64, i64* %21, i64 %23
store i64 %20, i64* %24
%25 = load i64, i64* %7
%26 = load i64, i64* %8
%27 = sub nsw i64 %25, %26
store i64 %27, i64* %8
%28 = load i64, i64* %8
%29 = icmp slt i64 %28, 0
br i1 %29, label %30, label %33
30: ; preds = %16
%31 = load i64, i64* %8
%32 = add nsw i64 %31, 1000000000
store i64 %32, i64* %8
br label %33
33: ; preds = %30, %16
%34 = load i64*, i64** %4
%35 = load i32, i32* %6
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds i64, i64* %34, i64 %36
%38 = load i64, i64* %37
store i64 %38, i64* %7
br label %39
39: ; preds = %33
%40 = load i32, i32* %5
%41 = add nsw i32 %40, 1
store i32 %41, i32* %5
br label %13
42: ; preds = %13
%43 = load i64*, i64** %4
%44 = getelementptr inbounds i64, i64* %43, i64 55
store i64 0, i64* %44
%45 = load i64*, i64** %4
%46 = getelementptr inbounds i64, i64* %45, i64 56
store i64 24, i64* %46
store i32 0, i32* %5
br label %47
47: ; preds = %53, %42
%48 = load i32, i32* %5
%49 = icmp slt i32 %48, 165
br i1 %49, label %50, label %56
50: ; preds = %47
%51 = load i64*, i64** %4
%52 = call i32 (i64*, ...) bitcast (i32 (...)* @lprand to i32 (i64*, ...)*)(i64* %51)
br label %53
53: ; preds = %50
%54 = load i32, i32* %5
%55 = add nsw i32 %54, 1
store i32 %55, i32* %5
br label %47
56: ; preds = %47
ret void
}

LLVM IR (-O2):
define void @sprand(i64 %0, i64* %1) {
store i64 %0, i64* %1
br label %3
3: ; preds = %3, %2
%4 = phi i32 [ 1, %2 ], [ %26, %3 ]
%5 = phi i64 [ 1, %2 ], [ %25, %3 ]
%6 = phi i64 [ %0, %2 ], [ %15, %3 ]
%7 = trunc i32 %4 to i16
%8 = mul nuw nsw i16 %7, 21
%9 = urem i16 %8, 55
%10 = zext i16 %9 to i64
%11 = getelementptr inbounds i64, i64* %1, i64 %10
store i64 %5, i64* %11
%12 = sub nsw i64 %6, %5
%13 = icmp slt i64 %12, 0
%14 = add nsw i64 %12, 1000000000
%15 = select i1 %13, i64 %14, i64 %12
%16 = trunc i32 %4 to i16
%17 = mul i16 %16, 21
%18 = add i16 %17, 21
%19 = urem i16 %18, 55
%20 = zext i16 %19 to i64
%21 = getelementptr inbounds i64, i64* %1, i64 %20
store i64 %15, i64* %21
%22 = sub nsw i64 %5, %15
%23 = icmp slt i64 %22, 0
%24 = add nsw i64 %22, 1000000000
%25 = select i1 %23, i64 %24, i64 %22
%26 = add nuw nsw i32 %4, 2
%27 = icmp eq i32 %26, 55
br i1 %27, label %28, label %3
28: ; preds = %3
%29 = getelementptr inbounds i64, i64* %1, i64 55
%30 = bitcast i64* %29 to <2 x i64>*
store <2 x i64> <i64 0, i64 24>, <2 x i64>* %30
br label %31
31: ; preds = %28, %31
%32 = phi i32 [ 0, %28 ], [ %34, %31 ]
%33 = tail call i32 (i64*, ...) bitcast (i32 (...)* @lprand to i32 (i64*, ...)*)(i64* %1)
%34 = add nuw nsw i32 %32, 1
%35 = icmp eq i32 %34, 165
br i1 %35, label %36, label %31
36: ; preds = %31
ret void
}