JP2001527673A

JP2001527673A - モントゴメリー乗算に基づくモジュラ乗算及び累乗の改善された装置と方法

Info

Publication number: JP2001527673A
Application number: JP54788398A
Authority: JP
Inventors: ハダッド，アイザック; アラジ，ベンジャミン; デビッドグレッセル，カーミ; ドロール，イタイ
Original assignee: フォートレスユーアンドティーリミティド
Priority date: 1997-05-04
Filing date: 1998-03-29
Publication date: 2001-12-25
Also published as: KR100684134B1; EP1008026B1; ATE372541T1; EP1008026A4; DE69838390T2; ES2293677T3; WO1998050851A1; AU6516498A; KR20010012222A; US6185596B1; JP2009230153A; DE69838390D1; EP1008026A1

Abstract

(57)【要約】本発明は、１つの桁上げ保存加算器（４１０）を有する１つのモジュラ乗算装置（４オペランド・モジュラ乗算器）を含む直並列論理演算装置（ＡＬＵ）を含むモジュラ乗算及び累乗の方法とシステムを開示する。

Description

【発明の詳細な説明】モントゴメリー乗算に基づくモジュラ乗算及び累乗の改善された装置と方法発明の分野本発明は、モジュラ乗算及び累乗と、直列整数除算のための装置と方法に関する。発明の背景大きな数についてモジュラ乗算及び累乗を行うコンパクトな小型電子装置は、その開示を引用によって本明細書の記載に援用する、本出願人の米国特許第５，５１３，１３３号で説明されている。その明細書で言及されているすべての出版物と本明細書に引用されている出版物の開示は、引用によって本明細書の記載に援用する。発明の概要本発明は、モジュラ乗算及び累乗と直列整数除算のための装置と方法を提供しようとする。従って、本発明の好適実施形態によれば、１つの桁上げ保存加算器を含む１つの乗算器と、好適には任意のビット長の被除数と任意のビット長の除数を受信し、商と剰余を計算するよう動作する直列除算装置とを含む直並列論理演算装置（ＡＬＵ）を含むモジュラ乗算及び累乗システムが提供される。さらに本発明の好適実施形態によれば、本システムは任意のビット長の少なくとも１対の整数入力を乗算するよう動作する。またさらに本発明の好適実施形態によれば、少なくとも１対の整数入力には２対の整数入力が含まれる。さらに追加して本発明の好適実施形態によれば、ＡＬＵは前もってゼロ強制モントゴメリー定数Ｊ₀を計算せずに、整数入力の積を生成し積の大きさを換算するよう動作する。また、本発明の別の好適実施形態によれば、任意のビット長の被除数と任意のビット長の除数を受信し商と剰余を計算するよう動作する直列除算装置を含む直列除算装置が提供される。さらに本発明の好適実施形態によれば、本装置は１対の整数入力を保存する１対のレジスタを含み、少なくとも１つがそれぞれのレジスタのビット長を越える１対の整数入力をインタリーブなしに乗算するよう動作する。また、本発明のまた別の好適実施形態によれば、１つだけの桁上げ保存累算器を有し１対の乗算を行いその結果を合計するよう動作する直並列乗算装置を含むモジュラ乗算及び累乗システムが提供される。さらに追加して、本発明のまた別の好適実施形態によれば、モジュラ乗算及び累乗方法が提供されるが、本方法には、１つの桁上げ保存加算器を有する１つのモジュラ乗算装置を含む直並列論理演算装置（ＡＬＵ）を提供するステップと、モジュラ乗算及び累乗を行うために直並列ＡＬＵを利用するステップとが含まれる。さらに、本発明のまた別の好適実施形態によれば、大きな整数の自然（非モジュラ）乗算の方法が提供されるが、本方法には、１つの桁上げ保存加算器を有する１つのモジュラ乗算装置を含む直並列論理演算装置（ＡＬＵ）を提供するステップと、大きな整数の自然（非モジュラ）乗算を行うために直並列ＡＬＵを利用するステップとが含まれる。さらに、本発明の好適実施形態によれば、利用ステップには、第１積を得るために任意のビット長の第１整数を任意のビット長の第２整数で乗算するステップと、第２積を得るために任意のビット長の第３整数を任意のビット長の第４整数で乗算するステップと、合計を得るために第１及び第２積を任意のビット長の第５整数と合計するステップとが含まれる。またさらに、本発明の好適実施形態によれば、利用ステップには、任意のビット長の被乗数、乗数及び法によるモジュラ乗算及び累乗を行うステップが含まれる。さらに追加して、本発明の好適実施形態によれば、本システムにはまた、１つの事前計算定数だけでモントゴメリー・モジュラ乗算を実行する二重被乗数事前計算システムが含まれる。さらに、本発明のまた別の好適実施形態によれば、利用ステップには、モントゴメリー乗算を行うステップが含まれ、そこには、乗数と被乗数を含む整数入力の積を生成するステップと、前もってモントゴメリー定数Ｊ₀を計算せずにモジュラ換算を実行するステップが含まれる。さらに、本発明の好適実施形態によれば、モントゴメリー定数Ｊ₀にはＮｍｏｄ２^kが含まれるが、ここでＮはモジュラ換算の法であり、ｋは被乗数のビット長である。またさらに、本発明の好適実施形態によれば、利用ステップには、一連のインタリーブド・モントゴメリー乗算演算を行うステップが含まれる。さらに追加して、本発明の好適実施形態によれば、各インタリーブド・モントゴメリー乗算演算は、少なくともｋの有意ゼロを伴う結果を出すために法を乗算演算の合同に合計しなければならない回数を前もって計算せずに行われる。またさらに、本発明の好適実施形態によれば、本システムにはまた、ｉ番目のインタリーブド・モントゴメリー乗算演算で生成される被乗数を収集及び直列合計し、それによって合計を生成し合計を（ｉ＋１）番目のモントゴメリー乗算演算に供給するよう動作するデータ・プロセッサが含まれる。さらに追加して、本発明の好適実施形態によれば、関数にはＮｍｏｄ２^kの乗法的逆元の加法的逆元が含まれる。さらに、本発明の好適実施形態によれば、本方法はまた、Ａ_i及びＢをゼロにリセットしＳ₀＝１を設定することでＪ₀を計算するステップを含む。本発明はまた、非常に大きな整数のモジュラ及び正規（整数の自然で負でない範囲）の乗算、除算、加算、減算及び累乗を行うコンパクトな小型電子論理演算装置ＡＬＵに関する。モントゴメリー法を使用するモジュラ乗算及び平方に言及する場合、装置の特定の部分はモジュラ演算コプロセッサと呼ばれ、ＭＡＰという頭字語が使用される。また、モントゴメリー乗算法はＭＭと呼ばれる。本発明はまた、大きな整数の演算処理に関する。こうした大きな数は、（負でない）整数の自然範囲、または素数及び複合素数の法のガロア域ＧＦ（ｐ）にある。すなわち、本発明は、公開鍵暗号認証及び暗号化プロトコルにとって不可欠な演算を行うのに適した大きな数のモジュラ乗算／累乗を実現できる装置に関するが、こうした演算は次第に大きなオペランドに対して行われるようになっており、現世代のモジュラ演算コプロセッサでは効率的に実行できず、ソフトウェアによる実現では確実に実行できないものである。本発明は任意の４ビットまたはそれ以上のプロセッサで動作し、現在のデジタル信号プロセッサを越える速度を達成している。本発明はまた、特に暗号化ソフトウェア指向システムで使用されることの多いインタリーブド・モントゴメリー・マルチプレシジョン・モジュラ乗算方法として知られる処理と、長いオペランドの整数に対する基本算術演算、すなわちＡ・Ｂ＋Ｃ・Ｄ＋Ｓ（ここでＡ、Ｂ、Ｃ、ＤまたはＳの大きさには理論的な制限はない）における固有の値との派生物の数値操作に関する大きなオペランドの整数演算のハードウェアによる実現に関する。さらに、本装置は特に、モジュラ乗算及び累乗を行うよう調整されている。基本装置は、モジュラ演算コプロセッサ（ＭＡＰ）に特に適しており、非常に大きな整数の除算を行う装置をも含むが、ここで除数は法レジスタＮと同じ長さのビット長を有し、被除数のビット長は２つの連結レジスタのビット長と同じ長さである。この装置は好適には、同じ次数の論理ゲートによって、半分未満の機械クロックサイクルで、米国特許第５、５１３，１３３号のすべての機能を行う。これは主として、同じ桁上げ保存累算器機構（従来の直並列乗算器の主要構成要素）を使用する２つの半分の大きさの乗算器の代わりに、１つだけの二重動作直並列乗算器が存在するためである。新しい論理演算装置ＡＬＵ、すなわちモジュラ演算コプロセッサＭＡＰは好適には本質的に新しい処理の前の別個の乗算処理を不要にする。この処理はまた２次モントゴメリー定数Ｊ₀を必要としたが、これも好適には不要になる。別言すれば、以前のモントゴメリー処置の２つの定数と、遭遇する遅延の代わりに、ここでは１つの定数だけが計算され、今では余分なＪ形乗算（後で説明する）によって発生する遅延は好適には除去される。さらに、ＣＰＵとこの周辺装置の間のデータ操作の制御向上によって、装置の自然レジスタ・サイズより長いオペランド上で行われるオペランドは、好適には少ない一時記憶メモリを使用して短縮された処理時間で行われる。モントゴメリーの方法論によってモジュラ乗算を行う３つの関連する方法が知られている。［Ｐ．Ｌ．Ｍｏｎｔｇｏｍｅｒｙ、「試行除算を行わないモジュラ乗算」、計算の数学、第４４巻、５１９〜５２１ページ、１９８５年］（以下「モントゴメリー」と呼ぶ）、［Ｓ．Ｒ．Ｄｕｓｓｅ及びＢ．Ｓ．ＫａｌｉｓｋｉＪｒ．、「モトローラＤＳＰ５６０００用暗号ライブラリ」、暗号化会議録 ’９０、Ｓｐｒｉｎｇｅｒ−Ｖｅｒｌａｇ、ベルリン、１９９０年］（以下「Ｄｕｓｓｅ」と呼ぶ）及びＭｉｙａｇｕｃｈｉへの米国特許第４，５１４，５９２号の方法、及びＥｖｅｎへの米国特許第５，１０１，４３１号の方法、及びＩｗａｍｕｒａへの米国特許第５，３２１，７５２号の方法、及びＡｒａｚｉへの米国特許第５，４４８，６３９号の方法、及びＧｒｅｓｓｅｌへの米国特許第５，５１３，１３３号の方法。好適なアーキテクチャは、ホスト・コンピュータのメモリにマップされた何らかの小型制御装置の設計に統合される機械のものであり、これはデータ供給機構との間で絶えずオペランドをスワップまたはやりとりする非常に長いコマンドについて、何らかの一般的な長さのモジュラ算術計算を可能にする制御装置と平行して動作するが、ここで操作のために必要なコプロセッサの揮発性メモリのサイズが最大オペランドの長さの３倍を越えることはまれである。このソリューションは、前の実現例で本質的に２つの乗算装置の役目を果たす１つの乗算装置だけを使用する。現在の一般的な技術を使用して、４×４．５× ０．２ｍｍの小型電子回路の上にメモリと共に小型制御装置を含むソリューション一式を集積することが可能である。本発明はまた、以下詳細に説明されるような、モントゴメリー及びＤｕｓｓｅによって発表された処理に関する計算、論理及びアーキテクチャの新しい特徴を有する従来のデジタル・プロセッサの周辺装置であることを企図するデジタル装置のアーキテクチャに向けられる。並列処理と独特のハードウェア・アーキテクチャか提供され、除算なしに、好適には従来の乗算／除算装置で行われるのと同じ数の演算でモジュラ累乗を行うが、その際従来の装置は各演算毎に乗算と除算の両方を行う。本発明の好適実施形態の特定の特徴は、基本的な大オペランド整数演算機能を考慮した、割込みのないリソースの効率的な使用による、無制限のオペランド長を考慮する装置によって行われる演算の同時性である。本発明の好適実施形態によって実現される利点は直列処理の同期シーケンスに起因するが、これはｎビット・オペランドに対して同時に（平行して）３つの乗算演算を達成するためにマージされ、（ｎ＋ｋ）有効クロックサイクルで１つの多重化ｋビット直並列乗算器を使用し、モントゴメリーによって規定されたように３つの乗算計算と同等の結果を達成する。次に使用されるオペランドの同期、実行中検出、実行中事前ロード及び同時追加によって、機械は決定論的に動作するが、その際すべての乗算及び累乗は所定の数のクロックサイクルで行われる。条件付き分岐はローカル検出及び補償装置によって置換され、それによって簡単な種類の制御機構の基礎を提供するが、これは、洗練された場合、通常一連の自励カスケード・カウンタを含む。ここで説明される基本演算は、スコットランドのＥａｓｔＫｉｌｂｒｉｄｅにあるＭｏｔｏｒｏｌａ（商品名ＳＣ−４９）と、フランスのＲｏｕｓｓｅｔにあるＳＧＳ −Ｔｈｏｍｓｏｎ（商品名ＳＴ１６−ＣＦ５４）の２社によって製造されるような、米国特許第５，５１３，１３３号で説明される装置を使用して決定論的時間で実行される。この機械は、演算の合計長さについてオペランドが装置にロードされ保存される際、大部分の演算の場合揮発性メモリに対する要求は特にわずかである。しかし、この機械は好適には機械が追加された先のＣＰＵを活用し、簡単なロードとアンロード、機械へのコマンドの順序付けを実行する一方で、機械は大きな数の演算を行う。累乗処理時間は、事実上それを制御するＣＰＵと無関係である。実際には、機械を何らかのＣＰＵに追加する際アーキテクチャの変更は不必要である。ハードウェア装置は自己完結型であり、任意のＣＰＵバスに追加できる。モジュラ乗算及び累乗処理を加速する装置が好適には提供されるが、これには必要な定数を事前計算する手段が含まれる。ここで説明される本発明の好適実施形態は、携帯型スマートカードの公開鍵暗号適用業務のためのモジュラ数学的演算子を提供するが、このスマートカードは通常、一般的な磁器ストライプ・クレジット及び銀行カードと形状と寸法が同一である。同様のスマートカード（米国特許第５，５１３，１３３号）が、新世代公開鍵暗号装置で、コンピュータ、データベース及び重要な施設へのアクセスの管理、商業用、軍事用及び家庭用のトランザクションにおけるデータフローの調整と保護、スクランブルされた有料テレビ番組のスクランブル解除等のために使用されている。こうした装置はまた、コンピュータ及びファックス端末、ドアロック、自動販売機等にも組み込まれることが認識されるだろう。説明されるハードウェアはρ演算子を新しい方法で適用することでモジュラ乗算及び累乗を実行する。さらに、それを互いに等しい被乗数と乗数に適用することで、同じ方法で平方が実行される。モジュラ累乗にはモジュラ乗算と平方の連続が含まれるので、上記の乗算、平方及び累乗法の反復される、適切に結合され方向付けされた応用を含む方法によって実行される。ＡＬＵの好適実施形態の動作を説明する際、有効クロックサイクルでの同期が説明されるが、これは実際のクロックサイクルとは異なって装置が算術演算を行っているサイクルのことであり、実際のクロックサイクルには、新しい演算段階の準備のために、ＡＬＵが停止し、乗算器、フリップフロップ及び他の装置の設定が変更されるアイドル・サイクルが含まれている。好適実施形態では、ＣＰＵの揮発性ＲＡＭまたはＳ_Aレジスタ１３０の何れかに保存される被乗数Ａ、７０と８０の連結であるＢレジスタ１０００の乗数Ｂ、及び２００と２１０の連結であるＮレジスタ１００５の法Ｎが各ｋビットのｍ文字を含み、被乗数と乗数は一般に法より大きくない、（平方と正規乗算に関連する）モントゴメリー・モジュラ乗算を実行する方法は、１）乗数Ｂと法Ｎをｎビット長のそれぞれのレジスタにロードするステップであって、ここでｎ＝ｍ・ｋであるステップと、｛正規フィールドで正、自然、整数を乗算し、Ｎは第２乗数である｝｛ｎがＢ、Ｎ及びＳレジスタより長い場合、値は通常反復の過程でこれらのレジスタにロード及びアンロードされ、機械が事実上あらゆる長さの法を操作できるようにする｝２）レジスタＳ_Bの出力をゼロに設定するステップであって、最初の反復についてＳ＊ｄフラッシュ（２５０）＝１であるステップと、３）無縁借り及び桁上げフラグをリセットするステップ（特許で規定されていない制御）と、４）ｍ回の反復を実行するステップであって、各反復が次の演算を含むステップとを含む。（０≦ｉ≦ｍ−１）ａ）被乗数Ａの次の文字Ａ_i-1を揮発性記憶装置からＡｉロード・バッファ２９０に転送する。ｂ）Ａｉロード・バッファのコンテントを回転する一方で、Ｎ₀（ＮのＬＳｋビット）を伴うＣｉロード・バッファを同時直列ロードすることで、直列加算器ＦＡ１、３３０によってＮ₀を伴うＡｊロード・バッファのコンテントを直列加算し、それによって合計Ｎ₀＋Ａ_i-1を伴うＡｉ＋Ｃｉロード・バッファを直列ロードする。事前ロード段階はここで終了する。この段階は通常、ＭＡＰが前の乗算反復を行っていた間に行われる。処理ａ）及びｂ）は同時に実行することができ、その際Ａ_i-1文字は対応するレジスタにロードされ、ＡｉストリームはＮ₀レジスタの回転、Ｒ２、３２０のロードに同期する。同時に、ＡｉストリームとＮ₀ストリームは合計され、Ｒ３レジスタ、３４０にロードされる。Ｂレジスタからの量の平方が実行されるが、そこでは初期化、ステップａ）及びｂ）で、Ｂ₀レジスタが回転される際に、Ｎ₀レジスタと同時に、Ｂｄの最初のｋビットがＲ１に挿入される。次のｋビットＢ_i文字列が、ＡＬＵに直列供給される際にＲ１レジスタに事前ロードされる。ｃ）機械が停止する。バッファＲ１、Ｒ２及びＲ３のオペランドはラッチＬ１、３６０、Ｌ２、３７０及びＬ３、３８０にラッチされる。Ｌ０−“０”ラッチは疑似ラッチであり、これは単に、３９０、乗算器の入力または出力の各ＡＮＤゲートに入るリテラル・コマンド信号である。ｄ）次のｋ回の有効クロックサイクルについて、ｉ）各有効クロックサイクルについて、Ｙ０ＳＥＮＳＥはＹ₀の次のビットを予想し、このビットをＭ３乗算器を通じてＣｉロード・バッファにロードする一方で、ＡｉビットをＲ１レジスタからシフトアウトすると同時にＣｉロード・バッファにＹ₀のｋビットをロードし、Ｒ１の出力にＹ₀を加算し、この値をＲ３バッファにロードする。ｉｉ）Ｍ＿Ｋ乗算器３９０によって４つのラッチＬ０、Ｌ１、Ｌ２またはＬ３の１つから望ましい値を論理的に選択することで、同時にＮ₀（Ｌ２、Ｃｉラッチ）を入力Ｙ₀ビットで乗算し、ＡｉをＢｄの次の入力ビットで乗算して、２つの結果を加算する。Ｙ₀ビットとＢビットのどちらも１でない場合、すべてのゼロ値はＣＳＡに多重化され、Ｎ₀ビットだけが１である場合、Ｎ₀だけがＣＳＡに多重化／加算され、Ｂビットだけが１である場合、Ａ_i-1がＣＳＡに加算され、ＢビットとＮ₀ビットの両方が１である場合、Ａ_i-1＋Ｎ₀がＣＳＡに加算される。ｉｉｉ）この合計は桁上げ保存ｋ＋１ビット累算器をビット毎に直列に出るので（Ｘストリーム）、直列加算器ＦＡ２、４６０を通じてＳｄの次の関連ビットをこの合計に加算する。ＭＭでは、Ｚストリームの最初のｋビットはゼロである。この第１段階で、Ｙ₀ ・Ｎ₀＋Ａ_i-1・Ｂ₀＋Ｓ₀の結果が計算され、ＬＳｋ全ゼロ・ビットがＺ^*ｏｕｔストリームに現れ、乗算装置のＭＳｋ＋１ビットがＣＳＡ桁上げ保存累算器で保存される。Ｒ１、Ｒ２及びＲ３事前ロード・バッファはそれぞれ値Ａ_i-1、Ｙ₀及びＹ₀＋Ａ_i-1を保持する。ｅ）最後の有効な（ｍ＋１）・ｋ番目のクロックサイクルで、機械は停止し、バッファＲ２及びＲ３はＬ２及びＬ３にラッチされる。Ｌ１の値は変化しない。次のｋ・（ｍ−１）有効クロックサイクルの初期および継続条件は次の通りである。乗数は、Ｂのｋ番目のビットから開始されるＢのビット・ストリームであり、Ｎからの残りのビット・ストリームもＮのｋ番目のビットから開始される。Ｌ１、Ｌ２及びＬ３の被乗数はＡ_i-1、Ｙ₀及びＹ₀＋Ａ_i-1であり、開始時にＣＳ加算器はｄ）で説明された値を含み、Ｓストリームは次のｋ・（ｍ−１）ビットをＦＡ２全加算器に供給する。次のｋ・ｍ有効クロックサイクルの間に、ユニット４７０のＮｄ、遅延ｋクロックサイクルが直列減算器４８０で、Ｚストリームから減算され、ＢまたはＳレジスタに進むその結果がＮより大きいかまたはそれ以上かを検知する（Ｚ／２^k ｍｏｄ２^k・m）。直列減算器４６０の検知結果に関わらず、｛（ｍ＋１）・ｋ｝番目の有効クロックサイクルでＣＳＡのＳＯ₁フリップフロップが１ならば、合計結果は確実にＮより大きく、Ｎは結果から減算され、その結果、部分的または最終的に、そのレジスタを出る。ｆ）次のｋ・（ｍ−１）有効クロックサイクルについて、Ｎ０レジスタ２１０は、入力Ａｉビットと同期してか、または他の適切なタイミングで回転し、ａ）及びｂ）で説明されたように、次の反復のためにＲ１、Ｒ２及びＲ３をロードする。このｋ・（ｍ−１）有効クロックサイクルで、Ｎの残りのＭＳビットがここでＹ₀を乗算し、残りのＭＳＢビットはＡ_i-1を乗算し続ける。ＮビットとＢビットがどちらも１でない場合、全ゼロ値がＣＳＡに多重化される。Ｎビットだけが１である場合、Ｙ₀だけがＣＳＡに多重化／加算される。Ｂビットだけが１である場合、Ａ_i-1がＣＳＡに加算される。ＢビットとＹ₀ビットの両方が１である場合、Ａ_i-1＋Ｙ₀がＣＳＡに加算される。同時に、ＣＳＡからの直列出力が、Ｚストリームを出力するＦＡ２加算器、ユニット４６０を通じて次のｋ・（ｍ−１）Ｓビットに加算される。Ｚ出力ストリームの関連部分は、Ｚの最初の非ゼロｋ・（ｍ−１）ビットである。Ｚストリームは、最初のｍ−１回の反復についてＳ_Bレジスタにスイッチされ、最後の反復について定義されたように、Ｓ_BまたはＢレジスタにスイッチされる。最後の反復で、Ｚストリームは、ＬＳｋゼロ・ビットを無視し、最後のＢ^*ストリームである。このストリームはＢレジスタに送られ、次の乗算及び平方で使用される際必要ならばＮだけ減らされる。最後の反復で、Ｎｄ、遅延ｋクロックサイクルが直列減算器によってＺストリームから減算され、Ｂに進むその結果がＮより大きいかまたはそれと等しいかを検知する。この段の終了時に、Ｎ、Ｂ及びＳ_BレジスタからのすべてのビットはＡＬＵに供給されており、結果の最終ｋ＋１ビットはＣＳＡにあり、フラッシュアウトの用意が整っている。ｇ）装置が停止する。Ｓフラッシュ２５０、Ｂフラッシュ２４０及びＮフラッシュ２６０が出力ゼロ文字列に設定され、次の段階で最後のｋ＋１最上位ビットが確実にＣＳＡからフラッシュアウトされるようにする（正規乗算では、Ｍ７ＭＵＸ４５０はＳの前の反復から最終桁上げを受け入れるよう設定される）。Ｓはｍ・ｋ＋１有効ビットを有するが、Ｓレジスタはこのデータを受信するｍｋセルだけを有する。この最終ビットは本質的にオーバフロー機構で保存される。ｅで説明されたように、４７０のＮｄ、遅延ｋクロックサイクルはＺストリームから減算され、Ｘからの有効出力と同期されて、ＢまたはＳレジスタに進むその結果がＮより大きいかまたはそれ以上かを判定する検知機構の微調整を提供する。４８０及び４９０は直列比較装置を備え、ここでモジュラ計算の場合の最後の借りコマンド・ビットだけと整数の自然範囲の正規乗算の場合の（ｋ・ｍ＋１）番目のビットが保存される。このオーバフロー／借りコマンドはｍ・ｋ番目の有効クロックサイクルで検出される。ｈ）装置は別のｋサイクルでクロックされ、ＣＳＡから完全にフラッシュアウトされるが、一方別のｋビットは定義された出力レジスタに出るＺである。次の出口ストリームでＮの減算を実行する関連フリップフロップ・コマンド直列減算器９０または５００への命令は、（Ｚ／２^k−Ｎ）≧Ｎ（Ｚはｍ・ｋ番目のＭＳビットを含む）が、次の信号、すなわち、ＣＳＡの第２最下位セルからのデータ出力ビットであるＳＯ₁ビットが１であるか、または、Ｘ＋Ｓ加算器４６０の内部キャリーアウトであるＣＯ₂ビットが１であるか、または、４８０センス減算器からの借りビットが設定されていない、によって検知されるならば、反復の最後の有効な（ｍ＋１）ｋ番目のクロックサイクルに設定される。この機構はＭｏｔｏｒｏｌａ及びＳＧＳ−Ｔｈｏｍｓｏｎの２社によって製造されるものとして米国特許第５，５１３，１３３号に現れる。自然数範囲の乗算の場合、ｍ・ｋ番目のＭＳビットが１であるならばスーパースカラ乗算器で発生し、米国特許第５、５１３，１３３号の機構では発生しないオーバフローを検出することが好適である。その後このオーバフローを次の反復で使用して、Ｓ（一時的結果）ストリームでＭＳ１を挿入することができる。ｊ）これは最後の反復か。いいえ、ならばｃ）に戻る。はい、ならばｍ）に進む。ｋ）結果の正確な値をここでＢまたはＳレジスタから出すことができる。Ｙ０ビットは、５つの決定論的量からＹ０Ｓ−Ｙ０ＳＥＮＳＥユニット４３０で次の方法で予想される。ｉＡ_i−Ｌ１のＬＳビットとＢｄストリームの次のビット、Ａ₀・Ｂ_d ｉｉ桁上げ保存累算器からのＬＳキャリーアウト・ビット、ＣＯ₀ ｉｉｉＣＳＡの第２ＬＳセルからのＳ_outビット、ＳＯ₁ ｉｖＳストリームからの次のビット、Ｓ_d ｖ４６０、全加算器からのキャリーアウト・ビット、ＣＯ₂ これらの５つの値は一緒に排他的論理和演算され、次のＹ₀ビット、Ｙ_0iを生成する。Ｙ_0iビットが１である場合、（必要な２の累乗を乗算された）同じランクの別のＮが通常加算されるが、そうでない場合、Ｎ（法）は通常加算されない。数の正規範囲での長い自然整数の乗算。この装置は、正規整数の乗算及び合計を効率的に行うのに適している。こうしたオペランドがすべてｋビット長より長くない場合、処理は好適にはインタリーブなしで実行されるが、そこでは２ｋ＋１ビットのＺストリームが最終記憶装置に送られる。ｋビットより長い整数の場合、処理は所定のインタリーブド・モジュラ演算処理と同様であるが、結果がここでは潜在的に最長オペランドの長さの２倍より１ビット長いことが異なっている。さらに、本発明の装置は、説明された装置で利用可能なリソースを使用して、被乗数Ａ（好適にはＲ１−Ａｉレジスタのセグメントにロードされる）×Ａの乗数Ｂ（好適には前もって指定されたＢレジスタにロードされる）＋第２乗数Ｎ（好適にはＮレジスタにロードされる） ×オペランドＣ（Ｒ２レジスタにロードされる）＋装置に入るビット・ストリームＳといった２つの別個の乗算を、最初の反復で、Ｓｄ信号ラインだけから、好適にはＳ_Aレジスタから同時に行うことができる。Ｙ０ＳＥＮＳＥ装置は使用されない。被乗数は、反復開始前にＲ３レジスタに合計される。反復開始時に、レジスタＲ１、Ｒ２及びＲ３は、反復終了までラッチＬ１、Ｌ２及びＬ３にコピーされる。一方、反復のｍｋ＋ｋ＋１有効クロックサイクル中、Ａ及びＣの次のセグメントが再び事前ロード及び合計されて、次の反復を準備する。各反復で、Ｚストリームに関する結果の最初のＬＳｋビットは、ここでは（ＭＭの場合のように、定義ゼロによってではなく）別個の記憶装置に向けられ、空にされて結果のＬＳ部分、ここでも適切にはＳ_Aレジスタを累算する。最上位ｍｋ＋１ビットは、次の反復の一時量Ｓ_Bを含む。最終段階では、ｇ、ｉ及びｊと同様に、ＣＳＡは累算値をフラッシュアウトする。乗数レジスタより長い数に関するＬＳ部分は、正規データアウト・レジスタとアンローダ（それぞれユニット６０及び３０）によって出ることができる。結果のＭＳ２ｍ番目のビットは、ＯＶＥＲＦＬＯＷ信号ラインを通じてＦＡ２、ユニット４６０のＬＡＳＴＣＡＲＲＹビットから読み出される。図面の簡単な説明本発明は、図面と共に行われる以下の詳細な説明から理解され認識されるだろう。図１Ａ〜図１Ｂは、一緒に取り上げられ、本発明の１つの実施形態によって構成され動作する直並列スーパースカラ論理演算装置（ＡＬＵ）の単純化された構成図を形成する。図２は、特に非常に長い数の直列整数除算適用業務について単独でも有益である、図１Ａの直列整数除算装置の好適実現例の単純化された構成図である。図３は、図１Ａ〜図１Ｂの直並列論理演算装置を含む、スマートカードまたは端末用の公開鍵暗号計算の単純化された構成図である。図４は、除数の有効ビット長が被除数の有効ビット長の半分である例について、図２の除算装置を使用する除数による被除数の除算の演算の段階を示す表である。図５は、除数の有効ビット長が被除数の有効ビット長の半分より少ない例について、図２の除算装置を使用する除数による被除数の除算の演算の段階を示す表である。好適実施形態の説明図１Ａ〜図１Ｂは、一緒に取り上げられ、本発明の好適実施形態によって構成され動作する直並列論理演算装置（ＡＬＵ）の単純化された構成図を形成する。図１Ａ〜図１Ｂの装置には、好適には次の構成要素が含まれる。単一多重装置−制御されたスイッチング素子であり、多数の信号入力から１つの信号またはビット・ストリームを選択し、この選択された信号を１つの出力に送る。多重装置はＭ１〜Ｍ１３の記号が付けられ、大きな要素の組み込み部分である。Ｍ＿Ｋ多重装置３９０はｋ＋１の単一多重装置のアレーであり、４つのｋまたはｋ＋１入力のどれをＣＳＡ４１０に加算するかを選択する。Ｂ（１０００）、Ｓ_A（１３０）、Ｓ_B（１８０）及びＮ（１００５）は好適実施形態における４つの主直列主レジスタである。Ｓ_Aは概念上及び実際上冗長であるが、非常に長い数の計算をかなり加速し、特に法の長さが２・ｋ・ｍビット長である場合揮発性メモリ・リソースを節約し、かつ長い除算計算を単純化する。直列加算器と直列減算器は２つの直列入力と１つの直列出力を有する論理要素であり、ビットの２つの長い文字列に対して合計または減算を行う。構成要素９０及び４８０は減算器であり、３３０及び４６０は直列加算器である。入力から出力への伝播時間は非常に小さい。直列減算器９０は、Ｂ^*がＮより大きいかまたはそれと等しいである場合、Ｂ^*をＢに換算する。直列減算器４８０が、比較器構成要素の一部として使用され、Ｂ^*がＮより大きいかまたはそれと等しいかを検出する。全加算器３３０は、ロード・バッファ３４０に、２９０及び３２０ロード・バッファの値の合計に等しい値を供給する２つのビット・ストリームを加算する。高速ローダ及びアンローダ（それぞれ１０及び２０、及び３０及び４０）は、ＣＰＵ制御装置からのデータフローを加速する装置である。これは、好適実施形態ではＤＭＡまたは他のハードウェア・アクセラレータから構成される。２０及び４０はデータ語を逆転するためのもので、図２の除算入力及び出力を調和させるために必要である。データ入力５０は並列入力直列出力装置であるが、本ＡＬＵ装置は直列供給シストリック・プロセッサであるので、データは並列に入力され、直列に出力される。データ出力６０は直列入力並列出力装置であり、コプロセッサから結果を出力する。商生成器は、除算機構の各反復で商ビットを生成する図２のその部分である。Ｂｄ、２４０、Ｓ^*ｄ、２５０及びＮｄ、２６０に対するフラッシュ信号は、最後のｋ＋１ビットが確実にＣＳＡからフラッシュアウトされるようにするために作成されているが、代替物は複雑化ｋ＋１ビット並列出力要素であり、累算器のＭＳｋ＋１ビットを検索する。ロード・バッファＲ１、２９０、Ｒ２、３２０及びＲ３、３４０は直列入力並列出力シフトレジスタであり、３つの可能なゼロより多い被乗数の組み合わせを受信することができる。ラッチＬ１、３６０、Ｌ２、３７０及びＬ３、３８０はロード・バッファからの出力を受信するよう作成されるので、ロード・バッファは、データが好適にもＬ２、Ｌ２及びＬ３にラッチされる前に、このデータの次の段階の処理が一時的にできるようになる。Ｙ０センス、４３０は、ＬＳゼロのｋビット文字列がモントゴメリー乗算及び平方でＺで出るようにするため、法が累算される回数を決定する論理装置である。１ビット遅延装置１００、２２０及び２３０はそれぞれのデータ・ストリームに挿入され、図１Ａのデータ準備装置と図１Ｂのデータ処理装置の間の同期の問題を調停する。ｋビット遅延、シフトレジスタ、４７０は、Ｚ／２^kがＮより大きいかまたはそれと等しければ、Ｚ／２^kとＮの比較が同期して確実になされるようにする。桁上げ保存累算器は直並列乗算器とほとんど同一であるが、普通直並列乗算器の入力にラッチされる１つの値の代わりに、３つの異なったゼロより大きい値を合計することができる点が異なっている。挿入最終桁上げ、４４０は、Ｓストリームのｍｋ＋１番目のビットを挿入するために使用されるが、Ｓレジスタはｍｋビット長しかない。借り／オーバフロー検出、４９０は、結果が（Ｎからの）法より大きいかまたはそれと等しいか、またはｍｋ番目のビットが１かの何れかを検出することができる。制御機構は示されないが、好適にはシストリック・データフローのためのスイッチセットを有する１組のカスケード計数装置であるものと理解される。数の素数及び複合素数範囲でのモジュラ乗算について、ここでＡ及びＢを被乗数及び乗数として定義し、Ｎを法として定義するが、Ｎは普通ＡまたはＢより大きい。Ｎはまた、法の値が保存されるレジスタをも示す。Ｎは、場合によっては、Ａより小さいことがある。Ａ、Ｂ及びＮはｍ・ｋ＝ｎビット長オペランドとして定義される。各ｋビットのグループは文字と呼ばれ、グループの大きさは乗算装置の大きさによって定義される。次に、Ａ、Ｂ及びＮは各々ｍ文字長である。以下のステップ毎に手順を追う説明を容易にするために、Ａ、Ｂ及びＮは５１２ビット長（ｎ＝５１２）であると想定し、ｋは６４ビット長であると想定するが、これはこうした乗算器の現在の費用効果的な長さと、簡単なＣＰＵのデータ操作速度のためである。またｍ＝８はオペランド中の文字数であり、また５１２ビット・オペランドの平方または乗算ループにおける反復の数でもある。すべてのオペランドは正の整数である。より一般的には、Ａ、Ｂ、Ｎ、ｎ、ｋ及びｍは任意の適切な値を想定する。非モジュラ関数では、Ｎ及びＳレジスタは、他の演算オペランドの一時記憶装置として使用できる。記号≡が使用され、例えば１６≡２ｍｏｄ７というように、モジュラ数の合同を示すが、ここで、１６を７で割った時剰余が２であるので１６は２モジュロ７に合同であるという。ＹｍｏｄＮ≡ＸｍｏｄＮと書く時、ＹとＸはどちらもＮより大きいが、正のＸ及びＹについて、剰余は同一である。また、負の整数Ｙの合同はＹ＋ｕＮであるが、ここでＮは法であり、Ｙの合同がＮより小さい場合、ｕは正の結果を与える最小の整数である。記号￥が使用され、より制限された意味での合同を示す。ここで説明される処理中、値は望ましい値、または望ましい値プラス法の何れかである。例えばＸ￥２ｍｏｄ７という場合、Ｘは２または９に等しい。Ｘは２ｍｏｄ７と制限的合同であるという。Ｘ＝ＡｍｏｄＮと書く場合、ＸはＡをＮで割った剰余として定義される。例えば、３＝４５ｍｏｄ７である。数論では、モジュラ乗算逆が基本概念である。例えば、Ｘのモジュラ乗算逆はＸ^-1と書かれるが、これはＸＸ^-1ｍｏｄＮ＝１によって定義される。Ｘ＝３でＮ＝１３の場合Ｘ^-1＝９である。すなわち、３・９を１３で割った剰余は１である。頭字語ＭＳ及びＬＳは、デジタル用語で慣用されているように、ビット、文字及びフルオペランド値に関連する場合、最上位及び最下位を示すために使用される。本明細書を通じて、Ｎは値Ｎと、Ｎを収容するシフトレジスタの名称の両方を示す。数値の上に付いた星印は、数値が現状では潜在的に不完全または変化する可能性があることを示す。Ａは累乗される数の値であり、ｎはＮオペランドのビット長である。初期化の後、ＡがＡ^*に「モントゴメリー正規化」される場合（Ａ^*＝２ⁿＡ、後で説明される）、Ａ^*とＮは累乗の中間ステップを通じて一定の値である。最初の反復中、累乗の初期化後、ＢはＡ^*に等しい。Ｂはまた、最終的に累乗の望ましい結果に等しい累算値が存在するレジスタの名前でもある。Ｓ^* は一時値を示し、Ｓ、Ｓ_A及びＳ_Bはまた、Ｓの単一ＭＳビットだけが保存されるレジスタ（単数または複数）を示す（このＭＳビットと連結されたＳ^*はＳと同一である）。Ｓ（ｉ−１）はｉ番目の反復の開始時のＳの値を示し、Ｓ₀は、Ｓ（ｉ）番目の値のＬＳ文字を示す。Ｐ範囲での乗算として処理ρ（Ａ・Ｂ）Ｎ（後で定義される）または、簡単に乗算演算が言及される。二重動作直並列乗算器を構成する基礎として標準構造の直並列乗算器を使用したため、我々は桁上げ保存累算に基づく乗算器の合計部分を（桁上げ先見加算器またはリップル加算器、前者はかなり複雑になり、後者は非常に低速である、と）区別して、それを桁上げ保存加算器または累算器と呼び、事前ロード機構及び多重装置及びラッチと別個に扱うが、これによってＡ×ＢとＣ×Ｄの乗算を同時に行い２つの結果を加算すること、すなわちＡ×Ｂ＋Ｃ×Ｄが可能になり、この累算器を非常に強力なエンジンに変換する。付加的論理がこの乗算器に追加され、モジュラ減算と直列合計のために必要な予想センス操作を提供するが、これは非常に大きな数の強力なモジュラ演算と普通の整数演算を提供する。モントゴメリー・モジュラ乗算モジュラ乗算を計算する従来のアプローチでは、Ａ・ＢｍｏｄＮ、すなわち積Ａ・Ｂの剰余が除算処理によって計算される。大きなオペランドの従来の除算を実現するのは、直並列乗算を行うより困難である。モントゴメリーのモジュラ換算法を使用すると、除算は本質的に２つの事前計算定数を使用する乗算によって置き換えられる。ここで示される手順では、法の関数である１つの事前計算定数だけが存在する。この定数はこのＡＬＵ装置を使用して計算される。この装置で使用されるモントゴメリー処理の単純化した表示がここで提供され、それに完全な好適説明が続く。奇数（ＬＳビット１）、例えば１０１０００１（＝８１₁₀）がある場合、それに奇数を補償する別のフィクシング、例えば１１１１（＝１５₁₀）を加算することでこの奇数を偶数（１つのＬＳビットがゼロ）に変換することができる、すなわち、１１１１＋１０１０００１＝１１０００００（９６₁₀）である。この特定の場合では、前もって全体の文字列８１を知っていたため５つのＬＳゼロを生じる数を発見し、８１に加算する二進数を容易に発見し、必要な数のＬＳゼロを有する新しい二進数を生成した。このフィクシング数は奇数であり、さもなければ結果のＬＳビットを累加する作用はない。行うべき処理がクロックされた直並列桁上げ保存処理で、ＬＳゼロの連続数を有することが望ましく、各クロックサイクルで次のビットを固定しなければならないだけだとすれば、次のビットが潜在的に１である場合フィクシングを加算し、潜在的ビットが０であった場合加算しないことで十分である。しかし、ビット間オーバフロー（二重キャリア）を発生しないために、このフィクシングは好適には前もって被乗数と合計され、関連乗数ビットが１であり、Ｙセンスも１を検出する時累算器に加算される。ここでは、モジュラ演算の場合のように、法で割った値の剰余だけに関心があり、法の任意の倍数をその値に加算しても、同じ剰余を有することを知っている。これは任意の整数にＹＮ＝Σ_yi２ⁱＮを加算しても同じ剰余を有するということを意味するが、ここでＹは必要なＬＳゼロを生成するために加算する倍数である。説明されたように、加算することのできる法は奇数だけである。偶数の法が、ｉが偶数中のＬＳゼロの数である時生じる奇数の２ⁱ倍として定義される方法が存在する。モントゴメリー・インタリーブド変化によって解決される問題は、数について有する制限された場所を低減し、乗算器の費用効果的な大きさを得ることを目指している。これは公開鍵暗号関数を行う場合特に有益であるが、そこでは絶えず大きな整数、たとえばｎ＝１０２４ビットと別の大きな整数の乗算が行われ、普通２倍の長さである２０４８ビットの整数を生成する。乗算（または平方）の処理中、十分な倍数のＮｓ（法）がＡ・Ｂ＝ＸまたはＡ・Ｂ＋Ｓ＝Ｘに加算されるので、ｎ個のＬＳゼロと、最大ｎ＋１個のＭＳビットを有する数Ｚが得られる。これらのゼロを無視することで望ましい結果を２ⁿで除算したことを思い出すならば、ＬＳｎビットを無視してこの数を使用し先に進むことができる。ＬＳｎビットが無視され、最上位ｎ（またはｎ＋１）ビットだけを使用する場合、結果は２^-n、すなわち２ⁿのモジュラ逆によって乗算された。次にこの結果に再び２ⁿｍｏｄＮ（または２ⁿ）を乗算するならば、（同じ剰余を有する）望ましい結果と合同な値をＡ・Ｂ＋ＳｍｏｄＮとして得ることができる。見られるように、ＭＭを使用すると、結果は好適には２²ⁿによって乗算され、ＭＭによって導入される２^-n寄生要素を克服する。例Ａ・Ｂ＋ＳｍｏｄＮ＝（１２・１１＋１０）ｍｏｄ１３＝（１１００・１０１１＋１０１０）₂ｍｏｄ１０１１₂ ｎＬＳビットにフィクシングが必要な場合常に２ⁱＮを加算する。Ｂ１０１１ ×Ａ１１００加算Ｓ１０１０加算Ａ（０）Ｂ００００ −−−−ＬＳビットの合計＝０Ｎを加算しない加算２⁰（Ｎ・０）００００合計０１０１ →０ＬＳビットは桁上げ保存加算器を出る加算Ａ（１）Ｂ００００ −−−− ＬＳビットの合計＝０ − Ｎを加算する加算２¹（Ｎ・１）１１０１合計１００１ →０ＬＳビットはＣＳ加算器を出る加算Ａ（２）Ｂ１０１１ −−−− ＬＳビットの合計＝０Ｎを加算しない加算２²（Ｎ・０）００００合計１０１０ →０ＬＳビットはＣＳ加算器を出る加算Ａ（３）Ｂ１０１１ −−−− 合計ＬＳビット＝１Ｎを加算する加算２³（Ｎ・１）１１０１合計１０００１ →０ＬＳビットはＣＳ加算器を出る結果は１０００１００００₂ｍｏｄ１３＝１７・２⁴ｍｏｄ１３である。１７は１３より大きいので１３を減算し、結果は次のようになる。１７・２⁴＝４・２⁴ｍｏｄ１３正式には２^-n（ＡＢ＋Ｓ）ｍｏｄＮ＝９（１２・１１＋１０）ｍｏｄ１３≡４となる。モントゴメリー演算ではＭＳ非ゼロ結果（４）だけを利用し、真の結果は２ｎで除算され、ｎゼロはＭＭ結果にされることを有効に記憶する。結果に２⁴ｍｏｄ１３≡３を有効に乗算した（８＋２）・１３＝１０・１３を加算した。実際には、余分のゼロを使用したので、Ａ・Ｂ＋Ｙ・Ｎ＋Ｓ−（１２・１１＋１０・１３＋１０）を１つの処理で行ったと言うことができるが、これは好適実施形態で可能であることが説明される。チェック−（１２・１１＋１０）ｍｏｄ１３＝１２，４・３＝１２要約すると、モントゴメリー乗算の結果は２^-nで乗算された望ましい結果である。前の結果から同じ乗算法を使用した望ましい結果を検索するために、前の結果を２²ⁿによってモントゴメリー乗算しなければならないが、各ＭＭは２^-nの寄生要素を残すので２²ⁿはＨと呼ばれる。モントゴメリー乗算関数ρ（Ａ・Ｂ）ＮはＰ範囲（上記の例で４を導出した）へのＡ・Ｂ積の乗算モジュロＮを行う。Ｐ範囲から正規モジュラ範囲への検索は、事前計算定数Ｈを使用してρ（Ａ・Ｂ）Ｎの結果にＰを規定することによって行われる。ここではＰ≡ρ（Ａ・Ｂ）Ｎならば、ρ（Ｐ・Ｈ）Ｎ≡Ａ・ＢｍｏｄＮということになり、２つのＰ範囲乗算で正規モジュラ乗算を行う。モントゴメリー・モジュラ減算は、ｎまたはｎ＋１ビット長であるオペランドに対して一連の乗算、加算及び減算を行うことによって、ｎ及び２ｎビット長であるオペランドに対する一連の乗算及び除算演算を回避する。処理全体はＮより小さいかまたはそれと等しい結果を生じる。所与のＡ、Ｂ及び奇数Ｎについて、Ａ・Ｂ＋Ｑ・Ｎが、ｎＬＳビットがゼロであるような数、またはＰ・２ⁿ＝Ａ・Ｂ＋Ｑ・ＮであるようなＱが常に存在する。これは、ｎＬＳビットがゼロである、２ｎビット長の式を有することを意味する。ここで、Ｉ・２ｎ≡１ｍｏｄＮ（Ｉはすべての奇数Ｎについて存在する）とする。前の等式の両側にＩを乗算することで次の合同が生じる。等式の左辺から、Ｐ・Ｉ・２ⁿ≡ＰｍｏｄＮ（Ｉ・２ⁿ≡１ｍｏｄＮを想起されたい）また右辺から、Ａ・Ｂ・Ｉ＋Ｑ・Ｎ・Ｉ≡ＡＢ・ＩｍｏｄＮ（Ｑ・Ｎ・Ｉ≡０ｍｏｄＮを想起されたい）従って、Ｐ≡Ａ・Ｂ・ＩｍｏｄＮこれはまた、Ｐ範囲乗算が行われる度毎に寄生要素Ｉ＝２^-nｍｏｄＮが導入されることを意味する。 ρ演算子は次のように定義されるが、Ｐ≡Ａ・Ｂ・ＩｍｏｄＮ≡ρ（Ａ・Ｂ）Ｎこれを「Ｐ範囲でのＡ×Ｂの乗算」またはモントゴメリー乗算と呼ぶ。Ｐ範囲からの検索はＰ・Ｈに対するρの演算によって計算され、次式を生じる。 ρ（Ｐ・Ｈ）Ｎ≡Ａ・ＢｍｏｄＮ前の合同でＰを減算することでＨの値が導出され、次式が発見される。 ρ（Ｐ・Ｈ）Ｎ≡（Ａ・Ｂ・Ｉ）（Ｈ）（Ｉ）ｍｏｄＮ（次を参照：Ａ・Ｂ・Ｉ←Ｐ、Ｈ←Ｈ、Ｉ←及びすべての乗算演算は寄生Ｉを導入する）ＨがＩ²の倍数の逆と合同な場合、合同は有効であり、すなわち次式である。Ｈ＝Ｉ^-2ｍｏｄＮ≡２²ⁿｍｏｄＮ（ＨはＮの関数であり、Ｈパラメータと呼ばれる）従来のモントゴメリー法では、Ａ・Ｂにρ演算子を定めるため、事前計算定数Ｊを使用する次の処理が利用される。１）Ｘ＝Ａ・Ｂ２）Ｙ＝（Ｘ・Ｊ）ｍｏｄ２ⁿ（ｎＬＳビットだけが必要である）３）Ｚ＝Ｘ＋Ｙ・Ｎ４）Ｓ＝Ｚ／２ⁿ（Ｊに対する要求はＺを２ｎで割り切れるものにするということである）５）Ｐ￥ＳｍｏｄＮ（Ｓ≧Ｎの場合、ＮはＳから減算される）最後に、ステップ５）で、Ｐ￥ρ（Ａ・Ｂ）Ｎ［Ｎの減算の後、必要ならば、Ｐ＝ρ（Ａ・Ｂ）Ｎ］上記に続いて、Ｙ＝Ａ・Ｂ・Ｊｍｏｄ２ⁿ（ｎＬＳビットだけを使用する）かつＺ＝Ａ・Ｂ＋（Ａ・Ｂ・Ｊｍｏｄ２ⁿ）・ＮＺが２ⁿで割り切れるようにするため（ＺのｎＬＳビットは好適にはゼロ）かつ次の合同が存在する。［Ａ・Ｂ＋（Ａ・Ｂ・Ｊｍｏｄ２ⁿ）・Ｎ］ｍｏｄ２ⁿ≡０この合同が存在するためには、Ｎ・Ｊｍｏｄ２ⁿが−１と合同であるか、またはＪ≡−Ｎ^-1ｍｏｄ２ⁿ であり、定数Ｊが発見された。従って、ＪはＮだけの関数である事前計算定数である。しかし、ＭＭ結果をビット毎に出力する機械では、ＬＳ文字列の出力ビットがＮｓを加算しなければゼロになる場合毎にＮｓを加算することを規定し、Ｊを事前計算し、次いでＹ＝Ａ・Ｂ・Ｊｍｏｄ２ⁿ（Ｙはハードワイヤード論理を使用してビット毎に検出される）を計算する必要を回避すべきである。また、この方法は奇数のＮｓに対してだけ使用できることも説明した。従って、明らかなように、説明された処理は、所与のＡ、Ｂ、Ｎ及び事前計算定数について３つの乗算、１つの合計及び最大１つの減算を使用してρ（Ａ・Ｂ）Ｎを得る。この結果と、同じ処理及び事前計算定数Ｈ（モジュロＮの関数）を使用して、ＡＢｍｏｄＮを求めることができる。ＡはＢと等しくてもよいので、この基本演算子は平方またはモジュラ演算での乗算のための装置としても使用することができる。インタリーブド・モントゴメリー・モジュラ乗算本節は、全てｎビット長でありその結果が２ｎ＋１ビットの記憶領域を必要とする乗算を含んだモジュラ乗算の方法を説明する。（Ｄｕｓｓｅによる上記の論文で説明された）モントゴメリーのインタリーブド減算を使用すると、より短いオペランド、レジスタ及びハードウェア乗算器で乗算演算を行うことができ、比較的少ない論理ゲートを有する電子装置の実現が可能になる。まず、インタリーブの各反復毎に、Ｊ₀定数を使用してＮを加算する回数を計算する装置の動作方法を説明する。その後、Ｙ₀のハードワイヤ導出を使用するインタリーブの方法を説明するが、これは各乗算のＪ₀₊段階を除去し｛次の例の２）｝、２つの別個の直列乗算器の機能を、新しい単一汎用乗算器に統合することができるが、これは同様のシリコン・リソースを使用して２倍を越える速度でＡ・Ｂ＋Ｃ・Ｎ＋Ｓを行うことができる。ｋビット乗数を使用すると、ｋビット長の文字を定義することが好都合である。ｎビットにはｍ文字が存在する、すなわちｍ・ｋ＝ｎである。Ｊ₀はＪのＬＳ文字である。従って、Ｊ₀≡−Ｎ₀ ^-1ｍｏｄ２^k（Ｊ₀はＮが奇数の時存在する）Ｊ及びＪ₀は、補償数であり、あり得る出力に対して規定される場合、最下位ゼロの所定の数を得るために何回法を加算すべきかを知らせることに注意されたい。後でこの直列装置に対する付加的な利点を説明するが、出力の次の直列ビットが用意に判定できるので、常に法（常に奇数）を次の中間結果に加算することができる。これは、この加算なしで、出力ビットである、ＣＳＡを出るＬＳ直列ビットが「１」であり、それによって法を前の偶数中間結果に加算し、それによって出力文字列で次のＬＳゼロを約束する場合そうである。法を何回結果に加算しようとも、剰余は一定であるので、合同が維持されることを想起されたい。モントゴメリーのインタリーブド減算の従来の使用では、ρ（Ａ・Ｂ）Ｎはステップ（１）〜（５）で説明されるように、ｍ回の反復で規定される。まずＳ（０）＝０（最初の反復開始時のＳの￥値）ｉ＝１，２．．．ｍについて、１）Ｘ＝Ｓ（ｉ−１）＋Ａ_i-1・Ｂ（Ａ_i-1はＡのｉ−１番目の文字、Ｓ（ｉ −１）はｉ番目の反復開始時のＳの値）２）Ｙ₀＝Ｘ₀・Ｊ₀ｍｏｄ２^k（Ｘ₀・Ｊ₀の積のＬＳｋビット）（この処理は、ｋＬＳビットだけ、例えば最下位６４ビットを使用し計算する）好適実現例では、直列機械ではＹ₀はビット毎に予想できるためこのステップは回避される。３）Ｚ＝Ｘ＋Ｙ₀・Ｎ４）Ｓ（ｉ）＝Ｚ／２^k（ＺのｋＬＳビットは常に０であるので、Ｚは常に２^kで割り切れる。ＺのＬＳｋビットはすべてゼロであるのでこの除算はｋビット右シフトと同等である、または回路に見られるように、ＺのＬＳｋビットは単純に無視される）（５）Ｓ（ｉ）＝Ｓ（ｉ）ｍｏｄＮ（Ｎは、Ｎより大きいＳ（ｉ）から減算される）最後に、最後の反復で（必要な場合Ｎの減算の後で）、Ｃ＝Ｓ（ｍ）＝ρ（Ａ・Ｂ）Ｎである。Ｐ＝Ａ・ＢｍｏｄＮを導出するために、Ｐ範囲計算、ρ（Ｃ・Ｈ）Ｎが行われる。好適実施形態では、すべてのＳ（ｉ）について、Ｓ（ｉ）が２Ｎより小さいことを知ることが望ましい。これはまた、最終結果（Ｓ（ｍ））が常に、最大１回のＮの減算で、Ｎより小さい量まで減らせることを意味する。この処理で使用されるオペランドについて以下が観察される。Ｓ（ｉ−１）＜２ⁿ⁺¹（一時レジスタはＢまたはＮレジスタより１ビット長いことがある）Ｂ＜Ｎ＜２ⁿかつＡ_i-1＜２^k 定義により、Ｓ（ｉ）＝Ｚ／２^k（あり得る減算の前の、処理終了時のＳの値）すべてのＺについて、Ｚ（ｉ）＜２^n+k+1 Ｘ_max＝Ｓ_max＋Ａ_i・Ｂ＜２ⁿ⁺¹−１＋（２^k−１）（２ⁿ−１）Ｑ_max＝Ｙ₀Ｎ＜（２^k−１）（２ⁿ−１）従って、Ｚ_max＜２^k+n+1−２^k-1＋１＜２^k+n+1−１であり、Ｚ_maxを２^kで除算すると、Ｓ（ｍ）＜２ⁿ⁺¹−２¹ Ｎ_min＞２ⁿ−２なので、Ｓ（ｍ）_maxは常に２・Ｎ_minより小さく、従って最終結果のために必要なのは１つの減算だけである。Ｓ（ｍ）_max−Ｎ_min＝（２ⁿ⁺¹−２¹−１）−（２ⁿ−１）＝２ⁿ−４＜Ｎ_min モントゴメリー・インタリーブド・モジュラ乗算の例１６進法形式での次の計算はインタリーブド法の意味を明らかにする。Ｎ＝ａ５９（モジュロ）、Ａ＝９９ｂ（乗数）、Ｂ＝５ｃ３（被乗数）、ｎ＝１２（Ｎのビット長）、ｋ＝４（乗数のビット単位のサイズ、かつ文字のサイズ）及びｎ＝ｋ・ｍなのでｍ＝３７．９≡−１ｍｏｄ１６及びＨ≡２^2・12ｍｏｄａ５９≡４４ｂなのでＪ₀＝７予想される結果は、Ｆ≡Ａ・ＢｍｏｄＮ≡９９ｂ・５ｃ３ｍｏｄａ５９≡３７５８１１ｍｏｄａ５９＝２２０₁₆ 初め、Ｓ（０）＝０ステップ１Ｘ＝Ｓ（０）＋Ａ₀・Ｂ＝０＋ｂ・５ｃ３＝３ｆ６１Ｙ₀＝Ｘ₀・Ｊ₀ｍｏｄ２^k＝７（Ｙ₀−新しいＭＡＰにおけるハードワイヤ予想）Ｚ＝Ｘ＋Ｙ₀・Ｎ＝３ｆ６１＋７・ａ５９＝８７ｄ０Ｓ（１）＝Ｚ／２^k＝８７ｄステップ２Ｘ＝Ｓ（１）＋Ａ₁・Ｂ＝８７ｄ＋９・５ｃ３＝３ｃ５８Ｙ₀＝Ｘ₀・Ｊ₀ｍｏｄ２^k＝８・７ｍｏｄ２⁴＝８（ハードワイヤ予想）Ｚ＝Ｘ＋Ｙ₀−Ｎ＝３ｃ５８＋５２ｃ８＝８ｆ２０Ｓ（２）＝Ｚ／２^k＝８ｆ２ステップ３Ｘ＝Ｓ（２）＋Ａ₂・Ｂ＝８ｆ２＋９／５ｃ３＝３ｃｃｄＹ₀＝ｄ・７ｍｏｄ２⁴＝ｂ（ハードワイヤ予想）Ｚ＝Ｘ＋Ｙ₀・Ｎ＝３ｃｃｄ＋ｂ・ａ５９＝ａｅａ０Ｓ（３）＝Ｚ／２^k＝ａｅａＳ（３）＞Ｎなので、Ｓ（ｍ）＝Ｓ（３）−Ｎ＝ａｅａ−ａ５９＝９１従ってＣ＝ρ（Ａ・Ｂ）Ｎ＝９１₁₆ Ｐ範囲からの検索がρ（Ｃ・Ｈ）Ｎを計算することによって行われる。再び初め、Ｓ（０）＝０ステップ１Ｘ＝Ｓ（０）＋Ｃ₀・Ｈ＝０＋１・４４ｂ＝４４ｂＹ₀＝ｄ（新しいＭＡＰにおけるハードワイヤ予想）Ｚ＝Ｘ＋Ｙ₀・Ｎ＝４４ｂ＋８６８５＝８ａｄ０Ｓ（１）＝Ｚ／２ｋ＝８ａｄステップ２Ｘ＝Ｓ（１）＋Ｃ₁・Ｈ＝８ａｄ＋９・４４ｂ＝２ｆ５０Ｙ₀＝０（新しいＭＡＰにおけるハードワイヤ予想）Ｚ＝Ｘ＋Ｙ₀・Ｎ＝２ｆ５０＋０＝２ｆ５０Ｓ（２）＝Ｚ／２^k＝２ｆ５ステップ３Ｘ＝Ｓ（２）＋Ｃ₂・Ｈ＝２ｆ５＋０・４４ｂ＝２ｆ５Ｙ₀＝３（新しいＭＡＰにおけるハードワイヤ予想）Ｚ＝Ｘ＋Ｙ₀・Ｎ＝２ｆ５＋３ａ・５９＝２２００Ｓ（３）＝Ｚ／２^k＝２２０₁₆ これは９９ｂ５ｃ３ｍｏｄａ５９の期待値である。各ステップでｋＬＳゼロを無視すると、本質的にｎＭＳビットを２^kで乗算することになる。同様に、各ステップで、乗数のｉ番目のセグメントも２^ikによって乗算された数であり、それにＳ（ｉ）と同じランクを与える。別の好適実施形態でも、それがＪ₀定数を知るためのある予想値のものであることに注意されたい。累乗次のシーケンスの導出［Ｄ．Ｋｎｕｔｈ、コンピュータ・プログラミングの技術、第２巻、準数値的アルゴリズム、Ａｄｄｉｓｏｎ−Ｗｅｓｌｅｙ、マサチューセッツ州レディング、１９８１年］（以下「Ｋｎｕｔｈ」と呼ぶ）はモジュラ累乗を実現する平方と乗算の順序を説明する。モントゴメリー定数Ｈ＝２²ⁿの事前計算後、この装置はＰ範囲での平方と乗算を行うことができるので、以下を計算する。Ｃ＝Ａ^EｍｏｄＮＥ（ｊ）が、インデックスがＩであるＭＳビットから始まりインデックスがｑであるＬＳビットで終わる指数Ｅのバイナリ表示でｊビットを示すとすると、奇数の指数について以下のように累乗を行うことができる。Ａ^*￥ρ（Ａ・Ｈ）ＮＡ^*はここではＡ・２ⁿに等しい。Ｂ＝Ａ^* ｊ＝２〜ｑ−１についてＢ￥ρ（Ｂ・Ｂ）ＮＥ（ｊ）＝１ならば、Ｂ￥ρ（Ｂ・Ａ^*）Ｎ終了Ｂ￥ρ（Ｂ・Ａ^*）ＮＥ（０）＝１：Ｂは２ⁿで乗算された最終的な望ましい結果であり、ＡはもとのＡである。Ｃ＝ＢＣ≧ＮならばＣ＝Ｃ−Ｎ最後の反復後、値Ｂは￥からＡ^EｍｏｄＮであり、Ｃは最終的な値である。これを明らかにするために、次の例を使用する。Ｅ＝１０１１→Ｅ（１）；Ｅ（２）＝０；Ｅ（３）＝１；Ｅ（４）＝１Ａ¹⁰¹¹ｍｏｄＮ；ｑ＝４を発見するために、Ａ^*＝ρ（Ａ・Ｈ）Ｎ＝ＡＩ^-2Ｉ＝ＡＩ^-1ｍｏｄＮＢ＝Ａ^* ｊ＝２〜ｑについて、Ｂ＝ρ（Ｂ・Ｂ）Ｎ、これはＡ²（Ｉ^-1）²Ｉ＝Ａ²・Ｉ^-1を生じる。Ｅ（２）＝０、Ｂ＝Ａ²・Ｉ^-1 ｊ＝３Ｂ＝ρ（Ｂ・Ｂ）Ｎ＝Ａ²（Ｉ^-1）²Ｉ＝Ａ⁴Ｉ^-1 Ｅ（３）＝１Ｂ＝ρ（Ｂ・Ａ^*）Ｎ＝（Ａ⁴・Ｉ^-1）（ＡＩ^-1）・Ｉ＝Ａ⁵・Ｉ^- ¹ ｊ＝４Ｂ＝ρ（Ｂ・Ｂ）Ｎ＝Ａ¹⁰・Ｉ^-2・Ｉ＝Ａ¹⁰・Ｉ^-1 Ｅ（４）は奇数だったので、最後の乗算はＡによって、寄生Ｉ^-1を除去するために行われる。Ｂ＝ρ（Ｂ・Ａ）Ｎ＝Ａ¹⁰・Ｉ^-1・Ａ・Ｉ＝Ａ¹¹ Ｃ＝Ｂ逆数処理によってＨパラメータを計算する方法は米国特許第５，５１３，１３３号で説明されている。長い除算処理図２は、図１Ａのプロセッサで利用可能なデータ操作装置を使用して長い除算を行う決定論的処理装置の好適実施形態を例示する。図４及び図５は、図２の装置の演算の例である。除数ｄがｄ、２^n-1＜ｄ＜２ⁿの範囲内で、被除数ＤがＤ、２^2n-1＜Ｄ＜２²ⁿの範囲内である除算処理では、装置は最も簡単に使用される。ｄが図２のＮ右シフト・レジスタ１００５（Ｎ０及びＮ１レジスタ、２１０及び２００の連結）に事前ロードされ、被除数ＤのＭＳｎビットがＢ右シフト・レジスタ（１０００）に事前ロードされた後、ＤのＬＳビットがＳ_A右シフト・レジスタ（１３０）に逆ロードされる。これは本質的に、手動の長い除算のために数字を配置するように、新しい試行減算毎に新しいＬＳビットがＳ_AからＳ_Bに供給されることを実現する。Ｓ_Bレジスタはすべてゼロで事前ロードされる。初期化反復では、図２のオーバフロー・フリップフロー１７０がまずリセットされる。ここでは除数ｄを含むＮレジスタがＤＣのＭＳビットであるＢから試行減算される一方、Ｂ及びＮが回転されるが、ここでそれらの出力は直列減算器９０に供給され、その出力はＢ−ｆＮであり、減算の場合ｆ＝１、減算しない場合ｆ＝０である。商生成器１２０はＢ≧Ｎかを判定する検出器であり、９０で次の反復の際ｄがＢから減算されるかを判定する「次の減算」を送信する。この信号は、ＢがＮより大きいかそれと等しい場合、かつその場合のみ１である。この「次の減算」ビットは成功を示し、商の最上位ビットであり、やはりＳ_Bレジスタ１８０にクロックされる。これがＳ_Bレジスタにクロックされているので、ゼロはＳ_BレジスタからＳ_Aレジスタにシフトされ、「新しいＬＳビット」をＳ_Aレジスタから出す。これらのレジスタは全て右シフト・レジスタなので、Ｓ_Bの商の値とＳ_Aの被除数の値はどちらも逆の順序で保持されることに注意されたい。次の反復では、最初の有効クロックサイクルで、Ｓ_Aからシフトされた「新しいＬＳビット」が９０の出力Ｒの前に置かれ、それによってＲに２を乗算し、「新しいＬＳビット」の値を２Ｒに加算する。この連結値は回転されてＢレジスタに戻され、１２０で試験されて、次のラウンドでｆ＝１かｆ＝０かを判定する。最後に、剰余がＢレジスタにあり、商のｎビットがＳ_Bレジスタにあり、商の最上位ビットがＳ_Aレジスタにある。被除数と除数両方の最上位は好適には、すべての大きさのＤとｄについて、ＢレジスタとＮレジスタの最上位ビット・セルにある。結果を得るために必要な反復の数は２^2n-1より小さいＤｓの場合減少し、２ｎ−１より小さいｄｓの場合増大する。装置はハードウェア駆動型であり、装置をアンロードする際オペランドをシフトするファームウェア補償が提供される。図３の不揮発性メモリにあるプログラムは、好適にはレジスタがロードされていることを確認し、制御レジスタのために除算処理を成功させるために必要な反復の数を定義する。逆データ出力アンローダ、ユニット６０及び４０を通じて処理される場合、商ビットはバイト毎に再逆転される。このプロセッサはＨパラメータの計算に有益な要素であり、また好適にはユークリッド関数の計算でも使用される。図１Ａの直列整数除算装置と、Ａ×Ｂ＋Ｃ×Ｄ＋Ｓを行う図１Ｂの二重動作乗算器は、通常同時に動作しない。図４の例では、被除数Ｂは１８７（１０１１１０１１）であり、除数Ｎは７（１１１）であり、除算が実行されると、商は２０（１０１００）となり余りは７（１１１）となる。図５の例では、被乗数Ｂは１７３（１０１０１１０１）、除数Ｎは５（１０１）であり、除算か実行されると、商は３４（１０００１０）となり余りは３（１１）となる。桁上げ保存累算器はＧｒｅｓｓｅｌに対する米国特許第５，５１３，１３３号の図５に例示されている。直列全加算器はＧｒｅｓｓｅｌに対する上記の米国特許の図７に例示されている。直列減算器はＧｒｅｓｓｅｌに対する上記の米国特許の図８に例示されている。除算装置はＧｒｅｓｓｅｌに対する上記の米国特許の図９に例示されている。「整数の正規範囲」という術語は負でない整数、例えば自然数のことである。本発明の好適実施形態によれば、ここで示され説明されたシステムはａ及びｂをゼロにリセットしＳ₀＝１を設定することでＪ₀を計算するよう動作する。この特許書類の開示に一部は著作権保護の対象となる資料を含んでいる。著作権所有者は特許書類または特許開示が特許商標局の特許ファイルまたは記録にある限り、その何者かによるファクシミリ複写に反対しないが、そうでない場合は何があろうともすべての著作権を留保する。本発明のソフトウェア構成要素は、望ましい場合、ＲＯＭ（読み出し専用記憶素子）の形態で実現される。ソフトウェア構成要素は、一般に、望ましい場合、従来の技術を使用してハードウェアで実現される。本発明の様々な特徴は、わかりやすくするために別個の実施形態の文脈で説明されたが、１つの実施形態に組み合わせて提供されることもあることが認識されるだろう。逆に、簡潔にするために１つの実施形態の文脈で説明された様々な特徴が個別または何らかの適切な形でさらに組み合わされて提供されることもある。当業技術分野に熟練した者には、本発明はここに示され説明されたものに制限されないことが認識されるだろう。むしろ、本発明の範囲は以下の請求項のみによって定義される。

