id Software (the method is older; its origins have been traced to Greg Walsh and Cleve Moler), 1999
Lighting a 3D scene means normalising millions of vectors a second, and each one needs one over the square root of a number. In 1999 that was expensive. This function, found in the Quake III Arena source when id Software released it under the GPL in 2005, does it in a handful of integer operations, to within a fifth of a percent.
The trick is to treat the bits of a floating-point number as an integer. Because of how IEEE 754 lays out exponent and mantissa, shifting that integer right by one and subtracting it from a magic constant gives a surprisingly good first guess at the inverse square root, and one step of Newton's method then polishes it. The constant 0x5f3759df was found by search, and the method has been traced back through Ardent Computer and SGI in the late 1980s. The comments are part of the work. "What the fuck?" records the reaction of whoever read it next, and it has been quoted ever since. Hardware does this directly now, so the function is kept for what it shows: that code can be fast because its author understood how numbers are stored.