SECRET BATCH APPROXIMATION SYSTEM, SECURE COMPUTATION DEVICE, SECRET BATCH APPROXIMATION METHOD, AND PROGRAM

Abstract

Calculation time is reduced without degrading approximation accuracy in calculation of a complicated function through secure computation. A secret batch approximation system calculates a concealed text [z] of an approximate value z for a function value y satisfying yj=f(xj) by using a concealed text [x] of a value x as input. g is defined as a polynomial for approximating each section of m sections into which the function f is divided. A parameter acquisition unit acquires a concealed text [a] of a parameter a corresponding to the value x for each integer j that is not less than 1 and not more than n, where aj is defined as a parameter pi corresponding to a section Ri including a value xj. A polynomial calculation unit calculates a polynomial g([x], [a]) by using the concealed text [x] of the value x as input based on the concealed text [a].

Description

TECHNICAL FIELD

The present invention relates to a technique for performing approximation calculation for a function through secure computation.

BACKGROUND ART

As a method for obtaining a specified operation result without reconstructing encrypted numerical values, there is a method called secure computation (see Non-patent Literature 1, for example). In the method described in Non-patent Literature 1, encryption for distributing shares of a numerical value over three secure computation devices is performed, and the three secure computation devices perform cooperative calculation and are accordingly allowed to hold results of addition-and-subtraction, constant addition, multiplication, constant multiplication, logical operations (a NOT, an AND, an OR, and an XOR), and data format conversion (an integer or a binary) in a state that the results are distributed over the three secure computation devices, that is, in an encrypted state, without reconstructing numerical values.

As a method for realizing calculation of a complicated function such as an exponential function and a logarithmic function on secure computation, calculation is performed by approximating the function to a polynomial which can be calculated through a combination of addition-and-subtraction and multiplication. Non-patent Literature 2, for example, describes a method for performing calculation by approximating an exponential function to a polynomial.

PRIOR ART LITERATURE

Non-Patent Literature

• Non-patent Literature 1: Koji Chida, Koki Hamada, Dai lkarashi, Katsumi Takahashi, “A Three-party Secure Function Evaluation with Lightweight Verifiability Revisited”, CSS, 2010
• Non-patent Literature 2: Liina Kamm and Jan Willemson, “Secure floating point arithmetic and private satellite collision analysis,” International Journal of Information Security, Vol. 14, No. 6, pp. 531-548, 2015.

SUMMARY OF THE INVENTION

Problems to be Solved by the Invention

However, in polynomial approximation of a complicated function, the order of the polynomial needs to be increased so as to improve approximation accuracy, in related art. For example, an approximate formula expressed by the fourth-order polynomial is calculated so as to realize calculation of an exponential function in Non-patent Literature 2. Further, a sigmoid function used in logistic regression and a neural network requires at least the tenth-order polynomial so as to secure sufficient accuracy. When the order n of a polynomial is increased, the number of times of addition-and-subtraction and multiplication is increased in O(n), increasing calculation time.

In view of the above-described technical problems, the present invention is aimed to reduce calculation time without degrading approximation accuracy when a complicated function is calculated through secure computation.

Means to Solve the Problems

In order to solve the above-described problem, a secret batch approximation system according to one aspect of the present invention includes a plurality of secure computation devices, in which a concealed text [z]:=([z₁], . . . , [z_n]) of an approximate value z:=(z₁, . . . , z_n) for a function value y (y₁, . . . , y_n) satisfying y_j=f(x_j) is calculated for each integer j, the integer j being not less than 1 and not more than n, by using a concealed text [x]:=([x₁], . . . , [x_n]) of n pieces of values x:=(x₁, . . . , x_n) as input when n is defined as an integer being not less than 1. m is defined as an integer being not less than 2, g is defined as a polynomial for approximating each section of m sections into which the function f is divided, i is defined as each integer being not less than 1 and not more than m, R_i is defined as the section, and p_i is defined as a parameter of the polynomial g corresponding to the section R_i. Each of the secure computation devices comprises: a parameter acquisition unit that acquires a concealed text [a]:=([a₁], . . . , [a_n]) of a parameter a:=(a₁, . . . , a_n) corresponding to the value x:=(x₁, . . . , x_n) for each integer j, where a_j is defined as the parameter p_i corresponding to the section R_i including the value x_j; and a polynomial calculation unit that calculates the polynomial g([x], [a]) by using the concealed text [x] of the value x as input based on the concealed text [a] of the parameter a so as to obtain the concealed text [z] of the approximate value z for the function value y.

