Traditional condensed matter physics is based on two theories: symmetry breaking theory for phases and phase transitions, and Fermi liquid theory for metals. Mean-field theory is a powerful method to describe symmetry breaking phases and phase transitions by assuming the ground state wavefunctions for many-body systems can be approximately described by direct product states. The Fermi liquid theory is another powerful method to study electron systems by assuming that the ground state wavefunctions for the electrons can be approximately described by Slater determinants. From the encoding point of view, both methods only use a polynomial amount of information to approximately encode many-body ground state wavefunctions which contain an exponentially large amount of information. Moreover, another nice property of both approaches is that all the physical quantities (energy, correlation functions, etc.) can be efficiently calculated (polynomially hard). In this talk, I'll introduce a new class of states: (Grassmann-number) tensor product states. These states only need polynomial amount of information to approximately encode many-body ground states. Many classes of states, such as Slater determinant states, projective states, string-net states and their generalizations, etc., are subclasses of (Grassmann-number) tensor prodcut states. Thus, (Grassmann-number) tensor prodcut states provide us a unified variational approach to describe both symmetry breaking and topologically ordered states. However, calculating the physical quantities for these state can be exponentially hard in general. To solve this difficulty, we develop the (Grassmann) Tensor-Entanglement Renormalization Group (TERG) method to efficiently calculate the physical quantities. We demonstrate our algorithm by studying several interesting boson/fermion models. Finally, we will show its non-trivial application to t-J model on a honeycomb lattice.