• <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>

            素?cái)?shù)篩法

            Posted on 2010-03-12 14:04 rikisand 閱讀(1544) 評(píng)論(0)  編輯 收藏 引用 所屬分類(lèi): 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 缺陷:合數(shù)可能被多次篩掉,例如 30被2,3,5篩掉了3次 然后 線性篩法限定了 任何一個(gè)合數(shù)只被它的最小質(zhì)因數(shù)篩掉一次,怎么做到這一點(diǎn)~~

            if(i%prime[t]==0) break; 如果此時(shí)篩掉的合數(shù)從小到大找到第一個(gè)可以整除的質(zhì)數(shù),那么顯然他找到了它的最小質(zhì)因數(shù),此時(shí)我們停止搜索質(zhì)數(shù)表,因?yàn)楹竺娴馁|(zhì)數(shù)比當(dāng)前的prime[t]要大,如果我們用prime[t+n]*i 篩掉了一個(gè)合數(shù),這個(gè)合數(shù)必然可以表述成為 prime[t]*someK  *prime[t+n] 也就是說(shuō)這個(gè)合數(shù)的最小質(zhì)因數(shù)也是prime[t],他應(yīng)該被 prime[t]篩掉-->當(dāng)程序運(yùn)行到 someK*prime[t+n] 的時(shí)候~~~~

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

            久久久久高潮综合影院| 久久被窝电影亚洲爽爽爽| 狠狠色丁香婷婷综合久久来| 久久久久久久久久久| 无夜精品久久久久久| 欧美麻豆久久久久久中文| 久久精品国产精品亚洲精品| 久久伊人精品一区二区三区| 无码人妻少妇久久中文字幕| 久久精品国产精品亚洲精品| 日本强好片久久久久久AAA | 久久午夜电影网| 久久精品aⅴ无码中文字字幕重口 久久精品a亚洲国产v高清不卡 | 国产精品99久久久精品无码| 伊人久久大香线蕉精品不卡| 噜噜噜色噜噜噜久久| 久久天天躁狠狠躁夜夜avapp| AV狠狠色丁香婷婷综合久久| 99久久夜色精品国产网站| 国产午夜福利精品久久| 日韩AV毛片精品久久久| 亚洲国产一成人久久精品| 久久精品国产亚洲AV麻豆网站 | 久久国产精品久久| 好久久免费视频高清| 久久国产欧美日韩精品免费| 色欲久久久天天天综合网精品| 日本精品久久久久中文字幕8| 亚洲精品无码久久久| 久久综合久久自在自线精品自| 久久久久久亚洲精品不卡 | 久久成人国产精品| 99热热久久这里只有精品68| 久久AV高潮AV无码AV| 91精品国产91久久久久久青草| 亚洲AV无码一区东京热久久| 国产精品毛片久久久久久久| 久久人人爽人人爽人人片av高请| 亚洲国产精品无码久久久久久曰| 久久国产精品久久国产精品| 亚洲欧美国产精品专区久久|