大家好,又见面了,我是你们的朋友全栈君。
这次轮到RSA加密算法了。RSA加密过程相对DES和MD5要简单很多,但作为现在还在使用的加密算法之一,它还是有需要认真思索的地方哒~
首先是密钥对的生成:
(1)选取两个大素数p和q(目前两个数的长度都接近512bit是安全的)
(2)计算乘积n=p*q,Φ(n)=(p-1)(q-1),其中Φ(n)为n的欧拉函数(因为两素数乘积的欧拉函数等于两数分别减一后的乘积)
(3)随机选取整数e(1<e<Φ(n))作为公钥d,要求满足e与Φ(n)的最大公约数为1,即两者互素
(4)用Euclid扩展算法计算私钥d,已满足d * e ≡ 1 (mod Φ(n)),即d ≡ e^(-1) (mod Φ(n))。则e与n是公钥,d是私钥
注意:e与n应公开,两个素数p和q不再需要,可销毁,但绝不可泄露。
加密过程:
将接收到的明文转换成特定的编码方式。如p=43,q=59,e=13,明文为cybergreatwall,按照英文字母表的顺序a=00,b=01,… ,z=25进行编码后为022401041706001922001111。
然后将转码后的字符串分块,分组要求:每个分组对应的十进制数小于0。这个要求是什么意思呢?我个人的理解通过举例向大家说明:上文字符串分组如下0224 0104 1706 0019 2200 1111。每一分组的数都小于n(2537),而2537能接受的最大的数为2525(也就是‘zz’的情况),所以是4位1组,即两字符一组。这样一来,m1=0224,m2=0104,… ,m6=1111
现在可以加密了~~加密算法就是这个式子—-ci ≡ mi^e (mod n),如第一分组 0224^13 ≡ mod 2537 ≡ 1692=c1 。这里有个隐藏的算法是需要了解的:
在RSA算法过程中容易出现天文数字(像上文的0224^13),而这些天文数字会为我们编程的过程造成一定的麻烦,更可恶的是会影响速度!!为了避免这种情况,快速取模指数算法可以很有效地算出c≡m^e mod n的准确结果且避免过程中出现天文数字~~
下面用伪代码为大家介绍下这种神奇的算法(个人感觉伪代码里的 ‘<-’ 就是平时用的 ‘=’ ):
t<-0;c<-1
for i<-k downto 0
do t<-2*t
c<-(c*c)mod n
if bi=1 then t<-t+1
c<-(c*m)mod n
return c
(p.s:e的二进制表示为bk bk-1 … b0,如e=13=(1101),所以k为3)
所以,用快速取模指数算法计算上文例子里的c1过程就该是这样子哒:
i 3 2 1 0
bi 1 1 0 1
t 1 3 6 13
ci (1*224)mod2537 (224*224*224)mod2537 (514*514)mod2537 (348*348*224)mod2537
=224 =514 =348 =1692
到这里RSA加密的算法就讲完了,下面附上代码
#include<stdio.h>
#include<stdlib.h>
/* 函数申明 */
int long_n(int n);
int shuru(char *arr, int k, char *wei, int is_first);
void jiami(char *arr, int k, int e, int n);
/* 输入函数,记录从键盘输入的明文*/
int shuru(char *arr, int k, char *wei, int is_first)
{
int i;
char ch;
/*判断是否为第一分组的输入,如果是则获取输入的字符,否则就将上一分组最后获取的字符作为这一分组的第一个字符*/
if (is_first == 1)
ch = getchar();
else
ch = *wei;
for (i = 0; (i < k) && (ch != '\n');i++) //获取字符直到获取到回车符为止
{
arr[i] = ch;
ch = getchar();
}
*wei = ch; //最后获取到的字符准备作为下一分组的第一个字符
for (i = i; i < k; i++)
arr[i] = 'a'; //输入不够一组个数的整数倍则补'a'(即为补零)
if (ch == '\n') //接收到回车符返回0,否则为1
return 0;
else
return 1;
}
/*加密函数*/
void jiami(char *arr, int k, int e, int n)
{
int m = 0,c=1, i, j,t=0, shu,temp,num=0;
int *array;
/*Mi赋值过程*/
for (i = 0; i < k; i++)
{
temp = 1;
for (j = 0; j < (k-i-1)*2; j++)
temp = temp * 10;
shu = (int)arr[i] - 97;
m = m + temp * shu;
}
temp = e;
/*获取e的二进制表达形式的位数*/
do{
temp = temp / 2;
num++;
} while (temp != 0);
array = (int *)malloc(sizeof(int)*k); //申请动态数组
temp = e;
/*动态数组存储e的二进制表达形式*/
for (i = 0; i < num; i++)
{
array[i] = temp % 2;
temp = temp / 2;
}
/*避免出现天文数字的算法,详情见上文文字说明*/
for (i = num - 1; i >= 0; i--)
{
t = t * 2;
temp = c*c;
if (temp > n)
{
for (j = 0; temp - n*j >= 0; j++);
j--;
c = temp - n*j;
}
else
c = temp;
if (array[i] == 1)
{
t = t + 1;
temp = c*m;
if (temp > n)
{
for (j = 0; temp - n*j >= 0; j++);
j--;
c = temp - n*j;
}
else
c = temp;
}
e = e / 2;
}
temp = c;
i = 0;
/*c的位数小于分组长度则在前补零*/
do{
temp = temp / 10;
i++;
} while (temp != 0);
for (i; i < num; i++)
printf("0");
printf("%d", c);
}
/*获取分组的长度*/
int long_n(int n)
{
int temp,i,j,k,shi,comp=0;
temp = n;
/*获取n的位数*/
for (i = 1; temp / 10 != 0; i++)
{
temp = temp / 10;
}
temp = i;
/*若n的位数为基数*/
if (i % 2 != 0)
{
i = i - 1;
return i;
}
/*若位数为偶数*/
else
{
for (j = 0; j < i/2; j++)
{
shi = 1;
for (k = 0; k < temp - 2; k++)
shi = shi * 10;
comp = comp + shi * 25;
temp = temp - 2;
}
if (comp <= n)
return i;
else
{
i = i - 2;
return i;
}
}
}
/*主函数*/
int main()
{
int p, q, e, d, n, fai_n, k, i,is_first=1;
char ch,*arr,wei='a';
printf("请输入p、q、e值,用空格间隔开\n");
scanf_s("%d%d%d", &p, &q, &e); //从键盘获取p、q、e值
n = p*q;
fai_n = (p-1)*(q-1); //Φ(n)
for (k = 0; (k*n + 1) % e != 0; k++);
if ((k*n + 1) % e == 0)
d = (k*n + 1) / e; //d * e ≡ 1 (mod Φ(n))
k = long_n(n);
k = k / 2; //分组的长度
ch = getchar(); //缓冲回车符
arr = (char *)malloc(sizeof(char)*k); //申请动态数组
printf("请输入明文\n");
while (1)
{
i=shuru(arr,k,&wei,is_first); //调用输入字符的函数,接收到回车符返回0,否则为1
is_first = 0; //第一分组录入结束设为0
jiami(arr,k,e,n); //调用加密函数
if (i == 0) //接收到返回值为0跳出循环
break;
}
printf("\n");
return 0;
}
运行结果:
如果有对算法和代码不理解或者有不同看法的请联系我哦,邮箱在我的信息里~~
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/132381.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...