【PAT】3-1 二分法求多项式单根-白红宇

【PAT】3-1 二分法求多项式单根

阅读量：6160 次

发布时间：2019-06-21

本文共 976 字，大约阅读时间需要 3 分钟。

二分法的步骤为：

检查区间长度，如果小于给定阈值，则停止，输出区间中点(a+b)/2；否则

如果f(a)f(b)<0，则计算中点的值f((a+b)/2)；

如果f((a+b)/2)正好为0，则(a+b)/2就是要求的根；否则

如果f((a+b)/2)与f(a)同号，则说明根在区间[(a+b)/2, b]，令a=(a+b)/2，重复循环；

如果f((a+b)/2)与f(b)同号，则说明根在区间[a, (a+b)/2]，令b=(a+b)/2，重复循环；

#include 
     
      #include 
      
       #include 
       
        #include 
        
         #include 
         
          #include 
          
           #include 
           
            #include 
            
             #include 
             
              using namespace std;#define read() freopen("data.in", "r", stdin)#define write() freopen("data.out", "w", stdout)#define clr( a , x ) memset ( a , x , sizeof a ) #define cpy( a , x ) memcpy ( a , x , sizeof a ) #define _max(a,b) ((a>b)?(a):(b))#define _min(a,b) ((a
              
               >a3>>a2>>a1>>a0>>a>>b; while(b-a >= 0.001)//因为精度要求是两位，那就保证最后的解在小数点后第三位精确就行 { mid = (a+b)/2; if (f(mid) == 0) { break; }else if (f(a)*f(mid) < 0) { b = mid; }else { a = mid; } } mid = (a+b)/2; printf("%.2f\n",mid ); return 0; }