标签:
题目:函数的渐进解,利用式子f(x+dx)= f(x)+ f`(x)dx,求x+dx;
其中x为f(x)处理求解的整数解,反推dx = (f(x+dx)- f(x))/f`(x)。
分析:分治。这里f(x) = x^3,直接利用二分求解x的最近整数解,然后接dx = (n-a^3)/(3a^2);
整理求解x+dx即可。
说明:这里不能使用库函数,否则精度会出现问题。
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
int bs(double x)
{
int l = 1,r = 100,mid;
while (l < r) {
mid = (r+l+1)/2;
if (x >= 1.0*mid*mid*mid)
l = mid;
else r = mid-1;
}
return r;
}
int main()
{
double x,y,dx;
while (~scanf("%lf",&x) && x) {
y = bs(x);
dx = (x-y*y*y)/3.0/y/y;
printf("%.4lf\n",y + dx);
}
return 0;
}
UVa 10509 - R U Kidding Mr. Feynman?
标签:
原文地址:http://blog.csdn.net/mobius_strip/article/details/43854815