Floor and ceiling functions

From testwiki
Jump to navigation Jump to search

Template:Short description Template:Use dmy dates Template:Multiple image

In mathematics, the floor function is the function that takes a real number Template:Mvar as input and returns the greatest integer less than or equal to Template:Mvar, written Template:Math or Template:Math. Similarly, the ceiling function returns the least integer greater than or equal to Template:Math, written Template:Math or Template:Math.[1]

For example, for floor: Template:Math, Template:Math, and for ceiling: Template:Math, and Template:Math.

The floor of Template:Mvar is also called the integral part, integer part, greatest integer, or entier of Template:Mvar, and was historically denoted Template:Math (among other notations).[2] However, the term "integer part" is ambiguous, as it can also mean truncation towards zero, which differs from the floor function for negative numbers.

For an integer Template:Mvar, Template:Math.

Although Template:Math and Template:Math are equal for non-integer values of Template:Mvar, and thus produce graphs that appear exactly alike, they differ when Template:Mvar is an integer. For example, when Template:Math, Template:Math. However, if Template:Math, then Template:Math but Template:Math.

Examples
x Floor Template:Math Ceiling Template:Math Fractional part Template:Math
2 2 2 0
2.0001 2 3 0.0001
[[e (mathematical constant)|Template:Mvar]] 2 3 0.7182...
2.9 2 3 0.9
2.999 2 3 0.999
−[[Pi|Template:Mvar]] −4 −3 0.8584...
−2 −2 −2 0

Notation

The integral part or integer part of a number (Template:Lang in the original) was first defined in 1798 by Adrien-Marie Legendre in his proof of the Legendre's formula.

Carl Friedrich Gauss introduced the square bracket notation Template:Math in his third proof of quadratic reciprocity (1808).[3] This remained the standard[4] in mathematics until Kenneth E. Iverson introduced, in his 1962 book A Programming Language, the names "floor" and "ceiling" and the corresponding notations Template:Math and Template:Math.[5][6] (Iverson used square brackets for a different purpose, the Iverson bracket notation.) Both notations are now used in mathematics, although Iverson's notation will be followed in this article.

In some sources, boldface or double brackets Template:Math are used for floor, and reversed brackets Template:Math or Template:Math for ceiling.[7][8]

The fractional part is the sawtooth function, denoted by Template:Math for real Template:Mvar and defined by the formula

Template:Math[9]

For all x,

Template:Math.

These characters are provided in Unicode and HTML entities:

In the LaTeX typesetting system, these symbols can be specified with the Template:Mono and Template:Mono commands in math mode. LaTeX has supported UTF-8 since 2018, so the Unicode characters can now be used directly.[10] Larger versions areTemplate:Mono and Template:Mono.

Definition and properties

Given real numbers Template:Mvar and Template:Mvar, integers Template:Mvar and Template:Mvar and the set of integers ℤ, floor and ceiling may be defined by the equations ⌊x⌋=max⁡{m∈ℤ∣m≤x}, ⌈x⌉=min⁡{n∈ℤ∣n≥x}.

Since there is exactly one integer in a half-open interval of length one, for any real number Template:Mvar, there are unique integers Template:Mvar and Template:Mvar satisfying the equation x−1<m≤x≤n<x+1, where ⌊x⌋=m and ⌈x⌉=n may also be taken as the definition of floor and ceiling.

Equivalences

These formulas can be used to simplify expressions involving floors and ceilings.[11]

⌊x⌋=m   if and only if m≤x<m+1,⌈x⌉=n if and only if   n−1<x≤n,⌊x⌋=m if and only if x−1<m≤x,⌈x⌉=n if and only if x≤n<x+1.

In the language of order theory, the floor function is a residuated mapping, that is, part of a Galois connection: it is the upper adjoint of the function that embeds the integers into the reals.

x<n if and only if ⌊x⌋<n,n<x if and only if n<⌈x⌉,x≤n if and only if ⌈x⌉≤n,n≤x if and only if n≤⌊x⌋.

