The recently discovered nickelate superconductor La_{3}Ni_{2}O_{7} has a high transition temperature near 80 K under pressure, providing an additional avenue for exploring unconventional superconductivity. Here, with state-of-the-art tensor-network methods, we study a bilayer t-J-J_{⊥} model for La_{3}Ni_{2}O_{7} and find a robust s-wave superconductive (SC) order mediated by interlayer magnetic couplings. Large-scale density matrix renormalization group calculations find algebraic pairing correlations with Luttinger parameter K_{SC}≲1. Infinite projected entangled-pair state method obtains a nonzero SC order directly in the thermodynamic limit, and estimates a strong pairing strength Δ[over ¯]_{z}∼O(0.1). Tangent-space tensor renormalization group simulations elucidate the temperature evolution of SC pairing and further determine a high SC temperature T_{c}^{*}/J∼O(0.1). Because of the intriguing orbital selective behaviors and strong Hund's rule coupling in the compound, t-J-J_{⊥} model has strong interlayer spin exchange (while negligible interlayer hopping), which greatly enhances the SC pairing in the bilayer system. Such a magnetically mediated pairing has also been observed recently in the optical lattice of ultracold atoms. Our accurate and comprehensive tensor-network calculations reveal a robust SC order in the bilayer t-J-J_{⊥} model and shed light on the pairing mechanism of the high-T_{c} nickelate superconductor.