NAF(Non-adjacent form) w-NAF及其在curve25519-dalek中scalar的实现

2023-11-01 12:20

本文主要是介绍NAF(Non-adjacent form) w-NAF及其在curve25519-dalek中scalar的实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

1. 引言

The non-adjacent form (NAF) of a number is a unique signed-digit representation. Like the name suggests, non-zero values cannot be adjacent. For example:
( 0 1 1 1 ) 2 = 4 + 2 + 1 = 7 (0\ 1\ 1\ 1)_2 = 4 + 2 + 1 = 7 (0 1 1 1)2=4+2+1=7
( 1 0 − 1 1 ) 2 = 8 − 2 + 1 = 7 (1\ 0\ −1\ 1)_2 = 8 − 2 + 1 = 7 (1 0 1 1)2=82+1=7
( 1 − 1 1 1 ) 2 = 8 − 4 + 2 + 1 = 7 (1\ −1\ 1\ 1)_2 = 8 − 4 + 2 + 1 = 7 (1 1 1 1)2=84+2+1=7
( 1 0 0 − 1 ) 2 = 8 − 1 = 7 (1\ 0\ 0\ −1)_2 = 8 − 1 = 7 (1 0 0 1)2=81=7

All are valid signed-digit representations of 7, but only the final representation, ( 1 0 0 − 1 ) 2 (1\ 0\ 0\ −1)_2 (1 0 0 1)2, is in NAF.

NAF即以一组有符号数字表示,且非零值不可相邻(即每个非零值的左右相邻位必须均为0)。NAF表示的数字,可保证Hamming weight值最小,且相比于普通的二进制表示,其非零值比率可控制在1/3以内。

2. NAF的优势

因以NAF表示的数字相比于以二进制形式表示的数字,其非零值的个数有效减少了(由1/2减为1/3)。非零值个数的减少,将提高某些算法的效率。比如密码学中应用中,可减少exponentiation幂算法中的乘法运算数量(该数量取决于非零值的位数)。
exponentiation幂运算通常表示由基数 b b b和指数 n 组 成 n组成 n
b n = b × b . . . × b b^n=b\times b...\times b bn=b×b...×b

NAF中的1即表示与基数 b b b的一次乘积运算,-1表示与基数的倒数 1 b \frac{1}{b} b1的乘法运算。

3. 转NAF算法

将普通二进制数字表示转换为NAF表示的算法如下:

 Input     E = (em − 1 em − 2 ··· e1 e0)2Output     Z = (zm zm − 1 ··· z1 z0)NAFi ← 0while E > 0 doif E is odd thenzi ← 2 − (E mod 4)E ← E − zielsezi ← 0E ← E/2i ← i + 1return z

特别的,对于以 w w w进制表示的NAF,通称为width-w NAF。具体定义见 书《Guide to Elliptic Curve Cryptography》中Definition 3.32:
在这里插入图片描述
具体的算法实现为:
在这里插入图片描述
上面算法中2.1步骤中的 k m o d s 2 w k\ mods\ 2^w k mods 2w,对应的结果区间在 [ − 2 w − 1 , 2 w − 1 − 1 ] [-2^{w-1}, 2^{w-1}-1] [2w1,2w11]

举例为,下面例子中上面有横杠的数字表示的即为负数:
在这里插入图片描述
观察可发现: w w w值越大,NAF中非零值的个数越少。

3. curve25519-dalek中scalar的NAF

以NAF形式表示的位数最多比以二进制形式表示的位数多一位。因此,需保证scalar以二进制表示时,其最高位保持为0值,而相应的NAF表示时,最多只有256位。

针对算法:在这里插入图片描述
上面算法中2.1步骤中的 k m o d s 2 w k\ mods\ 2^w k mods 2w相当于求 u u u,使得 u ≡ k ( m o d 2 w ) u\equiv k(mod\ 2^w) uk(mod 2w) u u u对应的结果区间在 [ − 2 w − 1 , 2 w − 1 ) [-2^{w-1}, 2^{w-1}) [2w1,2w1)

