We establish a factorization theorem for the Hardy space in terms of a multilinear fractional integral operator. As an application, a characterization of Lipschitz spaces via the boundedness of commutators of the multilinear fractional integral operator is obtained.