变形得到 $a^2 = (c + b) (c - b)$，枚举 $a^2$ 的因数即可。

#include <bits/stdc++.h>
using namespace std;
int main()
{
	long long n;
	cin >> n;
	n *= n;
	for (long long i = 1; i * i < n; i++)
		if (n % i == 0 && i % 2 == n / i % 2)
			cout << (n / i - i) / 2 << ' ' << (n / i + i) / 2 << endl;
	return 0;
}