We compute the next-to-leading-order (NLO) perturbative corrections to the exclusive decay Υ → J/ψ +η c (χ cJ ) (J = 0, 1, 2) within the nonrelativistic QCD(NRQCD) factorization framework. The amplitudes of four decay patterns are separately evaluated up to relative order α s . We analyze the impact of the NLO corrections to the (un)polarized decay widths as well as branching fractions, which largely reduce the renormalization scale dependence for J/ψ + χ cJ production, while slightly improve the renormalization scale dependence for J/ψ + η c . With high numerical accuracy, our theoretical prediction for the branching fraction of Υ → J/ψ + χ c1 yields to 3.65 +4.99+0.12 −2.02−1.20 × 10 −6 where the uncertainties are from the ambiguities of charm/bottom masses and renormalization scale, which agrees perfectly with Belle measurement (3.90 ± 1.21 ± 0.23) × 10 −6 . In addition, the theoretical predictions for other double charmonia production from Υ decay are consistent with the upper limit set by Belle measurement with uncertainties, except for the process Υ(2S) → J/ψ + χ c1 , where there exists some tension between theory and measurement. Since the predicted branching fractions are around 10 −6 , it is prospective to measure these processes in high-luminosity e + e − colliders as the super B factory.