离散随机过程在自然界中普遍存在，在物理、生物化学、流行病学、社会学和金融等领域有广泛应用。尽管解析解往往无法获得，现有的仿真框架能够生成与随机现象背后的动态规律兼容的随机轨迹。然而，大多数仿真算法假设系统动力学是无记忆的（马尔可夫假设），在这种假设下，未来的事件仅依赖于系统的当前状态，这使得通过Gillespie算法进行高效且精确的仿真成为可能。但是，许多真实世界系统本质上是非马尔可夫的，并表现出记忆效应。此类系统难以进行解析研究，而当前的数值方法往往计算成本高昂或受限于与实证数据相冲突的强大简化假设。为了解决这些局限性，我们引入了非马尔可夫反应拒绝型Gillespie算法（REGIR），这是一种通用且可扩展的框架，用于模拟具有任意事件间时间分布的非马尔可夫随机系统。REGIR在提供用户定义精度的同时保持了与经典Gillespie算法相同的渐进计算复杂度。我们推导了REGIR近似精度的下界，并展示了其在三类代表性非马尔可夫系统中的能力：（1）具有延迟的反应通道，（2）由个体反应物属性驱动的随机过程，以及（3）受节点活动支配的时间网络。在所有情况下，REGIR准确捕捉了依赖记忆的动力学，并在灵活性和效率方面优于现有方法。

Discrete stochastic processes are prevalent in natural systems, with applications in physics, biochemistry, epidemiology, sociology, and finance. While analytic solutions often cannot be derived, existing simulation frameworks can generate stochastic trajectories compatible with the dynamical laws underlying the random phenomena. Still, most simulation algorithms assume the system dynamics are memoryless (Markovian assumption), under which assumption, future occurrences only depend on the system's present state. This enables efficient and exact simulation via the Gillespie algorithm. However, many real-world systems are inherently non-Markovian and exhibit memory effects. Such systems are difficult to study analytically, and current numerical methods are often computationally expensive or limited by strong simplifying assumptions that conflict with empirical data. To address these limitations, we introduce the Rejection-based Gillespie algorithm for non-Markovian Reactions (REGIR), a general and scalable framework for simulating non-Markovian stochastic systems with arbitrary inter-event time distributions. REGIR provides user-defined accuracy while preserving the same asymptotic computational complexity as the classical Gillespie algorithm. We derive a lower bound on REGIR's approximation accuracy and demonstrate its capabilities across three representative classes of non-Markovian systems: (1) reaction channels with delays, (2) stochastic processes driven by individual reactant properties, and (3) temporal networks governed by node activity. In all cases, REGIR accurately captures memory-dependent dynamics and outperforms existing approaches in terms of flexibility and efficiency.