您的位置:首页 > 其它

用perl实现的des算法

2016-02-03 15:31 218 查看
#!/usr/bin/perl -w

use strict;

use warnings;

# final

my (@PC1,@PC_1,@PC2,@E,@S,@Y,@C,@H);

#variable

my ($une,$pwd,$magicnum,$iv,@pc,@keyset);

@Y=(56, 48, 40, 32, 24, 16, 8,

    0, 57, 49, 41,33, 25, 17,

    9, 1, 58, 50, 42, 34, 26,

    18, 10, 2, 59, 51, 43, 35,

    62, 54, 46, 38, 30, 22, 14,

    6, 61, 53, 45, 37, 29, 21,

    13, 5, 60,52, 44, 36, 28,

    20, 12, 4, 27, 19, 11, 3);

@PC1=(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,

      56,48,40,32,24,16,8,0,

      58,50,42,34,26,18,10,2,

      60,52,44,36,28,20,12,4,

      62,54,46,38,30,22,14,6);

@PC_1=(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,

       32,0,40,8,48,16,56,24);

@E=(31,0,1,2,3,4,

    3,4,5,6,7,8,

    7,8,9,10,11,12,

    11,12,13,14,15,16,

    15,16,17,18,19,20,

    19,20,21,22,23,24,

    23,24,25,26,27,28,

    27,28,29,30,31,0);

@C=(1, 2, 4, 6,

    8, 10, 12, 14,

    15, 17, 19,21,

    23, 25, 27, 28 );

@S=(##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]

   );

@H=(13, 16, 10, 23, 0, 4, 

    2, 27, 14, 5, 20,9,

    22, 18, 11, 3, 25, 7,

    15, 6, 26, 19, 12, 1,

    40, 51, 30, 36, 46,54,

    29, 39, 50, 44, 32, 47, 

    43, 48, 38, 55, 33, 52, 

    45, 41, 49, 35,28, 31 );

@PC2=(15,6,19,20,28,11,27,16,

      0,14,22,25,4,17,30,9,

      1,7,23,13,31,26,2,8,

      18,12,29,5,21,10,3,24);

     $une = "mail139";

     $pwd = "12345";

$magicnum = pack "H*","0123456789abcdef";

      $iv = pack "H*","0000000000000000";

  @keyset = ();

# pretreatemnt the text

sub preHandlePlain{

  my ($text,@param,@bytes)=(shift);

  $text=~s/(.)/chr(0).uc($1)/ge;

  @param = split //,$text;

  if(@param%8!=0){

     @param=(@param,(chr(0),chr(0),chr(0),chr(0),chr(0),chr(0),chr(0),chr(0))[0..7-@param%8])

  }

    foreach my $i(0..@param-1){

    foreach my $j(0..7){

      $bytes[$i*8+$j]=((ord($param[$i])&1<<7-$j)>0?1:0);

      #printf $bytes[$i*8+$j];

    }

  }

  return @bytes;

}

# translate hex string into bytes

