The exploration of neural network quantum states has become widespread in studies of complicated quantum many-body systems. However, achieving high precision remains challenging due to the exponential growth of Hilbert space size and the intricate sign structures. Utilizing symmetries of the physical system, we propose a method to evaluate and sample the variational Ansatz within a symmetric subspace. This approach isolates different symmetry sectors, reducing the relevant Hilbert space size by a factor approximately proportional to the size of the symmetry group. It is inspired by exact diagonalization techniques and the work of Choo et al., [Phys. Rev. Lett. 121, 167204 (2018)10.1103/PhysRevLett.121.167204]. We validate our method using the frustrated spin-12J1-J2 antiferromagnetic Heisenberg chain, and we compare its performance to the case without symmetrization. The results indicate that our symmetric subspace approach achieves a substantial improvement over the full Hilbert space on optimizing the Ansatz, reducing the energy error by orders of magnitude. We also compare the results on degenerate eigenstates with different quantum numbers, highlighting the advantage of operating within a smaller Hilbert subspace.