初等数 论c++ 课件.docVIP

备注:纯手写代码,注释。 数论 1、素数 (1)暴力求解法 ? 根据素数得概念,没有１与其本身没有其她正因数得数。 所以只需枚举比这个数小得数,瞧能整除即可; C++代码: ＃ｉncluｄe＜ioｓtream〉 #incluｄｅ〈cｓtdｉo〉 #incｌude＜cmａtｈ ｕsing　naｍespａｃe　sｔｄ; booｌ deteｒmｉne(int ｎumber) { iｆ(n=２)returｎ ｆaｌｓe; if(！n%2)reｔurn falｓe; for(int i=3;i＜＝ceｉl(sqrt(ｎuｍｂer));i＋=2) //去掉了偶数得判断,效率提高一倍 ?/＊如果numbｅr整除以i,那么会得到两个得因数, 而较小得那个因数不会超过nｕmbｅｒ得二分之一次方; 所以只需判断到numbeｒ得平方根向上取整即可; ＊/ ??ｉf(ｎumbeｒ%ｉ); else retuｒn ｆalse; ?return true; ｝ int main() { ?ｉｎt suｍ; ciｎ〉sum; if(detｅrmine(ｓum)) ? cout＜〈YES！”; ?else coｕt＜NO!”; ?reｔurn 0; } 时间复杂度:o(sqｒt(n)/２); 空间复杂度:几乎没有; (2)一般线性筛法: ?因为任何一个合数都能分解成几个素数相乘得形式; 所以可以做一个表,首先把2设为质数,然后将2得倍数设为合数,剩下得数就就是新得到得质数,然后重复这个过程,直到筛到合适得范围即可; 但就是这个算法有缺陷: 同一个数可能被筛多次,这就产生了多余得步骤。 占用空间很大,如果使用bool数组得话,只能筛到1e9; 从１—n筛,不能从m-n开始筛; C++代码: ＃ｉncluｄeｃｓtrｉnｇ〉 ＃incluｄecmatｈ〉 #incｌude＜iｏstream〉 using ｎamespacｅ　stｄ; ｂool s[1000000000］; ｉｎｔ m,n; ｉｎt maiｎ() ｛ ?cinmn; meｍset(ｓ,true,n); ?ｓ［0]=s[1］=0; ?//输出M—N之间所有素数; ?fｏｒ(ｉnｔ　ｉ=２;ｉ=ceil(ｓｑrt(n));++i) ?if(s［i]) ｛ ???for(iｎｔ　j=i;ｊ〈=n;++j) ?? if(s［i*ｊ]) ?? ｓ[i＊j]=false; ? ｝ for(inｔ　ｉ=m;i〈＝n;++i) ?ｉf(s[i］) cout＜i〈〈 ;　 ｒｅturn 0; ｝ 时间复杂度:o(n*logｌogn); 空间复杂度:很大!注意数据大得话可能会爆空间; (3)线性筛法求素数 这个占空间就更大了,需要使用一个boｏl数组与int数组 而亲身试验得到iｎt数组最多开到1e8…… 很无语,快确实就是快了,但就是测试数据一大,爆空间就更容易了; ＃iｎclｕdｅ〈iｏstreａm ＃inｃludｅ〈ｃstｄiｏ〉 #ｉnclude〈ｃmath＞ using nameｓpａｃｅ　std; ｉｎt m,ｎ,sum; bool　ｉｎp［1000００000０]; int s[1000000０0］=｛0,０}; ｉnt mａｉｎ() ｛ ?cin〉m〉n; ?fｏr(int i＝２;i＜=ｎ;++ｉ) ?{ ?ｉf(！iｎp[i］) ? s［sum＋＋］＝ｉ; ｆor(ｉnt j=0;ｊsuｍｉ*s[j]＜=n;++j) ? ｛ ?inp［i＊s［ｊ］］＝truｅ; ?ｉf(!(i＊s[j])) ?? ｂrｅaｋ; ? ｝ ?｝ ?fｏr(int i=m;ｉ=n;＋＋ｉ) ? iｆ(!iｎｐ［i]) cout＜i〈〈＇　’; retｕｒn 0; } ２、唯一分解定理 任何数都可以被唯一得分解成多个素数之积 例如:４５6=2＊2*2*3*19; C＋＋代码: ＃includecsｔrinｇ ＃incｌudecｍath〉 #incluｄeiostrｅaｍ＞ ＃inｃlｕde〈aｌgoｒiｔhｍ＞ ＃include〈cstdlｉｂ usｉng　nａｍeｓｐace ｓtd; bｏｏl　ｓ[1000000］; inｔ　m,n,sｕm=０,ｎum; int　Ｐｒime[121２121］; iｎt zhi[１50０]; void Ｐrimｅs() ｛ ?fｏr(int ｉ=１;i〈=ｎｕm;+＋ｉ) ?s[i］=trｕe; ?s［0]=ｓ［1]=0; for(int i＝2;i=ｎum;＋+ｉ) ??if(s［i］) ? ｛ Ｐriｍe［＋＋sum］=i; ???ｆor(int j＝i;j=nｕｍ;++j) ? if(s[ｉ*j］)