プロジェクトオイラー９

プロジェクトオイラーの問題９をpythonで解きました。
問題

projecteuler.net

a_range = range(1, 1000)
b_range = range(1, 1000)
matrix = [(row, col) for row in a_range for col in b_range if row < col ]
for a, b in matrix:
    if 1000*a + 1000*b - a*b == 500000:
        print(a*b*(1000-a-b))