DES加密算法

借鉴博客

https://blog.csdn.net/m0_37962600/article/details/79912654

#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];
    }
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值