Effects of the Invention

According to the secret batch approximation technique of the present invention, in calculation of a complicated function through secure computation, approximation can be performed by a lower-order polynomial than that of related art and calculation time can be accordingly reduced without degrading approximation accuracy.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram for explaining an approximation method of related art.

FIG. 2 is a diagram for explaining a processing flow of the related art.

FIG. 3 is a diagram for explaining an approximation method of the present invention.

FIG. 4 is a diagram for explaining a processing flow of the present invention.

FIG. 5 is a diagram illustrating a functional configuration of a secret batch approximation system.

FIG. 6 is a diagram illustrating a functional configuration of a secure computation device.

FIG. 7 is a diagram illustrating a processing procedure of a secret batch approximation method.

DETAILED DESCRIPTION OF THE EMBODIMENTS

Notation and definitions of terms in this specification are first described.

<Notation>

A value obtained by concealing a certain value a through encryption, secret sharing, or the like is referred to as a concealed text of a and expressed as [a]. Further, a is referred to as a plaintext of [a]. When concealment is performed through secret sharing, a set of shares of secret sharing, which are held by each secure computation device, is referred to by [a].

[a, b] (square bracket) in a domain of definition of a variable denotes a closed section and (a, b) (parentheses) denotes an open section. For example, i∈[a, b] represents that i is a value which is not less than a and not more than b. Further, i∈[a, b) represents that i is a value which is not less than a and less than b.

<Addition, Subtraction, Multiplication>

In each operation of addition, subtraction, and multiplication with respect to concealed texts, concealed texts [c₁], [c₂], and [c₃] of respective calculation results c₁, c₂, and c₃ of a+b, a−b, and ab are calculated by using concealed texts [a] and [b] of two values a and b as input. Execution of these operations are respectively expressed as the following formulas.

[c₁]←Add([a],[b])

[c₂]←Sub([a],[b])

[c₃]←Mul([a],[b])

When there is no possibility of misunderstanding, Add([a], [b]), Sub([a], [b]), and Mul([a], [b]) are respectively abbreviated as [a]+[b], [a]−[b], and [a]×[b].

<Batch Mapping>

In an operation of batch mapping, concealed texts ([y₁], . . . , [y_n]) of n pieces of values (y₁, . . . , y_n) (here, it is assumed that each y_i satisfies y_i=a_j where j satisfies x_i<u_j when j=1 and u_j-1≤x_i<u_j otherwise) are calculated by using concealed texts ([x₁], . . . , [x_n]) of n pieces of values (x₁, . . . , x_n) (here, it is assumed that each x_j satisfies x_j<u_m for later-described u_m) and tuples of m pieces of values (u₁, . . . , u_m) (here, u_i<u_i+1) and m pieces of values (a₁, . . . , a_m) as input. Execution of this operation is expressed as the following formula.

([y₁], . . . ,[y_n])←BatchMap(([x₁], . . . ,[x_n]),(u₁, . . . ,u_m),(a₁, . . . ,a_m))

Reference Literatures 1 and 2 below describe a secure computation batch mapping algorithm for efficiently calculating batch mapping through secure computation.

