高精度实现快速幂

因为涉及到乘法，所以我们最多只能压四位。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int M=10000,
		  P=4; 	
struct bignum 
{
	int n[5000],l; 	
	bignum(){l=1,memset(n,0,sizeof(n));}
//-----------------------------------------------
	void init(int q) 	
	{
		n[0]=q;
		l=1;
	}
	void print()
	{
		printf("%d",n[l-1]);
		for(int i=l-2;i>=0;i--)	
		printf("%0*d",P,n[i]);  
		printf("\n");
	}
//-----------------------------------------------
	bignum operator + (bignum x) const 	
	{
		bignum t=*this;
		if(x.l>t.l)t.l=x.l;
		for(int i=0;i<t.l;i++)
		{
			t.n[i]+=x.n[i];
			if(t.n[i]>=M)
			{
				t.n[i+1]+=t.n[i]/M;
				t.n[i]%=M;
			}
		}
		return t;		
	} 
//------------------------------------------------
	bignum operator *(bignum x) const 
	{
		bignum t=*this,tep;
		tep.l=t.l+x.l-1;
		for(int i=0;i<t.l;i++)
		for(int j=0;j<x.l;j++)
		{
			tep.n[i+j]+=t.n[i]*x.n[j];
			if(tep.n[i+j]>=M)
			{
				tep.n[i+j+1]+=tep.n[i+j]/M;
				tep.n[i+j]%=M;
			}
		}
		while(tep.n[tep.l])
		{
			tep.n[tep.l+1]+=tep.n[tep.l]/M;
			tep.n[tep.l++]%=M;
		}
		return tep;
	}
//-----------------------------------------
}a,r;

void power(int ta,int b)  //快速幂
{
	r.init(1);
	a.init(ta);
	while(b)
	{
		if(b&1)r=r*a;
		a=a*a;
		b>>=1;
	}
}

int main()
{
	power(3,200);
	r.print(); 
	cout<<r.l<<endl;
	return 0;
}