求X个的质数,嗯,这个算法好简洁!

求X个的质数,嗯,这个算法好简洁
 算法P:
   求X个不同的质数,PRIME[]为一个数组,存储所有求得的质数
   P1. [初始化]PRIME[1]<-2, N<-3, J<-1(N跑遍所有奇数,因为只有奇数可以成为质数
的候选者;J用来标记至今已找到多少个质数)
   P2. [N是质数]置J<-J+1, PRIME[J]<-N
   P3. [找到X个?]若J=X,则全部要求的质数找到,退出程序
   P4. [增加N]置N<-N+2
   P5. [K<-2]置K<-2(PRIME[K]将跑遍N的所有可能的质因子)
   P6. [K/PRIME[K]?]以PRIME[K]除K;设Q是商,R是余数,若R=0(因此N不是质数),转
向P4
   P7. [PRIME[K]足够大了?]若Q<=PRIME[K],则转P2(在这种情况下,N必然是质数)
   P8. [增加K]K加1,并转P6
 
算法来自The Art of Programming,的确是大师,窥一斑而见全豹,就这么个小东西都能看出来人家的厉害之处,我以为求N是不是质数就要一个一个挨个除到N的平方根呢,这个方法可快多了。
 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
素数$p$的本原元是指一个与$p$互质的数$g$,且满足$g$的阶为$p-1$。下面是一个用C++实现的一个素数本原元的算法: ```cpp #include <iostream> #include <vector> #include <cmath> using namespace std; int mod_pow(int x, int n, int m) { int res = 1; while (n > 0) { if (n & 1) { res = (res * x) % m; } x = (x * x) % m; n >>= 1; } return res; } bool is_prime(int n) { if (n <= 1) { return false; } for (int i = 2; i <= sqrt(n); i++) { if (n % i == 0) { return false; } } return true; } vector<int> prime_factors(int n) { vector<int> factors; for (int i = 2; i <= sqrt(n); i++) { while (n % i == 0) { factors.push_back(i); n /= i; } } if (n > 1) { factors.push_back(n); } return factors; } int primitive_root(int p) { if (!is_prime(p)) { return -1; } vector<int> factors = prime_factors(p - 1); for (int g = 2; g <= p - 1; g++) { bool flag = true; for (int i = 0; i < factors.size(); i++) { if (mod_pow(g, (p - 1) / factors[i], p) == 1) { flag = false; break; } } if (flag) { return g; } } return -1; } int main() { int p; cout << "Enter a prime number: "; cin >> p; int g = primitive_root(p); if (g == -1) { cout << "No primitive root found for " << p << endl; } else { cout << "A primitive root of " << p << " is " << g << endl; } return 0; } ``` 该算法的思路是:首先判断给定的$p$是否为素数,然后出$p-1$的所有质因子,并枚举$g$从$2$到$p-1$,判断$g$是否为$p$的本原元。判断$g$是否为$p$的本原元的方法是:对于$p-1$的每个质因子$f$,计算$g^{(p-1)/f} \mod p$的值,如果都不等于1,则$g$是$p$的本原元。如果枚举完所有的$g$都无法找到$p$的本原元,则返回-1。 在幂的过程中,可以使用快速幂算法来避免超时。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值