UVa 10428 The Roots,牛顿迭代法

UVa 10428 The Roots


题意:
给定一个一元n次方程,求出所有解。(由代数基本定理可知n次一元多项式总是有n个根)


题解:
牛顿迭代法
通过牛顿迭代法每次都可以求出多项式的一个根,即方程的一个解x',
然后由 n次一元多项式 =(x-x')*(n-1次一元多项式)
计算出这个 n-1次一元多项式,即计算其系数
令 bn+1 = 0
bn = bn+1 * x' + an

bn-1 = bn*x' + an-1

........

b2 = b3*x' + a2

b1 = b2*x' + a1

b0 = b1*x' + a0

得这个n-1次一元多项式的系数列表 {bn, bn-1,...,b1}。


然后重复上面的操作,就可以把一元n次方程的n个解都找出来。





#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdio>
#include<map>
#include<string>
using namespace std;

int N;
double a[10], b[10], fd[10];


//f(x) = gn*x^n + gn-1 * x^n-1 + ... + g0
double Function(double* g, double x, int n) {
    double ret = 0, u = 1;
    for(int i=0; i<=n; ++i) {
        ret += g[i]*u;
        u *= x;
    }
    return ret;
}

double Newton(int n) {
    for(int i=0; i<n; ++i)
        fd[i] = a[i+1] * (i+1);  //求导
    double x = -100;  //根的范围 abs(root)<25
    for(int i=0; i<100; ++i) //迭代
        x = x - Function(a, x, n) / Function(fd, x, n-1);
    return x;
}

void NewPoly(double x, int n) {
    b[n+1] = 0;
    for(int i=n; i>=0; --i)
        b[i] = b[i+1]*x + a[i];

    for(int i=0; i<n; ++i) {
        a[i] = b[i+1];
    }
}

void solve() {
    for(int i=0; i<N; ++i) {
        double x = Newton(N-i);
        printf(" %.4f", x);
        NewPoly(x, N-i);
    }
}
int main() {
    int cas = 1;
    while(~scanf("%d",&N)&&N) {
        for(int i=N; i>=0; --i)
            scanf("%lf", &a[i]);
        printf("Equation %d:",cas++);
        solve();
        printf("\n");
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值