In this work, we establish a relation between entanglement entropy and fractal dimension D of generic many-body wave functions, by generalizing the result of Page [Phys. Rev. Lett. 71, 1291 (1993)PRLTAO0031-900710.1103/PhysRevLett.71.1291] to the case of sparse random pure states (SRPS). These SRPS living in a Hilbert space of size N are defined as normalized vectors with only N^{D} (0≤D≤1) random nonzero elements. For D=1, these states used by Page represent ergodic states at an infinite temperature. However, for 0<D<1, the SRPS are nonergodic and fractal, as they are confined in a vanishing ratio N^{D}/N of the full Hilbert space. Both analytically and numerically, we show that the mean entanglement entropy S_{1}(A) of a subsystem A, with Hilbert space dimension N_{A}, scales as S_{1}[over ¯](A)∼DlnN for small fractal dimensions D, N^{D}<N_{A}. Remarkably, S_{1}[over ¯](A) saturates at its thermal (Page) value at an infinite temperature, S_{1}[over ¯](A)∼lnN_{A} at larger D. Consequently, we provide an example when the entanglement entropy takes an ergodic value even though the wave function is highly nonergodic. Finally, we generalize our results to Renyi entropies S_{q}(A) with q>1 and to genuine multifractal states and also show that their fluctuations have ergodic behavior in a narrower vicinity of the ergodic state D=1.