我尝试以 1e-25 甚至更低的顺序测量概率。我运行一个 cpp 代码，它基本上是一个数字 ODE 求解器（使用简单的欧拉步骤）。我有很长的整数来存储单元格编号（每一步我掷骰子，我使用 64 位精度的 Mersenne Twister 随机数生成器更新单元格编号，并使用 mersenne 引擎从泊松分布中选择数字）。只要 long int 和 64 位 mersenne 和 poisson dist 就足够了，直到大约。1e19，我需要 128 位精度。我正在寻找一种可靠的方法来生成 128 位伪随机数并将它们存储为 uint_128 等，尤其是生成 128 位随机泊松分布数。
我已经在该主题中进行了搜索，但没有找到任何有用的信息，也许我可以将这些碎片拼凑在一起，但我无法...
先感谢您！