信息学奥赛一本通 1170：计算2的N次方 | OpenJudge NOI 1.6 12:计算2的N次方

【题目链接】

ybt 1170：计算2的N次方
OpenJudge NOI 1.6 12:计算2的N次方

【题目考点】

1. 高精度

考察：高精乘低精
高精度计算讲解

2. 快速幂

【解题思路】

先估算结果的位数。指数最大为100，求$2^{100}$的位数。
已知：若正整数x的的位数为n，那么有：$n=\lfloor lgx\rfloor +1$
$2^{100}$的位数为$\lfloor lg{2}^{100}\rfloor +1=\lfloor 100\cdot lg2\rfloor +1\le \lfloor 100\cdot lg10\rfloor +1=101$
可知结果不会超过101位，数字数组长度可以设为105

解法1：累乘

由于每次乘的数字是2，那么可以使用高精对低精乘法，尽量减少算法复杂度，高精乘低精累乘n次即可。

解法2：快速幂

低精度快速幂如下：

int Pow(int a, int b)//求a^b
{
	int r = 1;
	while(b > 0)
	{
		if(b % 2 == 1)
			r *= a;
		a *= a;
		b /= 2;
	}
	return r;
}

本题要求$2^n$，$n\le 100$，n是个低精度数字。
即上述函数中的a、r替换为高精，b为低精。因此只需要实现高精度计算中的高精乘高精。

【题解代码】

解法1：累乘

• 使用数组与函数

#include <bits/stdc++.h>
using namespace std;
#define N 105
void setLen(int a[], int i)
{
    while(a[i] == 0 && i > 1)//去除多余的0
        i--;
    a[0] = i;
}
void Multiply(int a[], int b)//a *= b 高精乘低精
{
    int c = 0, i;
    for(i = 1; i <= a[0]; ++i)
    {
        a[i] = a[i]*b + c;
        c = a[i] / 10;
        a[i] %= 10; 
    }
    while(c > 0)
    {
        a[i] = c % 10;
        c /= 10;
        i++;
    }
    setLen(a, i);
}
void toNum(char s[], int a[])
{
    a[0] = strlen(s);
    for(int i = 1; i <= a[0]; ++i)
        a[i] = s[a[0] - i] - '0';
}
void showNum(int a[])
{
	for(int i = a[0];i >= 1;--i)
		cout << a[i];
	cout << endl;
}
int main()
{
    int a[N] = {1,1}, n;//高精度数字a初始化为1 
    cin >> n;
    for(int i = 1; i <= n; ++i)
        Multiply(a, 2);
    showNum(a);
    return 0;
}

• 使用高精度数字类

#include<bits/stdc++.h>
using namespace std;
#define N 105
struct HPN
{
   int a[N];
   HPN()
   {
       memset(a, 0, sizeof(a));
   }
   HPN(char s[])
   {
       memset(a, 0, sizeof(a));
       a[0] = strlen(s);
       for(int i = 1; i <= a[0]; ++i)
           a[i] = s[a[0] - i] - '0';
   }
   int& operator [] (int i)
   {
       return a[i];
   }
   void setLen(int i)//确定数字位数
   {
       while(a[i] == 0 && i > 1)
           i--;
       a[0] = i;
   }
   void operator *= (int b)//a *= b
   {
       int c = 0, i;
       for(i = 1; i <= a[0]; ++i)
       {
           a[i] = a[i]*b + c;
           c = a[i] / 10;
           a[i] %= 10; 
       }
       while(c > 0)
       {
           a[i] = c % 10;
           c /= 10;
           i++;
       }
       setLen(i); 
   }
   void show()
   {
       for(int i = a[0]; i >= 1; --i)
           cout << a[i];
   }
};
int main()
{
   int n;
   cin >> n;
   HPN a("1");//a初值为1 
   for(int i = 1; i <= n; ++i)
       a *= 2;
   a.show();
   return 0;
}

解法2：快速幂

• 使用数组与函数

#include <bits/stdc++.h>
using namespace std;
#define N 105
void setLen(int a[], int i)
{
    while(a[i] == 0 && i > 1)//去除多余的0
        i--;
    a[0] = i;
}
void numcpy(int a[], int b[])
{
    for(int i = 0; i <= b[0]; ++i)
        a[i] = b[i];
}
void Multiply(int a[], int b[])//a *= b 高精乘高精
{
    int r[N] = {};
    for(int i = 1; i <= a[0]; ++i)
    {
        int c = 0;
        for(int j = 1; j <= b[0]; ++j)
        {
            r[i+j-1] += a[i]*b[j] + c;
            c = r[i+j-1] / 10;
            r[i+j-1] %= 10;
        }
        r[i+b[0]] += c;
    }
    setLen(r, a[0] + b[0]);
	numcpy(a, r);
}
void toNum(char s[], int a[])
{
    a[0] = strlen(s);
    for(int i = 1; i <= a[0]; ++i)
        a[i] = s[a[0] - i] - '0';
}
void showNum(int a[])
{
	for(int i = a[0];i >= 1;--i)
		cout << a[i];
	cout << endl;
}
void fastPower(int a[], int b, int r[])//快速幂 r = a^b 传入前r初始化为0 
{
	int c[N] = {};//c:变化的底数 
	numcpy(c, a);
	r[0] = r[1] = 1;//将r设为1
	while(b > 0)
	{
		if(b % 2 == 1)
			Multiply(r, c);
		Multiply(c, c);
		b /= 2;
	}
}
int main()
{
    int a[N] = {1,2}, r[N] = {}, n;//a:初值为2， r初值为空 
    cin >> n;
    fastPower(a, n, r);
    showNum(r);
    return 0;
}

• 使用高精度计算类

#include<bits/stdc++.h>
using namespace std;
#define N 105
struct HPN
{
    int a[N];
    HPN()
    {
        memset(a, 0, sizeof(a));
    }
    HPN(char s[])
    {
        memset(a, 0, sizeof(a));
        a[0] = strlen(s);
        for(int i = 1; i <= a[0]; ++i)
            a[i] = s[a[0] - i] - '0';
    }
    int& operator [] (int i)
    {
        return a[i];
    }
    void setLen(int i)//确定数字位数
    {
        while(a[i] == 0 && i > 1)
            i--;
        a[0] = i;
    }
    HPN operator * (HPN b)//r = a * b
    {
        HPN r;
	    for(int i = 1; i <= a[0]; ++i)
	    {
	        int c = 0;
	        for(int j = 1; j <= b[0]; ++j)
	        {
	            r[i+j-1] += a[i]*b[j] + c;
	            c = r[i+j-1] / 10;
	            r[i+j-1] %= 10;
	        }
	        r[i+b[0]] += c;
	    }
	    r.setLen(a[0] + b[0]);
		return r;
    }
    HPN operator ^ (int b)//快速幂 求a^b 
    {
    	HPN c = *this, r("1"); 
		while(b > 0)
		{
			if(b % 2 == 1)
				r = r * c;
			c = c * c;
			b /= 2;
		}
		return r;
	}
    void show()
    {
        for(int i = a[0]; i >= 1; --i)
            cout << a[i];
    }
};

int main()
{
    int n;
    cin >> n;
    HPN a("2"), r;
    r = a ^ n;//快速幂求乘方 
    r.show();
    return 0;
}