繁体   English   中英

生成第n个素数的程序

[英]Program to generate nth prime number

#include<iostream.h>
#include<conio.h>
#include<stdio.h>
main()
{
clrscr();
int num,b,c,d,f,primeset[500],j=0;
cout<<"Enter Number (n): ";
cin>>num;
primeset[0]=2;
for(b=3;b<=1000;b++)
{

    for(c=2;c<b;c++)
    {
        d=b%c;
        if(d==0)
        {
            j=1;break;
        }
        else
        {
            j=0;
        }
    }

    if(j==0)
    {
        for(f=1;f<500;f++)
        {
            primeset[f]=b;
        }
    }
    else
    {
    //do nothing
    }
    }
    cout<<"n th Prime is "<<primeset[num-1];
    getch();
    return 0;
    }

这是我找到第n个素数的程序,但它不起作用,请帮帮我。 请告诉我我的错误。 请告诉我,为什么我的程序不起作用,它为 num=1 提供 2,为 num 的所有其他值提供 997。 告诉我错误和修复它们的代码。 PS:我使用涡轮增压 C++

编辑:你不喜欢 Eratosthenes 的筛子,所以这里是另一种使用动态规划的方法。 这意味着您使用已经创建的素数数组来决定下一个候选者是否是素数。

即,如果您生成了素数 [2,3,5,7],要确定第 4 个(基于 0 的)素数是您计算候选 8,9,10,11 是否可被当前素数列表整除。 你不断从 8 上升,直到一个候选人不能被一个已经存在的素数整除

代码:

void createNextPrime(int* primes, int nextPrimeIndex) {
  int primeCandidate=primes[nextPrimeIndex-1];
  do {
    primeCandidate++;
    bool found=true;
    for(int primeIndex=0;primeIndex<nextPrimeIndex;primeIndex++) {
      if (primeCandidate%primes[primeIndex]==0) {
        found=false;
        break;
      }
    }
    if (found) {
      primes[nextPrimeIndex]=primeCandidate;
      return;
    }
  }
  while(true);
}

void createNPrimes(int primeCountToCreate) {
  int primes[primeCountToCreate];
  primes[0]=2;

  for(int primeIndexToCreate=1;primeIndexToCreate<primeCountToCreate;primeIndexToCreate++) {
    createNextPrime(primes,primeIndexToCreate);
  }

  for(int i=0;i<primeCountToCreate;i++) {
    cout<< primes[i]<< " ";
  }
  cout<<endl;
}

编辑:

这是你的代码,固定工作:

int main()
{
clrscr();
int num,c,d,f,primeset[500],composite=0;
cout<<"Enter Number (n): ";
cin>>num;
primeset[0]=2;
f=1;
for(int candidate=3;candidate<=1000;candidate++)
{
    composite=0;
    for(c=2;c*c<=candidate;c++)
    {
        d=candidate%c;
        if(d==0)
        {
            composite=1;
            break;
        }
    }

    if(composite==0)
    {
      primeset[f++]=candidate;
      // Not sure what this did
        /*for(f=1;f<500;f++)
        {
            primeset[f]=b;
        } */
    }
  }
  cout<<"n th Prime is "<<primeset[num-1]<<endl;
    getch();
  return 0;
}

这是逻辑

number:-Enter Number 质数是 (2,3,5,7,11)

number=100;
nth= ? ;//enter the nth number prime u want
int a=2;
int count=0;
boolean status=true;
boolean success=false;
for(int i=a;i<=number;i++)
{
  for(int j=2;j<i;j++)
  {
     if(i/j == 0)
          status=false;
  }
  if(!status)
  {
   count++;
    if(count == nth){
        //print out nth prime number
          success=true;
         break;
    } 
  }
  if(!success){
   //nth prime is beyond the limit
  }
}

您描述的问题称为素数分布问题,并且没有公式,因此我们只能使用蛮力。

这意味着找出哪个素数的唯一方法是遍历所有值,尽管我们可以跳过所有偶数值,因为 2 是唯一的偶数素数,我们对此有一个特殊情况。

#include <iostream>

bool isPrime(unsigned int number) 
{
    for (unsigned int i = 3; (i*i) <= number; i += 2) 
    {
        if (number % i == 0) 
            return false;
    }
    return true;
}

unsigned int nthPrime(unsigned int n)
{
    // Special case for checking for 1st prime which is 2.
    if (n == 1)
        return 2;

    unsigned int counter = 1;
    unsigned int number = 1;

    while (counter < n)
    {
        number += 2;
        counter += isPrime(number) ? 1 : 0;
    }
    return number;
}

int main(int argc, char** argv)
{
    unsigned int num;
    std::cout << "Nth Prime to calculate : ";
    std::cin >> num;
    std::cout << "Prime is " << nthPrime(num) << '\n';
    return 0;
}

根据您希望使用的素数的大小,您可能会超出 32 位 unsigned int 的范围。 long long代替unsigned int会更进一步,但是,质数会很快变得很大,所以即使这样也可能不够。 最终,您需要使用没有固定/有限长度的 n 位整数。 这就是加密密钥的情况,它具有无法用 C++ 中的基本类型直接表示的巨大值。

n 位整数不在此问题的范围内,因此我将留给读者对此进行调查。

[编辑:] 请注意,您不能isPrime(...)依赖isPrime(...)函数,因为没有检查奇数/偶数或前一对素数的特殊情况。 它旨在与nthPrime(...)函数结合使用。

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM