###### Project Euler：Problem 58 Spiral primes

Starting with 1 and spiralling anticlockwise in the following way, a square spiral with side length 7 is formed.

37 36 35 34 33 32 31
38 17 16 15 14 13 30
39 18  5  4  3 12 29
40 19  6  1  2 11 28
41 20  7  8  9 10 27
42 21 22 23 24 25 26
43 44 45 46 47 48 49

It is interesting to note that the odd squares lie along the bottom right diagonal, but what is more interesting is that 8 out of the 13 numbers lying along both diagonals are prime; that is, a ratio of 8/13 ≈ 62%.

If one complete new layer is wrapped around the spiral above, a square spiral with side length 9 will be formed. If this process is continued, what is the side length of the square spiral for which the ratio of primes along both diagonals first falls below 10%?

#include <iostream>
#include <string>
using namespace std;

int cp[100000000];

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

void count_prime(unsigned long long n)
{
cp[n] = cp[n - 1];
int a[3];
a[0] = (2 * n + 1)*(2 * n + 1) - 4 * n;
a[1] = (2 * n + 1)*(2 * n + 1) - (2 * n + 1) + 1;
a[2] = (2 * n + 1)*(2 * n + 1) - 6 * n;
for (int i = 0; i < 3; i++)
{
if (isPrime(a[i]))
cp[n]++;
}
}

int main()
{
memset(cp, 0, sizeof(cp));
cp[0] = 0;
unsigned long long ans;
double a, b, res;
for (unsigned long long i = 1; i < 100000000; i++)
{
count_prime(i);
a = cp[i] * 1.0;
b = (4 * i + 1)*1.0;
res = a / b*1.0;
cout << res << endl;
if (res < 0.10)
{
ans = 2 * i + 1;
break;
}
}
cout << ans << endl;
system("pause");
return 0;
}


#### Project-Euler problem 1-50

2015-01-31 00:09:48

#### Project Euler Problem 27 Quadratic primes

2017-03-30 19:46:48

#### Project Euler Problem 71-80

2016-02-27 20:34:47

#### Project Euler:Problem 8

2015-05-30 12:16:18

#### Project Euler problem 66

2013-10-01 18:19:27

#### 欧拉工程第47题：Distinct primes factors

2015-06-03 19:30:55

#### project euler problem 19 数周末

2013-10-07 15:02:21

#### 【Project Euler】8 第八题

2015-02-07 12:26:00

#### project euler 102

2015-12-23 22:19:58

#### Project Euler Problem 58

2009-10-24 20:01:00

## 不良信息举报

Project Euler：Problem 58 Spiral primes