JP2014164144A - 秘密表除算装置及び方法 - Google Patents

Abstract

【課題】従来よりも短い計算時間で表の除算を実現する秘密表除算技術を提供する。
【解決手段】２つの表SとTについて、SとTを上下に並べた表を作り、Tの属性についてソートを行った後にSにしか含まれない属性について集計をすることにより、直積の演算を用いずに表の除算をするための処理を構成する。これにより、計算途中で行数がO(m+n)の表しか作る必要がないため、従来よりも計算時間を短くすることができる。具体的には、計算時間がO((m+n)log(m+n))である秘密計算上での表の除算を実現することができる。
【選択図】図１

Description

この発明は、暗号技術に関する。特に、秘密計算によって、表に含まれる情報を秘密にしたまま演算を行う技術に関する。

暗号化された数値を復元すること無く特定の演算結果を得る方法として、秘密計算と呼ばれる技術が知られている（例えば、非特許文献１参照。）。非特許文献１に記載された技術では、３つの秘密計算装置に数値の断片を分散させるという暗号化を行い、これらの３つの秘密計算装置が協調計算を行うことにより、数値を復元すること無く、加減算、定数和、乗算、定数倍、論理演算（否定、論理積、論理和、排他的論理和）、データ形式変換（整数、二進数）の結果を３つの秘密計算装置に分散された状態、すなわち暗号化されたまま保持させることができる。

また、秘密計算上で表の除算を実現した方法として、非特許文献２の技術が知られている。非特許文献２の技術では、表の除算を直積と差の演算の組み合わせにより実現している。

千田浩司, 濱田浩気, 五十嵐大, 高橋克巳, "軽量検証可能３パーティ秘匿関数計算の再考", In CSS, 2010. 志村正法, 遠藤つかさ, 宮崎邦彦, 吉浦裕, "安全で機能制限のないデータベースを実現するマルチパーティプロトコルを用いた関係代数演算", 情報処理学会研究報告. CSEC, Vol.2008, No.71, pp.187-193, 2008.

しかしながら、従来技術は直積の演算を要するため、入力の２つの表の行数をそれぞれm，nとすると計算時間がO(m²n)と膨大である。

この発明の目的は、従来よりも短い計算時間で表の除算を実現する秘密表除算装置及び方法を提供することである。

この発明の一態様による秘密表除算装置は、aを秘匿化している情報を<a>とし、表Sを表Tで除算した表をR₁をとして、<S>及び<T>に基づいて<R₁>を秘密計算する秘密表除算装置であって、各表は第一方向に配置された少なくとも１つのレコードから構成されているとし、各レコードは少なくとも１つの属性にそれぞれ対応する第二方向に配置された少なくとも１つの属性値から構成されているとし、表Sを構成するレコード数をmとし、表Tを構成するレコード数をnとして、上記秘密表除算装置は、少なくとも１つの秘密計算装置を含み、これら秘密計算装置の協調計算によって、レコード数が表Tのレコード数と同じであり、表Sの属性と同じ属性を持ち、表Sと表Tに共通の属性についてはその共通の属性に対応する表Tの属性値を取り、表Sのみに存在する属性については任意の値を取る表をT’とし、表S及び表T’を第一方向に結合させた表をUとして、<S>及び<T>に基づいて<U>を計算し、0を第一方向にm個並べたベクトルを0^mとし、1を第一方向にn個並べたベクトルを1ⁿとし、0^m及び1ⁿを第一方向に結合させたベクトルをtとして、<0^m>及び<1ⁿ>に基づいて<t>を計算し、表Uの中の表Tの属性に対応する属性値により構成される表をKとして、表KをキーとしてU,tを第一方向に安定ソートした表をそれぞれU₁,t₁として、<U>,<t>に基づいて<U₁>及び<t₁>を計算し、表U₁の中の表Tの属性に対応する属性値により構成される表をVとして、<U₁>に基づいて<V>を計算し、表Vを構成する第二方向ベクトルの中に、対応する、ベクトルt₁の要素が1であるものがあれば、その対応する要素が1である第二方向ベクトルと同じ要素で構成される全ての第二方向ベクトルに対応する、ベクトルt₁の要素を1にする処理を行ったベクトルをgとして、<V>及び<t₁>に基づいて<g>を計算し、π_rを所定のランダム置換として、U₁,t₁,gをπ_rで置換したものをU₂,t₂,g₁として、<U₁>,<t₁>,<g>に基づいて<U₂>,<t₂>,<g₁>を計算し、ベクトルt₂の中の値が１である要素に対応する、U₂,g₁の第二方向ベクトルを削除したものをU₃,g₂として、<U₂>,<t₂>,<g₁>に基づいて<U₃>,<g₂>を計算し、表Sに存在するが表Tに存在しない属性の属性値を表U₃から取り出した表をRとして、<U₃>,<S>,<T>,に基づいて<R>を計算し、Rをキーとしてg₂について集約関数総和を行った結果をcとして、<R>に基づいて<c>を計算し、cの中の値がnでない要素に対応する要素の値を1とし、cの中の値がnである要素に対応する要素の値を0とするベクトルをeとして、<c>に基づいて<e>を計算し、e及びRを第二方向に結合させた表をR₁として、<e>及び<R>に基づいて<R₁>を計算する。

