int e[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
		} ;
unsigned char m_bit1[8] = {0x0,0x0,0x0,0x0,0x0,0x0,0x0,0x0} ;
unsigned char m_byte[64] ;
unsigned char m_byte1[64] ;
unsigned char key_n[16][8] ;
unsigned char l_bit[17][8] ;
unsigned char r_bit[17][8] ;
unsigned char e_byte[64] ;
unsigned char e_byte1[64] ;
unsigned char r_byte[64] ;
unsigned char r_byte1[64] ;
int i = 0 , j = 0;

/*initialize*/
for(i=0;i<=15;i++)
	{
	for(j=0;j<=7;j++)
		{
		l_bit[i][j] = 0x0 ;
		r_bit[i][j] = 0x0 ;
		key_n[i][j] = 0x0 ;
		}
	}
for(i=0;i<=63;i++)
	{
	m_byte[i] = 0x0 ;
	m_byte1[i] = 0x0 ;
	r_byte[i] = 0x0 ;
	r_byte1[i] = 0x0 ;
	e_byte[i] = 0x0 ;
	e_byte1[i] = 0x0 ;
	}

keychange(k_bit,key_n) ;/*get the 48bit key x 16 (16rows x 6byte in key_n)*/
bit2byte(m_bit,m_byte) ;/*change to byte*/
for(i=0;i<=63;i++) m_byte1[i] = m_byte[ip[i] - 1] ;/*goto IP swap bit*/
byte2bit(m_byte1,m_bit1) ;/*re-change to bit*/
for(i=0;i<=3;i++) l_bit[0][i] = m_bit1[i] ;/*move left 32bit -> l0*/
for(i=4;i<=7;i++) r_bit[0][i - 4] = m_bit1[i] ;/*move right 32bit -> r0*/

for(i=1;i<=16;i++) /*16 layer*/
	{
	for(j=0;j<=3;j++) l_bit[i][j] = r_bit[i-1][j] ;/*L(n) = R(n-1)*/

	/*comput f(R(n-1),k)*/
	bit2byte(r_bit[i-1],r_byte) ;
	for(j=0;j<=47;j++) r_byte1[j] = r_byte[e[j] - 1] ;/*goto E swap bit*/
	byte2bit(r_byte1,r_bit[i-1]) ;/*now r_bit is 48bit*/

	/*xor 48bit key*/
	for(j=0;j<=5;j++) r_bit[i-1][j] = r_bit[i-1][j] ^ key_n[i-1][j] ;

	/*goto <S replace> and <P swap bit>*/
	s_replace(r_bit[i-1]) ;/*change 48bit r_bit[i-1]->32bit r_bit[i-1]*/
	for(j=0;j<=3;j++)/*get next r_bit*/
		{
		 r_bit[i][j] = l_bit[i-1][j] ^ r_bit[i-1][j] ;/*f(R(n-1),k)*/
		}
	}/*end of endes*/
for(i=0;i<=3;i++) e_bit[i] = r_bit[16][i] ;
for(i=4;i<=7;i++) e_bit[i] = l_bit[16][i - 4] ;
/*r_bit + l_bit -> e_bit(64bit)*/

bit2byte(e_bit,e_byte) ;/*change to byte for swap bit IP-1*/
for(i=0;i<=63;i++) e_byte1[i] = e_byte[ip_1[i] - 1] ;/*goto IP-1 swap bit*/
byte2bit(e_byte1,e_bit) ;/*got e_bit*/

}/*end of data encryption*/

