public class Atkin_Algo : IEnumerable<ulong>
{
private List<ulong> primes;
private ulong limit;
public Atkin_Algo(ulong _limit)
{
limit = _limit;
primes = new List<ulong>();
}
public IEnumerator<ulong> GetEnumerator()
{
if (!primes.Any())
Find_Primes();
foreach (var p in primes)
yield return p;
}
IEnumerator IEnumerable.GetEnumerator()
{
return GetEnumerator();
}
private void Find_Primes()
{
var is_prime = new bool[limit + 1];
var sqrt = Math.Sqrt(limit);
for (ulong x = 1; x <= sqrt; x++)
{
for (ulong y = 1; y <= sqrt; y++)
{
var n = 4 * x * x + y * y;
if (n <= limit && (n % 12 == 1 || n % 12 == 5))
{
is_prime[n] ^= true;
}
n = 3 * x * x + y * y;
if (n <= limit && n % 12 == 7)
{
is_prime[n] ^= true;
}
n = 3 * x * x - y * y;
if (x > y && n <= limit && n % 12 == 11)
{
is_prime[n] ^= true;
}
}
}
for (ulong n = 5; n <= sqrt; n++)
{
if (is_prime[n])
{
var s = n * n;
for (ulong k = s; k <= limit; k += s)
{
is_prime[k] = false;
}
}
}
primes.Add(2);
primes.Add(3);
for (ulong n = 5; n <= limit; n += 2)
{
if (is_prime[n])
{
primes.Add(n);
}
}
}
}
所以我的问题是,当我想生成一个LARGE列表时,我得到OutOfMemoryException这是预期的,但我希望能够解决这个问题。在得到任何建议之前,我还没有做过这样的事情。
我想避免这个问题的原因是很快就会使用BigInteger类并能够生成大型素数。
提前谢谢。
答案 0 :(得分:1)
关键的限制是必须保留大量的bool。您可能正在运行数组大小的限制;而不是一个正确的数组,尝试使用多个数组。您可以编写函数以使其像数组一样易于使用,但在幕后它由多个数组支持(可能一个用于0-1百万的数字,另一个用于1mio-2mio,依此类推)。
然后,当这些数组对于内存来说真的太大时,你可以将它们保存到磁盘并一次加载一个你需要的那个(显然这会比内存访问慢得多。幸运的是,很多访问应该与上次相同的数组,所以,如果你保持它的帮助)。
通过将结果素数保存在内存中的大列表中,您也会给自己带来麻烦。你最好把它们写到磁盘上(也许在屏幕上),就像你找到它们一样,然后不要把它们留在内存中。