従来よりも短い計算時間で表の除算を実現することができる。

秘密表除算装置の例の機能ブロック図。 表の例を説明するための図。 表の例を説明するための図。

以下、この発明の実施形態について説明する。

［定義］
aを暗号化や秘密分散などにより秘匿化している情報を<a>と表現する。<a>のことをaの秘匿文と呼ぶこともある。aは、スカラー値、ベクトル、行列である。aがベクトルである場合には、<a>はベクトルの各要素を秘匿化した情報である。aが行列である場合には、<a>は行列の各要素を秘匿化した情報である。<a>を復元するための情報は複数の断片に分割されており、複数の断片は協調計算を行う少なくとも１つの秘密計算装置により保有されているものとする。

aがベクトルである場合には、<a>[i]はベクトルaのi番目の要素を意味する。aが行列である場合には、<a>[i]は行列aのi番目の第二方向のベクトルを意味する。i番目の第二方向のベクトルとは、後述するように第一方向が行方向であり第二方向が列方向である場合には、i番目の行ベクトルのことである。

ベクトルxとyの所定の方向での連結をx||yとする。aがスカラー値である場合には秘匿文<a>を所定の方向にn個並べたベクトルを<a>ⁿと表現する。

除算は２つの表による演算である。表Xに対する表Yの除算の結果を表Zとする。

ここでは、各表は、第一方向に配置された少なくとも１つのレコードから構成されている。各レコードは少なくとも１つの属性にそれぞれ対応する第二方向に配置された少なくとも１つの属性値から構成されている。このように、各表は行列である。表を構成するレコードの各属性値に対応する属性を、その表の属性と表現する。

第一方向が行方向である場合には、第二方向は列方向である。一方、第一方向が列方向である場合には、第二方向は行方向である。以下、第一方向が行方向であり、第二方向が列方向である場合を例に挙げて説明する。

図２に例示した表は、行方向に配置された３個のレコードから構成されている。各レコードは、４個の属性にそれぞれ対応する列方向に配置された属性値から構成されている。この表の属性は、ID、年齢、性別、職業である。

表Xに対する表Yによる除算は、表Xの属性の集合αと表Yの属性の集合βとがβ⊂αという関係を満たす場合に定義される。例えば、表Xの属性の集合α={a₁,a₂,a₃,a₄}であり、表Yの属性の集合β={a₃,a₄}である場合には、β⊂αの関係を満たすため表Xに対する表Yによる除算が定義される。＼を差集合の演算を表現するとして、表Xに対する表Yの除算の結果である表Zの属性はγ=α＼βである。この例では、γ=α＼β={a₁,a₂,a₃,a₄}＼{a₃,a₄}={a₁,a₂}である。

表Xに対する表Yの除算の結果である表Zは、以下のように定義される。
Z=π_γ(X)-π_γ((π_γ(X)×Y)-X)

ここで、π_γ(X)は、表Xの中から属性γの属性値を抜き出した表π_γ(X)’において、重複する第二方向ベクトル（行ベクトル）を除いた表である。例えば、属性が(a₁,a₂,a₃,a₄)である表X及び属性が(a₃,a₄)である表Yが以下のように与えられた場合、属性が(a₁,a₂)であるπ_γ(X)’,π_γ(X),(π_γ(X)×Y),(π_γ(X)×Y)-X,π_γ((π_γ(X)×Y)-X),Zは以下のようになる。

