霍纳法则（Horner Rule）介绍及C语言实现

参考自：http://flynoi.blog.hexun.com/31272178_d.html

霍纳法则简介

假设有n＋2个实数a₀，a₁，…，a_n,和x的序列，要对多项式P_n(x)= a_nxⁿ＋a_n－１x^n－１＋…＋a_１x＋a_０求值，直接方法是对每一项分别求值，并把每一项求的值累加起来，这种方法十分低效，它需要进行n＋(n－1)＋…＋1＝n(n＋１)/2次乘法运算和n次加法运算。有没有更高效的算法呢?答案是肯定的。通过如下变换我们可以得到一种快得多的算法，即P_n(x)= a_nxⁿ ＋a_n－１x^n－１＋…＋a_１x＋a_０＝((…(((a_nx ＋a_n－１)x＋a_n－２)x+ a_n－3)…)x＋a₁)x＋a_０，这种求值的安排我们称为霍纳法则。

霍纳法则C语言实例

/* 霍纳算法实例 */

#include <stdio.h>

long int
horner(int coefficient[], int n, int x) //coefficient[]为待求多项式的系数数组，n为数组大小，x为多项式中未知数x的具体值
{                                       //注意：coefficient[0]存放系数a0，coefficient[1]存放系数a1，以此类推…
    int i;
    long int result;

    result = coefficient[n-1];
    for(i = 1; i <= n-1; i++)
    {
        result = result * x + coefficient[n-1-i];
    }

    return result;
}

int
main(void)
{
    long int p;
    int a[4] = {3, 2, 1, 1};
    p = horner(a, 4, 1);
    printf("polynomial x^3 + x^2 + 2x + 3 = %ld\n", p);
}

实例测试结果：