在 C++ 中处理巨大数组的最佳方法是什么？

我正在尝试模拟一些等离子体现象，这转化为模拟大量粒子的动力学。通常的方法是将许多粒子聚集成一些“宏观粒子”，这样程序就不必循环处理 N 个粒子，而是循环处理一些 n

然而，即使在某些情况下这样做，我也需要使用 1e22 大粒子。我需要将这些 1e22 粒子的位置存储在某种形式的数组中

std::vector<double> pos(Npart)

其中 Npart 是输入的一部分。由于没有 C++ 类型可以在 Npart 中存储如此大的整数，我想知道在这些问题中使用的常用策略是什么。也许用更少的粒子定义许多数组？这里有什么好的做法？