秘密表除算装置は、図１に例示するように、Ｎを１以上の整数として、Ｎ個の秘密計算装置を備えている。Ｎ個の秘密計算装置は、協調計算により、秘匿化、復元、加算、減算、乗算、否定、ランダム置換、安定ソート、等号判定、集約関数総和の演算を行うことができる。

秘匿化の演算は、入力されたaに基づいて、aの秘匿文<a>を計算する処理である。復元の演算は、秘匿文<a>に基づいてaを計算する処理である。秘匿化及び復号の演算の詳細については、例えば非特許文献１を参照のこと。

加算、減算、乗算の各演算は、Ｎ個の秘密計算装置による協調計算であり、a,bの秘匿文<a>,<b>を入力とし、それぞれa+b，a-b，abの計算結果であるcの秘匿文<c>を計算する処理である。加算、減算、乗算の演算の詳細については、例えば非特許文献１を参照のこと。

否定の演算は、Ｎ個の秘密計算装置による協調計算であり、aの秘匿文<a>を入力とし、aの否定である¬aの計算結果であるcの秘匿文<c>を計算する処理である。否定の演算の詳細については、例えば非特許文献１を参照のこと。

ランダム置換の演算は、ランダム置換は協調計算により複数の秘匿化されたベクトルに対して同一のランダム置換π_rを施す処理である。すなわち、Lを１以上の整数として、L個の大きさnのベクトルa⁽¹⁾,…,a^(L)の秘匿文<a⁽¹⁾>,…,<a^(L)>を入力とし、あるランダムな全単射π_r:{1,…,n}→{1,…,n}についてb_πr(j) ⁽ⁱ⁾=a_j ⁽ⁱ⁾(1≦i≦L,1≦j≦n)を満たすL個の大きさnのベクトルb⁽¹⁾,…,b^(L)の秘匿文<b⁽¹⁾>,…,<b^(L)>を、Ｎ個の秘密計算装置はどの秘密計算装置もランダム置換π_rを知ることなく協調計算する。ランダム置換の演算の詳細については、例えば参考文献１を参照のこと。

〔参考文献１〕Sven Laur, Jan Willemson, Bingsheng Zhang, “Round-efficient oblivious database manipulation”, In Xuejia Lai, Jianying Zhou, Hui Li, editors, ISC, Vol. 7001 of LNCS, pp.262-277, Springer, 2011.

安定ソートの演算は、秘匿化されたベクトル又は行列の要素を第一方向に所定の順序で並び替える処理である。第一方向とは例えば行方向であり、所定の順序とは例えば辞書順である。この場合、行列Xの安定ソートとは、行列Xを構成するn個の行ベクトルX[1],…,X[n]を辞書順に並び替えることを意味する。なお、辞書順とは、ある行ベクトルx=(x₁,…,x_m)とある行ベクトルy=(y₁,…,y_m)との順序を定める際、x<y iff (x₁<y₁)∨((x₁=y₁)∧((x₂<y₂)∨(…)))という基準に基づいてxとyとの順序を定めるものである。iffは、if and only ifの意味である。なお、安定ソートの際、ある行ベクトルとある行ベクトルが同じであり、これらの行ベクトルに辞書順の順序を定めることができない場合には、元の行列Xにおけるこれらの行ベクトルの順序に基づいてこれらの行ベクトルの順序が定まるとする。また、この場合、行列Xをキーとした、ベクトル又は行列であるYの行方向の安定ソートの演算とは、X,Yの行数をnとして、行列Xの行方向の安定ソートに対応する全単射の所定の置換π_s:{1,…,n}→{1,…,n}に基づいて、Yを構成する行ベクトルを置換する処理を意味する。安定ソートの演算の詳細については、例えば参考文献２を参照のこと。

〔参考文献２〕濱田浩気, 五十嵐大, 千田浩司, 高橋克巳, “秘匿関数計算上の線形時間ソート”, In SCIS, pp.1-7, 2011.

例えば、行列X,Yが以下のように与えられた場合、安定ソート後の行列X’、行列Xをキーとした安定ソート後のY’は以下のようになる。

等号判定の演算は、a,bの秘匿文<a>,<b>を入力とし、a=bの真偽値cの秘匿文<c>を計算する処理である。真偽値は真のとき1、偽のとき0とする。この演算の実行を(<a>=^?<b>)と表現する。等号判定の演算の詳細については、例えば参考文献３を参照のこと。

