In the context of tensor network states, we for first time reformulate corner transfer matrix renormalization group (CTMRG) method into a variational bilevel optimization algorithm. The solution problem corresponds to fixed-point environment pursued in conventional CTMRG method, from which partition function classical statistical model, represented by an infinite network, can be efficiently eva...