In this paper, a novel physics-informed neural networks (PINNs) called deep reduced-order least-square (ROLS) is proposed. The deep ROLS is a combination of reduced-order and least-square methods in the context of PINN methodology. The key idea of the deep ROLS is to convert a higher-order PDE to a system of lower-order PDEs by using primary and secondary variables, then the loss function constituted by integrals of corresponding squared residuals over the problem domain is minimized. Specifically, the squared residuals are established based on L2 norm of respective reduced-order PDE terms, and the Gauss–Legendre quadrature rule is applied to approximate their integrals. A multi-network structure is also designed, in which each sub-network serves as one field variable corresponding to each reduced-order PDE. Moreover, this work proposes a scheme to handle both essential and natural boundary conditions (BCs) directly. The deep ROLS demonstrates its better advantages in comparison to original PINN method in terms of solution accuracy, computational cost and data efficiency for several beam bending problems having continuous and discontinuous solutions.