〔参考文献３〕Takashi Nishide, Kazuo Ohta, “Multiparty computation for interval, equality, and comparison without bitdecomposition protocol”, In PKC, pp.343-360, 2007.

集約関数総和の演算は、行列Xの秘匿文<X>及びベクトルaの秘匿文<a>を入力とし、Xの互いに異なる行ベクトルごとに対応するaの要素の値の総和を計算する演算である。ただし、X,aの行数はnであるとする。この演算を、<X>をキーとする<a>についての集約関数総和の演算とも表現する。出力は、Xを構成する行ベクトルを全単射の所定の置換π_g:{1,…,n}→{1,…,n}で置換することにより並べ替えた行列X’の秘匿文<X’>、及び、同じ置換π_gでベクトルaの要素を置換することにより並び替えた行列a’の秘匿文<a’>である。X’の互いに異なる行ベクトルxごとに、<X’>[i]=xを満たす１つのiについてa’[i]=Σ_j∈C(x)a[j]とし、それ以外のa’[i]については0である。ただし、C(x)={j|X[j]=x}である。集約関数総和の演算の詳細については、例えば参考文献４を参照のこと。

〔参考文献４〕五十嵐大, 千田浩司, 濱田浩気, 高橋克巳, “軽量検証可能3 パーティ秘匿関数計算の効率化及びこれを用いたセキュアなデータベース処理”, In SCIS, pp.1-8, 2011.

例えば、行列X及びベクトルaが以下のように与えられた場合、集約関数総和演算後の行列X’及びベクトルa’は以下のようになる。

［第一実施形態］
秘密表除算装置は、図１に例示するように、Ｎを１以上の整数として、Ｎ個の秘密計算装置を備えている。表Sに対する表Tの除算の結果を表R₁として、Ｎ個の秘密計算装置は、協調計算により、<S>及び<T>に基づいて<R₁>を計算する。

第一実施形態の秘密表除算装置は、秘匿文<R>として出力されるRのレコード数を秘匿化する。

表Sを構成するレコード数をmとし、表Tを構成するレコード数をnとする。以下に、表S,Tの例を挙げる。表Sの属性は(a₁,a₂,a₃,a₄)であり、表Tの属性は(a₃,a₄)であるとする。

＜ステップＳ１＞
Ｎ個の秘密計算装置は、協調計算により、<S>及び<T>に基づいて<U>を計算する（ステップＳ１）。

ここで、レコード数が表Tのレコード数と同じであり、表Sの属性と同じ属性を持ち、表Sと表Tに共通の属性についてはその共通の属性に対応する表Tの属性値を取り、表Sのみに存在する属性については任意の値*を取る表をT’とし、表S及び表T’を第一方向（この例では行方向）に結合させた表をUとする。

表T’,Uは以下のようになる。

＜ステップＳ２＞
Ｎ個の秘密計算装置は、協調計算により、<0^m>及び<1ⁿ>に基づいて<t>を計算する（ステップＳ２）。

ここで、0を第一方向にm個並べたベクトルを0^mとし、1を第一方向（この例では行方向）にn個並べたベクトルを1ⁿとし、0^m及び1ⁿを第一方向（この例では行方向）に結合させたベクトルをtとする。

すなわち、Ｎ個の秘密計算装置は、協調計算により、<t>=<0>^m||<1>ⁿを計算する。

tは以下のようになる。

＜ステップＳ３＞
Ｎ個の秘密計算装置は、協調計算により、<U>,<t>に基づいて<U₁>及び<t₁>を計算する（ステップＳ３）。

ここで、表Uの中の表Tの属性に対応する属性値により構成される表をKとして、表KをキーとしてU,tを第一方向（この例では行方向）に安定ソートした表をそれぞれU₁,t₁とする。

表k,U₁,t₁は、以下のようになる。

＜ステップＳ４＞
Ｎ個の秘密計算装置は、協調計算により、<U₁>に基づいて<V>を計算する（ステップＳ４）。

ここで、表U₁の中の表Tの属性に対応する属性値により構成される表をVとする。

Vは以下のようになる。

