To hear the shape of a quantum drum is in principle not a mission impossible, where the shape means the conformal field theory (CFT) description of certain quantum phases and phase transitions in a (2 + 1)d manifold and the drum beat means some characteristic (such as spectral) measurements of the quantum system whose structure reveals the CFT information. In realistic quantum many-body lattice models, however, measurements such as the finite size scaling of the entanglement entropy with reliable data to extract the CFT characteristics, are rare due to the numerical difficulties in the precise measurement of entanglement with replicas of the partition functions. To overcome this problem, we design a generic computation protocol to obtain the Rényi entanglement entropy of the aforementioned systems with efficiency and precision. Our method, dubbed "Qiu Ku" algorithm, is based on the massive parallelization of the nonequilibrium increment processes in quantum Monte Carlo simulations. To demonstrate its power, we show the results on a few important yet difficult (2 + 1)d quantum lattice models, ranging from Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, quantum critical point with O(3) CFT to the toric code Z2 topological ordered state and the Kagome Z2 quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method reveals either the precise CFT data from the logarithmic correction or quantum dimension of the topological order, with unprecedentedly large system sizes, controlled errorbars and minimal computational costs. Our "Qiu Ku" algorithm therefore helps one to clearly hear the shapes of various difficult yet important quantum drums.