The order of the substitution, equal to the least common multiple of the lengths of its cycles, determines the period length of the sequence of substitution iterations, which is important for cryptographic analysis of systems that use substitution iterations, in particular for data encryption. In this paper, we study the Landau function—the highest order \(\mu \left( n \right)\) of substitution of degree \(n>1\) , where the maximum is reached among all partitions of the number n into cycle lengths. The lengths of substitution cycles are partially ordered with respect to the binary floor division relation, and the set of maximal (in this sense) lengths of substitution cycles uniquely determines its order. It is shown that the value \(\mu \left( n \right)\) , which does not decrease with increasing \(n\) , is achieved on substitutions whose set of maximal cycle lengths is uniquely determined and consists of pairwise coprime primal numbers. For \(n = \sigma^{\left( k \right)} ,k = 1, \ldots ,8\) , where \(\sigma^{\left( k \right)}\) is the sum of the \(k\) smallest prime numbers, the exact values of \(\mu \left( n \right)\) , equal to the product of the \(k\) smallest prime numbers, are obtained. For the remaining \(n\) , formulas for estimating upper and lower bounds of \(\mu \left( n \right)\) are obtained. Numerical values of upper and lower bounds of \(\mu \left( n \right)\) for some \(n \le \sigma^{{\left( {17} \right)}}\) are given.