• [Reference Literature 1] Koki Hamada, Dai Ikarashi, Koji Chida, “A Batch Mapping Algorithm for Secure Function Evaluation”, The transactions of the Institute of Electronics, Information and Communication Engineers. A, Vol. J96-A, No. 4, pp. 157-165, 2013
• [Reference Literature 2] Peeter Laud, “Parallel Oblivious Array Access for Secure Multiparty Computation and Privacy-Preserving Minimum Spanning Trees”, PoPETs 2015(2), pp. 188-205, 2015.

General Outline of Invention

The present invention utilizes the fact that even a complicated function can be approximated with sufficient accuracy even by a low-order polynomial if sections of the function are limited. Specifically, a function is divided into a plurality of sections, each of the sections is approximated by a lower-order polynomial, and the whole function is thus approximated by a plurality of low-order polynomials. In practical calculation, a section in which an input is included is specified and a low-order polynomial is calculated by using a parameter corresponding to the section. Accordingly, approximation exhibiting accuracy equivalent to accuracy of related art can be realized by using approximate formulas of a lower order when a number of calculation is performed in batch. Since calculation time of a polynomial is proportional to an order of the polynomial, entire calculation time can be reduced without degrading approximation accuracy.

FIG. 1 illustrates a specific example of calculating a complicated function through polynomial approximation according to related art. A function f(x) which is a calculation target has five peaks on top and bottom, so that at least the sixth-order polynomial is required to apply approximation by a polynomial g(x, a). In order to obtain sufficient accuracy, approximation is performed with the use of a higher-order polynomial in practice. For example, a parameter a of the polynomial g(x, a) is a:=(b, c, d, e, f, g, h) and the polynomial g(x, a) can be defined as g(x, a)=bx⁶+cx⁵+dx⁴+ex³+fx²+gx+h. In this case, the number of times of calculation required for obtaining one solution is 11 times for multiplication and six times for addition-and-subtraction.

FIG. 2 illustrates a processing flow of the related art. The parameter a of the polynomial g(x, a) for approximating the function f(x) is first calculated. Then, a polynomial [z_j]=g([x_j], a) is calculated for each integer j which is not less than 1 and not more than n by using concealed texts [x]:=([x₁], . . . , [x_n]) of input values x:=(x₁, . . . , x_n) and the parameter a so as to obtain concealed texts [z] ([z₁], . . . , [z_n]) of approximate values z:=(z₁, . . . , z_n) of function values y (y₁, . . . , y_n).

FIG. 3 illustrates a specific example for calculating a complicated function through polynomial approximation according to the present invention. A domain of definition of a function f(x) which is a calculation target is divided into m pieces of sections and each of the sections is approximated by a quadratic polynomial (that is, a quadratic function). Here, a parameter a of a polynomial g(x, a) is a:=(b, c, d) and the polynomial g(x, a) can be defined as g(x, a)=bx²+cx+d. In this case, the number of times of calculation required for obtaining one solution can be reduced to three times for multiplication and twice for addition-and-subtraction.

FIG. 4 illustrates a processing flow of the present invention. A parameter p_i of the polynomial g(x, a) for approximating each section R_i (i=1, . . . , m) of the function f(x) is first calculated. Then, concealed texts [a]:=([a₁], . . . , [a_n]) of n pieces of parameters are calculated by using the concealed texts [x]:=([x₁], . . . , [x_n]) of input values x:=(x₁, . . . , x_n) and tuples of a section (R₁, . . . , R_m) and a parameter (p₁, . . . , p_m). Here, [a_j] is a concealed text of the parameter p_i corresponding to the section R_i to which x_j belongs. At the end, polynomials [z_j]=g([x_j], [a_j]) are calculated for respective integers j which are not less than 1 and not more than n by using the concealed texts [x]:=([x₁], . . . , [x_n]) of the input values x:=(x₁, . . . , x_n) and the concealed texts [a]:=([a₁], . . . , [a_n]) of the parameters a:=(a₁, . . . , a_n) so as to obtain concealed texts [z]:=([z₁], . . . , [z_n]) of approximate values z:=(z₁, . . . , z_n) for the function values y (y₁, . . . , y_n).