sub hexstrToBytes{

  my (@param,@bytes)=(split //,shift);

  foreach my $i(0..@param-1){

    foreach my $j(0..7){

      $bytes[$i*8+$j]=((ord($param[$i])&1<<7-$j)>0?1:0);

      #printf $bytes[$i*8+$j];

    }

  }

  return @bytes;

}

# deal with each block

sub eachBlock{

  my ($flg,@tempkey,@byte1,@left,@right,@temp,@result)=(shift,@_);

  ##PC-1 subsitute

  @byte1=map{$tempkey[$PC1[$_]]}(0..63);

  ##each round

  @left=@byte1[0..31];

  @right=@byte1[32..63];

  for my $ii (0..15){

    my $rnd = ($flg==0?$ii:15-$ii);

    my @b1=@right;

    my @b2=funF($rnd,@right);

    @right=map{$left[$_]^$b2[$_]}(0..31);

    @left=@b1;

  }

  @temp=(@right,@left);

 

  ##PC-1 subsitute

  $result[$_]=$temp[$PC_1[$_]]for(0..63);

  return @result;

}

sub genKey{

  my @tempkey = @_;

  my @byte1;

  my @byte2;

  for my $i (0..55){

    $byte1[$i]=$tempkey[$Y[$i]];

  }

  for my $i (0..15){

    my $mov = $C[$i];

    for my $j (0..47){

      my $temp=$H[$j];

      ## $temp = 26

      if($temp<28){

        ## $temp = 26 , $mov =1

        if($temp+$mov<28){

          $byte2[$j]=$byte1[$temp+$mov];

        }

        ## $temp = 26 , $mov =3

        else{

          $byte2[$j]=$byte1[$temp+$mov-28];

        }

      }

       ## $temp = 53

      else{

         ## $temp = 53,$mov =2

        if($temp+$mov<56){

          $byte2[$j]=$byte1[$temp+$mov];

        }

          ## $temp = 53,$mov =5

        else{

          $byte2[$j]=$byte1[$temp+$mov-28];

        }

      }

    }

    $keyset[$i]=[@byte2];

  }

}

sub funF{

  my($rnd,@right,@temp,@byte1,@byte2)=(shift,@_);

  ## expand the right key and xor with the subkey

  @temp=map{$right[$E[$_]]^$keyset[$rnd][$_] }(0..47);

  ##S box substitute

    for(0..7){

      my $row=$temp[$_*6]*2+$temp[$_*6+5];

      my $col=$temp[$_*6+1]*8+$temp[$_*6+2]*4+$temp[$_*6+3]*2+$temp[$_*6+4];

      $byte1[$_*4] = ($S[$_][$row*16+$col]&8)>0?1:0;

      $byte1[$_*4+1] = ($S[$_][$row*16+$col]&4)>0?1:0;

      $byte1[$_*4+2] = ($S[$_][$row*16+$col]&2)>0?1:0;

      $byte1[$_*4+3] = ($S[$_][$row*16+$col]&1)>0?1:0;

    }

    ## p box substitute

    $byte2[$_]=$byte1[$PC2[$_]]for(0..31);

    return @byte2;

}

#encrypt

sub desEncrypt{

 my ($flg,@temp,@b0,@b2)=(shift,@_);

 @b0=hexstrToBytes($iv);

 for(0..@temp/64-1){

   my @b1 = @temp[$_*64..($_*64+63)];

   @b0 = map{$b0[$_]^$b1[$_]}(0..63);

   ## 0 encrypt ,1 decrypt

   @b0=eachBlock(0,@b0);

   if($flg==0){

    @b2=(@b2,@b0);

   }

 }

 if($flg==1){

  return @b0;

 }

 else{

  return @b2;

 }

}

#decrypt

sub desDecrypt{

  my ($flg,@temp,@b0,@b2,@b3)=(shift,@_);

   @b0=hexstrToBytes($iv);

   for(0..@temp/64-1){

   my @b1 = @temp[$_*64..($_*64+63)];

   @b3=eachBlock(1,@b1);

   @b3 = map{$b0[$_]^$b3[$_]}(0..63);

   ## 0 encrypt ,1 decrypt

   @b0=@b1; 

   if($flg==0){

    @b2=(@b2,@b3);

   }

 }

 if($flg==1){

  return @b3;

 }

 else{

  return @b2;

 }

}

sub getPwdHash{

#16 rounds of generate key

genKey(hexstrToBytes($magicnum));

my @temp1=preHandlePlain($une.$pwd);

my @result = desE
4000
ncrypt(0,@temp1);

genKey(@result[@result-64..@result-64+63]);

my $str;

($str.=$result[$_])for(0..@result-1);

$str=~s/([01]{4})/unpack("H",pack("B4",$1))/ge;

printf "getPwdHash -> %s\n",$str;

$str="";

my @result1 = desEncrypt(0,@temp1);

($str.=$result1[$_])for(0..@result1-1);

$str=~s/([01]{4})/unpack("H",pack("B4",$1))/ge;

printf "getPwdHash -> %s\n",$str;  

}

sub getDec{

  genKey(hexstrToBytes(pack("H*","a60228af1437c0b5")));  

my @result = desDecrypt(0,hexstrToBytes(pack("H*","005c0b27184dff90355eed8d814784cf35b80dfeb685e590")));

my $str;

($str.=$result[$_])for(0..@result-1);

#my $str="1111000011011001";

$str=~s/([01]{4})/unpack("H",pack("B4",$1))/ge;

printf "desdecrypt -> %s\n",$str;  

}

getPwdHash();

getDec();
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: