素数

比赛中常用的素数，这里总结一下各方面 ，先贴一个搞笑的打印素数

貌似是线性级别的高效。。

/*遇到素数需要打表时，先估算素数的个数：num = n / lnx;num为大概数字，越大误差越小（只是估计，用于估算素数表数组大小）这个打表法效率貌似很高，网上说几乎达到了线性时间（不知道是真是假=。=）*/#include<iostream>#include<cstring>#include<cstdio>#include<ctime>#include<algorithm>using namespace std;int n;bool visit[10100000];int prime[10000000];void init_prim(){memset(visit, true, sizeof(visit));int num = 0;for (int i = 2; i <= n; ++i){if (visit[i] == true){num++;prime[num] = i;}for (int j = 1; ((j <= num) && (i * prime[j] <= n));  ++j){visit[i * prime[j]] = false;if (i % prime[j] == 0) break; //点睛之笔}}}int main(){memset(prime, 0, sizeof(prime));int count = 0;cin>>n;init_prim();for(int i = 0; i <= n; ++i)if(prime[i]){cout<<prime[i]<<" ";count++;}cout<<endl;cout<<"素数个数为："<<count<<endl;}        

 

判断素数

bool isPrime(int num){    if (num == 2 || num == 3)    {        return true;    }    if (num % 6 != 1 && num % 6 != 5)    {        return false;    }    for (int i = 5; i*i <= num; i += 6)    {        if (num % i == 0 || num % (i+2) == 0)        {            return false;        }    }    return true;}

 

 

最大公约数

int kgcd(int a, int b){    if (a == 0) return b;    if (b == 0) return a;    if (!(a & 1) && !(b & 1)) return kgcd(a>>1, b>>1) << 1;    else if (!(b & 1)) return kgcd(a, b>>1);    else if (!(a & 1)) return kgcd(a>>1, b);    else return kgcd(abs(a - b), min(a, b));