Quantum chemical investigations of reaction pathways connecting a reactant to its product, which provides detailed microscopic information on numerous intermediates and transition states that are challenging to detect experimentally, play crucial roles in predicting chemical reactions and elucidating their mechanisms. With the recent rapid advancement of computational technology, large-scale calculations of various chemical reactions have become feasible. However, the calculations of complex chemical reaction systems remain a challenge. Compared to the computations of low energy and stable structures of reactant and product, those of transition states are difficult because of their inherent instability, making the determination of reaction pathways difficult. In this study, we focus on finding kinetically feasible reaction pathways between a reactant and its product when one structure is known in advance for each reactant and product. This is called the reaction pathway finding problem and is one of the important problem settings in this field. As double-ended computational methods for determining reaction paths between known reactant and product-typically involving a single transition state-the nudged elastic band method and string method 9,10 are well established. For complex systems, where multiple local minima and transition states exist between the reactant and the product, an effective approach is to iteratively apply double-ended methods to pairs of local minima, obtained using structural sampling, which progressively extends the pathway. In contrast, singleended methods compute reaction paths and transition states from one local minimum to another. Among these, the artificial-force-induced reaction (AFIR) method utilizes bond formation or dissociation by applying artificial forces between atom pairs, enabling the effective identification of reaction paths and transition states. A well-established approach involves applying the AFIR method iteratively to local minima discovered along a reaction path, thus expanding the reaction path network. This approach is especially suited to examine chemical reactions and has been applied to various systems. More recently, a method for solving the reaction pathway finding problem has been investigated by exploiting the information of structural similarity between each local minimum and the target product to guide the direction of iterative AFIR applications. In this context, we investigate the following problem: Given two local minima corresponding to a reactant and its product, we initiate the search from the reactant and sample one reaction path at each iteration step, guided by the information of the product. Here, "reaction path" means not the path connecting the reactant and the product, but a reaction path from a local minimum structure to another via one transition state. The goal is to identify a kinetically feasible reaction pathway from the reactant to the product, one that can occur within a given time scale, e.g., 1 day. Because the calculation of all the possible reaction paths is practically impossible, a key challenge is to achieve this goal with a minimal number of computational trials.