We report the theoretical discovery of a systematic scheme to produce topological flat bands (TFBs) with arbitrary Chern numbers. We find that generically a multi-orbital high Chern number TFB model can be constructed by considering multi-layer Chern number C = 1 TFB models with enhanced translational symmetry. A series of models are presented as examples, including a twoband model on a triangular lattice with a Chern number C = 3 and an N -band square lattice model with C = N for an arbitrary integer N . In all these models, the flatness ratio for the TFBs is larger than 30 and increases with increasing Chern number. In the presence of appropriate inter-particle interactions, these models are likely to lead to the formation of novel Abelian and Non-Abelian fractional Chern insulators. As a simple example, we test the C = 2 model with hardcore bosons at 1/3 filling and an intriguing fractional quantum Hall (FQH) state is observed. Introduction -The experimental fractional quantum Hall effect (FQHE) arises from the highly degenerate Landau levels of continuum 2D electron systems, and is described by variational wave functions, first proposed by Laughlin for the primary FQHE states [1] and later generalized by Jain for composite fermion states [2], which are analytic functions of the 2D spatial coordinates. Many important properties of the FQHE, e.g., the hierarchy structures and fractionalized excitations [3,4], can be understood within this framework. It even leads to the predictions of intriguing non-Abelian FQH states at certain filling fractions [5][6][7][8]. Moreover, based on a classification of the pattern of zeros of symmetric (analytic) polynomials, a systematic way to classify FQH states [9,10] has been proposed. Thus, our current theoretical knowledge of FQHE is based on the analytic structure the 2D Landau level (LL) Hilbert space.