Recently, fractional Chern insulators (FCIs), also called fractional quantum anomalous Hall (FQAH) states, have been theoretically established in lattice systems with topological flat bands. These systems exhibit similar fractionalization phenomena to the conventional fractional quantum Hall (FQH) systems. Using the mapping relationship between the FQH states and the FCI/FQAH states, we construct the many-body wave functions of the fermionic FCI/FQAH states in disk geometry with the aid of the generalized Pauli principle (GPP) and Jack polynomials. Compared with the ground state by the exact diagonalization method, the wave-function overlap is higher than 0.97, even when the Hilbert space dimension is as large as 3×10 6 . We also use the GPP and the Jack polynomials to construct edge excitations for the fermionic FCI/FQAH states. The quasi-degeneracy sequences of fermionic FCI/FQAH systems reproduce the prediction of the chiral Luttinger liquid theory, complementing the exact diagonalization results with larger lattice sizes and more particles. OPEN ACCESS RECEIVED exclusion principle corresponds to the Fermi statistics, the GPP corresponds to the fractional statistics of anyons.Recently, finding the suitable trial wave functions of FCI/FQAH states has become an outstanding problem. It is conjectured that there is a one-to-one mapping between IQH and CI/QAH states [10,[34][35][36][37]. In other words, the single-particle orbitals of an LL can be mapped to the ones in a TFB. Based on one-dimensional (1D) maximally localized Wannier functions, Qi et al have used a mapping between the FQH and the FCI/FQAH states to construct generic wave functions of FCI/FQAH states in cylinder geometry with these 1D Wannier functions [10,[34][35][36][37]. In terms of this mapping relationship, the Haldane pseudo-potential [26] for these FCI/ FQAH states can be constructed [38,39] through a suitable gauge choice for the Wannier functions. An improved prescription has been adopted to construct variational wave functions of FCI/FQAH states in torus geometry by utilizing gauge-fixed (non-maximally) localized Wannier states [40][41][42]. From another aspect, conventional FQH states can also be obtained for 2D lattices analytically by using conformal field theory [43][44][45]. In contrast to the above analytical or semi-analytical approaches, we here pursue a very direct yet effective purely numerical prescription to construct FCI/FQAH states in disk geometry, without any variational parameter or adjustable gauge freedom, but just utilizing the powerful GPP and the Jacks structure of FQH states, and the information from exact numerical single-particle orbitals of TFB models.In this paper, we exploit the single-particle wave functions of CI/QAH states (in the Kagomé model [5, 46] and the Haldane [3, 9] model) with TFB parameters in disk geometry, and explore the polynomial structure of the continuum Laughlin wave functions [25] to establish the many-body wave functions of the FCI/FQAH states. We firstly construct finite-size lattices in...