Embodiments according to the present invention are described in detail below. It is to be noted that components mutually having the same function are identified with the same reference character in the drawings and duplicate description thereof are omitted.

First Embodiment

A configuration example of a secret batch approximation system 100 according to a first embodiment will be described with reference to FIG. 5. The secret batch approximation system 100 includes K (≥2) pieces of secure computation devices 1₁, . . . , 1_k. In the present embodiment, the secure computation devices 1₁, . . . , 1_K are each connected to a communication network 9. The communication network 9 is a circuit-switched or packet-switched communication network configured to allow mutual communication of devices connected thereto, and Internet, local area network (LAN), and wide area network (WAN), for example, can be used as the communication network 9. Here, each device is not necessarily communicable via the communication network 9 online. For example, the configuration may be employed in which information to be inputted into the secure computation devices 1₁, . . . , 1_K is stored in a portable recording medium such as a magnetic tape and a USB memory and the information is inputted into the secure computation devices 1₁, . . . , 1_K from the portable recording medium offline.

A configuration example of the secure computation device 1_k (k=1, . . . , K) included in the secret batch approximation system 100 according to the present embodiment will be described with reference to FIG. 6. The secure computation device 1_k includes a storage 10, an input unit 11, a parameter acquisition unit 12, a polynomial calculation unit 13, and an output unit 14, for example, as illustrated in FIG. 6. The secure computation device 1_k (k=1, . . . , K) performs processing of each step described later in cooperation with other secure computation devices 1_k′ (k′=1, . . . , K, where k≠k′), realizing the secret batch approximation method according to the present embodiment.

The secure computation device 1_k is a special device configured in a manner that a special program is read in a known or dedicated computer including a central processing unit (CPU) and a main storage unit (random access memory: RAM), for example. The secure computation device 1_k executes each processing under the control of the central processing unit, for example. Data inputted into the secure computation device 1_k and data obtained through each processing are stored, for example, in the main storage unit and the data stored in the main storage unit is read onto the central processing unit as needed and used for other processing. At least part of processing units of the secure computation device 1_k may be composed of hardware such as an integrated circuit. Each storage included in the secure computation device 1_k may be composed of a main storage unit such as a random access memory (RAM), an auxiliary storage unit composed of a hard disk, an optical disk, or a semiconductor memory element such as a flash memory, or middleware such as relational database and a key-value store, for example.

A processing procedure of the secret batch approximation method executed by the secret batch approximation system 100 according to the present embodiment will be described with reference to FIG. 7.

The polynomial g(x, a) and a sequence ((R₁, p₁), . . . , (R_m, p_m)) of tuples of a section and a parameter are stored in the storage 10 of each secure computation device 1_k. The polynomial g(x, a) is a polynomial for approximating each of sections obtained by dividing the function f(x), which is a calculation target, into predetermined m sections and is a polynomial of lower order than a polynomial for approximating the whole of the function f(x). a denotes a parameter for defining the polynomial g(x, a) and is an array of coefficients of respective terms, for example. R_i (i=1, . . . , m) denotes information indicating each of the sections obtained by dividing the function f(x) into m sections. p_i (i=1, . . . , m) denotes a parameter for approximating the section R_i of the function f(x) by the polynomial g(x, a).

In step S1, the input unit 11 of each secure computation device 1_k receives the concealed texts [x]:=([x₁], . . . , [x_n]) of n pieces of values x:=(x₁, . . . , x_n) being calculation targets, as input. The input unit 11 outputs the concealed texts [x] of the values x to the parameter acquisition unit 12.

