这篇讲一下如何使用中国剩余定理CRT来对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
d P = e − 1 m o d ( p − 1 ) dP = e^{-1} mod \ (p-1) dP=e−1mod (p−1)
d Q = e − 1 m o d ( q − 1 ) dQ=e^{-1} mod \ (q-1) dQ=e−1mod (q−1)
q I n v = q − 1 m o d p qInv = q^{-1} mod \ p qInv=q−1mod p
e − 1 e^{-1} e−1表示模逆,表达式 x = e − 1 m o d N x=e^{-1} mod \ N x=e−1mod 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) x⋅e≡1(mod N)。 N = n = p q N=n=pq N=n=pq。
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⋅(m1−m2)modp
m = m 2 + h ⋅ q m = m_2 + h \cdot q m=m2+h⋅q
我们把 ( 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)。
欧拉定理是费马小定理(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)是Euler′stotient函数,求小于n的正整数中与n互质的个数
一个质数p, ϕ ( n ) = p − 1 \phi(n)=p-1 ϕ(n)=p−1
我们需要计算 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+h⋅q
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=((x1−x2)(q−1 mod p)) mod p
CRT系数 q I n v = q − 1 m o d p qInv = q^{-1} \ mod \ p qInv=q−1 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使用相同的算法。
d m o d ( p − 1 ) = e − 1 m o d ( p − 1 ) , d \ mod \ (p-1)=e^{-1} \ mod \ (p-1), d mod (p−1)=e−1 mod (p−1),
d m o d ( q − 1 ) = e − 1 m o d ( q − 1 ) . d \ mod \ (q-1) = e^{-1} \ mod \ (q-1). d mod (q−1)=e−1 mod (q−1).
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=e−1 mod (p−1)=d mod (p−1)
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=e−1 mod (q−1)=d mod (q−1)
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=q−1 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⋅(m1−m2) mod p
m = m 2 + h ⋅ q m=m_2+h \cdot q m=m2+h⋅q