These formulas show how adding an integer Template:Mvar to the arguments affects the functions: ⌊x+n⌋=⌊x⌋+n,⌈x+n⌉=⌈x⌉+n,{x+n}={x}.

The above are never true if Template:Mvar is not an integer; however, for every Template:Mvar and Template:Mvar, the following inequalities hold: ⌊x⌋+⌊y⌋≤⌊x+y⌋≤⌊x⌋+⌊y⌋+1,⌈x⌉+⌈y⌉−1≤⌈x+y⌉≤⌈x⌉+⌈y⌉.

Monotonicity

Both floor and ceiling functions are monotonically non-decreasing functions: x1≤x2⇒⌊x1⌋≤⌊x2⌋,x1≤x2⇒⌈x1⌉≤⌈x2⌉.

Relations among the functions

It is clear from the definitions that ⌊x⌋≤⌈x⌉, with equality if and only if x is an integer, i.e. ⌈x⌉−⌊x⌋={0 if x∈ℤ1 if x∉ℤ

In fact, for integers Template:Mvar, both floor and ceiling functions are the identity: ⌊n⌋=⌈n⌉=n.

Negating the argument switches floor and ceiling and changes the sign: ⌊x⌋+⌈−x⌉=0−⌊x⌋=⌈−x⌉−⌈x⌉=⌊−x⌋ and: ⌊x⌋+⌊−x⌋={0if x∈ℤ−1if x∉ℤ, ⌈x⌉+⌈−x⌉={0if x∈ℤ1if x∉ℤ.

Negating the argument complements the fractional part: {x}+{−x}={0if x∈ℤ1if x∉ℤ.

The floor, ceiling, and fractional part functions are idempotent: ⌊⌊x⌋⌋=⌊x⌋,⌈⌈x⌉⌉=⌈x⌉,{{x}}={x}.

The result of nested floor or ceiling functions is the innermost function: ⌊⌈x⌉⌋=⌈x⌉,⌈⌊x⌋⌉=⌊x⌋ due to the identity property for integers.

Quotients

If Template:Mvar and Template:Mvar are integers and Template:Math, 0≤{mn}≤1−1|n|.

If Template:Mvar is positive[12] ⌊x+mn⌋=⌊⌊x⌋+mn⌋, ⌈x+mn⌉=⌈⌈x⌉+mn⌉.

If Template:Mvar is positive[13] n=⌈n1m⌉+⌈n−1m⌉+…+⌈n−m+1m⌉, n=⌊n1m⌋+⌊n+1m⌋+…+⌊n+m−1m⌋.

For Template:Math these imply n=⌊n12⌋+⌈n12⌉.

More generally,[14] for positive Template:Mvar (See Hermite's identity) ⌈mx⌉=⌈x⌉+⌈x−1m⌉+…+⌈x−m−1m⌉, ⌊mx⌋=⌊x⌋+⌊x+1m⌋+…+⌊x+m−1m⌋.

The following can be used to convert floors to ceilings and vice versa (with Template:Mvar being positive)[15] ⌈n1m⌉=⌊n+m−1m⌋=⌊n−1m⌋+1, ⌊n1m⌋=⌈n−m+1m⌉=⌈n+1m⌉−1,

For all Template:Mvar and Template:Mvar strictly positive integers:[16] ∑k=1n−1⌊kmn⌋=(m−1)(n−1)+gcd⁡(m,n)−12, which, for positive and coprime Template:Mvar and Template:Mvar, reduces to ∑k=1n−1⌊kmn⌋=12(m−1)(n−1), and similarly for the ceiling and fractional part functions (still for positive and coprime Template:Mvar and Template:Mvar), ∑k=1n−1⌈kmn⌉=12(m+1)(n−1), ∑k=1n−1{kmn}=12(n−1).


Since the right-hand side of the general case is symmetrical in Template:Mvar and Template:Mvar, this implies that ⌊m1n⌋+⌊2mn⌋+…+⌊(n−1)mn⌋=⌊n1m⌋+⌊2nm⌋+…+⌊(m−1)nm⌋.

More generally, if Template:Mvar and Template:Mvar are positive, ⌊x1n⌋+⌊m+xn⌋+⌊2m+xn⌋+…+⌊(n−1)m+xn⌋=⌊x1m⌋+⌊n+xm⌋+⌊2n+xm⌋+⋯+⌊(m−1)n+xm⌋.

This is sometimes called a reciprocity law.[17]

Division by positive integers gives rise to an interesting and sometimes useful property. Assuming m,n>0, m≤⌊xn⌋⟺n≤⌊xm⌋⟺n≤⌊x⌋m.

Similarly, m≥⌈xn⌉⟺n≥⌈xm⌉⟺n≥⌈x⌉m.

Indeed, m≤⌊xn⌋⟹m≤xn⟹n≤xm⟹n≤⌊xm⌋⟹…⟹m≤⌊xn⌋, keeping in mind that ⌊xn⌋=⌊⌊x⌋n⌋. The second equivalence involving the ceiling function can be proved similarly.

For Template:Mvar being a positive integer with Template:Mvar greater than Template:Mvar. Then [18] {xd}={x}+rd and ⌊xd⌋=⌊x⌋−rd, where 0≤r≤d−1 is the remainder of dividing ⌊x⌋ by Template:Mvar

Nested divisions

For a positive integer Template:Mvar, and arbitrary real numbers Template:Mvar and Template:Mvar:[19] ⌊⌊xm⌋n⌋=⌊xmn⌋⌈⌈xm⌉n⌉=⌈xmn⌉.

Continuity and series expansions

None of the functions discussed in this article are continuous, but all are piecewise linear: the functions ⌊x⌋, ⌈x⌉, and {x} have discontinuities at the integers.

⌊x⌋ is upper semi-continuous and ⌈x⌉ and {x} are lower semi-continuous.

Since none of the functions discussed in this article are continuous, none of them have a power series expansion. Since floor and ceiling are not periodic, they do not have uniformly convergent Fourier series expansions. The fractional part function has Fourier series expansion[20] {x}=12−1π∑k=1∞sin⁡(2πkx)k for Template:Mvar not an integer.

At points of discontinuity, a Fourier series converges to a value that is the average of its limits on the left and the right, unlike the floor, ceiling and fractional part functions: for Template:Mvar fixed and Template:Mvar a multiple of Template:Mvar the Fourier series given converges to Template:Math, rather than to Template:Math. At points of continuity the series converges to the true value.

Using the formula ⌊x⌋=x−{x} gives ⌊x⌋=x−12+1π∑k=1∞sin⁡(2πkx)k for Template:Mvar not an integer.

Applications

Mod operator

For an integer Template:Mvar and a positive integer Template:Mvar, the modulo operation, denoted by Template:Math, gives the value of the remainder when Template:Mvar is divided by Template:Mvar. This definition can be extended to real Template:Mvar and Template:Mvar, Template:Math, by the formula xmody=x−y⌊xy⌋.

Then it follows from the definition of floor function that this extended operation satisfies many natural properties. Notably, Template:Math is always between 0 and Template:Mvar, i.e.,

if Template:Mvar is positive, 0≤xmody<y, and if Template:Mvar is negative, 0≥xmody>y.

Quadratic reciprocity

Gauss's third proof of quadratic reciprocity, as modified by Eisenstein, has two basic steps.[21][22]

Let Template:Mvar and Template:Mvar be distinct positive odd prime numbers, and let m=12(p−1), n=12(q−1).

First, Gauss's lemma is used to show that the Legendre symbols are given by (qp)=(−1)⌊qp⌋+⌊2qp⌋+…+⌊mqp⌋,(pq)=(−1)⌊pq⌋+⌊2pq⌋+…+⌊npq⌋.

The second step is to use a geometric argument to show that ⌊qp⌋+⌊2qp⌋+…+⌊mqp⌋+⌊pq⌋+⌊2pq⌋+…+⌊npq⌋=mn.

Combining these formulas gives quadratic reciprocity in the form (pq)(qp)=(−1)mn=(−1)p−12q−12.

There are formulas that use floor to express the quadratic character of small numbers mod odd primes Template:Mvar:[23] (2p)=(−1)⌊p+14⌋,(3p)=(−1)⌊p+16⌋.

Rounding

For an arbitrary real number x, rounding x to the nearest integer with tie breaking towards positive infinity is given by rpi(x)=⌊x+12⌋=⌈12⌊2x⌋⌉; rounding towards negative infinity is given as rni(x)=⌈x−12⌉=⌊12⌈2x⌉⌋.

If tie-breaking is away from 0, then the rounding function is ri(x)=sgn⁡(x)⌊|x|+12⌋ (where sgn is the sign function), and rounding towards even can be expressed with the more cumbersome ⌊x⌉=⌊x+12⌋+⌈14(2x−1)⌉−⌊14(2x−1)⌋−1, which is the above expression for rounding towards positive infinity rpi(x) minus an integrality indicator for 14(2x−1).

Rounding a real number x to the nearest integer value forms a very basic type of quantizer – a uniform one. A typical (mid-tread) uniform quantizer with a quantization step size equal to some value Δ can be expressed as Q(x)=Δ⋅⌊xΔ+12⌋,

Number of digits

The number of digits in base Template:Mvar of a positive integer Template:Mvar is ⌊logbk⌋+1=⌈logb(k+1)⌉.

Number of strings without repeated characters

The number of possible strings of arbitrary length that doesn't use any character twice is given by[24]Template:Better source needed (n)0+⋯+(n)n=⌊en!⌋ where:

For Template:Math, this comes out to 1096259850353149530222034277.

Factors of factorials

Let Template:Mvar be a positive integer and Template:Mvar a positive prime number. The exponent of the highest power of Template:Mvar that divides Template:Math is given by a version of Legendre's formula[25] ⌊np⌋+⌊np2⌋+⌊np3⌋+…=n−∑kakp−1 where n=∑kakpk is the way of writing Template:Mvar in base Template:Mvar. This is a finite sum, since the floors are zero when Template:Math.

Beatty sequence

The Beatty sequence shows how every positive irrational number gives rise to a partition of the natural numbers into two sequences via the floor function.[26]

Euler's constant (γ)

There are formulas for Euler's constant Template:Math that involve the floor and ceiling, e.g.[27] γ=∫1∞(1⌊x⌋−1x)dx, γ=limn→∞1n∑k=1n(⌈nk⌉−nk), and γ=∑k=2∞(−1)k⌊log2k⌋k=12−13+2(14−15+16−17)+3(18−⋯−115)+⋯

Riemann zeta function (ζ)

The fractional part function also shows up in integral representations of the Riemann zeta function. It is straightforward to prove (using integration by parts)[28] that if φ(x) is any function with a continuous derivative in the closed interval Template:Math, ∑a<n≤bφ(n)=∫abφ(x)dx+∫ab({x}−12)φ′(x)dx+({a}−12)φ(a)−({b}−12)φ(b).

Letting φ(n)=n−s for real part of Template:Mvar greater than 1 and letting Template:Mvar and Template:Mvar be integers, and letting Template:Mvar approach infinity gives ζ(s)=s∫1∞12−{x}xs+1dx+1s−1+12.

This formula is valid for all Template:Mvar with real part greater than −1, (except Template:Math, where there is a pole) and combined with the Fourier expansion for Template:Math can be used to extend the zeta function to the entire complex plane and to prove its functional equation.[29]

For Template:Math in the critical strip Template:Math, ζ(s)=s∫−∞∞e−σω(⌊eω⌋−eω)e−itωdω.

In 1947 van der Pol used this representation to construct an analogue computer for finding roots of the zeta function.[30]

Formulas for prime numbers

The floor function appears in several formulas characterizing prime numbers. For example, since ⌊nm⌋−⌊n−1m⌋={1if m divides n0otherwise, it follows that a positive integer Template:Mvar is a prime if and only if[31] ∑m=1∞(⌊nm⌋−⌊n−1m⌋)=2.

One may also give formulas for producing the prime numbers. For example, let Template:Math be the Template:Mvar-th prime, and for any integer Template:Math, define the real number Template:Mvar by the sum α=∑m=1∞pmr−m2.

Then[32] pn=⌊rn2α⌋−r2n−1⌊r(n−1)2α⌋.

A similar result is that there is a number Template:Math (Mills' constant) with the property that ⌊θ3⌋,⌊θ9⌋,⌊θ27⌋,… are all prime.[33]

There is also a number Template:Math with the property that ⌊2ω⌋,⌊22ω⌋,⌊222ω⌋,… are all prime.[33]

Let Template:Math be [[Prime-counting function|the number of primes less than or equal to Template:Mvar]]. It is a straightforward deduction from Wilson's theorem that[34] π(n)=∑j=2n⌊(j−1)!+1j−⌊(j−1)!j⌋⌋.

Also, if Template:Math,[35] π(n)=∑j=2n⌊1∑k=2j⌊⌊jk⌋kj⌋⌋.

None of the formulas in this section are of any practical use.[36][37]

Solved problems

Ramanujan submitted these problems to the Journal of the Indian Mathematical Society.[38]

If Template:Mvar is a positive integer, prove that

  1. ⌊n3⌋+⌊n+26⌋+⌊n+46⌋=⌊n2⌋+⌊n+36⌋,
  2. ⌊12+n+12⌋=⌊12+n+14⌋,
  3. ⌊n+n+1⌋=⌊4n+2⌋.

Some generalizations to the above floor function identities have been proven.[39]

Unsolved problem

The study of Waring's problem has led to an unsolved problem:

Are there any positive integers Template:Math such that[40] 3k−2k⌊(32)k⌋>2k−⌊(32)k⌋−2 ?

Mahler has proved there can only be a finite number of such Template:Mvar; none are known.[41]

Computer implementations

int function from floating-point conversion in C

In most programming languages, the simplest method to convert a floating point number to an integer does not do floor or ceiling, but truncation. The reason for this is historical, as the first machines used ones' complement and truncation was simpler to implement (floor is simpler in two's complement). FORTRAN was defined to require this behavior and thus almost all processors implement conversion this way. Some consider this to be an unfortunate historical design decision that has led to bugs handling negative offsets and graphics on the negative side of the origin.Template:Cn

An arithmetic right-shift of a signed integer x by n is the same as ⌊x2n⌋. Division by a power of 2 is often written as a right-shift, not for optimization as might be assumed, but because the floor of negative results is required. Assuming such shifts are "premature optimization" and replacing them with division can break software.[42]

Many programming languages (including C, C++,[43][44] C#,[45][46] Java,[47][48] Julia,[49] PHP,[50][51] R,[52] and Python[53]) provide standard functions for floor and ceiling, usually called floor and ceil, or less commonly ceiling.[54] The language APL uses ⌊x for floor. The J Programming Language, a follow-on to APL that is designed to use standard keyboard symbols, uses <. for floor and >. for ceiling.[55] ALGOL usesentier for floor.

In Microsoft Excel the function INT rounds down rather than toward zero,[56] while FLOOR rounds toward zero, the opposite of what "int" and "floor" do in other languages. Since 2010 FLOOR has been changed to error if the number is negative.[57] The OpenDocument file format, as used by OpenOffice.org, Libreoffice and others, INT[58] and FLOOR both do floor, and FLOOR has a third argument to reproduce Excel's earlier behavior.[59]

See also

Citations

Template:Reflist

References

Template:Commons category

  1. ↑ Graham, Knuth, & Patashnik, Ch. 3.1
  2. ↑ 1) Luke Heaton, A Brief History of Mathematical Thought, 2015, Template:Isbn (n.p.)
    2) Albert A. Blank et al., Calculus: Differential Calculus, 1968, p. 259
    3) John W. Warris, Horst Stocker, Handbook of mathematics and computational science, 1998, Template:Isbn, p. 151
  3. ↑ Lemmermeyer, pp. 10, 23.
  4. ↑ e.g. Cassels, Hardy & Wright, and Ribenboim use Gauss's notation. Graham, Knuth & Patashnik, and Crandall & Pomerance use Iverson's.
  5. ↑ Iverson, p. 12.
  6. ↑ Higham, p. 25.
  7. ↑ Mathwords: Floor Function.
  8. ↑ Mathwords: Ceiling Function
  9. ↑ Graham, Knuth, & Patashnik, p. 70.
  10. ↑ Template:Cite web
  11. ↑ Graham, Knuth, & Patashink, Ch. 3
  12. ↑ Graham, Knuth, & Patashnik, p. 73
  13. ↑ Graham, Knuth, & Patashnik, p. 85
  14. ↑ Graham, Knuth, & Patashnik, p. 85 and Ex. 3.15
  15. ↑ Graham, Knuth, & Patashnik, Ex. 3.12
  16. ↑ Graham, Knuth, & Patashnik, p. 94.
  17. ↑ Graham, Knuth, & Patashnik, p. 94
  18. ↑ Problems and Solutions, The College Mathematics Journal, 56:4
  19. ↑ Graham, Knuth, & Patashnik, p. 71, apply theorem 3.10 with Template:Sfrac as input and the division by n as function
  20. ↑ Titchmarsh, p. 15, Eq. 2.1.7
  21. ↑ Lemmermeyer, § 1.4, Ex. 1.32–1.33
  22. ↑ Hardy & Wright, §§ 6.11–6.13
  23. ↑ Lemmermeyer, p. 25
  24. ↑ Template:OEIS el (See Formulas.)
  25. ↑ Hardy & Wright, Th. 416
  26. ↑ Graham, Knuth, & Patashnik, pp. 77–78
  27. ↑ These formulas are from the Wikipedia article Euler's constant, which has many more.
  28. ↑ Titchmarsh, p. 13
  29. ↑ Titchmarsh, pp.14–15
  30. ↑ Crandall & Pomerance, p. 391
  31. ↑ Crandall & Pomerance, Ex. 1.3, p. 46. The infinite upper limit of the sum can be replaced with n. An equivalent condition is n > 1 is prime if and only if ∑m=1⌊n⌋(⌊nm⌋−⌊n−1m⌋)=1.
  32. ↑ Hardy & Wright, § 22.3
  33. ↑ 33.0 33.1 Ribenboim, p. 186
  34. ↑ Ribenboim, p. 181
  35. ↑ Crandall & Pomerance, Ex. 1.4, p. 46
  36. ↑ Ribenboim, p. 180 says that "Despite the nil practical value of the formulas ... [they] may have some relevance to logicians who wish to understand clearly how various parts of arithmetic may be deduced from different axiomatzations ... "
  37. ↑ Hardy & Wright, pp. 344—345 "Any one of these formulas (or any similar one) would attain a different status if the exact value of the number α ... could be expressed independently of the primes. There seems no likelihood of this, but it cannot be ruled out as entirely impossible."
  38. ↑ Ramanujan, Question 723, Papers p. 332
  39. ↑ Template:Cite journal
  40. ↑ Hardy & Wright, p. 337
  41. ↑ Template:Cite journal
  42. ↑ Sutter, Herb; Alexandrescu, Andrei. C++ Coding Standards: 101 Rules, Guidelines, and Best Practices. Addison-Wesley 2004. ISBN 978-0-321-11358-0. ch. 8, 9
  43. ↑ Template:Cite web
  44. ↑ Template:Cite web
  45. ↑ Template:Cite web
  46. ↑ Template:Cite web
  47. ↑ Template:Cite web
  48. ↑ Template:Cite web
  49. ↑ Template:Cite web
  50. ↑ Template:Cite web
  51. ↑ Template:Cite web
  52. ↑ Template:Cite web
  53. ↑ Template:Cite web
  54. ↑ Sullivan, p. 86.
  55. ↑ Template:Cite web
  56. ↑ Template:Cite web
  57. ↑ Template:Cite web
  58. ↑ Template:Cite web
  59. ↑ Template:Cite web