//使用dev-cpp
#include <cstdlib>
#include <iostream>

using namespace std;

int main(int argc, char *argv[])
{
    int p,q,r;
    long long h,s,g;
    //1=1/p+1/q+1/r+1/s
    for(p=1;p<=1000;p++)
    {
        for(q=1;q<=1000;q++)
        {
            for(r=1;r<=1000;r++)
            {
                if((p!=q)&&(p!=r)&&(q!=r))
                {
                h=p*q*r-q*r-p*r-p*q;g=p*q*r;
                if((g%h==0)&&(h!=0))
                {
                    s=g/h;
                    cout<<p<<"---"<<q<<"---"<<r<<"---"<<s<<endl;
                }   
                }
               
            }
        }
    }
   
    cout<<endl;
    system("pause");
    return 0;
}

这段程序是干啥的  目测没测出来啥

就是一道数学练习题

终于找到原因了，谢谢大家的热情帮助

啥原因  贴出来看看

求 1/p + 1/q + 1/r + 1/s = 1 的正整数解
因为加法交换律的存在，所以只要求 p<=q<=r<=r 条件下的解就可以了

因为p最小，所以1/p最大，所以1/p应该大于等于1/4，即 2<=p<=4
当p为最小值2时，1/q+1/r+1/s=1-1/2。因为1/q最大，所以 1/q >= 1/3×(1-1/2)，即 p<=q<=6
当p为最小值2，q为最小值3时，1/r+1/s=1-1/2-1/3。因为1/r最大，所以 1/r >= 1/2×(1-1/2-1/3)，即 q<=r<=12

代码如下：

#include <iostream>
using namespace std;

int main()
{
    // 1=1/p+1/q+1/r+1/s  ===>>> 1/s = (abc-ab-bc-ca) / abc
    for( int p=2; p<=4; ++p )
    {
        for( int q=p; q<=6; ++q )
        {
            for( int r=q; r<=12; ++r )
            {
                int h = p*q*r-q*r-p*r-p*q;
                int g = p*q*r;

                if( h>0 && g%h==0 && g/h>=r )
                {
                    int s = g/h;
                    cout<<p<<"---"<<q<<"---"<<r<<"---"<<s<<'\n';
                }
            }
        }
    }
    cout << endl;

    return 0;
}