我认为下面这段代码中的全局范围和局部范围存在一些问题。如果有人请帮助我,为什么这段代码不起作用!这是一个寻找第 25 个神奇数字的程序。神奇数字意味着它除了 2,3,5 之外没有素因数。我从 16 开始,这是第 12 个神奇的数字。然后检查这些数字是否可以与 7 中的任何素数整除。所以我在全局数组中列出素数。并检查数组中那些素数的数字。如果检查完成,程序会尝试使用 next_prime 函数找到下一个素数并将它们列出到全局数组中。但我得到了错误的结果。它总是给出不期望的结果 input+=4 。
/*a program to find out 25th magical number. magical number
means it has no prime factor except 2,3,5*/
#include<stdio.h>
#include<math.h>
double primes[10000]={2,3,5};
int serial=3;
double next_prime(double f)//function, when called, returns the next prime number
{
int j=0,loop_breaker=0;
int count;
while(j==0){
for(count=1;primes[count]<sqrt(f)+1 &&count<serial;count++){
if(fmod(f,primes[count])==0){
f+=2;
loop_breaker=1;
break;
}
}
if(loop_breaker==0){
primes[serial]=f;
serial++;
j=1;
}
}
return f;
}
int main()
{
double f=7,prime_divisor,magic_serial=12,magic_number=16;
int c,loop_breaker,is_magic;
prime_divisor=next_prime(f);
f+=2;//this line always comes after the previous line so that everytime the value of f gets changed
while(magic_serial!=25){
is_magic=0;
loop_breaker=0;
for(c=3;c<serial;c++){
if(fmod(magic_number,primes[c])==0){
loop_breaker=1;
break;
}
}
if(loop_breaker==0){
while(prime_divisor<sqrt(magic_number)+1){
prime_divisor=next_prime(f);
f+=2;
if(fmod(magic_number,prime_divisor)==0){
is_magic=1;
break;
}
}
}
if(is_magic==0){
magic_serial++;
}
magic_number++;
}
printf("%lf",magic_number);
return 0;
}