モジュラ逆数

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:Expand English テンプレート:Refimprove 合同算術において、 テンプレート:Math に関する整数 テンプレート:Mvarモジュラ逆数(モジュラぎゃくすう、テンプレート:Lang-en-short)とは,

axxa1(modm)

という関係にある整数 テンプレート:Mvar の属する合同類(あるいはその標準的な代表元)のことで、これを テンプレート:Math で表す。これは、整数の法 テンプレート:Mvar に関する合同類環 /m における乗法逆元である。ある種の応用においては、モジュラ逆数 テンプレート:Math/m に属さないような場合を考えることもある。

テンプレート:Mvarテンプレート:Mvar を法とする逆数が存在するための必要十分条件テンプレート:Mvarテンプレート:Mvar とが互いに素(即ち、最大公約数 テンプレート:Mathテンプレート:Math)となることである。法 テンプレート:Mvar に関する テンプレート:Mvar のモジュラ逆数が存在するならば、テンプレート:Mvar を法とした テンプレート:Mvar による除法(「余り付き除法」ではない)を、モジュラ逆数を掛けることとして定義することができる。

整数 テンプレート:Math の法 テンプレート:Math に関するモジュラ逆数 テンプレート:Mvar を求める。これは

31x(mod11)

を満たすものを計算するということだが、その意味は

3x1(mod11)

なる x を求めることに他ならない。/11 においてこの合同式の解 テンプレート:Mvarテンプレート:Math であることは

3×4=121(mod11)

から明らかであり、従って テンプレート:Math を法として テンプレート:Math の逆数は テンプレート:Math である。こうして /11 において逆数が求まったが、上記の合同式を満たす テンプレート:Mvarテンプレート:Math のみではない。残りの解は、得られた テンプレート:Mathテンプレート:Math の倍数を加えたものとして得られる。一般にこの例において テンプレート:Mvar となり得る整数は

4+11z,z

の形をしており、これは テンプレート:Math という集合を成す。

逆数計算

拡張ユークリッド互除法

テンプレート:Mvar に関する テンプレート:Mvar のモジュラ逆数は、拡張ユークリッド互除法を用いて計算することができる。互除法のアルゴリズムはベズーの等式

ax+by=gcd(a,b)

の解を求めるものである。ただし、テンプレート:Math が既知で、整数 テンプレート:Math および テンプレート:Math がこのアルゴリズムから求まる。

さて、モジュラ逆数は合同式

ax1(modm)

の解であったから、合同式の定義により テンプレート:Mathテンプレート:Mvarテンプレート:Math約数)、即ち適当な整数 テンプレート:Mvar を用いて

ax1=qm

となる。これを整理した

axqm=1

は、テンプレート:Mvarテンプレート:Mvar が既知で、逆数となる テンプレート:Mvar と捨てられる テンプレート:Mvar は未知であるから、拡張ユークリッド互除法が解く等式とまったく同じ形をしている。違う点といえば テンプレート:Math が「求まる」のではなくて「所与」であること程度で、それゆえ テンプレート:Mvarテンプレート:Mvar と互いに素でなければならず、さもなくば逆元 テンプレート:Mvar は存在しない。

このアルゴリズムの計算量は テンプレート:Math と仮定すると テンプレート:Math で、一般に冪乗の計算よりも効率的である。

オイラーの定理の利用

拡張ユークリッド互除法の代わりにオイラーの定理をモジュラ逆数の計算に利用することもできる[1]

オイラーの定理によれば、テンプレート:Mvarテンプレート:Mvar互いに素ならば、オイラー函数 テンプレート:Math に対して

aφ(m)1(modm)

が成り立つ。これは テンプレート:Mvar が法 テンプレート:Mvar に関する既約合同類群 (/m)× に入るための必要十分条件テンプレート:Mvarテンプレート:Mvar と互いに素であることから従う。従って、モジュラ逆数が直接的に

aφ(m)1a1(modm)

と求まる。これから、テンプレート:Mvar が素数である特別の場合には、テンプレート:Math ゆえ

a1am2(modm)

となる。この方法は一般には拡張ユークリッド互除法よりも遅いが、冪剰余の実装が使える場合には利用されることがある。この方法が不利な点としては、

応用

モジュラ逆数を見つけるという操作は、モジュラ計算を使ったアルゴリズムに多く応用されている。例えば、暗号技術ではモジュラ演算を用いることで、演算をより早く、より少ない記憶容量で行う部分と、演算がより困難な部分がある。[2]これら二つの性質は共に利点として利用できる。特にRSA暗号では、暗号化・復号には適切に選ばれたモジュラ逆数となる2つの数が使われる。2つのうち1つは公開し、暗号化するために使用される。もう1つは安全に保管しておき、復号に使用する。ここまでの処理は高速だが、1つ目の公開鍵から2つ目の秘密鍵を作るには非現実的な計算コストは避けられないと考えられており、このことがRSA暗号の仕組みを支えている。[3]

ほかの変わった例を一つ挙げる。ワードのサイズの、kの倍数で奇数である整数を要素に持つ(長い)リストがあったとする。ここで、これらの数をすべてkで割りたい。実装の一つとして次のようなものがある。

  1. k1(mod2w)拡張ユークリッド互除法を用いて計算する(wはワードのbit数、通常32)。gcd(2w,k)=1であるから、これは必ず存在する。
  2. リストの各要素にk1を掛け、下位 w bit に相当する部分を取り出す(つまり2wで割った余りをとる)。

除算をハードウェアで実装していない処理系にとって、除算は掛け算に比べて遅く、この方法によって格段に処理速度を速めることができる。手順1.には時間はかかるが1回で済むため、扱うリストのサイズが大きければそれだけ効率的である。

このほかモジュラ逆数を用いて、中国の剰余定理によって存在が保証されている合同方程式の解を求めることができる。

例えば、合同方程式

X4(mod5)
X4(mod7)
X6(mod11)

は 5, 7, 11 が互いに素であり解をもつ。解は

X=t1(7×11)×4+t2(5×11)×4+t3(5×7)×6

であり

t1=(7×11)1(mod5)=3
t2=(5×11)1(mod7)=6
t3=(5×7)1(mod11)=6

である。上記の合同方程式を満たしていることは容易にわかる。よって

X350439(mod385)

385 は 5, 7, 11 の最小公倍数である。

また、モジュラ逆数はテンプレート:仮リンクの定義で重要な位置を占めている。

関連項目

参考文献

テンプレート:Reflist