rand.sagews

def pseudo_rand(seed, n):
    results = [Integer(seed)]
    for i in range(1, 31):
        results.append((16807*results[-1]) % 2147483647)
    results.append(results[0])
    results.append(results[1])
    results.append(results[2])
    for i in range(34, 344):
        results.append((results[-31]+results[-3]) % 2**32)
    output = []
    for i in range(n):
        results.append((results[-31]+results[-3]) % 2**32)
        output.append(results[-1] >> 1)
    return output