【算法】寻找1000000000（十亿）内素数并统计个数

最新推荐文章于 2023-07-08 19:42:00 发布

Rand Tsui

最新推荐文章于 2023-07-08 19:42:00 发布

阅读量9.7k

点赞数 1

分类专栏：数据结构+算法文章标签：算法质数素数十亿

本博客所有公开的原创博文皆允许转载，无须再征求作者本人（也就是我）同意

本文链接：https://blog.csdn.net/for_cxc/article/details/51838811

版权

数据结构+算法专栏收录该内容

19 篇文章 0 订阅

订阅专栏

已优化到 22 秒内，不知道世界最快算法在同样环境下运行速度是多少，下面贴出代码，C#

using System;

namespace prime
{
    class Program
    {
        static void Main(string[] args)
        {
            DateTime dt = DateTime.Now;
            int n = 1000000000;
            Console.WriteLine(Get_Prime_Count(n));
            DateTime dt2 = DateTime.Now;
            Console.WriteLine(dt2-dt);
        }
        static int Get_Prime_Count(int n)
        {
            int count = 0;  //后面用于统计素数个数
            bool[] big = new bool[n + 1]; //用于标记是否为素数

            int sl = 2 * 3 * 5 * 7 * 11 * 13 * 17 * 19;
            bool[] small = new bool[sl + 1]; //用于筛选掉13及以内质数的倍数

            int i;

            for (i = 2; i <= sl; i += 2)
                small[i] = true;
            for (i = 3; i <= sl; i += 3)
                small[i] = true;
            for (i = 5; i <= sl; i += 5)
                small[i] = true;
            for (i = 7; i <= sl; i += 7)
                small[i] = true;
            for (i = 11; i <= sl; i += 11)
                small[i] = true;
            for (i = 13; i <= sl; i += 13)
                small[i] = true;
            for (i = 17; i <= sl; i += 17)
                small[i] = true;
            for (i = 19; i <= sl; i += 19)
                small[i] = true;

            //以下通过拷贝去掉big中所有13及以内质数的倍数
            int remainder = n % sl; //余数
            int p = 1;
            int q = p + n - remainder;
            while (p < q)
            {
                Array.Copy(small, 1, big, p, sl);
                p += sl;
            }
            if (remainder > 0)
            {
                Array.Copy(small, 1, big, p, remainder);
            }

            //下面这几个不能漏掉
            big[2] = false;
            big[3] = false;
            big[5] = false;
            big[7] = false;
            big[11] = false;
            big[13] = false;
            big[17] = false;
            big[19] = false;

            int sqrt = (int)Math.Sqrt(n);//平方根
            int step,j;
            sqrt -= 2;
            for (i = 23; i <= sqrt; i += 6)
            {
                if (!big[i])
                {
                    step = i * 6;
                    for (j = i * i; j <= n; j += step)
                    {
                        big[j] = true;
                    }
                    for (j = i * i + 2 * i; j <= n; j += step)
                    {
                        big[j] = true;
                    }
                }
                int mm = i + 2;
                if (!big[mm])
                {
                    step = mm * 6;
                    for (j = mm * mm; j <= n; j += step)
                    {
                        big[j] = true;
                    }
                    for (j = mm * mm + 4 * mm; j <= n; j += step)
                    {
                        big[j] = true;
                    }
                }
            }

            //统计素数个数
            for (i = 7; i <= n-4; i += 6)   //100010  从7开始
            {
                if (!big[i]) count++;
                if (!big[i+4]) count++;
            }
            count+=3; //还有2、3、5
            return count;
        }
    }
}