牛顿迭代法-朋友Z与方程

题目链接：http://115.231.222.240:8081/JudgeOnline/problem.php?id=1470；

1470: 朋友Z与方程

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 1131   Solved: 246
[ Submit][ Status][ Web Board]

Description

在朋友Z的信中她问了L先生一道题目，可惜L先生太笨，你能帮助他解决吗？

已知等式：6*x^3 + 5*x^2 + 7*x + 9 == Y 给出Y的值请求出x( 0≤x≤10）

Input

输入有多组测试数据，处理到文件结束（测试数据数量<=100000），每组测试数据给出一个浮点数Y(fabs(Y)<10 5 )。

Output

每行输出一个答案。

对于每个测试例如果存在满足条件的x，那么就输出x的值，并精确到小数点后4位，否则输出No solution! 。

Sample Input

100
-10

Sample Output

2.0877
No solution!

HINT

eps<=10-8

思路：牛顿迭代法；

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define eps 0.000001

double f(double x)
{
    return 6*x*x*x+5*x*x+7*x+9;
}

int main()
{
    double y;
    while(~scanf("%lf",&y))
    {
        double l=0,r=10;
        if(f(0)>y+eps||f(10)<y-eps)
            puts("No solution!");
        else
        {
            for(int i=0; i<50; i++)
            {
                double mid=(l+r)/2;
                if(f(mid)>y)
                    r=mid;
                else
                    l=mid;
            }

            printf("%.4f\n",l);
        }
    }
    return 0;
}

关于牛顿迭代法：

参考：http://www.matongxue.com/madocs/205.html#/madoc；

上面的参考说的很详细；

以下是我根据参考的理解...

给一个函数然后求函数的零点；

在函数上随便选取一个点(xn,f(xn))，则这一点上的切线方程为y=f'(xn) (x-xn) + f(n);那么切线与x轴的交点就是x(n+1)=xn-f(xn)/f'(xn) ;

在根据(xn+1,f(xn+1))这个点找下一个近似根，那么到什么时候为止呢，就是前一个近似根与这个近似根的差接近为0时！

代码可以根据方程具体实现；

本题时给你了确定的y然后求x，可以用二分来求得x；

根据已知得方程可知这个函数时递增得，所以判断无解得条件是：

if(f(0)>y+eps||f(10)<y-eps)