使用中国剩余定理CRT对RSA运算进行加速

这篇讲一下如何使用中国剩余定理CRT来对RSA加密运算进行加速。

RSA运算

当我们使用RSA私钥(n,d)对密文c进行解密(或者计算数字签名时),我们需要计算模幂 m = c d m o d   n m=c^d mod \ n m=cdmod n。私钥指数 d d d并不像公钥指数 e e e那样方便。一个k比特的模n,对应的私钥指数d差不多跟它一样长。计算的工作量同长度k成正比,所以对于RSA私钥的运算,有更多的计算量。

我们可以使用CRT模式更有效的计算 m = c d m=c^d m=cd

    1. 使用 p , q , p > q p,q,p \gt q p,qp>q提前计算以下值:

d P = e − 1 m o d   ( p − 1 ) dP = e^{-1} mod \ (p-1) dP=e1mod (p1)
d Q = e − 1 m o d   ( q − 1 ) dQ=e^{-1} mod \ (q-1) dQ=e1mod (q1)
q I n v = q − 1 m o d   p qInv = q^{-1} mod \ p qInv=q1mod p

e − 1 e^{-1} e1表示模逆,表达式 x = e − 1 m o d   N x=e^{-1} mod \ N x=e1mod Ny也会写成 x = ( 1 / e ) m o d   N x=(1/e) mod \ N x=1/e)mod N。x是任意整数满足 x ⋅ e ≡ 1 ( m o d   N ) x \cdot e \equiv 1 (mod \ N) xe1(mod N) N = n = p q N=n=pq N=n=pq

    1. 使用密文c计算明文消息m

m 1 = c d P m o d   p m_1 = c^{dP} mod \ p m1=cdPmod p
m 2 = c d Q m o d   q m_2 = c^{dQ} mod \ q m2=cdQmod q
h = q I n v ⋅ ( m 1 − m 2 ) m o d p h = qInv \cdot (m_1-m_2) mod p h=qInv(m1m2)modp
m = m 2 + h ⋅ q m = m_2 + h \cdot q m=m2+hq

我们把 ( p , q , d P , d Q , q I n v ) (p,q,dP,dQ,qInv) (p,q,dP,dQ,qInv)作为私钥保存。

下面需要了解两个数论的原理,分别是中国剩余定义的一个特殊情况和欧拉定理。

中国剩余定理-特殊情况

中国剩余定理的特殊情况可以表述如下:

$p和q是不相同的素数,n=p \cdot q.对于任意的一对(x_1,x_2),0 \leq x_1 \lt p 且 0 \leq x_2 \lt q,存在唯一的数x,0 \leq x \lt n $
x 1 = x   m o d   p , 且 x 2 = x   m o d   q x_1=x \ mod \ p, 且 x_2 =x \ mod \ q x1=x mod p,x2=x mod q
所以任意整数x都可以使用CRT表示方法唯一的表示成 ( x 1 , x 2 ) (x_1,x_2) (x1,x2)

欧拉定理 Euler’s Theorem

欧拉定理是费马小定理(Fermat’s Little Theorem)的推广,也称作欧拉-费马定理(Euler-Fermat Theorem)。

如果n是一个正整数,a是任意整数,且 g c d ( a , n ) = 1 gcd(a,n)=1 gcd(a,n)=1,那么 a ϕ ( n ) ≡ 1   ( m o d   n ) , ϕ ( n ) 是 E u l e r ′ s t o t i e n t 函 数 , 求 小 于 n 的 正 整 数 中 与 n 互 质 的 个 数 a^{\phi(n)} \equiv 1 \ (mod \ n),\phi(n)是Euler's totient函数,求小于n的正整数中与n互质的个数 aϕ(n)1 (mod n),ϕ(n)Eulerstotientnn

一个质数p, ϕ ( n ) = p − 1 \phi(n)=p-1 ϕ(n)=p1

CRT表示法中的运算

我们需要计算 m = c d   m o d   n m=c^d \ mod \ n m=cd mod n。如果我们知道 ( c d   m o d   p , c d   m o d q ) (c^d \ mod \ p, c^d \ mod q) (cd mod p,cd modq)那么CRT告诉我们存在唯一的值 c d   m o d   n c^d \ mod \ n cd mod n在范围[0,n-1]。

使用CRT表示方法 ( x 1 , x 2 ) (x_1,x_2) (x1,x2)恢复出x,我们使用Garner’s方程式。

x = x 2 + h ⋅ q x=x_2 + h \cdot q x=x2+hq
h = ( ( x 1 − x 2 ) ( q − 1   m o d   p ) )   m o d   p h=((x_1-x_2)(q^{-1} \ mod \ p)) \ mod \ p h=((x1x2)(q1 mod p)) mod p
CRT系数 q I n v = q − 1   m o d   p qInv = q^{-1} \ mod \ p qInv=q1 mod p可以提前计算。模幂的运算量随着模的比特数k的立方增加而增加。所以做两次幂运算mod p和mod q,比做一次幂运算mod n效率要高。

计算 c d   m o d   p c^d \ mod \ p cd mod p,可以使用欧拉定理来减少指数d modulo (p-1):

$c^d \ mod \ p = c^{d \ mod \ \phi§} \ mod \ p = c^{d \ mod \ (p-1)} \ mod \ p $

对于q使用相同的算法。

RSA运算

d   m o d   ( p − 1 ) = e − 1   m o d   ( p − 1 ) , d \ mod \ (p-1)=e^{-1} \ mod \ (p-1), d mod (p1)=e1 mod (p1),
d   m o d   ( q − 1 ) = e − 1   m o d   ( q − 1 ) . d \ mod \ (q-1) = e^{-1} \ mod \ (q-1). d mod (q1)=e1 mod (q1).

d P = e − 1   m o d   ( p − 1 ) = d   m o d   ( p − 1 ) dP = e^{-1} \ mod \ (p-1) = d \ mod \ (p-1) dP=e1 mod (p1)=d mod (p1)
d Q = e − 1   m o d   ( q − 1 ) = d   m o d   ( q − 1 ) dQ = e^{-1} \ mod \ (q-1) = d \ mod \ (q-1) dQ=e1 mod (q1)=d mod (q1)
m 1 = c d P   m o d   p m_1 = c^{dP} \ mod \ p m1=cdP mod p
m 2 = c d Q   m o d   q m_2 = c^{dQ} \ mod \ q m2=cdQ mod q

q I n v = q − 1   m o d   p qInv = q^{-1} \ mod \ p qInv=q1 mod p
h = q I n v ⋅ ( m 1 − m 2 )   m o d   p h=qInv \cdot (m_1-m_2) \ mod \ p h=qInv(m1m2) mod p
m = m 2 + h ⋅ q m=m_2+h \cdot q m=m2+hq

你可能感兴趣的:(密码,智能卡,程序员,RSA,密码学,加密算法,数论)