RSA算法是优秀的公钥密码算法，使用十分广泛。

RSA可用于加密、签名等，应用于PKI。

1）初始化：随机选择两个大素数p、q， n = p * q

2) 密钥生成：随机选择整数e , 1<e<(p-1)(q-1),且 gcd(e,(p-1)(q-1))=1,求e*d=1 mod (p-1)(q-1)

3)公钥（e，n），私钥 （ d ）

加密 明文 m<n ，密文 c=m^e mod n

解密 明文 m = c^d mod n

小素数情况下的实现：

 

 
  
  1. #include<stdio.h> 
  2. #include<stdlib.h> 
  3. #include<time.h> 
  4. typedef long long ULL; 
  5.  
  6. ULL extended_eculid(ULL a , ULL b, ULL &s, ULL &t)//扩展Eculid定理 
  7. { 
  8.     if(b == 0) 
  9.     { 
  10.         s = 1; 
  11.         t = 0; 
  12.         return a; 
  13.     } 
  14.     ULL d, x, y; 
  15.     d = extended_eculid(b,a % b,x,y); 
  16.     s = y; 
  17.     t = x - a/b*y; 
  18.     return d; 
  19. } 
  20.  
  21. ULL Squ_Mul(ULL a, ULL b, ULL n)//模重复平方法求(a^b)mod n 
  22. { 
  23.     ULL ans; 
  24.     ans = 1; 
  25.     while(b!=0) 
  26.     { 
  27.         if(b&1==1) 
  28.             ans = (ans*a)%n; 
  29.         a = (a*a)%n; 
  30.         b>>=1; 
  31.     } 
  32.     return ans; 
  33. } 
  34.  
  35. int Miller_Rabin(ULL n,int k)//Miller_Rabin素性测试 
  36. { 
  37.     ULL m = n - 1; 
  38.     int i; 
  39.     while((m%2)!=1) 
  40.         m>>=1; 
  41.     ULL a; 
  42.     while(1) 
  43.     { 
  44.         time_t random = time(NULL);  
  45.         a = random%(n-1); 
  46.         if(a!=0 && a!=1) 
  47.             break; 
  48.     } 
  49.     ULL b = Squ_Mul(a,m,n); 
  50.     if(b == 1) 
  51.         return 1; 
  52.     for(i = 0; i < k; i++) 
  53.     { 
  54.         if( b == n -1) 
  55.         { 
  56.             return 1; 
  57.         } 
  58.         else 
  59.         { 
  60.             b = ((b%n)*(b%n))%n; 
  61.         } 
  62.     } 
  63.     return 0; 
  64. } 
  65.  
  66. int menu()//菜单 
  67. { 
  68.     int select=1; 
  69.     printf("***********RSA**********\n"); 
  70.     printf("1 生成素数对和密钥对\n"); 
  71.     printf("2 加密\n"); 
  72.     printf("3 解密\n"); 
  73.     printf("0 退出\n"); 
  74.     printf("***********************\n"); 
  75.     while(1) 
  76.     { 
  77.         printf("请选择:"); 
  78.         scanf("%d",&select); 
  79.         if(select == 0 || select == 1 || select == 2 || select == 3) 
  80.             return select; 
  81.     } 
  82. } 
  83.  
  84.  
  85. ULL GenPrime( ULL &p, int bitNum,int k)//生成素数 
  86. { 
  87. //  p = 1; 
  88. //  p |= (1<<(bitNum -1)); 
  89.     time_t random = time(NULL);  
  90.     int count = 100; 
  91.     ULL con = 1<<bitNum; 
  92.     if(p%2==0) 
  93.         p = (p+1)%con; 
  94.     while(count) 
  95.     { 
  96.         if(Miller_Rabin(p,k)) 
  97.             return 1;        
  98.         p = (p+time(NULL))%con; 
  99.         p |= (1<<(bitNum -1)); 
  100.         if(p%2==0) 
  101.             p = (p+1)%con; 
  102. //      p = (p+2)%con; 
  103.         if(p == 0 || p == 1) 
  104.             p = 3; 
  105.         count--; 
  106.     } 
  107.     return 0; 
  108.  
  109. } 
  110. void GenPQ(ULL &p, ULL &q) 
  111. { 
  112.     int BitNumP = 10,BitNumQ = 10; 
  113.     printf("请输入要生成的素数 P 的bit位数(4<=bit<31):"); 
  114.     scanf("%d",&BitNumP); 
  115.     p = q = 3; 
  116.     p|=1<<(BitNumP-1); 
  117.     if(GenPrime(p,BitNumP,5) == 0) 
  118.     { 
  119.         printf("没有找到素数p\n"); 
  120.         exit(0); 
  121.     } 
  122.     printf("请输入要生成的素数 Q 的bit位数(4<=bit<31):"); 
  123.     scanf("%d",&BitNumQ); 
  124.     q|=1<<(BitNumQ-1); 
  125.     int count = 100; 
  126.     while(count) 
  127.     { 
  128.         if(GenPrime(q,BitNumQ,5)==0) 
  129.         { 
  130.             printf("没有找到素数q\n"); 
  131.             exit(0); 
  132.         } 
  133.         if(q!=p) 
  134.             break; 
  135.          else 
  136.              q = q + 2; 
  137.         count--; 
  138.     } 
  139.     if(count == 0) 
  140.     { 
  141.         printf("没有找到素数q\n"); 
  142.         exit(0); 
  143.     } 
  144. } 
  145.  
  146. int main() 
  147. { 
  148.     int select = menu(); 
  149.     ULL m,n,c; 
  150.     ULL Euler_N; 
  151.     ULL e,d=0; 
  152.     while(select) 
  153.     { 
  154.     switch (select) 
  155.     { 
  156.     case 1: 
  157.     { 
  158.         ULL p,q; 
  159.         p = q = 0; 
  160.         GenPQ(p,q);//生成P Q 
  161.         n = p*q; 
  162.          
  163.         Euler_N = (p-1)*(q-1);//计算欧拉函数值 
  164.         time_t random = time(NULL);  
  165.         int count = 100; 
  166.         ULL t = 0; 
  167.         for(e = 3; e < Euler_N; e+=2) 
  168.         {                             
  169.             if(extended_eculid(e,Euler_N,d,t) == 1) 
  170.                 break; 
  171.         } 
  172.         d = (d+Euler_N)%Euler_N; 
  173.         if(e>=Euler_N) 
  174.         { 
  175.             printf("没有合适的e\n"); 
  176.             exit(0); 
  177.         } 
  178.         printf(" p= %lld q= %lld n= %lld e= %lld d= %lld \n",p,q,n,e,d); 
  179. /*      while(count) 
  180.         { 
  181.             e = time(0); 
  182.             e|= 1; 
  183.             e = e%Euler_N; 
  184.             if(e!=0 && e!=1 && extended_eculid(e,Euler_N,d,t) == 1) 
  185.             { 
  186.                 break; 
  187.             } 
  188.             count--; 
  189.         }  
  190.         if(count==0) 
  191.         { 
  192.             printf("没有e\n"); 
  193.             exit(0); 
  194.         }*/ 
  195.         break; 
  196.         } 
  197.         case 2:                         //用公钥加密 
  198.         { 
  199.             printf("导入公钥 N :"); 
  200.             scanf("%lld",&n); 
  201.             printf(" n= %lld ",n); 
  202.             printf("导入公钥 E :"); 
  203.             scanf("%lld",&e); 
  204.             printf(" e= %lld \n",e); 
  205.             ULL leftLen = n/(1<<30);//此处主要是控制m不能太小，否则得到的C太大，计算时会溢出，出错 
  206.             ULL rightLen = (n<(1<<30))?n:(1<<30);//m也不能太大，否则计算m^e mod n时，m*m会有问题 
  207.             printf("请输入明文(%lld<m<%lld) :",leftLen,rightLen); 
  208.             scanf("%lld",&m); 
  209.             if(m>=rightLen || m<=leftLen) 
  210.             { 
  211.                 printf("明文分组不符合要求！\n"); 
  212.                 break; 
  213.             } 
  214.          
  215.              
  216.             printf(" m=%lld ",m); 
  217.             c = Squ_Mul(m,e,n); 
  218.             printf(" c= %lld \n",c); 
  219.  
  220.             break; 
  221.         } 
  222.         case 3:                             //用私钥解密 
  223.         { 
  224.             printf("导入公钥 N :"); 
  225.             scanf("%lld",&n); 
  226.             printf(" n= %lld ",n); 
  227.              
  228.             printf("导入私钥 D :"); 
  229.             scanf("%lld",&d); 
  230.             printf(" d= %lld \n",d); 
  231.              
  232.             printf("请输入密文 :");       
  233.             scanf("%lld",&c); 
  234.          
  235.             printf(" c=%lld ",c); 
  236.             m = Squ_Mul(c,d,n); 
  237.             printf(" m= %lld \n",m); 
  238.             break; 
  239.         } 
  240.         default:break; 
  241.     } 
  242.     select = menu(); 
  243.     } 
  244.     return 0; 
  245. }