───────────────────────────────────────────────────── フロントページの続き (81)指定国ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＤＥ，ＤＫ，ＥＳ，ＦＩ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，ＬＵ，ＭＣ，ＮＬ，ＰＴ，ＳＥ)，ＯＡ(ＢＦ，ＢＪ，ＣＦ，ＣＧ，ＣＩ，ＣＭ，ＧＡ，ＧＮ，ＭＬ，ＭＲ，ＮＥ，ＳＮ，ＴＤ，ＴＧ)，ＡＰ(ＧＨ，ＧＭ，ＫＥ，ＬＳ，ＭＷ，ＳＤ，ＳＺ，ＵＧ，ＺＷ)，ＥＡ(ＡＭ，ＡＺ，ＢＹ，ＫＧ，ＫＺ，ＭＤ，ＲＵ，ＴＪ，ＴＭ)，ＡＬ，ＡＭ，ＡＴ，ＡＵ，ＡＺ，ＢＡ，ＢＢ，ＢＧ，ＢＲ，ＢＹ，ＣＡ，ＣＨ，ＣＮ，ＣＵ，ＣＺ，ＤＥ，ＤＫ，ＥＥ，ＥＳ，ＦＩ，ＧＢ，ＧＥ，ＧＨ，ＧＭ，ＧＷ，ＨＵ，ＩＤ，ＩＬ，ＩＳ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ，ＬＣ，ＬＫ，ＬＲ，ＬＳ，ＬＴ，ＬＵ，ＬＶ，ＭＤ，ＭＧ，ＭＫ，ＭＮ，ＭＷ，ＭＸ，ＮＯ，ＮＺ，ＰＬ，ＰＴ，ＲＯ，ＲＵ，ＳＤ，ＳＥ，ＳＧ，ＳＩ，ＳＫ，ＳＬ，ＴＪ，ＴＭ，ＴＲ，ＴＴ，ＵＡ，ＵＧ，ＵＳ，ＵＺ，ＶＮ，ＹＵ，ＺＷ (72)発明者グレッセル，カーミデビッドイスラエル国，85530 モービルポストネジェブ，キブーツユーリム (72)発明者ドロール，イタイイスラエル国，84451 ビアシェバ，ミブツァナクソンストリート 78／32

