百練 迭代法解方程 2697

2697:迭代法解方程

• 查看
• 提交
• 統計
• 提示
• 提問

總時間限制:

2000ms

記憶體限制:

65536kB

描述

對函數y = f(x) = x*x*x + x + a（其中a大於0）。要直接求y=0時x的取值並不容易，但因為y是單調遞增的函數，我們可以利用這個性質求x的近似值。取x1＝-a，顯然f(-a)小於0，取x2=0，顯然f(a)大於0。我們規定近似誤差不超過b（比如b=0.001）。這時我們取x=(x2-x1)/2，再判斷(x2-x1)是否小於b，如果小於b，x就是我們要求的值，否則我們需要判斷f(x)的取值，如果f(x)大於0，令x2=x，重新求解；如果f(x)小於0，令x1=x，也重新求解；如果f(x)恰好等於0，則以此x應作為最終x的取值。我們把這種求解的過程稱為迭代法。

輸入

第一行為資料數量n。其餘個行每行包含兩個浮點數a和b。

輸出

輸出為n行，每行輸出對應的x和y的值。

範例輸入

323 0.0123 0.00123 0.0001

範例輸出

-2.726196 0.012314-2.726547 0.004137-2.726767 -0.000974

#include<stdio.h>#include<stdlib.h>int main(){int n,i;double x,b,a;scanf("%d",&n);for(i=0; i<n; i++){scanf("%lf %lf",&a,&b);if(a>0){x = (-a-0)/2;}if(x<0){while(x<0){x=(a-0)/2;}}if(x>0){while(x>0){x=(-x);}}if(x==0 || x<b){printf("%lf %lf\n",x,x*x*x+x+a);}}return 0;}