在curve25519-dalek的实际代码实现中,具体的思路为已知 k = ( k m k m − 1 . . . k w + 1 k w k w − 1 . . . k 1 k 0 ) 2 k=(k_mk_{m-1}...k_{w+1}k_wk_{w-1}...k_1k_0)_2 k=(kmkm1...kw+1kwkw1...k1k0)2,求相应的 N A F w ( k ) = ( n m . . . . n 1 n 0 ) w NAF_w(k)=(n_m....n_1n_0)_w NAFw(k)=(nm....n1n0)w
1)将k以二进制表示为 k = ( k m k m − 1 . . . k w + 1 k w k w − 1 . . . k 1 k 0 ) 2 = ∑ i = 0 w − 1 k i 2 i + 2 w ∗ ∑ i = 0 k i + w 2 i k=(k_mk_{m-1}...k_{w+1}k_wk_{w-1}...k_1k_0)_2=\sum_{i=0}^{w-1}k_i2^i+2^w*\sum_{i=0}k_{i+w}2^i k=(kmkm1...kw+1kwkw1...k1k0)2=i=0w1ki2i+2wi=0ki+w2i,其中 ∑ i = 0 w − 1 k i 2 i ≡ k m o d 2 w \sum_{i=0}^{w-1}k_i2^i\equiv k\ mod\ 2^w i=0w1ki2ik mod 2w
2)若 k m o d 2 w k\ mod\ 2^w k mod 2w为奇数,且 k m o d 2 w < 2 w − 1 k\ mod\ 2^w < 2^{w-1} k mod 2w<2w1时, n 0 = k m o d 2 w n_0=k\ mod\ 2^w n0=k mod 2w,对于其中的 k = k − n 0 k=k-n_0 k=kn0计算,其实即为 k = k > > w k=k>>w k=k>>w
3)若 k m o d 2 w k\ mod\ 2^w k mod 2w为奇数,且 k m o d 2 w ≥ 2 w − 1 k\ mod\ 2^w \ge 2^{w-1} k mod 2w2w1时, n 0 = k m o d 2 w − 2 w n_0=k\ mod\ 2^w-2^w n0=k mod 2w2w,对于其中的 k = k − n 0 k=k-n_0 k=kn0计算,其实即为 k = k − n 0 = ∑ i = 0 w − 1 k i 2 i + 2 w ∗ ∑ i = 0 k i + 2 2 i − ( k m o d 2 w − 2 w ) ≡ k m o d 2 w + 2 w ∗ ∑ i = 0 k i + w 2 i − ( k m o d 2 w − 2 w ) ≡ 2 w + 2 w ∗ ∑ i = 0 k i + w 2 i k=k-n_0=\sum_{i=0}^{w-1}k_i2^i+2^w*\sum_{i=0}k_{i+2}2^i-(k\ mod\ 2^w-2^w)\equiv k\ mod\ 2^w+2^w*\sum_{i=0}k_{i+w}2^i-(k\ mod\ 2^w-2^w)\equiv 2^w+2^w*\sum_{i=0}k_{i+w}2^i k=kn0=i=0w1ki2i+2wi=0ki+22i(k mod 2w2w)k mod 2w+2wi=0ki+w2i(k mod 2w2w)2w+2wi=0ki+w2i,最终可有 k = k > > w , c a r r y = 1 k=k>>w, carry=1 k=k>>w,carry=1
4)若 k m o d 2 w k\ mod\ 2^w k mod 2w为偶数,则有 n 0 = 0 n_0=0 n0=0 k = k > > 1 k=k>>1 k=k>>1

