#include<stdio.h>
int main(){
    long long int p=1,q=1,count=0;
    while(1){
        while(p*p<(2*q*q+1)){
            p++;
        }
        if(p*p==(2*q*q+1)){
            count++;
            printf("(%d,%d) ",p,q);
        }
        q++;
        if(count==10) break;
    }
    return 0;
}
我只会这么写，但这么写根本输出不了那么多。。。而且也不是递归。。。。

你用 count==10 限制了只处理10对数据， 当然不会出现50个

printf("(%d,%d) ",p,q);
// 改为
printf("(%lld,%lld) ",p,q);

[此贴子已经被作者于2020-5-30 21:10编辑过]

编译你的代码报错：
printf("(%d,%d) ",p,q);
warning: format '%d' expects argument of type 'int', but argument 2 has type 'long long int'

另外，你题目中说50个，怎么代码中却是 if(count==10) break ？

您好，改完之后还是出现不了50个

我现在有几个问题

1. 你的这个方程式的解，到后面数大概多大，你自己知道么？（找到第12个解，就已经9位数了，那按照你写的这个方法，要循环多少次才能找到50个答案，你需要等多久？）
2. 以你目前的方法，根本不要考虑改成递归，压栈会死掉的

所以我的目的就是来论坛问问大佬们，用什么方法做，我知道我这个算法不行的，所以来请教下

路过，纯吃瓜群众，我只是想看看 long long int 能不能容得下前50个解。

int main( void )
{
    typedef unsigned long long TYPE;
   
    TYPE d = 2;
    TYPE x1=3, y1=2;
    TYPE xn=x1, yn=y1;
   
    for( size_t i=0; i!=50; ++i )
    {
        printf( "%zu: %llu, %llu\n", i+1, xn, yn );

        if( xn*x1/x1!=xn || d*yn*y1/y1/yn!=d || xn*x1+d*yn*y1<xn*x1 )
        {
            puts( "overflow" );
            break;
        }

        TYPE x = xn*x1 + d*yn*y1;
        TYPE y = xn*y1 + yn*x1;
        xn=x, yn=y;
    }
}

结果是只能容纳前25组解，输出结果是

然后我将 unsigned long long 换成 非标准的 __uint128_t，正巧能容纳下第50组，接下来就溢出了。

p只可能是奇数，所以p可以+=2
话说1    0算一对解么

玩数论的时候一般对0有歧视，只玩正整数，所以不算。

了解了