In step S2, the parameter acquisition unit 12 of each secure computation device 1_k receives the concealed texts [x] of the values x from the input unit 11, and acquires the concealed texts [a]:=([a₁], . . . , [a_n]) of n pieces of parameters a:=(a₁, . . . , a_n) satisfying ∀j∈[1, n] and ∃i s.t. a_j=p_i, x_j∈R_i from the sequence ((R₁, p₁), . . . , (R_m, p_m)) of tuples of a section and a parameter which are stored in the storage 10, for respective integers j which are not less than 1 and not more than n. That is, the concealed texts [a₁], . . . , [a_n] of the parameters a₁, . . . , a_n respectively corresponding to the values x₁, . . . , x_n are generated for respective integers j which are not less than 1 and not more than n, where a_j is defined as the parameter p_i of the section R_i corresponding to the value x_j. The parameter acquisition unit 12 outputs the concealed texts [x] of the values x and the concealed texts [a] of the parameters a to the polynomial calculation unit 13.

In step S3, the polynomial calculation unit 13 of each secure computation device 1_k receives the concealed texts [x] of the values x and the concealed texts [a] of the parameters a from the parameter acquisition unit 12 and calculates [z_j]=g([x_j], [a_j]) in accordance with the polynomial g(x, a) stored in the storage 10 for respective integers j which are not less than 1 and not more than n. The polynomial calculation unit 13 outputs the concealed texts [z]:=([z₁], . . . , [z_n]) of the approximate values z:=(z₁, . . . , z_n) for the function values y (y₁, . . . , y_n) to the output unit 14.

In step S4, the output unit 14 of each secure computation device 1_k receives the concealed texts [z] of the approximate values z for the function values y from the polynomial calculation unit 13 and sets the concealed texts [z] as output of the secure computation device 1_k.

Second Embodiment

A second embodiment describes a more specific example in which the following sigmoid function, a domain of definition of which is expressed as X:=[0, 1000), is approximated by a quadratic function.

$f\left(x\right)=\sigma \left(x\right)=\frac{1}{1+e^{-x}}$

Differences from the first embodiment will be mainly described below.

A polynomial g(x, a) and a sequence ((R₁, p₁), . . . , (R_m, p_m)) of tuples of a section and a parameter are stored in the storage 10 of each secure computation device 1_k. A parameter of the polynomial g(x, a) of the present embodiment is a:=(b, c, d) and the polynomial g(x, a) is defined as g(x, a)=bx²+cx+d. The section R_i of the present embodiment is defined as R_i:=[1_i, u_i) for i∈[1, m] (here, 1_i=0, u_m=1000, u_i=1_i+i, and 1_i≤u_i for i∈[1, m)). The parameter p_i of the present embodiment is defined as p_i:=(b_i, c_i, d_i).

The parameter acquisition unit 12 of each secure computation device 1_k according to the present embodiment executes a secure computation batch mapping algorithm as the following formula by using the concealed texts [x], (u₁, . . . , u_m), and ((b₁, c₁, d₁), . . . , (b_m, c_m, d_m)) as input, acquiring concealed texts [a]:=([a₁], . . . , [a_n]) of the parameters a (a₁, . . . , a_n).

([a₁], . . . ,[a_n])←BatchMap([x],(u₁, . . . ,u_m),((b₁,c₁,d₁), . . . ,(b_m,c_m,d_m)))

Here, a_j satisfies a_j=(b_i, c_i, d_i) for certain i satisfying x_j∈R_i.

The polynomial calculation unit 13 of each secure computation device 1_k according to the present embodiment calculates the following formula based on [a_j]:=([b′_j], [c′_j], [d′_j]) for each integer j which is not less than 1 and not more than n, obtaining a concealed text [z_j] of an approximate value z_j.

[z_j]←[b′_i]×[x_j]×[x_j]+[c′_i]×[x_j]+[d′_i]

