解法のひとつです。 幾何数列に着目するってもの考えたんだけどなあ。 いまいち無理でしたｗ 幾何数列「a, ca, c^2 a」における比率を「c」とします。 gcd(p, q) == 1 の場合、c は p/q と書くことができます。 したがって、n^2 == c^3*a^2 + a がわかります…

メインから。 時間を計るので、ちょいゴタゴタしてます。 1兆までなのでかなり時間がかかります。 いえ、力技コードなので処理時間は気にしてないですｗ #include <iostream> #include <chrono> int main() { long long limit = 1000000000000; std::chrono::system_clock::tim</chrono></iostream>…