11.1 介绍
大数一般指的是位数很多的数。计算机表示的数的大小是有限的,精度也是有限的,它不能支持大数运算。密码学中采用了很多大数计算,为了让计算机实现大数运算,用户需要定义自己的大数表示方式并及实现各种大数运算。Openssl为我们提供了这些功能,主要用于非对称算法。
11.2 openssl大数表示
crypto/bn.h中定义了大数的表示方式,如下:
struct bignum_st
{
BN_ULONG *d;
int top;
int dmax;
int neg;
int flags;
};
各项意义如下:
d:BN_ULONG(应系统而异,win32下为4个字节)数组指针首地址,大数就存放在这里面,不过是倒放的。比如,用户要存放的大数为12345678000(通过BN_bin2bn放入),则d的内容如下:0x30 0x30 0x30 0x38 0x37 0x36 0x35 0x34 0x33 0x32 0x31;
top:用来指明大数占多少个BN_ULONG空间,上例中top为3。
dmax:d数组的大小。
neg:是否为负数,如果为1,则是负数,为0,则为正数。
flags:用于存放一些标记,比如flags含有BN_FLG_STATIC_DATA时,表明d的内存是静态分配的;含有BN_FLG_MALLOCED时,d的内存是动态分配的。
11.3 大数函数
大数函数一般都能根据函数名字知道其实现的功能。下面简单介绍了几个函数。
1) BN_rand/BN_pseudo_rand
2) BN_rand_range/BN_pseudo_rand_range
3) BN_dup
大数复制。
4) BN_generate_prime
生成素数。
5) int BN_add_word(BIGNUM *a,BN_ULONG w)
给大数a加上w,如果成功,返回1。
示例:
#include <openssl/bn.h>
int main()
int ret;
BIGNUM *a;
BN_ULONG w;
a=BN_new();
BN_one(a);
w=2685550010;
ret=BN_add_word(a,w);
if(ret!=1)
printf("a+=w err!\n");
BN_free(a);
return -1;
}
BN_free(a);
return 0;
6) BIGNUM *BN_bin2bn(const unsigned char *s,int len,BIGNUM *ret)
将内存中的数据转换为大数,为内存地址,len为数据长度,ret为返回值。
int main()
BIGNUM *ret1,*ret2;
ret1=BN_new();
ret1=BN_bin2bn("242424ab",8,ret1);
ret2=BN_bin2bn("242424ab",NULL);
BN_free(ret1);
BN_free(ret2);
return 0;
注意:输入参数“242424ab”是asc码,对应的大数值为16进制的0x3234323432346162
7) int BN_bn2bin(const BIGNUM *a,unsigned char *to)
将大数转换为内存形式。输入参数为大数a,to为输出缓冲区地址,缓冲区需要预先分配,返回值为缓冲区的长度。
int main()
BIGNUM *ret1=NULL;
char bin[50],*buf=NULL;
int len;
ret1=BN_bin2bn("242424ab",u5b8bu4f53; font-size:16px; text-indent:2em"> len=BN_bn2bin(ret1,bin);
len=BN_num_bytes(ret1);
buf=malloc(len);
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',buf);
free(buf);
BN_free(ret1);
return 0;
本例的缓冲区分配有两种方法:静态分配和动态分配。动态分配时,先调用
8) char *BN_bn2dec(const BIGNUM *a)
将大数转换成整数字符串。返回值中存放整数字符串,它由内部分配空间,用户必须在外部用OPENSSL_free函数释放该空间。
#include <openssl/crypto.h>
BIGNUM *ret1=NULL;
char *p=NULL;
int len=0;
p=BN_bn2dec(ret1);
printf("%s\n",p); /* 3617571600447332706 */
OPENSSL_free(p);
getchar();
9) char *BN_bn2hex(const BIGNUM *a)
将大数转换为十六进制字符串。返回值为生成的十六进制字符串,外部需要用OPENSSL_free函数释放
p=BN_bn2hex(ret1);
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',p);
OPENSSL_free(p);
输出的结果为:323432346162
10) BN_cmp
比较两个大数。
11)BIGNUM *BN_mod_inverse(BIGNUM *in, const BIGNUM *a,
const BIGNUM *n,BN_CTX *ctx)
计算ax=1(mod n)。
用户使用openssl函数编程时,一般用不着进行大数运算。BN_bin2bn、BN_hex2bn、BN_dec2bn、BN_bin2bn、BN_bn2bin、BN_bn2hex和BN_bn2dec比较常用。比如给定RSA密钥的内存形式,用户可以调用BN_bin2bn来构造RSA密钥的大数元素来进行RSA运算,或者已经生成了RSA密钥,用户调用BN_bn2bin将RSA各个元素导出到内存中再写入密钥文件。
11.4 使用示例
1)示例1
#include <string.h>
#include <openssl/bio.h>
int main()
{
BIGNUM *bn;
BIO *b;
char a[20];
int ret;
bn=BN_new();
strcpy(a,"32");
ret=BN_hex2bn(&bn,a);
b=BIO_new(BIO_s_file());
ret=BIO_set_fp(b,stdout,BIO_NOCLOSE);
BIO_write(b,"aaa",3);
BN_print(b,bn);
BN_free(bn);
return 0;
}
2)示例2
加法运算
#include <openssl/bn.h>
#include <string.h>
#include <openssl/bio.h>
int main()
{
BIGNUM *a,*b,*add;
BIO *out;
char c[20],d[20];
a=BN_new();
strcpy(c,u5b8bu4f53; font-size:16px; text-indent:2em"> ret=BN_hex2bn(&a,c);
b=BN_new();
strcpy(d,"100");
ret=BN_hex2bn(&b,d);
out=BIO_new(BIO_s_file());
ret=BIO_set_fp(out,u5b8bu4f53; font-size:16px; text-indent:2em"> add=BN_new();
ret=BN_add(add,a,b);
if(ret!=1)
{
printf("err.\n");
return -1;
}
BIO_puts(out,"bn 0x32 + 0x100 = 0x");
BN_print(out,add);
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',"\n");
BN_free(a);
BN_free(b);
BN_free(add);
BIO_free(out);
}
3) 示例3
减法运算
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',*sub;
BIO *out;
sub=BN_new();
ret=BN_sub(sub,"bn 0x100 - 0x32 = 0x");
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',sub);
BN_free(sub);
4)示例4
乘法运算
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',*mul;
BN_CTX *ctx;
ctx=BN_CTX_new();
mul=BN_new();
ret=BN_mul(mul,b,ctx);
}
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',"bn 0x32 * 0x100 = 0x");
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',mul);
BN_free(mul);
BN_CTX_free(ctx);
}
5)示例5
除法运算
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',*div,*rem;
ctx=BN_CTX_new();
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',"17");
div=BN_new();
rem=BN_new();
ret=BN_div(div,rem,"bn 0x100 / 0x17 =0x");
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',div);
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',"bn 0x100 % 0x17 =0x");
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',rem);
BN_free(div);
BN_free(rem);
6)示例6
平方运算
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',*sqr;
char c[20];
sqr=BN_new();
ret=BN_sqr(sqr,"bn 0x100 sqr =0x");
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',sqr);
BN_free(sqr);
}
7)示例7
次方运算
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',*exp,*b;
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',"3");
exp=BN_new();
ret=BN_exp(exp,"bn 0x100 exp 0x3 =0x");
ottom:10px; padding-top:0px; padding-bottom:0px; font-family:'Hiragino Sans GB W3',exp);
BN_free(exp);
}
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。