The present invention utilizes the fact that even a complicated function can be approximated with sufficient accuracy even by a low-order polynomial if sections of the function are limited. Specifically, a function is divided into the predetermined number of sections to obtain parameters for approximating respective sections by low-order polynomials in advance, and low-order polynomial approximation is performed through acquiring a parameter corresponding to an input value depending on the section including the input value. Accordingly, approximation exhibiting accuracy equivalent to accuracy of related art can be realized by using lower-order polynomials when a number of calculation is performed in batch. Especially, polynomial approximation can be more efficiently performed through acquiring parameters corresponding to respective input values of a plurality of input values by using a secure computation batch mapping algorithm for the input values. For example, it is reported that at least tenth order is experimentally required so as to obtain sufficient accuracy in logistic regression for a sigmoid function which is heavily required in each iterative calculation in learning of logistic regression (see Reference Literature 3). However, if the number of sections is set to 5000, for example, approximation can be realized by a cubic polynomial with almost double accuracy and the number of times of multiplication and addition for real numbers can be reduced by 70 percent.

• [Reference Literature 3] Payman Mohassel and Yupeng Zhang, “SecureML: A system for scalable privacy-preserving machine learning,” In 2017 IEEE Symposium on Security and Privacy, S P 2017, San Jose, Calif., USA, May 22-26, 2017, pp. 19-38. IEEE Computer Society, 2017.

While the embodiments of the present invention have been described, specific configurations are not limited to these embodiments, but design modifications and the like within a range not departing from the spirit of the invention are encompassed in the scope of the invention, of course. The various processes described in the embodiments may be executed in parallel or separately depending on the processing ability of a device executing the process or on any necessity, rather than being executed in time series in accordance with the described order.

[Program and Recording Medium]

When various types of processing functions in the devices described in the above embodiments are implemented on a computer, the contents of processing function to be contained in each device is written by a program. With this program executed on the computer, various types of processing functions in the above-described devices are implemented on the computer.

This program in which the contents of processing are written can be recorded in a computer-readable recording medium. The computer-readable recording medium may be any medium such as a magnetic recording device, an optical disk, a magneto-optical recording medium, and a semiconductor memory.

Distribution of this program is implemented by sales, transfer, rental, and other transactions of a portable recording medium such as a DVD and a CD-ROM on which the program is recorded, for example. Furthermore, this program may be stored in a storage unit of a server computer and transferred from the server computer to other computers via a network so as to be distributed.

A computer which executes such program first stores the program recorded in a portable recording medium or transferred from a server computer once in a storage unit thereof, for example. When the processing is performed, the computer reads out the program stored in the storage unit thereof and performs processing in accordance with the program thus read out. As another execution form of this program, the computer may directly read out the program from a portable recording medium and perform processing in accordance with the program. Furthermore, each time the program is transferred to the computer from the server computer, the computer may sequentially perform processing in accordance with the received program. Alternatively, a configuration may be adopted in which the transfer of a program to the computer from the server computer is not performed and the above-described processing is executed by so-called application service provider (ASP)-type service by which the processing functions are implemented only by an instruction for execution thereof and result acquisition. It should be noted that a program in this form includes information which is provided for processing performed by electronic calculation equipment and which is equivalent to a program (such as data which is not a direct instruction to the computer but has a property specifying the processing performed by the computer).

In this form, the present device is configured with a predetermined program executed on a computer. However, the present device may be configured with at least part of these processing contents realized in a hardware manner.

Claims

