关于Diffie-Hellman算法

//-----------------NO1 求X^n次方------------------------

//纯C跑就小地方改改就可以跑了

#include <iostream>

#include <algorithm>



using namespace std;



int tvalue;

int effectbit;

#define  N 100

#define  ENDPOS N-1

int temp[N];

int result[N];

void moveonebit()

{

	int newpos,endpos,prevpos;

	endpos = ENDPOS;

	newpos = endpos-temp[0]; //此位置位进位后的数开始位置 

	prevpos = newpos+1; //此位置为原temp中数开始的位置

	for(int loopcnt = temp[0]; loopcnt>0; loopcnt--,newpos++,prevpos++)

	{

		temp[newpos] = temp[prevpos];

	}

	temp[endpos] = 0;

	temp[0] += 1;

}



void multiply(int bitvalue)//处理乘法进位

{

	int bitlen = temp[0];

	if(bitlen>ENDPOS) exit(-1); //进位到temp[0]的标志位时，说明长度不够，出现异常



	int newcnt = 0;    //从新统计目标数组的有效位数

	int cur;

	int bkvalue=0;



	for(int i=ENDPOS; bitlen--; i--)//此处bitlen该为temp的有效位数

	{	

		cur = i;

		result[i] += temp[i] * bitvalue;



		while(result[cur]/10 && cur>0)

		{

			result[cur-1] += result[cur] / 10;

			result[cur] = result[cur] % 10;

			cur --; //从尾部向前移

		}

	}

	newcnt = ENDPOS-cur+1;

	result[0] = newcnt;

}



int* Power(int arr[], int n)

{	

	if(1 == n)

	{

		int effectnum = effectbit;

		int elen = ENDPOS;

		while(effectnum-- > 0)

		{

			temp[elen] = arr[elen];

			elen --;

		}

		return temp;

	}

	else

	{   

		Power(arr,--n);

		int  len = ENDPOS;

		

		int elen,ilen;

		elen = ilen = len;

		int  loopcnt = effectbit;

		

		while(loopcnt--)//此循环次数位输入数组的有效位数

		{               //需要循环每一位让其与目标数组的每一位进行乘法进位

			multiply(arr[len]); //将这值

			if(loopcnt > 0)

			{

				moveonebit(); //temp 进一位

			}

			len --;	

		}

		memcpy(temp,result,N*4); //更新下一次递归的temp值

		memset(result,0,N*4);    //将result清空,不清空则会与旧值累加，变成就阶乘的情况

		return temp;		

	}

	

}



int main()

{

	int array[N],n;

	cout << "calculate X^n!\n";

	cout << "input X :";

	char cNumber='\0';



	scanf("%c",&cNumber);

	int icnt = 0;

	int end = N-1;

	while(cNumber != '\n' && icnt++ < N)

	{

		array[end--] = cNumber - '0'; 

		scanf("%c",&cNumber);

	}

	if(icnt > 1) //处理 输入123 保存是从后保存的，所以变成321，现要在数组中改成实际输入状况

	{

		int tswap=0;

		int ix , jx;

		int wid = icnt;

		end = ENDPOS;

		for(ix=jx= ENDPOS-wid+1; (wid)/2 && ix<=ENDPOS ; ix++,end--,wid-=2)

		{

			tswap = array[ix];

			array[ix] = array[end];

			array[end] = tswap;

		}

	}



	cout << "input n :";

	cin >> n;

	if(n<=0)

	{

		cout << "n value must greate 0" <<endl;

		exit(-1);

	}



	int emp = N - icnt;

	for(int i=0; emp--; i++) //对空位进行标识

	{

		array[i] = -1;

	}

	

	for(i=1; i<N; i++)

	{

		temp[i] = 0;

	}



	array[0] = icnt; 

	temp[0]  = icnt;

        effectbit = icnt;

	

	Power(array,n);



	cout << endl << "X^n result is:" ;

	memcpy(result,temp,N*4);

	copy(result+(ENDPOS-result[0]+1),result+N,ostream_iterator<int>(cout," "));

	cout << endl;

	system("pause");

	return 0;

}