HDU - 2866 Special Prime（数论+数学证明）

HDU - 2866 Special Prime（数论+数学证明）

题目大意：给定一个区间，判断形如$n^3+n^{2}p=m^3$这样的非线性丢番图方程是否存在正整数解$(n_0,m_0)$，其中$p$为区间中的任一质数。输出这个区间中使得上述方程有解的质数的数量。
题解：若$p\mid n$，令$k=n/p$，那么有方程$k^3{p}^3+k^2{p}^3=m^3\Rightarrow p^3(k^3+k^2)=m^3$。显然，当且仅当$k^3+k^2$为立方数时，上述方程成立。又$k^3+k^2=k^2(k+1)$且$gcd(k,k+1)=1$。我们可以分两种情况讨论。
$(1):$ 对于$\forall k\ge 2$，$k^2(k+1)$可质因数分解为$p_1^{2s_1}{p}_2^{2s_2}\cdots p_r^{2s_r}{q}_1^{t_1}{q}_2^{t_2}\cdots q_w^{t_w}$，其中$k=p_1^{s_1}{p}_2^{s_2}\cdots p_r^{s_r}$，$k+1=q_1^{t_1}{q}_2^{t_2}\cdots q_w^{t_w}$，且 { p 1 , p 2 , ⋯   , p r } ∩ { q 1 , q 2 , ⋯   , q w } = ∅ \{p_1,p_2,\cdots ,p_r\} \cap \{q_1,q_2,\cdots ,q_w\}=\varnothing {p1​,p2​,⋯,pr​}∩{q1​,q2​,⋯,qw​}=∅，此时当且仅当$s_1\%3=0,s_2\%3=0,\cdots ,s_r\%3=0,t_1\%3=0,t_2\%3=0,\cdots ,t_w\%3=0$成立时，$k^3+k^2$才为立方数，上述等式的成立与$k$和$k+1$同时为立方数为等价命题，而这显然是不可能的。
$(2):$当$k=1$时，$k^3+k^2=1+1=2$，$2$不是立方数。
综上所述，当$p\mid n$时，方程$n^3+n^{2}p=m^3$不存在正整数解。

若$p\nmid n$，可以对方程左边的$n^2(n+p)$的$n^2$进行质因数分解。有$p_1^{2s_1}{p}_2^{2s_2}\cdots p_r^{2s_r}(n+p)$，若 ∃ s i ∈ { s 1 , s 2 , ⋯   , s r } , s . t .   3 ∤ s i \exist s_i \in \{s_1,s_2,\cdots,s_r\},s.t.\space 3\nmid s_i ∃si​∈{s1​,s2​,⋯,sr​},s.t. 3∤si​，那么$n+p$的质因数分解中，必有质因子$p_i^{3-(2s_{i}mod3)+x}(x\ge 0)$，否则方程右边不满足立方数的定义。令$n+p=y{p}_i^{3-(2s_{i}mod3)+x}$，那么我们有：
$$\begin{gathered} p=n+p-n \\ p=y{p}_i^{3-(2s_{i}mod3)+x}-p_1^{s_1}{p}_2^{s_2}\cdots p_r^{s_r} \\ p=p_i(y{p}_i^{2-(2s_{i}mod3)+x}-p_1^{s_1}{p}_2^{s_2}\cdots p_i^{s_i-1}\cdots p_r^{s_r}) \end{gathered}$$
说明$y{p}_i^{2-(2s_{i}mod3)+x}-p_1^{s_1}{p}_2^{s_2}\cdots p_i^{s_i-1}\cdots p_r^{s_r}$必须要等于$1$，此时有$p_i=p\Rightarrow p\mid n$，与题设矛盾。
由上述可知，$n^2$只能为一立方数，同时$n+p$也为一立方数。设$n^2={\beta }^{2\ast 3},n+p={\alpha }^3$，则有：
$$\begin{gathered} p=n+p-n \\ p={\alpha }^3-{\beta }^3 \\ p=(\alpha -\beta )({\alpha }^2+\alpha \beta +{\beta }^2) \end{gathered}$$
此时$\alpha -\beta$必须等于$1$，否则与$p$为质数的定义矛盾。
综上所述，当且仅当$p$可以表为两个相邻正整数的立方数之差时，方程$n^3+n^{2}p=m^3$有唯一一组正整数解$(n_0,m_0)$，$n_0$为那两个相邻正整数中较小的那一个。
那么问题便转化为：对于一给定区间，求区间中有多少质数能表为两个相邻正整数的立方数之差。
下面是AC代码。

#include<bits/stdc++.h>
using namespace std;
struct fastio{fastio(){ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);}}fio;
#define int long long
#define endl '\n'
#define all(x) x.begin(),x.end()
const int N=1e6;
vector<int> prime,res(N+1),judge(N+1,1);
void init(){
    judge[0]=judge[1]=0;
    for(int i=2;i<=N;++i){
        if(judge[i]==1) prime.push_back(i);
        for(int j=0;i*prime[j]<=N;++j){
            judge[i*prime[j]]=0;
            if(i%prime[j]==0) break;
        }
    }
    for(int i=1;(i+1)*(i+1)*(i+1)-i*i*i<=N;++i){
        if(judge[(i+1)*(i+1)*(i+1)-i*i*i]==1) {
            ++judge[(i+1)*(i+1)*(i+1)-i*i*i];
        }
    }
    res[0]=res[1]=0;
    for(int i=2;i<=N;++i){
        res[i]=res[i-1]+(judge[i]==2?1:0);
    }
}
signed main(){
    int l;
    init();
    while(cin>>l){
        if(res[l]==0) cout<<"No Special Prime!"<<endl;
        else cout<<res[l]<<endl;
    }
}