我正在编写一个找到素数总和的程序。它必须使用重定向输入。我已经编写了它,以便它找到输入的最大数字,然后将其用作第 n 个素数。然后它使用第 n 个素数来设置数组的大小。它一直有效,直到我尝试打印总和。我不知道为什么我在所有地方都出现了段错误。我想我已经用 malloc 正确分配了数组。为什么在我使用数组时不会在 printf 上发生故障?也欢迎对我的代码提出任何建议。
EDIT 使用了 2000 表格 1 到 2000 的测试输入并且它有效,但是 10000 表格 1 到 10000 崩溃的完整测试文件仍在调查原因。我猜我没有分配足够的空间
编辑我的问题是在我的筛子上我没有使用 sqrt(nthprime) 所以它找到了更多的素数然后数组可以容纳
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
int nprime (int max);
void sieve_sum ( int *primes, int nthprime, int tests,int *input);
int main(void)
{
int i=0;
int max=0; //largest input
int tests; //number of tests
int nthprime; //estimated nth prime
int *primes; //array of primes
int *input; // numbers to put in to P(n), where p(n) is the summation of primes
scanf("%d",&tests); //gets number of tests
input = malloc(sizeof(int)*tests);
//test values
for(i=0; i<=tests-1; i++)
scanf("%d",&input[i]);
//finds max test value
i=0;
for (i = 0; i < tests; i++ )
{
if ( input[i] > max-1 )
max = input[i];
}
// calls nprime places value in n
nthprime = nprime(max);
primes = malloc(sizeof(int)*nthprime);
// calls sieve_sum
sieve_sum( primes, nthprime, tests, input);
//free memory
free(input);
free(primes);
return 0;
}
//finds Primes and their sum
void sieve_sum ( int *primes, int nthprime, int tests,int *input)
{
int i;
int j;
//fills in arrays with 1's
for(i=2; i<=nthprime; i++)
primes[i] = 1;
//replaces non primes with 0's
i=0;
for(i=2; i<=sqrt(nthprime); i++)
{
if(primes[i] == 1)
{
for(j=i; (i*j)<=(nthprime); j++)
primes[(i*j)] = 0;
}
}
//rewrites array with only primes
j=1;
i=0;
for(i=2; i<=nthprime; i++)
{
if(primes[i] == 1)
{
primes[j] = i;
j++;
}
}
//sums
i=0;
for ( i=1; i<=tests; i++ )
{
int sum=0;//sum of primes
j=0;
for(j=1; j<=input[i-1]; j++)
{
sum = primes[j] + sum;
}
printf("%d\n", sum );
}
return 0;
}
//finds the Nth number prime
int nprime (int max)
{
//aproximization of pi(n) (the nth prime) times 2 ensures correct allocation of memory
max = ceil( max*((log (max)) + log ((log (max)))))*2;
return (max);
}
示例输入文件:
20
1
2
3
4
5
6
7
8
9
10
10
9
8
7
6
5
4
3
2
1
示例输出应为:
2
5
10
17
28
41
58
77
100
129
129
100
77
58
41
28
17
10
5
2