Modulo

From testwiki
Revision as of 10:00, 25 January 2025 by imported>Infiternity (Gauss didn't introduce modular arithmetic, he was behind the modern approach to it.)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:About Template:Short description In computing, the modulo operation returns the remainder or signed remainder of a division, after one number is divided by another, called the modulus of the operation.

Given two positive numbers Template:Math and Template:Math, Template:Math modulo Template:Math (often abbreviated as Template:Math) is the remainder of the Euclidean division of Template:Math by Template:Math, where Template:Math is the dividend and Template:Math is the divisor.[1]

For example, the expression "5 mod 2" evaluates to 1, because 5 divided by 2 has a quotient of 2 and a remainder of 1, while "9 mod 3" would evaluate to 0, because 9 divided by 3 has a quotient of 3 and a remainder of 0.

Although typically performed with Template:Math and Template:Math both being integers, many computing systems now allow other types of numeric operands. The range of values for an integer modulo operation of Template:Math is 0 to Template:Math. Template:Math mod 1 is always 0.

When exactly one of Template:Math or Template:Math is negative, the basic definition breaks down, and programming languages differ in how these values are defined.

Variants of the definition

In mathematics, the result of the modulo operation is an equivalence class, and any member of the class may be chosen as representative; however, the usual representative is the least positive residue, the smallest non-negative integer that belongs to that class (i.e., the remainder of the Euclidean division).[2] However, other conventions are possible. Computers and calculators have various ways of storing and representing numbers; thus their definition of the modulo operation depends on the programming language or the underlying hardware.

In nearly all computing systems, the quotient Template:Math and the remainder Template:Math of Template:Math divided by Template:Math satisfy the following conditions: Template:NumBlk

This still leaves a sign ambiguity if the remainder is non-zero: two possible choices for the remainder occur, one negative and the other positive; that choice determines which of the two consecutive quotients must be used to satisfy equation (1). In number theory, the positive remainder is always chosen, but in computing, programming languages choose depending on the language and the signs of Template:Math or Template:Math.Template:Efn Standard Pascal and ALGOL 68, for example, give a positive remainder (or 0) even for negative divisors, and some programming languages, such as C90, leave it to the implementation when either of Template:Math or Template:Math is negative (see the table under Template:Section link for details). Some systems leave Template:Math modulo 0 undefined, though others define it as Template:Math.

Template:Bulleted list

If both the dividend and divisor are positive, then the truncated, floored, and Euclidean definitions agree. If the dividend is positive and the divisor is negative, then the truncated and Euclidean definitions agree. If the dividend is negative and the divisor is positive, then the floored and Euclidean definitions agree. If both the dividend and divisor are negative, then the truncated and floored definitions agree.

As described by Leijen, Template:Quote

However, truncated division satisfies the identity (a)/b=(a/b)=a/(b).[3]

Notation

Template:About

Some calculators have a Template:Math function button, and many programming languages have a similar function, expressed as Template:Math, for example. Some also support expressions that use "%", "mod", or "Mod" as a modulo or remainder operator, such as Template:Code or Template:Code.

For environments lacking a similar function, any of the three definitions above can be used.

Common pitfalls

When the result of a modulo operation has the sign of the dividend (truncated definition), it can lead to surprising mistakes.

For example, to test if an integer is odd, one might be inclined to test if the remainder by 2 is equal to 1:

bool is_odd(int n) {
    return n % 2 == 1;
}

But in a language where modulo has the sign of the dividend, that is incorrect, because when Template:Math (the dividend) is negative and odd, Template:Math mod 2 returns −1, and the function returns false.

One correct alternative is to test that the remainder is not 0 (because remainder 0 is the same regardless of the signs):

bool is_odd(int n) {
    return n % 2 != 0;
}

Or with the binary arithmetic:

bool is_odd(int n) {
    return n & 1;
}

Performance issues

Modulo operations might be implemented such that a division with a remainder is calculated each time. For special cases, on some hardware, faster alternatives exist. For example, the modulo of powers of 2 can alternatively be expressed as a bitwise AND operation (assuming Template:Math is a positive integer, or using a non-truncating definition):

x % 2n == x & (2n - 1)

Examples:

Template:Code
Template:Code
Template:Code

In devices and software that implement bitwise operations more efficiently than modulo, these alternative forms can result in faster calculations.[4]

Compiler optimizations may recognize expressions of the form Template:Code where Template:Code is a power of two and automatically implement them as Template:Code, allowing the programmer to write clearer code without compromising performance. This simple optimization is not possible for languages in which the result of the modulo operation has the sign of the dividend (including C), unless the dividend is of an unsigned integer type. This is because, if the dividend is negative, the modulo will be negative, whereas Template:Code will always be positive. For these languages, the equivalence x % 2n == x < 0 ? x | ~(2n - 1) : x & (2n - 1) has to be used instead, expressed using bitwise OR, NOT and AND operations.

Optimizations for general constant-modulus operations also exist by calculating the division first using the constant-divisor optimization.

Properties (identities)

Template:See also Some modulo operations can be factored or expanded similarly to other mathematical operations. This may be useful in cryptography proofs, such as the Diffie–Hellman key exchange. The properties involving multiplication, division, and exponentiation generally require that Template:Math and Template:Math are integers.

In programming languages

Modulo operators in various programming languages
Language Operator Integer Floating-point Definition
ABAP Template:Code Template:Yes Template:Yes Euclidean
ActionScript Template:Code Template:Yes Template:No Truncated
Ada Template:Code Template:Yes Template:No Floored[5]
Template:Code Template:Yes Template:No Truncated[5]
ALGOL 68 Template:Code, Template:Code Template:Yes Template:No Euclidean
AMPL Template:Code Template:Yes Template:No Truncated
APL |Template:Efn Template:Yes Template:Yes Floored
AppleScript Template:Code Template:Yes Template:No Truncated
AutoLISP Template:Code Template:Yes Template:No Truncated
AWK Template:Code Template:Yes Template:No Truncated
bash Template:Code Template:Yes Template:No Truncated
BASIC Template:Code Template:Yes Template:No Varies by implementation
bc Template:Code Template:Yes Template:No Truncated
CTemplate:BreakC++ Template:Code, Template:Code Template:Yes Template:No TruncatedTemplate:Efn
Template:Code (C)Template:BreakTemplate:Code (C++) Template:No Template:Yes Truncated[6]
Template:Code (C)Template:BreakTemplate:Code (C++) Template:No Template:Yes Rounded
C# Template:Code Template:Yes Template:Yes Truncated
Template:Code Template:No Template:Yes Rounded[7]
Clarion Template:Code Template:Yes Template:No Truncated
Clean Template:Code Template:Yes Template:No Truncated
Clojure Template:Code Template:Yes Template:No Floored[8]
Template:Code Template:Yes Template:No Truncated[9]
COBOL Template:Code Template:Yes Template:No Floored[10]
Template:Code Template:Yes Template:Yes Truncated[10]
CoffeeScript Template:Code Template:Yes Template:No Truncated
Template:Code Template:Yes Template:No Floored[11]
ColdFusion Template:Code, Template:Code Template:Yes Template:No Truncated
Common Intermediate Language Template:Code (signed) Template:Yes Template:Yes Truncated[12]
Template:Code (unsigned) Template:Yes Template:No Template:N/A
Common Lisp Template:Code Template:Yes Template:Yes Floored
Template:Code Template:Yes Template:Yes Truncated
Crystal Template:Code, Template:Code Template:Yes Template:Yes Floored
Template:Code Template:Yes Template:Yes Truncated
CSS Template:Code Template:Yes Template:Yes Floored[13]
Template:Code Template:Yes Template:Yes Truncated[14]
D Template:Code Template:Yes Template:Yes Truncated[15]
Dart Template:Code Template:Yes Template:Yes Euclidean[16]
Template:Code Template:Yes Template:Yes Truncated[17]
Eiffel Template:Code Template:Yes Template:No Truncated
Elixir Template:Code Template:Yes Template:No Truncated[18]
Template:Code Template:Yes Template:No Floored[19]
Elm Template:Code Template:Yes Template:No Floored[20]
Template:Code Template:Yes Template:No Truncated[21]
Erlang Template:Code Template:Yes Template:No Truncated
Template:Code Template:No Template:Yes Truncated (same as C)[22]
Euphoria Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
F# Template:Code Template:Yes Template:Yes Truncated
Template:Code Template:No Template:Yes Rounded[7]
Factor Template:Code Template:Yes Template:No Truncated
FileMaker Template:Code Template:Yes Template:No Floored
Forth Template:Code Template:Yes Template:No Implementation defined
Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Fortran Template:Code Template:Yes Template:Yes Truncated
Template:Code Template:Yes Template:Yes Floored
Frink Template:Code Template:Yes Template:No Floored
Full BASIC Template:Code Template:Yes Template:Yes Floored[23]
Template:Code Template:Yes Template:Yes Truncated[24]
GLSL Template:Code Template:Yes Template:No Undefined[25]
Template:Code Template:No Template:Yes Floored[26]
GameMaker Studio (GML) Template:Code, Template:Code Template:Yes Template:No Truncated
GDScript (Godot) Template:Code Template:Yes Template:No Truncated
Template:Code Template:No Template:Yes Truncated
Template:Code Template:Yes Template:No Euclidean
Template:Code Template:No Template:Yes Euclidean
Go Template:Code Template:Yes Template:No Truncated[27]
Template:Code Template:No Template:Yes Truncated[28]
Template:Code Template:Yes Template:No Euclidean[29]
Template:Code Template:Yes Template:No Truncated[30]
Groovy Template:Code Template:Yes Template:No Truncated
Haskell Template:Code Template:Yes Template:No Floored[31]
Template:Code Template:Yes Template:No Truncated[31]
Template:Code (GHC) Template:No Template:Yes Floored
Haxe Template:Code Template:Yes Template:No Truncated
HLSL Template:Code Template:Yes Template:Yes Undefined[32]
J |Template:Efn Template:Yes Template:No Floored
Java Template:Code Template:Yes Template:Yes Truncated
Template:Code Template:Yes Template:No Floored
JavaScriptTemplate:BreakTypeScript Template:Code Template:Yes Template:Yes Truncated
Julia Template:Code Template:Yes Template:Yes Floored[33]
Template:Code, Template:Code Template:Yes Template:Yes Truncated[34]
Kotlin Template:Code, Template:Code Template:Yes Template:Yes Truncated[35]
Template:Code Template:Yes Template:Yes Floored[36]
ksh Template:Code Template:Yes Template:No Truncated (same as POSIX sh)
Template:Code Template:No Template:Yes Truncated
LabVIEW Template:Code Template:Yes Template:Yes Truncated
LibreOffice Template:Code Template:Yes Template:No Floored
Logo Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Lua 5 Template:Code Template:Yes Template:Yes Floored
Lua 4 Template:Code Template:Yes Template:Yes Truncated
Liberty BASIC Template:Code Template:Yes Template:No Truncated
Mathcad Template:Code Template:Yes Template:No Floored
Maple Template:Code (by default), Template:Code Template:Yes Template:No Euclidean
Template:Code Template:Yes Template:No Rounded
Template:Code Template:Yes Template:Yes Rounded
Mathematica Template:Code Template:Yes Template:No Floored
MATLAB Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Maxima Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Maya Embedded Language Template:Code Template:Yes Template:No Truncated
Microsoft Excel Template:Code Template:Yes Template:Yes Floored
Minitab Template:Code Template:Yes Template:No Floored
Modula-2 Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
MUMPS Template:Code Template:Yes Template:No Floored
Netwide Assembler (NASM, NASMX) Template:Code, Template:Code (unsigned) Template:Yes Template:No Template:N/A
Template:Code (signed) Template:Yes Template:No Implementation-defined[37]
Nim Template:Code Template:Yes Template:No Truncated
Oberon Template:Code Template:Yes Template:No Floored-likeTemplate:Efn
Objective-C Template:Code Template:Yes Template:No Truncated (same as C99)
Object Pascal, Delphi Template:Code Template:Yes Template:No Truncated
OCaml Template:Code Template:Yes Template:No Truncated[38]
Template:Code Template:No Template:Yes Truncated[39]
Occam Template:Code Template:Yes Template:No Truncated
Pascal (ISO-7185 and -10206) Template:Code Template:Yes Template:No Euclidean-likeTemplate:Efn
Perl Template:Code Template:Yes Template:No FlooredTemplate:Efn
Template:Code Template:No Template:Yes Truncated
Phix Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
PHP Template:Code Template:Yes Template:No Truncated[40]
Template:Code Template:No Template:Yes Truncated[41]
PIC BASIC Pro Template:Code Template:Yes Template:No Truncated
PL/I Template:Code Template:Yes Template:No Floored (ANSI PL/I)
PowerShell Template:Code Template:Yes Template:No Truncated
Programming Code (PRC) Template:Code Template:Yes Template:No Undefined
Progress Template:Code Template:Yes Template:No Truncated
Prolog (ISO 1995) Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
PureBasic Template:Code, Template:Code Template:Yes Template:No Truncated
PureScript Template:Code Template:Yes Template:No Euclidean[42]
Pure Data Template:Code Template:Yes Template:No Truncated (same as C)
Template:Code Template:Yes Template:No Floored
Python Template:Code Template:Yes Template:Yes Floored
Template:Code Template:No Template:Yes Truncated
Template:Code Template:No Template:Yes Rounded
Q# Template:Code Template:Yes Template:No Truncated[43]
R Template:Code Template:Yes Template:Yes Floored[44]
Racket Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Raku Template:Code Template:No Template:Yes Floored
RealBasic Template:Code Template:Yes Template:No Truncated
Reason Template:Code Template:Yes Template:No Truncated
Rexx Template:Code Template:Yes Template:Yes Truncated
RPG Template:Code Template:Yes Template:No Truncated
Ruby Template:Code, Template:Code Template:Yes Template:Yes Floored
Template:Code Template:Yes Template:Yes Truncated
Rust Template:Code Template:Yes Template:Yes Truncated
Template:Code Template:Yes Template:Yes Euclidean[45]
SAS Template:Code Template:Yes Template:No Truncated
Scala Template:Code Template:Yes Template:Yes Truncated
Scheme Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Scheme R6RS Template:Code Template:Yes Template:No Euclidean[46]
Template:Code Template:Yes Template:No Rounded[46]
Template:Code Template:No Template:Yes Euclidean
Template:Code Template:No Template:Yes Rounded
Scratch Template:Code Template:Yes Template:Yes Floored
Seed7 Template:Code Template:Yes Template:Yes Floored
Template:Code Template:Yes Template:Yes Truncated
SenseTalk Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
[[POSIX shell|Template:Code (POSIX)]] (includes bash, mksh, &c.) Template:Code Template:Yes Template:No Truncated (same as C)[47]
Smalltalk Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Snap! Template:Code Template:Yes Template:No Floored
Spin Template:Code Template:Yes Template:No Floored
Solidity Template:Code Template:Yes Template:No Truncated[48]
SQL (SQL:1999) Template:Code Template:Yes Template:No Truncated
SQL (SQL:2011) Template:Code Template:Yes Template:No Truncated
Standard ML Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
Template:Code Template:No Template:Yes Truncated
Stata Template:Code Template:Yes Template:No Euclidean
Swift Template:Code Template:Yes Template:No Truncated[49]
Template:Code Template:No Template:Yes Rounded[50]
Template:Code Template:No Template:Yes Truncated[51]
Tcl Template:Code Template:Yes Template:No Floored
Template:Code Template:No Template:Yes Truncated (as C)
tcsh Template:Code Template:Yes Template:No Truncated
Torque Template:Code Template:Yes Template:No Truncated
Turing Template:Code Template:Yes Template:No Floored
Verilog (2001) Template:Code Template:Yes Template:No Truncated
VHDL Template:Code Template:Yes Template:No Floored
Template:Code Template:Yes Template:No Truncated
VimL Template:Code Template:Yes Template:No Truncated
Visual Basic Template:Code Template:Yes Template:No Truncated
WebAssembly Template:Code, Template:Code (unsigned) Template:Yes Template:No Template:N/A[52]
Template:Code, Template:Code (signed) Template:Yes Template:No Truncated[52]
x86 assembly Template:Code Template:Yes Template:No Truncated
XBase++ Template:Code Template:Yes Template:Yes Truncated
Template:Code Template:Yes Template:Yes Floored
Zig Template:Code,

Template:Code, Template:Code

Template:Yes Template:Yes Truncated[53]
Z3 theorem prover Template:Code, Template:Code Template:Yes Template:No Euclidean

In addition, many computer systems provide a Template:Code functionality, which produces the quotient and the remainder at the same time. Examples include the x86 architecture's Template:Code instruction, the C programming language's Template:Code function, and Python's Template:Code function.

Generalizations

Modulo with offset

Sometimes it is useful for the result of Template:Mvar modulo Template:Mvar to lie not between 0 and Template:Math, but between some number Template:Mvar and Template:Math. In that case, Template:Mvar is called an offset and Template:Math is particularly common.

There does not seem to be a standard notation for this operation, so let us tentatively use Template:Math. We thus have the following definition:[54] Template:Math just in case Template:Math and Template:Math. Clearly, the usual modulo operation corresponds to zero offset: Template:Math.

The operation of modulo with offset is related to the floor function as follows:

amoddn=anadn.

To see this, let x=anadn. We first show that Template:Math. It is in general true that Template:Math for all integers Template:Mvar; thus, this is true also in the particular case when b=adn; but that means that xmodn=(anadn)modn=amodn, which is what we wanted to prove. It remains to be shown that Template:Math. Let Template:Mvar and Template:Mvar be the integers such that Template:Math with Template:Math (see Euclidean division). Then adn=k, thus x=anadn=ank=d+r. Now take Template:Math and add Template:Mvar to both sides, obtaining Template:Math. But we've seen that Template:Math, so we are done.

The modulo with offset Template:Math is implemented in Mathematica as Template:Code .[54]

Implementing other modulo definitions using truncation

Despite the mathematical elegance of Knuth's floored division and Euclidean division, it is generally much more common to find a truncated division-based modulo in programming languages. Leijen provides the following algorithms for calculating the two divisions given a truncated integer division:[55]

/* Euclidean and Floored divmod, in the style of C's ldiv() */
typedef struct {
  /* This structure is part of the C stdlib.h, but is reproduced here for clarity */
  long int quot;
  long int rem;
} ldiv_t;

/* Euclidean division */
inline ldiv_t ldivE(long numer, long denom) {
  /* The C99 and C++11 languages define both of these as truncating. */
  long q = numer / denom;
  long r = numer % denom;
  if (r < 0) {
    if (denom > 0) {
      q = q - 1;
      r = r + denom;
    } else {
      q = q + 1;
      r = r - denom;
    }
  }
  return (ldiv_t){.quot = q, .rem = r};
}

/* Floored division */
inline ldiv_t ldivF(long numer, long denom) {
  long q = numer / denom;
  long r = numer % denom;
  if ((r > 0 && denom < 0) || (r < 0 && denom > 0)) {
    q = q - 1;
    r = r + denom;
  }
  return (ldiv_t){.quot = q, .rem = r};
}

For both cases, the remainder can be calculated independently of the quotient, but not vice versa. The operations are combined here to save screen space, as the logical branches are the same.

See also

Notes

Template:Notelist

References

Template:Reflist

de:Division mit Rest#Modulo