We calculate the form factors for the baryon number violation process of a heavy-flavor baryon decaying into a pseudoscalar meson and a lepton. In the framework of the Standard Model effective field theory, the leptoquark operators at the bottom quark scale, whose matrix elements define the form factors, are derived by integrating out the high energy physics. Under the QCD factorization approach, the form factors of the baryon number violation process at leading power can be factorized into the convolution of the long-distance hadron wave function as well as the short-distance hard and jet functions representing the hard scale and hard-collinear scale effects, separately. Based on measurements of the baryon number violation process by LHCb, we further impose constraints on the new physics constants of leptoquark operators.