如何在Eratosthenes的分段筛中跳过偶数？

Question

我正在使用Eratosthenes的Sieve，它使用适合内存的段，而另一个段则跳过偶数。 我很难将两者结合起来。 贝娄是我的代码。 它按原样产生正确的输出，但不会跳过2的倍数。 如何跳过倍数？ 我知道我需要用i/2替换索引（例如i ），但是当我尝试获取分段错误时。 如果任何人都有关于该概念的一般信息，即使没有代码，也请分享（因此使用语言不可知标记）。

std::vector<unsigned long long> segmentedSieveOfEratosthenes(const unsigned long long limit)
{
    const unsigned long long sqrtOflimit = floor((sqrtl(limit))/2)+1;
    /*assume primesLtSqrt holds all primes from 2 to square root of limit*/
    std::vector<unsigned long long> primesLtSqrt = interpretBits(sieveOfEratosthenes(floor(sqrt(limit))+1));
    std::vector<unsigned long long> primes(primesLtSqrt);

    /*divide sieve into segments*/
    unsigned long long segBegin = sqrtOflimit; 
    unsigned long long segEnd = 2 * sqrtOflimit;

    while(segBegin < limit)
    {
        if(segEnd >= limit)
        {
            segEnd = limit;
        }

        std::vector<unsigned char> currentPrimes(sqrtOflimit + 1, true);

        for(unsigned long long i = 0; i < primesLtSqrt.size(); i++)
        {
            /*find first composite segment*/
            unsigned long long firstCompositeInSeg = floor(segBegin/static_cast<unsigned long long>(primesLtSqrt[i])) * static_cast<unsigned long long>(primesLtSqrt[i]);
            if(firstCompositeInSeg < segBegin)
            {
                firstCompositeInSeg += primesLtSqrt[i];
            }


            /*mark multiples of primesLtSqrt[i] in current segment*/
            for(unsigned long long j = firstCompositeInSeg; j < segEnd; j += primesLtSqrt[i])
            {
                currentPrimes[(j-segBegin)] = false;
            }
        }

        for(unsigned long long i = segBegin; i < segEnd; i++)
        {
            if(currentPrimes[(i - segBegin)] == true)
            {
                primes.push_back(i);
            }
        }

        segBegin = segBegin + sqrtOflimit;
        segEnd = segEnd + sqrtOflimit;
    }

    return primes;
}

我愿意接受其他改进。 随时提出建议。