Claims

【特許請求の範囲】１．モジュラ乗算及び累乗システムであって、１つの桁上げ保存加算器を有する１つのモジュラ乗算装置を含む直並列論理演算装置（ＡＬＵ）、を備えるシステム。２．任意のビット長の少なくとも１対の整数入力を乗算するよう動作する、請求項１に記載のシステム。３．前記少なくとも１対の整数入力が２対の整数入力を備える、請求項２に記載のシステム。４．前記ＡＬＵが、前もってゼロ強制モントゴメリー定数Ｊ₀を計算せずに、整数入力の積を生成し、前記積の大きさを低減する、先行する請求項の何れかに記載のシステム。５．直列整数除算システムであって、任意のビット長の被除数と任意のビット長の除数を受信し商と剰余を計算するよう動作する直列除算装置、を備えるシステム。６．１対の整数入力を保存する１対のレジスタを更に備え、前記システムが、その少なくとも１つがその対応するレジスタのビット長を越える対応する対の整数入力を、インタリーブなしで乗算するよう動作する、請求項１から請求項３までのいずれか１項に記載のシステム。７．任意のビット長の被除数と任意のビット長の除数を受信し商と剰余を計算するよう動作する直列除算装置を更に備える、請求項１から請求項３までのいずれか１項に記載のシステム。８．モジュラ乗算及び累乗システムであって、１つだけの桁上げ保存累算器を有し、１対の乗算を行いその結果を合計するよう動作する直並列乗算装置、を備えるシステム。９．モジュラ乗算及び累乗方法であって、１つの桁上げ保存加算器を有する１つのモジュラ乗算装置を含む直並列論理演算装置（ＡＬＵ）を提供するステップと、モジュラ乗算及び累乗を行うために前記直並列ＡＬＵを利用するステップと、を含む方法。１０．大きな整数の自然（非モジュラ）乗算の方法であって、１つの桁上げ保存加算器を有する１つのモジュラ乗算装置を含む直並列論理演算装置（ＡＬＵ）を提供するステップと、大きな整数の自然（非モジュラ）乗算を行うために前記直並列ＡＬＵを利用するステップと、を含む方法。１１．前記利用ステップが、第１積を得るために任意のビット長の第１整数を任意のビット長の第２整数によって乗算するステップと、第２積を得るために任意のビット長の第３整数を任意のビット長の第４整数によって乗算するステップと、合計を得るために前記第１及び第２積と任意のビット長の第５整数を合計するステップと、を含む、請求項１０に記載の方法。１２．前記利用ステップが任意のビット長の被乗数、乗数及び法によってモジュラ乗算及び累乗を行うステップを含む、請求項９に記載の方法。１３．１つだけの事前計算定数でモントゴメリー・モジュラ乗算を実行する二重被乗数事前計算システムを備える、請求項８に記載のシステム。１４．前記利用ステップが、乗数と被乗数を含む整数入力の積を生成するステップと、前もってモントゴメリー定数Ｊ₀を計算せずに、モジュラ換算を実行するステップとを含むモントゴメリー乗算を行うステップと、を含む、請求項９に記載の方法。１５．前記モントゴメリー定数Ｊ₀が、Ｎが前記モジュラ換算の法でありｋが前記被乗数のビット長である、Ｎｍｏｄ２^kの関数を含む、請求項１４に記載の方法。１６．前記利用ステップが一連のインタリーブド・モントゴメリー乗算演算を行うステップを含む、請求項９に記載の方法。１７．前記インタリーブド・モントゴメリー乗算演算が、少なくともｋの有意ゼロを伴う結果を得るために法を前記乗算演算の合同に合計しなければならない回数を前もって計算せずに、前記インタリーブド・モントゴメリー乗算演算が行われる、請求項１６に記載の方法。１８．ｉ番目のインタリーブド・モントゴメリー乗算演算で生成される被乗数を収集し直列合計するよう動作し、それによって合計を生成し前記合計を（ｉ＋１）番目のモントゴメリー乗算演算に供給するデータ・プロセッサを備える、請求項８に記載のシステム。１９．前記関数がＮｍｏｄ２^kの乗法的逆元の加法的逆元を含む、請求項１５に記載の方法。２０．ａ及びｂをゼロにリセットしＳ₀＝１を設定することでＪ₀を計算するステップを更に含む、請求項１５に記載の方法。２１．ａ及びｂをゼロにリセットしＳ₀＝１を設定することでＪ₀を計算するステップを更に含む、請求項１５に記載の方法。２２．前もって計算せずに実行する前記ステップが、整数入力の前記積を生成する過程で生成される乗法的合計に法を加算しなければならないか否かを予想するステップを含む、請求項１４に記載の方法。２３．前記ＡＬＵが、前記積が、前記法乗算装置が動作する法より小さいか否かを判定し、それによって前記積の大きさを換算すべきか否かを判定する手段を備える、請求項４に記載のシステム。２４．１対の整数入力を保存する１対のレジスタを更に備え、前記システムが、その少なくとも１つが対応するレジスタのビット長を越える対応する対の整数入力を、インタリーブなしで乗算するよう動作する、請求項４に記載のシステム。２５．任意のビット長の被除数と任意のビット長の除数とを受信し商と剰余とを計算するよう動作する直列除算装置を更に備える、請求項４に記載のシステム。２６．任意のビット長の被除数と任意のビット長の除数とを受信し商と剰余とを計算するよう動作する直列除算装置を更に備える、請求項２４に記載のシステム。