关闭

神奇的线性筛法

196人阅读 评论(0) 收藏 举报
分类:

关于线筛

所谓线性筛法,顾名思义它是线性的筛素数的方法,所以是O(N)

神奇的东西

int getprime(int lim)
{
    int i,j,num=0;
    fo(i,2,lim) 
    {
        if (!bz[i])
        {
            prime[++num]=i;
            for(j=1;i*j<=lim;j++) bz[i*prime[j]]=1;
        }   
    }
}

这是普通筛法

int getprime(int lim)
{
    int i,j,num=0;
    fo(i,2,lim) 
    {
        if (!bz[i]) prime[++num]=i;
        for(j=1;i*prime[j]<=lim&&j<=num;j++)
        {
            bz[i*prime[j]]=1;
            if (!(i%prime[j])) break;
        } 
    }
}

这是线性筛法

实测

某神犇已经测过,我直接上结果

  • 筛 [0, 100000) 范围内的素数
    第一种素数筛法 0 毫秒
    第二种素数筛法 0 毫秒

  • 筛 [0, 200000) 范围内的素数
    第一种素数筛法 15 毫秒
    第二种素数筛法 0 毫秒

  • 筛 [0, 400000) 范围内的素数
    第一种素数筛法 16 毫秒
    第二种素数筛法 15 毫秒

  • 筛 [0, 800000) 范围内的素数
    第一种素数筛法 47 毫秒
    第二种素数筛法 16 毫秒

  • 筛 [0, 1600000) 范围内的素数
    第一种素数筛法 62 毫秒
    第二种素数筛法 63 毫秒

  • 筛 [0, 3200000) 范围内的素数
    第一种素数筛法 297 毫秒
    第二种素数筛法 109 毫秒

  • 筛 [0, 6400000) 范围内的素数
    第一种素数筛法 922 毫秒
    第二种素数筛法 266 毫秒

  • 筛 [0, 12800000) 范围内的素数
    第一种素数筛法 2187 毫秒
    第二种素数筛法 563 毫秒

  • 筛 [0, 25600000) 范围内的素数
    第一种素数筛法 4828 毫秒
    第二种素数筛法 1187 毫秒

为什么是线性的呢?

回到程序

可以发现,普通筛法将每个质数的倍数都筛掉,每个合数都被筛了它的质因子次
That’s too slow!!

我们再看线性筛法
筛的时候,把当前到的这一个数(不论质数合数)乘上筛出来的每个质数筛一遍

关键是这一句

if (!(i%prime[j])) break;

这样的情况,意思是prime[j]这一个质数在prime[j]×i中指数已经大于等于2了。

并且,这个合数的任何倍数以后一定会被prime[j]筛到
所以后面的都不必筛了。

这样,每个合数都只被它的最小的质因子筛过一遍,所以是线性的复杂度。

1
0
查看评论

求素数: 一般线性筛法 + 快速线性筛法

From: http://blog.csdn.net/dinosoft/article/details/5829550 素数总是一个比较常涉及到的内容,掌握求素数的方法是一项基本功。 基本原则就是题目如果只需要判断少量数字是否为素数,直接枚举因子2 。。N^(0.5) ,看看能否整除N。 如...
  • JoeBlackzqq
  • JoeBlackzqq
  • 2015-03-11 14:53
  • 1139

几种线性筛法

几种线性筛法 本文出自blog.csdn.net/svitter 写在之前 maker 关于线性筛素数 之前一直没有正视线性筛素数的问题。今天特意来写一个伪证明。 如果当前的i不是素数,那么必然被之前的某个素数筛掉了。i × prime[j]。 一个合数必然可以写成几个素数的乘积...
  • svitter
  • svitter
  • 2014-08-05 15:46
  • 1115

bzoj1416&1498: [NOI2006]神奇的口袋

传送门 傻逼题。 直接计算方案数,除上总方案数就可以了。 注意要用高精度。 为啥oj上有两道题?#include<cmath> #include<cstdio> #include<cstring> #include<iostream> #inc...
  • zhouyuyang233
  • zhouyuyang233
  • 2017-04-29 15:42
  • 249

一段神奇的代码

做一个优雅的程序员,好像是很多程序员的一大目标,何为优雅,如何优雅?下面是一段神奇的代码,什么编码风格,什么编程规范,统统丢到一边吧,化腐朽为神奇才是更加难能可贵的! #define M 002354l #define A 000644l #define G 00...
  • u011836061
  • u011836061
  • 2015-06-04 12:55
  • 1099

浅析线性筛法

今天晚上,在家里自学了一下线性筛法,还是比较容易理解的,并没有想象的那么难。 所谓线性就是O(n)的效率咯,线性筛法较快的一个本质就是他不会干重复的事,例如12——普通的筛法在素数2时候筛一次,素数3时筛一次,很明显这些是不必要的选择,我们可以舍去。 所以,在线性筛法里的优化就是...
  • John_pascal
  • John_pascal
  • 2016-08-05 22:00
  • 238

记一次问题解决历程——脚本的神奇功能

今天一个测试的同事问我,你对adb熟悉吗?我说简单的会。然后说adb reboot(手机重启指令)知道吧?我说知道的。然后她的问题是:能不能一直让手机执行重启命令:我想了下,用脚本应该是可以的。之前接触过shell脚本,感受到了shell脚本的神奇之处,但是在Dos窗口下,对于bat脚本我是一颗大白...
  • Simon_Crystin
  • Simon_Crystin
  • 2017-04-01 20:42
  • 605

{算法}乱谈线性筛法

指在线性时间O(n)用筛选的方法求一些东西的算法。 伪目录 1. 求质数表 2. 筛出欧拉函数φ As we all know......
  • FarmerJohnOfZS
  • FarmerJohnOfZS
  • 2016-06-17 20:31
  • 2474

百练2755 神奇的口袋 【深搜】or【动规】or【普通递归】or【递推】

总Time Limit:  10000ms  Memory Limit:  65536kB 有一个神奇的口袋,总的容积是40,用这个口袋可以变出一些物品,这些物品的总体积必须是40。John现在有n个想要得到的物品,每个物品的体积分别是a1,a2...
  • u012846486
  • u012846486
  • 2014-07-19 18:39
  • 1181

noi2006神奇口袋

神奇的口袋,考数学甚于考算法,【考虑无自由取球情况,那么,设A(X,Y)为第x次,取出颜色为Y的球,P{A(1,Y1),A(2,Y2)…}=P{A(1,Y1),A(2,Y1),A(3,Y1)}因为首先分母是确定的,这毋庸置疑的,而分母取值只与在它之前取过多少同颜色球有关,与其他颜色球无关,所以可以进...
  • huyuncong
  • huyuncong
  • 2012-02-15 18:51
  • 1638

Poj 2755:神奇的口袋

总时间限制: 10000ms 内存限制: 65536kB 描述有一个神奇的口袋,总的容积是40,用这个口袋可以变出一些物品,这些物品的总体积必须是40。John现在有n个想要得到的物品,每个物品的体积分别是a1,a2……an。John可以从这些物品中选择一些,...
  • qq_34125999
  • qq_34125999
  • 2016-05-19 21:34
  • 570
    个人资料
    • 访问:72888次
    • 积分:3281
    • 等级:
    • 排名:第12254名
    • 原创:254篇
    • 转载:2篇
    • 译文:0篇
    • 评论:21条