#include "rand48_impl.h"

#include <stdint.h>

uint64_t __rand48_step(unsigned short* xi, unsigned short* lc) {
    uint64_t a, x;
    x = xi[0] | (unsigned)xi[1] << 16 | (unsigned long long)xi[2] << 32;
    a = lc[0] | (unsigned)lc[1] << 16 | (unsigned long long)lc[2] << 32;
    x = a * x + lc[3];
    xi[0] = x;
    xi[1] = x >> 16;
    xi[2] = x >> 32;
    return x & 0xffffffffffffull;
}