[en] Exact algebraic solutions of a hard-core Bose-Hubbard model for the description of nearest orbit pairing interaction in well-deformed even-even nuclei were derived based on a simple algebraic routine. Excitation energies and the corresponding wavefunctions of the model were obtained. Further extensions to arbitrary spin cases and the corresponding Fermi-Hubbard model were also presented