借鉴博客
#include <stdio.h>
#include <stdlib.h>
#include<string.h>
#define mmax 64
int kn[16][48]={0};//十六个子密钥
int cn[17][28]={0};
int dn[17][28]={0};
int NMW[64]={0};//使用IP表的明文
int sh[32]={0};//s盒置换的结果
int RL[64]={0};//最后的二进制密文
int R0[32]={0};
int L0[32]={0};
int NR0[48]={0};//扩展后的R0
int E[48]={0};//E盒异或后的值
int P[32];//P盒转换的值
int t1[56]={57,49,41,33,25,17,9,1,58,50,42,34,26,18,10,2,59,51,43,35,27,19,11,3,60,52,44,36,63,55,47,39,31,23,15,7,62,54,46,38,30,22,14,6,61,53,45,37,29,21,13,5,28,20,12,4};
int t0[16]={1,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1};
int T2[48]={14,17,11,24,1,5,3,28,15,6,21,10,23,19,12,4,26,8,16,7,27,20,13,2,41,52,31,37,47,55,30,40,51,45,33,48,44,49,39,56,34,53,46,42,50,36,29,32};
int T3[64]={58,50,42,34,26,18,10,2,60,52,44,36,28,20,12,4,62,54,46,38,30,22,14,6,64,56,48,40,32,24,16,8,57,49,41,33,25,17,9,1,59,51,43,35,37,19,11,3,61,53,45,37,29,21,13,5,63,55,47,39,31,23,15,7};
int T4[48]={32,1,2,3,4,5,4,5,6,7,8,9,8,9,10,11,12,13,12,13,14,15,16,17,16,17,18,19,20,21,20,21,22,23,24,25,24,25,26,27,28,29,28,29,30,31,32,1};
int t5[32]={16,7,20,21,29,12,28,17,1,15,23,26,5,18,31,10,2,8,24,14,32,27,3,9,19,13,30,6,22,11,4,25};//P盒
int T6[64]={40,8,48,16,56,24,64,32,39,7,47,15,55,23,63,31,38,6,46,14,54,22,62,30,37,5,45,13,53,21,61,29,36,4,44,12,52,20,60,28,35,3,43,11,51,19,59,27,34,2,42,10,50,18,58,26,33,1,41,9,49,17,57,25};
int sbox[8][4][16]={
//S1
14, 4,13, 1, 2,15,11, 8, 3,10, 6,12, 5, 9, 0, 7,
0,15, 7, 4,14, 2,13, 1,10, 6,12,11, 9, 5, 3, 8,
4, 1,14, 8,13, 6, 2,11,15,12, 9, 7, 3,10, 5, 0,
15,12, 8, 2, 4, 9, 1, 7, 5,11, 3,14,10, 0, 6,13,
//S2
15, 1, 8,14, 6,11, 3, 4, 9, 7, 2,13,12, 0, 5,10,
3,13, 4, 7,15, 2, 8,14,12, 0, 1,10, 6, 9,11, 5,
0,14, 7,11,10, 4,13, 1, 5, 8,12, 6, 9, 3, 2,15,
13, 8,10, 1, 3,15, 4, 2,11, 6, 7,12, 0, 5,14, 9,
//S3
10, 0, 9,14, 6, 3,15, 5, 1,13,12, 7,11, 4, 2, 8,
13, 7, 0, 9, 3, 4, 6,10, 2, 8, 5,14,12,11,15, 1,
13, 6, 4, 9, 8,15, 3, 0,11, 1, 2,12, 5,10,14, 7,
1,10,13, 0, 6, 9, 8, 7, 4,15,14, 3,11, 5, 2,12,
//S4
7,13,14, 3, 0, 6, 9,10, 1, 2, 8, 5,11,12, 4,15,
13, 8,11, 5, 6,15, 0, 3, 4, 7, 2,12, 1,10,14, 9,
10, 6, 9, 0,12,11, 7,13,15, 1, 3,14, 5, 2, 8, 4,
3,15, 0, 6,10, 1,13, 8, 9, 4, 5,11,12, 7, 2,14,
//S5
2,12, 4, 1, 7,10,11, 6, 8, 5, 3,15,13, 0,14, 9,
14,11, 2,12, 4, 7,13, 1, 5, 0,15,10, 3, 9, 8, 6,
4, 2, 1,11,10,13, 7, 8,15, 9,12, 5, 6, 3, 0,14,
11, 8,12, 7, 1,14, 2,13, 6,15, 0, 9,10, 4, 5, 3,
//S6
12, 1,10,15, 9, 2, 6, 8, 0,13, 3, 4,14, 7, 5,11,
10,15, 4, 2, 7,12, 0, 5, 6, 1,13,14, 0,11, 3, 8,
9,14,15, 5, 2, 8,12, 3, 7, 0, 4,10, 1,13,11, 6,
4, 3, 2,12, 9, 5,15,10,11,14, 1, 7, 6, 0, 8,13,
//S7
4,11, 2,14,15, 0, 8,13, 3,12, 9, 7, 5,10, 6, 1,
13, 0,11, 7, 4, 0, 1,10,14, 3, 5,12, 2,15, 8, 6,
1, 4,11,13,12, 3, 7,14,10,15, 6, 8, 0, 5, 9, 2,
6,11,13, 8, 1, 4,10, 7, 9, 5, 0,15,14, 2, 3,12,
//S8
13, 2, 8, 4, 6,15,11, 1,10, 9, 3,14, 5, 0,12, 7,
1,15,13, 8,10, 3, 7, 4,12, 5, 6,11, 0,14, 9, 2,
7,11, 4, 1, 9,12,14, 2, 0, 6,10,13,15, 3, 5, 8,
2, 1,14, 7, 4,10, 8,13,15,12, 9, 0, 3, 5, 6,11
};
void DecToBit(int kk[] , int k[] , int n);//十进制转二进制
void HexToBit(char sstring[] , int k[]); //十六进制转二进制
void BitToHex(int bit[] , int n);//二进制转为十六进制
void ORfac(int k[] , int kk[] , int n , int m[]);//异或
void Change1(int k[]);
void Chang2(int k[] , int cn[][28] , int dn[][28]);
void Chang3(int kn[][48] , int cn[][28] , int dn[][28]);
void MWdes(int MW[],int NMW[]);
void R0L0(int M[64],int L0[],int R0[]);//得到R0L0
void KZ(int R0[],int NR0[]);//扩展成48位
void YH(int NR0[],int kn[16][48],int a,int E[]);
void XY(int E[]);
void Pchang(int R[],int p[]);
void IPChange(int a[]);
void EnChange(int a[] , int b[] , int n);//数组交换
int main()
{
printf("请输入需要加密的明文M(16位十六位数):\n");
char sstring1[16] , sstring2[16];
int Mkey[64];//0123456789ABCDEF
int i , j;
int k[64];//密钥133457799BBCDFF1
int pen[mmax];//p盒异或
gets(sstring1);
HexToBit(sstring1 , Mkey);
printf("请输入密钥K(16位十六位数):\n");
gets(sstring2);
HexToBit(sstring2 , k);
Change1(k);
Chang2(k,*cn,*dn);
Chang3(*kn,*cn,*dn);
//ChangeEnd(Mkey);
MWdes(Mkey , NMW);
R0L0(NMW , L0 , R0);
for(i=0 ; i<16 ; i++){
KZ(R0 , NR0);
ORfac(NR0 , kn[i] , 48 , E);
XY(E);
Pchang(sh,P);
ORfac(P,L0,32,pen);
EnChange(R0,L0,32);
EnChange(pen,R0,32);
}
for(i=0 ; i<64 ; i++){
if(i<32){
RL[i]=R0[i];
}
else{
RL[i]=L0[i-32];
}
}
IPChange(RL);
printf("\n\n加密结果为:");
BitToHex(RL , 64);
return 0;
}
void ORfac(int k[] , int kk[] , int n, int m[])//n为k和kk的长度,异或结果存于m[]
{
int i=0;
for(i=0 ; i<n ; i++){
m[i]=k[i]^kk[i];
// printf("%d " , m[i]);
}
}
void Change1(int k[])//密钥k的第一次转化
{
int k1[mmax];
int i;
for(i=0 ; i<64 ; i++){
k1[i]=k[i];
}
int j=0;
for(i=0 ; i<63 ; i++){
if((i+1)%8==0)i++;
k[j]=k1[t1[j]-1];
// printf("%d" , k[j]);
j++;
}
}
void Chang2(int k[] , int cn[][28] , int dn[][28])
{
int i=0 , j=0 , f=0;
for(i=0 ; i<56 ; i++){
if(i<28){
cn[0][j]=k[i];
j++;
}
else{
dn[0][f]=k[i];
f++;
}
}
int p=t0[0];//记录cn相对c0移动的位数
for(i=1 ; i<17 ; i++){
for(j=0 ; j<28-p ; j++){
cn[i][j]=cn[0][j+p];//1,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1
//printf("%d",cn[i][j]);
}
for(f=0 ; f<p;f++){
cn[i][j]=cn[0][f];
// printf("%d",cn[i][j]);
j++;
}
// printf("\n");
p+=t0[i];
}
p=t0[0];
for(i=1 ; i<17 ; i++){
for(j=0 ; j<28-p ; j++){
dn[i][j]=dn[0][j+p];//1,1,2,2,2,2,2,2,1,2,2,2,2,2,2,1
// printf("%d",dn[i][j]);
}
for(f=0 ; f<p;f++){
dn[i][j]=dn[0][f];
// printf("%d",dn[i][j]);
j++;
}
// printf("\n");
p+=t0[i];
}
}
void Chang3(int kn[][48] , int cn[][28] , int dn[][28])
{
int i , j;
int pn[16][56];
int f=0;
for(i=1 ; i<17 ; i++){
for(j=0 ; j<56 ; j++){
if(j<28)
pn[i-1][j]=cn[i][j];
else{
pn[i-1][j]=dn[i][j-28];
}
}
}
for(i=0 ; i<16 ; i++){
f=0;
for(j=0 ; j<55 ; j++){
if((j+1)%7==0)j++;
kn[i][f]=pn[i][T2[f]-1];
// printf("%d" , kn[i][f]);
f++;
}
// printf("\n");
}
}
void DecToBit(int kk[] , int k[] , int n)//十进制转二进制
{
int i=0 , j=0 , p;
for(i=0 ; i<n ; i++){
p=3;
for(; p>=0 ; p--){
if(kk[i]!=0){
k[j+p]=kk[i]%2;
kk[i]/=2;
}
else{
k[j+p]=0;
}
}
j+=4;
}
}
void HexToBit(char sstring[] , int k[]) //十六进制转二进制
{
int i=0 , j=0;
int kk[mmax];
int l=strlen(sstring);
for(i=0 ; i<l ; i++){
if(isdigit(sstring[i])!=0){
kk[i]=(sstring[i]-'0')*1;
}
else{
kk[i]=(sstring[i]-'A'+10)*1;
}
}
DecToBit(kk , k , l);
}
void BitToHex(int bit[] , int n)//二进制数组及其长度
{
int i=0;
int l=n/4;
int sum[mmax]={0};//存二进制转化的十进制数
for(i=0 ; i<l ; i++){
int j=(i+1)*4-1;
int p;
for(p=0 ; p<4 ; p++){
sum[i]+=bit[j]*pow(2,p);
j--;
}
//printf("%d " , sum[i]);
}
for(i=0 ; i<l ; i++){
if(sum[i]<10){
printf("%d" , sum[i]);
}
else{
char c;
c=sum[i]+'A'-10;
printf("%c" , c);
}
}
}
void MWdes(int MW[],int NMW[]){//对明文加密,扩展到64位
int i;
for(i=0;i<64;i++){
NMW[i]=MW[T3[i]-1];
// printf("%d %d\n" ,NMW[i] , T3[i]);
}
//R0(NMW);
}
void R0L0(int M[64],int L0[],int R0[]){//得到R0L0
int i,j=0,k=0;
//printf("%d\n",M[63]);
for(i=0;i<32;i++){
L0[k]=M[i];
k++;
}
for(i=32;i<64;i++){
R0[j]=M[i];
//printf("%d\n",R0[j]);
j++;
}
//KZ(R0);
}
void KZ(int R0[],int NR0[]){//扩展成48位
//int NR0[48];
int i,j,k;
for(i=0;i<48;i++){
NR0[i]=R0[T4[i]-1];
}
//YH(NR0,k1);
}
void YH(int NR0[],int kn[16][48],int a,int E[]){
int i=0;
for(i=0;i<48;i++){
if(NR0[i]==kn[a][i]){
E[i]=0;
}else{
E[i]=1;
}
}
}
void XY(int E[]){
int i=0 , j=0;
int x[9],y[9];
for(i=0;i<48;){
x[j]=E[i]*2+E[i+5];
y[j]=E[i+1]*8+E[i+2]*4+E[i+3]*2+E[i+4];
i+=6;
j++;
}
int a=31;
int dec[mmax];
for(i=0 ; i<8 ; i++){
dec[i]=sbox[i][x[i]][y[i]];
// printf("%d %d %d\n" , dec[i] , x[i] , y[i]);
}
DecToBit(dec,sh,32);
}
void Pchang(int R[],int P[])
{
int i=0;
//int p[32];
for(i=0 ; i<32 ; i++){
P[i]=R[t5[i]-1];
}
}
/*void ECG(int P[],int R0[]){
int R1[32],L1[32];
*R1=*P;
*L1=*R0;
KZ(R0);
}*/
void EnChange(int a[] , int b[] , int n)//a=>b
{
int i , j;
for(i=0 ; i<n ; i++){
b[i]=a[i];
}
}
void IPChange(int a[])
{
int i=0;
int b[64];
for(i=0 ; i<64 ; i++){
b[i]=a[i];
}
for(i=0 ; i<64 ; i++){
a[i]=b[T6[i]-1];
}
}