快速幂，不多说了。

　　#include

　　#include

　　int a[15001],b[15001],c[15001],n,p;

　　void copy(int c[],int a[])

　　{

　　for(int i=0;i<=c[0];i++)

　　a[i]=c[i];

　　}

　　void mutil(int a[],int b[],int c[])

　　{

　　for(int i=0;i<=a[0]+b[0];i++)

　　c[i]=0;

　　for(int i=1;i<=a[0];i++)

　　{

　　for(int j=1;j<=b[0];j++)

　　{

　　c[i+j-1]+=a[i]*b[j];

　　}

　　}

　　c[0]=a[0]+b[0]-1;

　　for(int i=1;i<=c[0];i++)

　　{

　　c[i+1]+=c[i]/10;

　　c[i]%=10;

　　}

　　if(c[c[0]+1]>0)

　　c[0]++;

　　}

　　int main()

　　{

　　freopen("cruel1.in","r",stdin);

　　freopen("cruel1.out","w",stdout);

　　scanf("%d%d",&n,&p);

　　c[0]=1;

　　c =n;

　　while(c[c[0]]>=10)

　　{

　　c[0]++;

　　c[c[0]]=c[c[0]-1]/10;

　　c[c[0]-1]%=10;

　　}

　　copy(c,a);

　　b[0]=1;

　　b =1;

　　while(p>0)

　　{

　　if(p%2!=0)

　　{

　　mutil(a,b,c);

　　copy(c,b);

　　}

　　mutil(a,a,c);

　　copy(c,a);

　　p/=2;

　　}

　　int counter=0;

　　for(int i=b[0];i>=1;i--)

　　{

　　if(counter==70)

　　{

　　printf("\n");

　　counter=0;

　　}

　　counter++;

　　printf("%d",b[i]);

　　}

　　fclose(stdin);

　　fclose(stdout);

　　return 0;

　　}