c++ 如何有效地从其根部找到多项式的系数?

栏目: C++ · 发布时间: 7年前

内容简介:http://stackoverflow.com/questions/33067755/how-to-efficiently-find-coefficients-of-a-polynomial-from-its-roots

参见英文答案 > Sum of multiplication of all combination of m element from an array of n elements 3

给定的是主导系数为1的多项式的n根.

如何有效地找出这个多项式的系数?

在数学上,

我知道如果第一个系数是1,那么一次取k产品根的总和将是多项式的第k个第一个系数.

我的代码是基于这种方法.

换句话说,如何从一次列出的列表中最佳地找到数字乘积的总和.

int main()
{

    int n, sum;
    cin >> n;
    int a[n];
    for (int i=0; i<n; i++) cin >> a[i];
    //for the second coefficient
    sum=0;
    for (int i=0; i<n; i++)
    {
        sum+=a[i];
    }
    cout << sum << endl;
    //for the third coefficient
    sum=0;
    for (int i=0; i<n; i++)
    {
        for (int j=i+1; j<n; j++)
        {
            sum+=a[i]*a[j];
        }
    }
    cout << sum << endl;
}

我想到的是,为了更高系数的目的,我是否已经把它们加入到产品中,但是没有编写代码,因为如果多项式的程度很大,实际上是没有用的.

你需要计算线性因子的乘积

(x-z1)·(x-z2)·…·(x-zn)

这可以通过将多项式与线性因子重复相乘来进行感应地实现

(a[0]+a[1]·x+…+a[m-1]·x^(m-1))·(x-zm) 
= (-a[0]·zm) + (a[0]-a[1]·zm)·x + … + (a[m-2]-a[m-1]·zm) ·x^(m-1) + a[m-1]·x^m

就这样可以实现为循环

a[m] = a[m-1]
for k = m-1 downto 1
    a[k] = a[k-1] - a[k]*zm
end
a[0] = -a[0]*zm

这给出了对所有n个线性因子的乘法的总和n 2/2乘法和相似数量的减法.

http://stackoverflow.com/questions/33067755/how-to-efficiently-find-coefficients-of-a-polynomial-from-its-roots


以上就是本文的全部内容,希望本文的内容对大家的学习或者工作能带来一定的帮助,也希望大家多多支持 码农网

查看所有标签

猜你喜欢:

本站部分资源来源于网络,本站转载出于传递更多信息之目的,版权归原作者或者来源机构所有,如转载稿涉及版权问题,请联系我们

谋局者

谋局者

何常在 / 北京联合出版公司 / 2017-1 / 39.80

★商战版《官场笔记》!全面超越《问鼎》《交手》!商战小说*大神何常在迄今为止至为满意之作! ★以马云、马化腾、李彦宏、雷军、刘强东、张朝阳等大佬为原型,写透高手们的大智慧、大手腕、大谋略! ★善谋者胜,善算者赢!内含大量阳谋诡计、商业运作、商业谈判、事件营销等可以读以致用的知识!是商界人士必看读物! ★全景再现互联网三大帝国七大诸侯从无到有从有到强从强到吞并一切的成长和并购史! ......一起来看看 《谋局者》 这本书的介绍吧!

随机密码生成器
随机密码生成器

多种字符组合密码

MD5 加密
MD5 加密

MD5 加密工具

html转js在线工具
html转js在线工具

html转js在线工具