具体代码实现为:

    pub(crate) fn non_adjacent_form(&self, w: usize) -> [i8; 256] {// required by the NAF definitiondebug_assert!( w >= 2 );// required so that the NAF digits fit in i8debug_assert!( w <= 8 );use byteorder::{ByteOrder, LittleEndian};let mut naf = [0i8; 256];let mut x_u64 = [0u64; 5];LittleEndian::read_u64_into(&self.bytes, &mut x_u64[0..4]);let width = 1 << w;let window_mask = width - 1;let mut pos = 0;let mut carry = 0;while pos < 256 {// Construct a buffer of bits of the scalar, starting at bit `pos`let u64_idx = pos / 64;let bit_idx = pos % 64;let bit_buf: u64;if bit_idx < 64 - w {// This window's bits are contained in a single u64bit_buf = x_u64[u64_idx] >> bit_idx;} else {// Combine the current u64's bits with the bits from the next u64bit_buf = (x_u64[u64_idx] >> bit_idx) | (x_u64[1+u64_idx] << (64 - bit_idx));}// Add the carry into the current windowlet window = carry + (bit_buf & window_mask);if window & 1 == 0 {// If the window value is even, preserve the carry and continue.// Why is the carry preserved?// If carry == 0 and window & 1 == 0, then the next carry should be 0// If carry == 1 and window & 1 == 0, then bit_buf & 1 == 1 so the next carry should be 1pos += 1;continue;}if window < width/2 {carry = 0;naf[pos] = window as i8;} else {carry = 1;naf[pos] = (window as i8).wrapping_sub(width as i8);}pos += w;}naf}

参考资料:
[1] https://en.wikipedia.org/wiki/Non-adjacent_form
[2] 《Guide to Elliptic Curve Cryptography》

这篇关于NAF(Non-adjacent form) w-NAF及其在curve25519-dalek中scalar的实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/322937

相关文章

Oracle查询优化之高效实现仅查询前10条记录的方法与实践

《Oracle查询优化之高效实现仅查询前10条记录的方法与实践》:本文主要介绍Oracle查询优化之高效实现仅查询前10条记录的相关资料,包括使用ROWNUM、ROW_NUMBER()函数、FET... 目录1. 使用 ROWNUM 查询2. 使用 ROW_NUMBER() 函数3. 使用 FETCH FI

Python脚本实现自动删除C盘临时文件夹

《Python脚本实现自动删除C盘临时文件夹》在日常使用电脑的过程中,临时文件夹往往会积累大量的无用数据,占用宝贵的磁盘空间,下面我们就来看看Python如何通过脚本实现自动删除C盘临时文件夹吧... 目录一、准备工作二、python脚本编写三、脚本解析四、运行脚本五、案例演示六、注意事项七、总结在日常使用

Java实现Excel与HTML互转

《Java实现Excel与HTML互转》Excel是一种电子表格格式,而HTM则是一种用于创建网页的标记语言,虽然两者在用途上存在差异,但有时我们需要将数据从一种格式转换为另一种格式,下面我们就来看看... Excel是一种电子表格格式,广泛用于数据处理和分析,而HTM则是一种用于创建网页的标记语言。虽然两

Java中Springboot集成Kafka实现消息发送和接收功能

《Java中Springboot集成Kafka实现消息发送和接收功能》Kafka是一个高吞吐量的分布式发布-订阅消息系统,主要用于处理大规模数据流,它由生产者、消费者、主题、分区和代理等组件构成,Ka... 目录一、Kafka 简介二、Kafka 功能三、POM依赖四、配置文件五、生产者六、消费者一、Kaf

使用Python实现在Word中添加或删除超链接

《使用Python实现在Word中添加或删除超链接》在Word文档中,超链接是一种将文本或图像连接到其他文档、网页或同一文档中不同部分的功能,本文将为大家介绍一下Python如何实现在Word中添加或... 在Word文档中,超链接是一种将文本或图像连接到其他文档、网页或同一文档中不同部分的功能。通过添加超

windos server2022里的DFS配置的实现

《windosserver2022里的DFS配置的实现》DFS是WindowsServer操作系统提供的一种功能,用于在多台服务器上集中管理共享文件夹和文件的分布式存储解决方案,本文就来介绍一下wi... 目录什么是DFS?优势:应用场景:DFS配置步骤什么是DFS?DFS指的是分布式文件系统(Distr

NFS实现多服务器文件的共享的方法步骤

《NFS实现多服务器文件的共享的方法步骤》NFS允许网络中的计算机之间共享资源,客户端可以透明地读写远端NFS服务器上的文件,本文就来介绍一下NFS实现多服务器文件的共享的方法步骤,感兴趣的可以了解一... 目录一、简介二、部署1、准备1、服务端和客户端:安装nfs-utils2、服务端:创建共享目录3、服

C#使用yield关键字实现提升迭代性能与效率

《C#使用yield关键字实现提升迭代性能与效率》yield关键字在C#中简化了数据迭代的方式,实现了按需生成数据,自动维护迭代状态,本文主要来聊聊如何使用yield关键字实现提升迭代性能与效率,感兴... 目录前言传统迭代和yield迭代方式对比yield延迟加载按需获取数据yield break显式示迭

Python实现高效地读写大型文件

《Python实现高效地读写大型文件》Python如何读写的是大型文件,有没有什么方法来提高效率呢,这篇文章就来和大家聊聊如何在Python中高效地读写大型文件,需要的可以了解下... 目录一、逐行读取大型文件二、分块读取大型文件三、使用 mmap 模块进行内存映射文件操作(适用于大文件)四、使用 pand

python实现pdf转word和excel的示例代码

《python实现pdf转word和excel的示例代码》本文主要介绍了python实现pdf转word和excel的示例代码,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价... 目录一、引言二、python编程1,PDF转Word2,PDF转Excel三、前端页面效果展示总结一