yingying.cpp

加密解密现代密码学的内容c语言编写uoto

#include<stdio.h>
#include<malloc.h>
#define MAX_NUM 40

static long (*sp)[64];
/*二维数组sp[8][64]中存放轮函数中经过"S盒代换/选择和P置换"后的所有可能输出结果*/
/*注:具体说明详见函数spinit,函数spinit先于DES加密过程进行*/

static unsigned char (*kn)[8];
/*当用户输入初始密钥之后，二维数组kn[16][8]中将存放由初始密钥生成的16轮子密钥*/
/*关于kn的操作函数setkey独立于DES加密过程进行*/

static char (*iperm)[16][8];    /*用于初始置换*/
static char (*fperm)[16][8];    /*用于逆初始置换*/


static int bytebit[]={0200,0100,040,020,010,04,02,01};
/*任意一个小于256(二进制表示为100000000)的数与bytebit的第i个数作按位与运算,将*/
/*可得到这个数二进制表示中的第i+1比特(数组下标从0开始).若二进制表示不足八位,*/
/*程序默认高位自动填0*/
/*例如14(即00001110)这个数与bytebit中的各个数作按位于运算,依次得到0,0,0,0,1,1,*/
/*1,0. 具体使用详见函数setkey*/

static int nibblebit[]={010,04,02,01}; /*用于一些函数中的辅助计算*/

static int permute(char *,char [16][16][8],char *);  /*见此函数体里的说明*/
static void round(int ,unsigned long *);             /*见此函数体里的说明*/
static void spinit();          /*见此函数体里的说明*/
static void perminit(char [16][16][8],char [64]);    /*见此函数体里的说明*/
static long f(unsigned long, unsigned char [8]);     /*见此函数体里的说明*/


/*明文所要进行的初始置换的比特变换顺序*/
static char ip[]={
    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,27,19,11,3,
    61,53,45,37,29,21,13,5,
    63,55,47,39,31,23,15,7
    };

/*逆初始置换的比特变换顺序*/
static char fp[]={
    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
};


#ifdef notdef
/*轮函数中选择扩展运算比特存放顺序*/
  static char ei[]={
    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
  };
#endif


/*密钥编排的初始置换的比特变换顺序*/
static char pc1[]={
    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,
};


/*存放密钥编排的16轮置换中左循环移位的位数*/
static char totrot[]={
    1,2,4,6,8,10,12,14,15,17,19,21,23,25,27,28
}
/*注:totrot[i]实际存放的是当前轮与之前所有轮循环左移位数的总和,具体应用详见函*/
/*数setkey*/;


/*密钥编排的16轮置换中所用到的比特变换顺序*/
static char pc2[]={
    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
};

/*S盒*/
static char si[8][64]={
    /*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,9,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,9,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
};


/*轮换数中置换预算的比特变换顺序*/
static char p32i[]={
    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
};


int desinit(){

    /*对DES操作中会用的一些变量进行初始化*/

    if(sp!=NULL)   return 0;

    if((sp=(long(*)[64])malloc(sizeof(long)*8*64))==NULL)    
        return -1;
    /*初始化二维数组sp[8][64]*/

    spinit();  /*功能详见spinit函数定义处*/

    kn=(unsigned char (*)[8])malloc(sizeof(char)*8*16);
    /*初始化二维数组kn[16][8]*/

    if(kn==NULL){
       free((char *)sp);
       return -1;
    }

    iperm=(char(*)[16][8])malloc(sizeof(char)*16*16*8);
    /*初始化三维数组iperm[16][16][8]*/
    if(iperm==NULL){
        free((char *)sp);
        free((char *)kn);
        return -1;
    }
    perminit(iperm,ip);

    fperm=(char(*)[16][8])malloc(sizeof(char)*16*16*8);
    /*初始化三维数组fperm[16][16][8]*/
    if(fperm==NULL){
       free((char *)sp);
       free((char *)kn);
       free((char *)iperm);
       return -1;
    }
    perminit(fperm,fp);
    return 0;
}

void desdone(){
    
    /*此函数作用为回收DES加密过程中所占用的空间*/
    
    if(sp==NULL)
        return;      

    free((char *)sp);
    free((char *)kn);

    if(iperm!=NULL)    free((char *)iperm);
    if(fperm!=NULL)    free((char *)fperm);

    sp=NULL;
    iperm=NULL;
    fperm=NULL;
    kn=NULL;
    return;
}


void setkey(char *key){
    /* 用户输入64bit初始密钥，只将使用56bit*/
    /*当用户输入初始密钥之后，经此函数二维数组kn[16][8]中将存放由初始密钥生成的*/
    /*16轮子密钥*/
    /*此程序先于DES加密过程(函数endes)进行*/

    char pclm[56]; /*存放初始密钥经置换表pc1[64]置换后每一比特*/
    char pcr[56]; /*存放当前轮左移后的结果*/
    register int i,j,l;
    int m;

    /* 初始化kn */
    for(i=0;i<16;i++)
       for(j=0;j<8;j++)
          kn[i][j]=0;

    for(j=0;j<56;j++){         /* 把key[]按置换表pc1[]转成比特存储*/
       l=pc1[j]-1;              /* 用户输入的key[]的下标从0开始，要减1*/
       m=l&07;
       pclm[j]=(key[l>>3]&bytebit[m])?1:0;
    }

    /*l按比特共6位.l的高3比特位位决定在key[]中第几字节即在key[]中的下标，*/
    /*l的低3位决定在key[]某个字节(即下标)的某个比特*/

    for(i=0;i<16;i++){
        for(j=0;j<56;j++)
          pcr[j]=pclm[(l=j+totrot[i])<(j<28?28:56)?l:l-28];
          /*循环左移*/
       for(j=0;j<48;j++) {

          if(pcr[pc2[j]-1]){

             l=j%6;
             kn[i][j/6]|=bytebit[l]>>2;
            /*将移位后的结果经pc-2置换表的48bit存入本轮密钥存储器kn中*/

          }
       }
    }
}


void endes(char *block){

    /*此函数实现DES加密的过程(在函数desinit,setkey被调用后,才会执行此函数)*/

    register int i;
    unsigned long work[2];   /*存放轮函数迭代中当前轮的左右分组各32bit*/
    long tmp;

    permute(block,iperm,(char *)work);   /*初始置换*/

/*    printf("分别加密的16轮每轮所得结果为:\n");
     // 以下进行16轮迭代
    for(i=0;i<16;i++)
        round(i,work);
    getch();                     */
    /* 16轮迭代后,还要将左右各32bit交换一下 */
    tmp=work[0];
    work[0]=work[1];
    work[1]=tmp;

    permute((char *)work,fperm,block);   /*逆初始置换*/
    printf("\n\n");
}



void  dedes(char *block){

    /*此函数实现DES解密 */
    register int i;