＜ステップＳ５＞
Ｎ個の秘密計算装置は、協調計算により、<V>及び<t₁>に基づいて<g>を計算する（ステップＳ５）。

ここで、表Vを構成する第二方向ベクトル（この例では行ベクトル）の中に、対応する、ベクトルt₁の要素が1であるものがあれば、その対応する要素が1である第二方向ベクトル（この例では行方向）と同じ要素で構成される全ての第二方向ベクトル（この例では行方向）に対応する、ベクトルt₁の要素を1にする処理を行ったベクトルをgとする。

ベクトルgは以下のようになる。

Ｎ個の秘密計算装置は、協調計算により、例えば以下の処理を行うことにより<g>を計算する。

まず、<f>=<t₁>とする。すなわち、f=t₁として、<t₁>に基づいて<f>を計算する。また、d=1とする。そして、d<nの間、以下の(a)から(c)の処理を繰り返し行う。

(a) 各i∈[1,m+n]について、以下の式により定義される<p>を計算する。

すなわち、i+d≦nであればf[i+d]*(V[i]=^?V[i+d])をp[i]とし、i+d≦nでなければ0をp[i]として、このように定義されるp[i]に基づいて定まるpを秘匿化した<p>を<f>,<V>に基づいて計算する。

(b) [f]=[f]+[p]とする。つまり、[f],[p]に基づいて、f=f+pのように定義される[f]を計算する。

(c) d=2dとする。

＜ステップＳ６＞
Ｎ個の秘密計算装置は、協調計算により、<U₁>,<t₁>,<g>に基づいて<U₂>,<t₂>,<g₁>を計算する（ステップＳ６）。

ここで、π_rを所定のランダム置換として、U₁,t₁,gをπ_rで置換したものをU₂,t₂,g₁とする。

U₂,t₂,g₁は例えば以下のようになる。

＜ステップＳ７＞
Ｎ個の秘密計算装置は、協調計算により、<U₂>,<t₂>,<g₁>に基づいて<U₃>,<g₂>を計算する（ステップＳ７）。

ここで、ベクトルt₂の中の値が１である要素に対応する、U₂,g₁の第二方向ベクトル（この例では行ベクトル）を削除したものをU₃,g₂とする。

U₃,g₂は以下のようになる。

＜ステップＳ８＞
Ｎ個の秘密計算装置は、協調計算により、<U₃>,<S>,<T>,に基づいて<R>を計算する（ステップＳ８）。

ここで、表Sに存在するが表Tに存在しない属性の属性値を表U₃から取り出した表をRとする。

Rは以下のようになる。

＜ステップＳ９＞
Ｎ個の秘密計算装置は、協調計算により、<R>に基づいて<c>を計算する（ステップＳ９）。

ここで、Rをキーとしてg₂について集約関数総和を行った結果をcとする。

cは以下のようになる。

＜ステップＳ１０＞
Ｎ個の秘密計算装置は、協調計算により、<c>に基づいて<e>を計算する（ステップＳ１０）。

ここで、cの中の値がnでない要素に対応する要素の値を1とし、cの中の値がnである要素に対応する要素の値を0とするベクトルをeとする。

eは以下のようになる。

＜ステップＳ１１＞
Ｎ個の秘密計算装置は、協調計算により、<e>及び<R>に基づいて<R₁>を計算する（ステップＳ１１）。

ここで、e及びRを第二方向（この例では列方向）に結合させた表をR₁とする。

R₁は以下のようになる。

R₁の３列目が0の行の要素うち属性がa₁,a₂である要素を抜き出した行列R₂は、表Sに対する表Tの除算の結果となる。このように、e[i]が1であればその行は空ベクトルであり、e[i]が0であればその行は空ベクトルではない。

したがって、R₁は表の大きさ（行数）を隠しつつ表Sに対する表Tの除算の結果を表している。

従来技術は直積の演算を要するため、入力の2つの表の行数をそれぞれm,nとすると、計算途中で行数がmnの表を作る必要があった。この発明は、直積の演算を用いずに構成されており計算途中で行数がO(m+n)の表しか作る必要がないため、従来よりも計算時間が短い。具体的には、計算時間がO((m+n)log(m+n))である秘密計算上での表の除算を実現することができる。

［第二実施形態］
秘密表除算装置は、第一実施形態を、入力される表S,Tが秘匿化された表である場合に拡張したものである。

