相同算法的C ++和Python版本给出不同的结果

Question

以下代码是一种算法，用于确定最大三角形的边小于等于MAX的整数三角形的数量，这些三角形的中位数为整数。 Python版本可以使用，但是对于较大的N来说太慢了，而C ++版本则快了很多，但没有给出正确的结果。

当MAX为10时，C ++和Python都返回3。

当MAX为100时，Python返回835，而C ++返回836。

当MAX为200时，Python返回4088，而C ++返回4102。

当MAX为500时，Python返回32251，而C ++返回32296。

当MAX为1000时，Python返回149869，而C ++返回150002。

这是C ++版本：

#include <cstdio>
#include <math.h>

const int MAX = 1000;

int main()
{
    long long int x = 0;
    for (int b = MAX; b > 4; b--)
    {
        printf("%lld\n", b);
        for (int a = b; a > 4; a -= 2){
            for (int c = floor(b/2); c < floor(MAX/2); c+=1)
            {
                if (a+b > 2*c){
                    int d = 2*(pow(a,2)+pow(b,2)-2*pow(c,2));
                    if (sqrt(d)/2==floor(sqrt(d)/2))
                        x+=1;
                }
            }
            }
    }
    printf("Done: ");       
    printf("%lld\n", x);
}

这是原始的Python版本：

import math

def sumofSquares(n):
    f = 0
    for b in range(n,4,-1):
        print(b)
        for a in range(b,4,-2):
            for C in range(math.ceil(b/2),n//2+1):
                if a+b>2*C:
                    D = 2*(a**2+b**2-2*C**2)
                    if (math.sqrt(D)/2).is_integer():
                        f += 1
    return f

a = int(input())
print(sumofSquares(a))
print('Done')

我对C ++不太熟悉，所以我不知道是什么原因导致的（也许是溢出错误？）。

当然，欢迎对该算法进行任何优化！

Answer 1

问题是您的c （在python中为C ）变量的范围不匹配。 为了使它们等效于您现有的C ++范围，您可以将python循环更改为：

for C in range(int(math.floor(b/2)), int(math.floor(n/2))):
    ...

为了使它们等效于您现有的python范围，您可以将C ++循环更改为：

for (int c = ceil(b/2.0); c < MAX/2 + 1; c++) {
    ...
}

根据最初正确的循环，这将使结果匹配。

Answer 2

它似乎在这里带来一些麻烦：

(sqrt(d)==floor(sqrt(d)))