Algorithms

// NOTE: The order of the bits is inconsequential, as long as each string of 4 bits comes from a different 64-bit block
output = (bit[ 0] << 255) | (bit[ 64] << 254) | (bit[128] << 253) | (bit[192] << 252)
       | (bit[ 1] << 251) | (bit[ 65] << 250) | (bit[129] << 249) | (bit[193] << 248)
       | ...
       | (bit[63] <<   3) | (bit[127] <<   2) | (bit[191] <<   1) | (bit[255])

// Interleave
__uint16 *data   = new __uint16[32 / sizeof(__uint16)]; // Data to interleave
__uint64 *result = new __uint64[32 / sizeof(__uint64)]();

for (int resIdx = 0; resIdx < 4; ++resIdx) {
    for (int dataIdx = resIdx; dataIdx < 16; dataIdx += 4) {
        __uint64 temp = data[dataIdx];
        temp |= temp << 24;
        temp |= temp << 12;
        temp &= 0x000F000F000F000FULL;
        temp |= temp <<  6;
        temp |= temp <<  3;
        temp &= 0x1111111111111111ULL;
        result[resIdx] |= temp << (dataIdx & 3);
    }
}

// Deinterleave
__uint64 *data   = new __uint64[32 / sizeof(__uint64)]; // Data to deinterleave
__uint16 *result = new __uint16[32 / sizeof(__uint16)];
for (int dataIdx = 0; dataIdx < 4; ++dataIdx) {
    for (int resIdx = dataIdx; resIdx < 16; resIdx += 4) {
        __uint64 temp = data[dataIdx] >> (resIdx >> 2);
        temp &= 0x1111111111111111ULL;
        temp |= temp >>  3;
        temp |= temp >>  6;
        temp &= 0x000F000F000F000FULL;
        temp |= temp >> 12;
        temp |= temp >> 24;
        result[resIdx] = (__uint16)temp;
    }
}

Imports System.Numerics

Public Class BlumBlumShub

    Private _x As BigInteger
    Private _p As BigInteger
    Private _q As BigInteger
    Private _m As BigInteger

    Private _pow As New BigInteger(2)

    Public Function NextNumber() As BigInteger
        _x = BigInteger.ModPow(_x, _pow, _m)
        Return _x
    End Function

    Public Sub New(ByVal seed As BigInteger)
        _p = BigInteger.Parse("32416190071")
        _q = BigInteger.Parse("32416185031")
        _x = seed
        _m = BigInteger.Multiply(_p, _q)
    End Sub

End Class

Public Class RandomBitStream

    Private _b As BlumBlumShub

    Public Function ReadByte() As Byte
        Dim num As Byte = 0
        For i = 1 To 8
            num += Math.Pow(i, 2) * If(_b.NextNumber().IsEven, 1, 0)
        Next
        Return num
    End Function

    Public Sub New(ByVal seed As BigInteger)
        _b = New BlumBlumShub(seed)
    End Sub

End Class

Public Class BlumXor

    Private _bitSrc As RandomBitStream
    Private _key As BigInteger

    Public Sub Cipher(ByRef message As Byte())
        _bitSrc = New RandomBitStream(_key)
        For i = 0 To message.Length - 1
            message(i) = _bitSrc.ReadByte() Xor message(i)
        Next
    End Sub

    Public Function ByteToStr(ByVal inByte As Byte) As String
        Return inByte.ToString().PadLeft(3, "0")
    End Function

    Public Sub New(ByVal keyStr As String)
        Dim encoder As New System.Text.ASCIIEncoding()
        Dim keyBytes As Byte() = encoder.GetBytes(keyStr)
        _key = BigInteger.Parse(String.Join("", Array.ConvertAll(Of Byte, String)(keyBytes, New System.Converter(Of Byte, String)(AddressOf ByteToStr))))
    End Sub

End Class