以下、第一実施形態と異なる部分を中心に説明し、第一実施形態と同様の部分については重複説明を省略する。

表S,Tは、図３に例示するように、第一方向に配置されたダミーのレコードを含む少なくとも１つのレコードから構成されている。各レコードは少なくとも１つの属性にそれぞれ対応する第二方向に配置された少なくとも１つの属性値及びその各レコードがダミーであるかどうかを表す空フラグから構成されている。例えば、空フラグが1であればその各レコードはダミーを表し空フラグが0であればその各レコードはダミーではないことを表す。

以下、第一方向が行方向であり、第二方向が列方向である場合を例に挙げて説明する。

図３に例示した表は、行方向に配置されたダミーのレコードを含む３個のレコードから構成されている。この例では、３行目のレコードがダミーである。各レコードは、４個の属性にそれぞれ対応する列方向に配置された属性値及び空フラグから構成されている。この例では、空フラグは第１列目に配置されている。この表の属性は、ID、年齢、性別、職業である。このように、表にダミーのレコードを含めることにより、表の真の大きさを秘匿化することができる。

このダミーのレコードの追加は、秘密表除算装置により事前に行われる。

ダミーのレコードを追加した表S,Tの例を以下に挙げる。この例では、空フラグは第１列目に配置されている。

ステップＳ５以外のステップの処理は、第一実施形態と同様である。

第二実施形態のステップＳ５では、<f>=<t₁>とする代わりに、<f>=<t₁>*(¬<L>)とする。すなわち、f=t₁*(¬L)として、<t₁>,<L>に基づいて<f>を計算する。そして、また、d=1とする。そして、d<nの間、上述の(a)から(c)の処理を繰り返し行う。

ここで、LはU₁の空フラグを表す第一方向ベクトル（この例では列ベクトル）である。

この場合、ステップＳ２の後に、U,t₁は以下のようになる。

ステップＳ５の後に、U₁,V,t₁,gは以下のようになる。

また、ステップＳ７の後にU₃,g₂は以下のようになる。

ステップＳ１０の後にR,c,eは以下のようになる。

［変形例］
上記装置及び方法において説明した処理は、記載の順にしたがって時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。

また、秘密表除算装置における各処理をコンピュータによって実現する場合、秘密表除算装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、秘密表除算装置における処理手段がコンピュータ上で実現される。

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。

また、各処理手段は、コンピュータ上で所定のプログラムを実行させることにより構成することにしてもよいし、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

その他、この発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。

Claims (2)

1. aを秘匿化している情報を<a>とし、表Sを表Tで除算した表をR₁をとして、<S>及び<T>に基づいて<R₁>を秘密計算する秘密表除算装置であって、
   各表は第一方向に配置された少なくとも１つのレコードから構成されているとし、各レコードは少なくとも１つの属性にそれぞれ対応する第二方向に配置された少なくとも１つの属性値から構成されているとし、表Sを構成するレコード数をmとし、表Tを構成するレコード数をnとして、
   上記秘密表除算装置は、少なくとも１つの秘密計算装置を含み、これら秘密計算装置の協調計算によって、
   レコード数が表Tのレコード数と同じであり、表Sの属性と同じ属性を持ち、表Sと表Tに共通の属性についてはその共通の属性に対応する表Tの属性値を取り、表Sのみに存在する属性については任意の値を取る表をT’とし、表S及び表T’を第一方向に結合させた表をUとして、<S>及び<T>に基づいて<U>を計算し、
   0を第一方向にm個並べたベクトルを0^mとし、1を第一方向にn個並べたベクトルを1ⁿとし、0^m及び1ⁿを第一方向に結合させたベクトルをtとして、<0^m>及び<1ⁿ>に基づいて<t>を計算し、
   表Uの中の表Tの属性に対応する属性値により構成される表をKとして、表KをキーとしてU,tを第一方向に安定ソートした表をそれぞれU₁,t₁として、<U>,<t>に基づいて<U₁>及び<t₁>を計算し、
   表U₁の中の表Tの属性に対応する属性値により構成される表をVとして、<U₁>に基づいて<V>を計算し、
   表Vを構成する第二方向ベクトルの中に、対応する、ベクトルt₁の要素が1であるものがあれば、その対応する要素が1である第二方向ベクトルと同じ要素で構成される全ての第二方向ベクトルに対応する、ベクトルt₁の要素を1にする処理を行ったベクトルをgとして、<V>及び<t₁>に基づいて<g>を計算し、
   π_rを所定のランダム置換として、U₁,t₁,gをπ_rで置換したものをU₂,t₂,g₁として、<U₁>,<t₁>,<g>に基づいて<U₂>,<t₂>,<g₁>を計算し、
   ベクトルt₂の中の値が１である要素に対応する、U₂,g₁の第二方向ベクトルを削除したものをU₃,g₂として、<U₂>,<t₂>,<g₁>に基づいて<U₃>,<g₂>を計算し、
   表Sに存在するが表Tに存在しない属性の属性値を表U₃から取り出した表をRとして、<U₃>,<S>,<T>,に基づいて<R>を計算し、
   Rをキーとしてg₂について集約関数総和を行った結果をcとして、<R>に基づいて<c>を計算し、
   cの中の値がnでない要素に対応する要素の値を1とし、cの中の値がnである要素に対応する要素の値を0とするベクトルをeとして、<c>に基づいて<e>を計算し、
   e及びRを第二方向に結合させた表をR₁として、<e>及び<R>に基づいて<R₁>を計算する、
   秘密表除算装置。
