Factorial of 401

Counting with 401!

There are 2.5677 × 10^871 different ways to arrange 401 different books on a shelf: 401 choices for the first place, 400 for the second, and so on down to one choice for the last. If you only fill the first three places, as with gold, silver and bronze among 401 runners, there are 401 × 400 × 399 = 63,999,600 possible podiums. If the order does not matter, divide by 3! = 6 to get 10,666,600 ways to choose a group of three. Choosing a pair from 401 people can be done in 80,200 ways.

In probability, if 401 names are drawn out of a hat one at a time, the chance of drawing them in one particular order is 1 in 2.5677 × 10^871. That is why shuffled decks and passwords with many symbols are so hard to guess.

How big is 401!?

401! has 872 digits, so it is larger than the estimated number of atoms in the observable universe (about 10^80). The digits add up to 3339, a multiple of 9, as they are for every factorial from 6! on. Stirling's formula, √(2πn) × (n/e)^n, gives a close estimate without multiplying everything out, and gets more accurate as n grows.

Prime factors of 401!

Every prime up to 401 divides 401!. Legendre's formula counts how often: add ⌊n/p⌋ + ⌊n/p²⌋ + ⌊n/p³⌋ + … for each prime p. For 401! that gives 2^397 × 3^196 × 5^99 × 7^66 × 11^39 × 13^32 × 17^24 × 19^22 × 23^17 × 29^13 × … (79 primes in all). The 99 trailing zeros come from pairing the 99 factors of 5 with factors of 2.

Shuffles where nothing stays put

A related count is the number of derangements of 401 items: orders in which no item ends up in its original place, as in a Secret Santa where nobody draws their own name. The count is the whole number nearest to 401! ÷ e, so about 36.8% of all 2.5677 × 10^871 orders are derangements. From about four items on, that share stays close to 1 ÷ e, roughly 37%, however many items there are.

Computing 401! yourself

In Excel or Google Sheets, type =FACT(401). Spreadsheets keep only 15 significant digits, so from 18! on the result is rounded, and FACT stops at 170!, the largest factorial that fits in a standard floating-point number. In Python, math.factorial(401) returns the exact integer with all 872 digits, because Python integers have no size limit. The gamma function extends factorials to fractions: Γ(n + 1) = n!, so (1/2)! = √π ÷ 2 ≈ 0.8862.

Why is 0! equal to 1?

There is exactly one way to arrange nothing, and the rule n! = n × (n − 1)! only works for n = 1 if 0! = 1. Defining it this way keeps formulas for combinations and series correct at the edges.