数论——Special Prime

题目链接

题解：
没写出来，看了题解豁然开朗。
我们先化简式子得到$n^2\ast (n+p)=m^3$假设n和p不互质,则$n=k\ast p$带入式子可以得到$p^3\ast (k^3+k^2)=m^3$显然开三方过后$(k^3+k^2)$不是一个整数，所以假设不成立，n于p互质。然后我们令$n=a^3,n+k=b^3$。因为p为质数，因式分解过后$p=(a-b)(a²+ab+b²)$所以$a-b=1$所以带入得到$p=3a^2+3a+1$所以只要满足这个式子的质数就是满足的。预处理一下就可以了。

#pragma GCC optimize(2)
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+7;
int cnt,prime[N],vis[N],f[N];
void init()
{
    for(int i=2;i<N;i++){
        if(!vis[i]) prime[cnt++]=i;
        for(int j=0;j<cnt&&prime[j]*i<N;j++){
            vis[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
    for (int i = 0; prime[i] <= N; i++) {
        for (int j = 1; j <= sqrt(prime[i]); j++) {
            if (3 * j * j + 3 * j + 1 == prime[i]) f[prime[i]]++;
        }
    }
    for(int i=3;i<N;i++) f[i]+=f[i-1];
}
int main()
{
    init();
    int n;
    while(scanf("%d",&n)==1) {
        ll res = 0;
        res=f[n];
        if (!res) puts("No Special Prime!");
        else printf("%lld\n", res);
    }
}