/*--------------------------------------------------*/
/*                  data uncryption                 */
/*--------------------------------------------------*/
undes(unsigned char m_bit[8] , unsigned char k_bit[8] , unsigned char e_bit[8])
/*NOTE: in fact , m_bit is encryption data , e_bit is uncryption*/
{
int ip[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,27,19,11,3,
		61,53,45,37,29,21,13,5,
		63,55,47,39,31,23,15,7
		} ;
int ip_1[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 e[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
		} ;
unsigned char m_bit1[8] = {0x0,0x0,0x0,0x0,0x0,0x0,0x0,0x0} ;
unsigned char m_byte[64] ;
unsigned char m_byte1[64] ;
unsigned char key_n[16][8] ;
unsigned char l_bit[17][8] ;
unsigned char r_bit[17][8] ;
unsigned char e_byte[64] ;
unsigned char e_byte1[64] ;
unsigned char l_byte[64] ;
unsigned char l_byte1[64] ;
int i = 0 , j = 0 ;

/*initialize*/
for(i=0;i<=15;i++)
	{
	for(j=0;j<=7;j++)
		{
		l_bit[i][j] = 0x0 ;
		r_bit[i][j] = 0x0 ;
		key_n[i][j] = 0x0 ;
		}
	}
for(i=0;i<=63;i++)
	{
	m_byte[i] = 0x0 ;
	m_byte1[i] = 0x0 ;
	l_byte[i] = 0x0 ;
	l_byte1[i] = 0x0 ;
	e_byte[i] = 0x0 ;
	e_byte1[i] = 0x0 ;
	}

keychange(k_bit,key_n) ;/*get the 48bit key x 16 (16rows x 6byte in key_n)*/

bit2byte(m_bit,m_byte) ;/*change to byte*/
for(i=0;i<=63;i++) m_byte1[i] = m_byte[ip[i] - 1] ;/*goto IP swap bit*/
byte2bit(m_byte1,m_bit1) ;/*re-change to bit*/
for(i=0;i<=3;i++) r_bit[16][i] = m_bit1[i] ;/*move left 32bit -> r16*/
for(i=4;i<=7;i++) l_bit[16][i - 4] = m_bit1[i] ;/*move right 32bit -> l16*/

for(i=16;i>=1;i--) /*->(16) layer from 16 -> 1*/
	{
	for(j=0;j<=3;j++) r_bit[i-1][j] = l_bit[i][j] ;/*R(n-1) = L(n)*/

	/*comput f(L(n),k)*/
	bit2byte(l_bit[i],l_byte) ;
	for(j=0;j<=47;j++) l_byte1[j] = l_byte[e[j] - 1] ;/*goto E swap bit*/
	byte2bit(l_byte1,l_bit[i]) ;/*now r_bit is 48bit*/
	/*xor 48bit key*/
	for(j=0;j<=5;j++) l_bit[i][j] = l_bit[i][j] ^ key_n[i-1][j] ;
	/*goto <S replace> and <P swap bit>*/
	s_replace(l_bit[i]) ;/*change 48bit l_bit[i]->32bit l_bit[i]*/
	for(j=0;j<=3;j++)/*get PREV l_bit*/
		{
		 l_bit[i-1][j] = r_bit[i][j] ^ l_bit[i][j] ;/*f(L(n),k)*/
		}
	}/*end of undes*/
for(i=0;i<=3;i++) e_bit[i] = l_bit[0][i] ;
for(i=4;i<=7;i++) e_bit[i] = r_bit[0][i - 4] ;
/*r_bit + l_bit -> e_bit(64bit)*/

bit2byte(e_bit,e_byte) ;/*change to byte for swap bit IP-1*/
for(i=0;i<=63;i++) e_byte1[i] = e_byte[ip_1[i] - 1] ;/*goto IP-1 swap bit*/
byte2bit(e_byte1,e_bit) ;/*got e_bit*/
/*now ! in fact , the e_bit is the uncryption data.*/

}/*end of uncryption data*/

/*--------------------------------------------------*/
/*                 S_replace                        */
/*--------------------------------------------------*/
s_replace(unsigned char s_bit[8])
{
int p[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
		} ;
unsigned char s1[4][16] = {	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
				} ;
unsigned char s2[4][16] = {	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
				};
unsigned char s3[4][16] = {	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
				};
unsigned char s4[4][16] = {	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
				};
unsigned char s5[4][16] = {	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,
				};
unsigned char s6[4][16] = {	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
				};
unsigned char s7[4][16] = {	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
				};
unsigned char s8[4][16] = {	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
				};
int i=0 , j=0 ;
unsigned char s_byte[64] ;
unsigned char s_byte1[64] ;
unsigned char s_bit_temp[8] = {0x0,0x0,0x0,0x0,0x0,0x0,0x0,0x0} ;
unsigned char row = 0 , col = 0 ;
unsigned char s_out_bit[8] = {0x0,0x0,0x0,0x0,0x0,0x0,0x0,0x0} ;
unsigned char s_out_bit1[8] = {0x0,0x0,0x0,0x0,0x0,0x0,0x0,0x0} ;

/*initialize*/
for(i=0;i<=63;i++) s_byte[i] = 0x0 ;
for(i=0;i<=63;i++) s_byte1[i] = 0x0 ;

/*change 48bit(8bit x 6byte) to 48bit(6bit x 8byte)*/
bit2byte(s_bit,s_byte) ;
for(i=0;i<=7;i++)
	{
        for(j=0;j<=63;j++) s_byte1[j] = 0x0 ;/*clear temp*/

	s_byte1[6] = s_byte[i*6] ;  /*get bit 0 in 6 bit*/
	s_byte1[7] = s_byte[i*6+5] ;/*get bit 5 in 6 bit*/
	byte2bit(s_byte1,s_bit_temp) ;/* 0000 00?? */
	row = s_bit_temp[0] ;/*get row[i]*/

	for(j=0;j<=63;j++) s_byte1[j] = 0x0 ;/*clear temp*/

	s_byte1[4] = s_byte[i*6+1] ;/*0000 ????*/
	s_byte1[5] = s_byte[i*6+2] ;
	s_byte1[6] = s_byte[i*6+3] ;
	s_byte1[7] = s_byte[i*6+4] ;
	byte2bit(s_byte1,s_bit_temp) ;
	col = s_bit_temp[0] ;/*get column in S table*/

	/*get number from S table with row and col*/
	switch (i) {
		case 0 :
			s_out_bit[i] = s1[row][col] ;
			break ;
		case 1 :
			s_out_bit[i] = s2[row][col] ;
			break ;
		case 2 :
			s_out_bit[i] = s3[row][col] ;
			break ;
		case 3 :
			s_out_bit[i] = s4[row][col] ;
			break ;
		case 4 :
			s_out_bit[i] = s5[row][col] ;
			break ;
		case 5 :
			s_out_bit[i] = s6[row][col] ;
			break ;
		case 6 :
			s_out_bit[i] = s7[row][col] ;
			break ;
		case 7 :
			s_out_bit[i] = s8[row][col] ;
			break ;
		} ;


	} /*s_out_bit[0-7]:0000???? 0000???? 0000???? 0000????...0000???? */

	/*change 64bit to 32bit : clean 0000(high 4bit)*/
	s_out_bit1[0] = (s_out_bit[0]<<4) + s_out_bit[1] ;
	s_out_bit1[1] = (s_out_bit[2]<<4) + s_out_bit[3] ;
	s_out_bit1[2] = (s_out_bit[4]<<4) + s_out_bit[5] ;
	s_out_bit1[3] = (s_out_bit[6]<<4) + s_out_bit[7] ;
	/*now s_out_bit1[0-7] = ???????? ???????? ???????? ???????? 0000..*/
        for(i=0;i<=63;i++) s_byte[i] = 0x0 ;
	for(i=0;i<=63;i++) s_byte1[i] = 0x0 ;
	bit2byte(s_out_bit1 , s_byte) ;/*change byte for P swap bit*/
	for(i=0;i<=31;i++) s_byte1[i] = s_byte[p[i] - 1] ;/*goto P swap bit*/
	for(i=0;i<=7;i++) s_bit[i] = 0x0 ;
	byte2bit(s_byte1 , s_bit) ;
	/*now ! we got 32bit f(R,K)*/
}/*end of S_replace*/