• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            素數篩法

            Posted on 2010-03-12 14:04 rikisand 閱讀(1555) 評論(0)  編輯 收藏 引用 所屬分類: C/C++

            #include <iostream> 
            #include <vector>
            #include "time.h"
            using namespace std; 
            void sieve(int n){
                vector<bool> isprime(n,true);
                vector<int> prime;
                int cnt=0;
                for(int i=2;i<n;i++){
                    if(isprime[i])cnt++,prime.push_back(i);
                    for(int t=0;t<cnt&&i*prime[t]<n;t++){
                        isprime[i*prime[t]]=false;
                        if(i%prime[t]==0)break;
                    }
                }
                /*for(int i=0;i<cnt;i++)
                    cout<<prime[i]<<" ";*/
            }
            void oldseive(int n){
                vector<bool> isprime(n,true);
                vector<int> prime;
                for(int i=2;i<n;i++){
                    if(isprime[i]){
                        prime.push_back(i);
                        for(int j=i*2;j<n;j+=i)
                            isprime[j]=false;
                    }
                }
                /*for(int i=0;i<prime.size();i++)
                    cout<<prime[i]<<" ";*/
            }
            int main(){
                clock_t start,end;
                start = clock();
                 sieve(2000000);
                 //oldseive(2000000);
                end  = clock();
                double time = double(end-start)/CLOCKS_PER_SEC;
                cout<<endl<< time<<endl;
            } 

            線性篩法sieve 1.546s oldsieve 2.875s 快了將近一倍

            old sieve 缺陷:合數可能被多次篩掉,例如 30被2,3,5篩掉了3次 然后 線性篩法限定了 任何一個合數只被它的最小質因數篩掉一次,怎么做到這一點~~

            if(i%prime[t]==0) break; 如果此時篩掉的合數從小到大找到第一個可以整除的質數,那么顯然他找到了它的最小質因數,此時我們停止搜索質數表,因為后面的質數比當前的prime[t]要大,如果我們用prime[t+n]*i 篩掉了一個合數,這個合數必然可以表述成為 prime[t]*someK  *prime[t+n] 也就是說這個合數的最小質因數也是prime[t],他應該被 prime[t]篩掉-->當程序運行到 someK*prime[t+n] 的時候~~~~

            over--------------------------------------------------------------------

            久久99热这里只有精品国产| 中文字幕亚洲综合久久2| 久久久久久久波多野结衣高潮| 久久精品国产日本波多野结衣| 99久久精品午夜一区二区| 国产综合免费精品久久久| 老男人久久青草av高清| 国产情侣久久久久aⅴ免费| 品成人欧美大片久久国产欧美| 狠狠色综合网站久久久久久久高清| 成人久久久观看免费毛片| 色播久久人人爽人人爽人人片aV | 三级韩国一区久久二区综合| 久久精品亚洲一区二区三区浴池 | 国产91色综合久久免费分享| 亚洲国产成人乱码精品女人久久久不卡 | 99久久国产综合精品女同图片| 久久se精品一区二区| 无码日韩人妻精品久久蜜桃| 久久久久国产精品麻豆AR影院| avtt天堂网久久精品| 亚洲国产欧美国产综合久久| 久久精品二区| 狠狠久久综合| 国产—久久香蕉国产线看观看| 久久久精品国产sm调教网站| 国产精品99久久久精品无码| 亚洲一区精品伊人久久伊人| 久久精品国产一区二区三区不卡| 久久人人爽人人爽人人AV东京热 | 99久久中文字幕| 精品人妻久久久久久888| 乱亲女H秽乱长久久久| 亚洲中文字幕久久精品无码APP| 欧美亚洲国产精品久久久久| 国产女人aaa级久久久级| 精品国产乱码久久久久久浪潮| 久久最近最新中文字幕大全 | 国产精品99久久久精品无码 | 亚洲国产精品一区二区久久hs| 久久精品国产色蜜蜜麻豆|