这道题简单来做的话就是找规律,虽然我没找出来。。。。。。。。
还有就是严格的数学推到,虽然我不会。。。。。。。
下面是代码
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll MAX = 1e6;
ll prime[MAX + 10] = { 0 };
void init()
{
prime[0] = prime[1] = 1;
for (ll i = 2; i <= MAX; i++)
{
if (!prime[i])
{
for (int j = i * 2; j <= MAX; j += i)
prime[j] = j;
}
}
}
ll n;
int main()
{
init();
while (cin >> n)
{
int cnt = 0;
for (ll i = 1; i * 3 + 3 * i * i + 1 <= n; i++)
{
if (!prime[i * 3 + 3 * i * i + 1]) cnt++;
}
if (cnt)
cout << cnt << endl;
else cout << "No Special Prime!\n";
}
return 0;
}