Within the framework of relativistic independent quark model (RIQM) we studied the exclusive two-body nonleptonic decays of \(B_{d(s)} \rightarrow M_{c\bar{c}}(1S,2S) K_s\) , induced by the \(b\rightarrow c\bar{c}s(d)\) transition. The RIQM based on a flavor-independent interaction potential in scalar-vector harmonic form, by which we calculate the weak form factors from the overlapping integrals of meson wave function obtained in this model. Using the factorization approximation, we predict the branching fractions(BF) and the BF ratio, which can be compared with future theoretical predictions. The branching fractions for \(B_{d(s)}\rightarrow \psi (1S,2S) K_s\) decays are found to be within the detection ability of the LHCb and CMS Collaborations, respectively. These results indicate that the present approach works well in the description of exclusive nonleptonic \(B_{d(s)}\) decays within the framework of the RIQ model.