HDU 2199 牛顿迭代

HDU 2199

题目链接：

http://acm.hdu.edu.cn/showproblem.php?pid=2199

题意：

求一个给定高次方程的根。

思路：

练习牛顿迭代。至于为什么要遍历所有整数点来取根暂且不知，精度的取法和退出方式好像也有讲究。

源码：

#include <cstdio>

#include <cstring>

#include <cmath>

#include <cstdlib>

#include <algorithm>

#include <iostream>

#include <queue>

using namespace std;

const double eps = 1e-6;

double Newton_iterative(double a, double b, double c, double d, double e, double x0)

{

    int cnt = 0;

    double x = x0;

    while(1){

        x0 = x;

        long double f1 = ((((a * x + b) * x + c) * x + d) * x + e);

        long double f2 = (((4 * a * x + 3 * b) * x + 2 * c) * x + d);

        x = x0 - f1 / f2;

        cnt++;

        if(fabs(f1) < eps)

            break;

        if( cnt > 50)

            return -1;

    }

    return x;

}

int main()

{

    double a, b, c, d, e;

    a = 8, b = 7, c = 2, d = 3;

    int t;

    scanf("%d", &t);

    while(t--){

        double lv;

        scanf("%lf", &lv);

        e = 6 - lv;

        int ok = 0;

        double u = 0;

        double x;

        for(u = 0.0 ; u < 100.0 ; u++){

            x = Newton_iterative(a, b, c, d, e, u);

            if(x >= 0 && x <= 100){

//              printf("u = %f\n", u);

                ok = 1;

                break;

            }

        }

        if(ok)

            printf("%.4f\n", x);

        else

            printf("No solution!\n");

    }

    return 0;

}