1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:

Function BM_Add(Num1, Num2: TBigNumber): TBigNumber;
Function BM_Sub(Num1, Num2: TBigNumber): TBigNumber;
Function BM_Divide(Num1, Num2: TBigNumber): TBigNumber;
Function BM_Modulo(Num1, Num2: TBigNumber): TBigNumber;
Function BM_DivMod(Var Num1: TBigNumber; Num2: TBigNumber): TBigNumber; //Returns Int as Num1, Mod as Result
Function BM_Multiply(Num1, Num2: TBigNumber): TBigNumber;
Function BM_MultiplyMod(Num1, Num2, Module: TBigNumber): TBigNumber;
Function BM_Power(Base, Exp: TBigNumber): TBigNumber;
Function BM_PowerMod(Base, Exp, Module: TBigNumber): TBigNumber;
Function BM_Log(Power, Base: TBigNumber): TBigNumber;
Function BM_SquareRoot(Num: TBigNumber): TBigNumber;
Function BM_XRoot(Num, Exp: TBigNumber): TBigNumber;

//Compare Funktionen nach CPU-Flag-Notation
Function BM_CompareC(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1<Num2
Function BM_CompareNC(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1>=Num2
Function BM_CompareNZ(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1<>Num2
Function BM_CompareZ(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1=Num2

//Compare-Funktionen nach Arithmetischer Notation
Function BM_CompareE(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1=Num2
Function BM_CompareNE(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1<>Num2
Function BM_CompareL(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1<Num2
Function BM_CompareLE(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1<=Num2
Function BM_CompareG(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1>Num2
Function BM_CompareGE(Num1, Num2: TBigNumber): Boolean; //Returns True if Num1>=Num2

//ggT und kgV
Function BM_GCD(Num1, Num2: TBigNumber): TBigNumber;
Function BM_LCM(Num1, Num2: TBigNumber): TBigNumber;

Procedure BM_Assign(Var Result: TBigNumber; Const Number: TBigNumber);
Procedure BM_Optimize(Var Result: TBigNumber);
Procedure BM_GCDOptimize(Var Num1, Num2: TBigNumber);

Function BMD_StrToBigNum(Str: String; HexInput: Boolean): TBigNumber;
Function BMD_BigNumToStr(Num: TBigNumber; HexOutput: Boolean): String;
Function BMD_FractionToStr(Num1, Num2: TBigNumber): String;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:

Function BM_PowerMod(Base, Exp, Module: TBigNumber): TBigNumber;
Const
    C: Array[0..7] Of Byte = ($01, $02, $04, $08, $10, $20, $40, $80);
Var
    i, j: Integer;
Begin
    Result := BMD_StrToBigNum('1', false);
    For i := High(Exp) Downto 0 Do
    Begin
        For j := 7 Downto 0 Do
        Begin
            Result := BM_Multiply(Result, Result);
            If (Exp[i] And C[j]) > 0 Then
                Result := BM_Multiply(Result, Base);
            Result := BM_Modulo(Result, Module);
        End;
    End;
End;