削除された内容 追加された内容
GonJii3 (会話 | 投稿記録)
m →‎解法: ちょっとしたTIPを、追記。
GonJii3 (会話 | 投稿記録)
m →‎解法: 表記の曖昧さを、説明文の追加でカバー。
10行目:
は全てペル方程式の解になる。また逆にペル方程式の全ての解は最小解の[[冪乗]]になることが知られている。
 
最小解を得る法としては、[[連分数]]展開からの近似分数を利用する方法が良く使わ用いられる。
 
具体的には、{{math|{{sqrt|''n''}}}} の[[連分数]]展開を、{{math|{{sqrt|''n''}}}} = {{math|''A''}} = {{math|[''a''<sub>0</sub>; ''a''<sub>1</sub>, ''a''<sub>2</sub>, ..., ''a<sub>m</sub>'']}} (と置き、近似分数{{math|''aP''<sub>0</sub>''Q''}} が[[整数部分]]、{{math|''aP''<sub>1</sub>, ''aQ''<sub>2</sub>,}} ...,= {{math|''a<sub>m</sub>B''}} が循環節)とし、= {{math|[''a''<sub>0</sub>; ''a''<sub>1</sub>, ''a''<sub>2</sub>, ..., ''a''<sub>''m'' &minus; 1</sub>]}} ({{mvar|a{{sub|m}}}} が除かれている) から近似分数 {{math|''P''/''Q''}} が得られとすると、{{math|1=(''x'', ''y'') = (''P'', ''Q'')}} が解になる。但し、周期 {{mvar|m}} が奇数の場合は、右辺 = −1 の解が得られるので、1 の解を得るには上記の式で二乗する必要がある。ここで、''A'' は{{math|''a''<sub>0</sub>}} を[[整数部分]]、{{math|''a''<sub>1</sub>, ''a''<sub>2</sub>, ..., ''a<sub>m</sub>''}} を循環節とする無限連分数で、''B'' は循環節を一周期だけ採り、最後の数{{mvar|a{{sub|m}}}}を除いた、有限連分数である。ちなみに、{{mvar|a{{sub|m}}}} = {{mvar|2a{{sub|0}}}} が常に成立する。
 
例えば {{mvar|n}} が {{math|7}} ならば、{{math|1={{sqrt|7}} = [2; 1, 1, 1, 4]}} (周期は 4 で偶数) なので、{{math|[2; 1, 1, 1]}} から近似分数 {{math|8/3}} が得られ、{{math|1=(''x'', ''y'') = (8, 3)}} が最小解となる。{{mvar|n}} が {{math|61}} の場合は 、<math>\sqrt{61} = [7; 1, 4, 3, 1, 2, 2, 1, 3, 4, 1, 14]</math> (周期は 11 で奇数) なので近似分数 29718/3805 が得られ、右辺 = −1 の最小解は <math>(x_1, y_1) = (29718, 3805)</math> となる。右辺 = 1 の最小解は、<math>x + y\sqrt{61} = (x_1 + y_1\sqrt{61})^2\,</math> から {{math|1=(''x'', ''y'') = (1766319049, 226153980)}} となる。