A study of B-c(+) -> J/psi D-s(+) and B-c(+)-> J/psi D-s*(+) decays using 139 fb(-1) of in- tegrated luminosity collected with the ATLAS detector from root s = 13 TeV pp collisions at the LHC is presented. The ratios of the branching fractions of the two decays to the branching fraction of the B-c(+) -> J/psi pi(+) decay are measured: B(B-c(+) -> J/psi D-s(+))/B(B-c(+) -> J/psi pi(+)) = 2.76 +/- 0.47 and B(B-c(+)-> J/psi D-s*(+))/B(B-c(+) -> J/psi pi(+)) = 5.33 +/- 0.96. The ratio of the branching fractions of the two decays is found to be B(B-c(+)-> J/psi D-s*(+))/B(B-c(+) -> J/psi D-s(+)) = 1.93 +/- 0.26. For the B-c(+)-> J/psi D-s*(+) decay, the transverse polarization fraction, Gamma(+/-+/-)/Gamma, is measured to be 0.70 +/- 0.11. The reported uncertainties include both the statistical and systematic components added in quadrature. The precision of the measurements exceeds that in all previous studies of these decays. These results supersede those obtained in the earlier ATLAS study of the same decays with root s = 7 and 8 TeV pp collision data. A comparison with available theoretical predictions for the measured quantities is presented.