1. A secret batch approximation system including a plurality of secure computation devices for calculating a concealed text [z]:=([z1], . . . , [zn]) of an approximate value z:=(z1, . . . , zn) for a function value y:=(y1, . . . , yn) satisfying yj=f(xj) for each integer j, the integer j being not less than 1 and not more than n, by using a concealed text [x]:=([x1], . . . , [xn]) of n pieces of values x:=(x1, . . . , xn) as input when n is defined as an integer being not less than 1, wherein m is defined as an integer being not less than 2, g is defined as a polynomial for approximating each section of m sections into which the function f is divided, i is defined as each integer being not less than 1 and not more than m, Ri is defined as the section, and pi is defined as a parameter of the polynomial g corresponding to the section Ri,each of the secure computation devices comprises processing circuitry configured to: acquire a concealed text [a]:=([a1], . . . , [an]) of a parameter a:=(a1, . . . , an) corresponding to the value x:=(x1, . . . , xn) for each integer j, where aj is defined as the parameter pi corresponding to the section Ri including the value xj; andcalculate the polynomial g([x], [a]) by using the concealed text [x] of the value x as input based on the concealed text [a] of the parameter a so as to obtain the concealed text [z] of the approximate value z for the function value y.

2. The secret batch approximation system according to claim 1, wherein the processing circuitry acquires the concealed text [a]:=([a1], . . . , [an]) of the parameter a:=(a1, . . . , an) based on a formula: ([a1], . . . ,[an])←BatchMap([x],(u1, . . . ,um),(p1, . . . ,pm))

3. The secret batch approximation system according to claim 2, wherein the function f is a sigmoid function with a domain of definition expressed as X:=[0, 1000), m≤5000 holds, and the polynomial g is a cubic polynomial.

4. A secure computation device that is included in a secret batch approximation system, the secret batch approximation system calculating a concealed text [z]:=([z1], . . . , [zn]) of an approximate value z:=(z1, . . . , zn) for a function value y:=(y1, . . . , yn) satisfying yj=f(xj) for each integer j, the integer j being not less than 1 and not more than n, by using a concealed text [x] ([x1], . . . , [xn]) of n pieces of values x:=(x1, . . . , xn) as input when n is defined as an integer being not less than 1, wherein m is defined as an integer being not less than 2, g is defined as a polynomial for approximating each section of m sections into which the function f is divided, i is defined as each integer being not less than 1 and not more than m, Ri is defined as the section, and pi is defined as a parameter of the polynomial g corresponding to the section Ri,the secure computation device comprises processing circuitry configured to: acquire a concealed text [a]:=([a1], . . . , [an]) of a parameter a:=(a1, . . . , an) corresponding to the value x:=(x1, . . . , xn) for each integer j, where aj is defined as the parameter pi corresponding to the section Ri including the value xj; andcalculate the polynomial g([x], [a]) by using the concealed text [x] of the value x as input based on the concealed text [a] of the parameter a so as to obtain the concealed text [z] of the approximate value z for the function value y.

5. A secret batch approximation method that is executed by a secret batch approximation system, the secret batch approximation system including a plurality of secure computation devices for calculating a concealed text [z]: ([z1], . . . , [zn]) of an approximate value z:=(z1, . . . , zn) for a function value y:=(y1, . . . , yn) satisfying yj=f(xj) for each integer j, the integer j being not less than 1 and not more than n, by using a concealed text [x]:=([x1], . . . , [xn]) of n pieces of values x:=(x1, . . . , xn) as input when n is defined as an integer being not less than 1, wherein m is defined as an integer being not less than 2, g is defined as a polynomial for approximating each section of m sections into which the function f is divided, i is defined as each integer being not less than 1 and not more than m, Ri is defined as the section, and pi is defined as a parameter of the polynomial g corresponding to the section Ri,the secret batch approximation method comprising:acquiring, by processing circuitry of each of the secure computation devices, a concealed text [a]:=([a1], . . . , [an]) of a parameter a:=(a1, . . . , an) corresponding to the value x:=(x1, . . . , xn) for each integer j, where aj is defined as the parameter pi corresponding to the section Ri including the value xj, andcalculating, by the processing circuitry of each of the secure computation devices, the polynomial g([x], [a]) by using the concealed text [x] of the value x as input based on the concealed text [a] of the parameter a so as to obtain the concealed text [z] of the approximate value z for the function value y.

6. A non-transitory computer readable medium having a program recorded thereon for making a computer function as the secure computation device according to claim 4.