2. aを秘匿化している情報を<a>とし、表Sを表Tで除算した表をR₁をとして、<S>及び<T>に基づいて<R₁>を秘密計算する秘密表除算方法であって、
   各表は第一方向に配置された少なくとも１つのレコードから構成されているとし、各レコードは少なくとも１つの属性にそれぞれ対応する第二方向に配置された少なくとも１つの属性値から構成されているとし、表Sを構成するレコード数をmとし、表Tを構成するレコード数をnとして、
   少なくとも１つの秘密計算装置が、協調計算によって、
   レコード数が表Tのレコード数と同じであり、表Sの属性と同じ属性を持ち、表Sと表Tに共通の属性についてはその共通の属性に対応する表Tの属性値を取り、表Sのみに存在する属性については任意の値を取る表をT’とし、表S及び表T’を第一方向に結合させた表をUとして、<S>及び<T>に基づいて<U>を計算するステップと、
   0を第一方向にm個並べたベクトルを0^mとし、1を第一方向にn個並べたベクトルを1ⁿとし、0^m及び1ⁿを第一方向に結合させたベクトルをtとして、<0^m>及び<1ⁿ>に基づいて<t>を計算するステップと、
   表Uの中の表Tの属性に対応する属性値により構成される表をKとして、表KをキーとしてU,tを第一方向に安定ソートした表をそれぞれU₁,t₁として、<U>,<t>に基づいて<U₁>及び<t₁>を計算するステップと、
   表U₁の中の表Tの属性に対応する属性値により構成される表をVとして、<U₁>に基づいて<V>を計算するステップと、
   表Vを構成する第二方向ベクトルの中に、対応する、ベクトルt₁の要素が1であるものがあれば、その対応する要素が1である第二方向ベクトルと同じ要素で構成される全ての第二方向ベクトルに対応する、ベクトルt₁の要素を1にする処理を行ったベクトルをgとして、<V>及び<t₁>に基づいて<g>を計算するステップと、
   π_rを所定のランダム置換として、U₁,t₁,gをπ_rで置換したものをU₂,t₂,g₁として、<U₁>,<t₁>,<g>に基づいて<U₂>,<t₂>,<g₁>を計算するステップと、
   ベクトルt₂の中の値が１である要素に対応する、U₂,g₁の第二方向ベクトルを削除したものをU₃,g₂として、<U₂>,<t₂>,<g₁>に基づいて<U₃>,<g₂>を計算するステップと、
   表Sに存在するが表Tに存在しない属性の属性値を表U₃から取り出した表をRとして、<U₃>,<S>,<T>,に基づいて<R>を計算するステップと、
   Rをキーとしてg₂について集約関数総和を行った結果をcとして、<R>に基づいて<c>を計算するステップと、
   cの中の値がnでない要素に対応する要素の値を1とし、cの中の値がnである要素に対応する要素の値を0とするベクトルをeとして、<c>に基づいて<e>を計算するステップと、
   e及びRを第二方向に結合させた表をR₁として、<e>及び<R>に基づいて<R₁>を計算するステップと、
   を含む秘密表除算方法。

Images