update

Author: GiggleLiu

src/longlonguint.jl

@@ -191,6 +191,17 @@ function Base.div(x::LongLongUInt{C}, y::LongLongUInt{C}) where {C}
 
     return quotient
 end
+function Base.rem(x::LongLongUInt{C}, y::LongLongUInt{C}) where {C}
+    y == zero(LongLongUInt{C}) && throw(DivideError())
+    x < y && return x
+    x == y && return zero(LongLongUInt{C})
+    
+    # Calculate remainder using the relationship: rem(x, y) = x - div(x, y) * y
+    quotient = div(x, y)
+    remainder = x - quotient * y
+    
+    return remainder
+end
 
 Base.count_ones(x::LongLongUInt) = sum(count_ones, x.content)
 Base.bitstring(x::LongLongUInt) = join(bitstring.(x.content), "")

test/longlonguint.jl

@@ -245,4 +245,11 @@ end
     b3 = bmask(LongLongUInt{1}, 3)
 
     @test hash((b1, b2, b3)) == hash((b1, b2, b3))
+end
+
+@testset "LongLongUInt rem and mod" begin
+    x = LongLongUInt((1234, 213))
+    y = LongLongUInt((0, 17))
+    @test rem(x, y) == rem(BigInt(x), BigInt(y))
+    @test mod(x, y) == mod(BigInt(x), BigInt(y))
 end