信息学奥赛一本通 1093：计算多项式的值 | OpenJudge NOI 1.5 36:计算多项式的值

【题目链接】

ybt 1093：计算多项式的值
OpenJudge NOI 1.5 36:计算多项式的值

【题目考点】

1. 调用乘方函数pow()（存在于<cmath>中）

double pow(double a, double b); 求$a^b$
pow函数的时间复杂度为$O(logn)$

2. 快速幂

一般循环求$a^n$时间复杂度为$O(n)$，快速幂求$a^n$时间复杂度为$O(logn)$
以下为快速幂算法：
设底数为a，指数为b，求$a^b$

• 设结果变量r初始值为1：int r = 1;，
• 此时r需要乘以$a^b$才能得到最终结果
  • 如果此时b为奇数，那么r *= a;，此时r还需要乘以b-1个a，使b = b - 1。如果b是偶数，则跳过这一步。（由于b是奇数，整除运算 b/2 的值与 (b-1)/2的值是相等的，所以这一步中b = b - 1可以省略）
  • 此时r需要乘以$a^b$，已知$a^b=(a^2{)}^{\frac{b}{2}}$，所以使底数a *= a，指数b /= 2;后，r需要乘的仍然是$a^b$
• 重复上述步骤，直到b为0为止。r即为最终结果

3.秦九韶算法

秦九韶算法可以做到O(n)复杂度求多项式的值

【解题思路】

注意：n能取到的最大值是$10^6$，若使用循环求幂算法，求出每一项的值，然后加和，该算法的复杂度为$O(n^2)$，是无法在1s内解决规模为$10^6$的问题的。因此只能选择其他算法。

【题解代码】

解法1：使用pow函数

该程序的时间复杂度为$O(nlogn)$

#include <bits/stdc++.h>
using namespace std;
int main()
{
    double x, n, s = 1;//s:加和 
    cin >> x >> n;
    for(int i = 1; i <= n; ++i)
        s += pow(x, i);
    cout << fixed << setprecision(2) << s;
    return 0;
}

解法2：使用快速幂

该程序的时间复杂度为$O(nlogn)$

#include <bits/stdc++.h>
using namespace std;
double fastPow(double a, int b)//快速幂求a^b 
{
    double r = 1;
    while(b > 0)
    {
        if(b % 2 == 1)
            r *= a;
        a *= a;
        b /= 2;
    }
    return r;
} 
int main()
{
    double x, n, s = 1;//s:加和 
    cin >> x >> n;
    for(int i = 1; i <= n; ++i)
        s += fastPow(x, i);
    cout << fixed << setprecision(2) << s;
    return 0;
}

解法3：使用秦九韶算法

该程序的时间复杂度为$O(n)$

#include <bits/stdc++.h>
using namespace std;
int main()
{
    double x, s = 1;
    int n;
    cin >> x >> n;
    for(int i = 0; i < n; ++i)
        s = x * s + 1;
    cout << fixed << setprecision(2) << s;
    return 0;
}