The aim of this work is to develop the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) theory based on the point-coupling density functionals and extend it to provide a unified description for all even-even nuclei in the nuclear chart by overcoming all possible challenges. The nuclear superfluidity is considered via Bogoliubov transformation. Densities and potentials are expanded in terms of Legendre polynomials to include the axial deformation degrees of freedom. Sophisticated relativistic Hartree-Bogoliubov equations in coordinate space are solved in the DiracWoods-Saxon basis to consider the continuum effects. Numerical checks are performed from light nuclei to heavy nuclei. The techniques to construct the DRHBc mass table for even-even nuclei are explored. The DRHBc theory is extended to study heavier nuclei beyond magnesium isotopes. Taking Nd isotopes as examples, the experimental binding energies, two-neutron separation energies, quadrupole deformations, and charge radii are reproduced rather well. The deformation and continuum play essential roles in the description of nuclear masses and prediction of drip-line nuclei. By examining the single-particle levels in the canonical basis and their contributions to the total density, the thickness of the neutron skin, the particles number in continuum, and the Coulomb barrier, the exotic structures including the neutron skin and the proton radioactivity are predicted.