信息学奥赛一本通：1241：二分法求函数的零点

最新推荐文章于 2024-01-12 17:04:16 发布

阿耨多羅

最新推荐文章于 2024-01-12 17:04:16 发布

阅读量260

点赞数

文章标签：算法 c++ c语言

本文链接：https://blog.csdn.net/ffsdfg/article/details/129944163

版权

1241：二分法求函数的零点

时间限制: 1000 ms 内存限制: 65536 KB
提交数: 10342 通过数: 6391

【题目描述】

有函数：f(x)=x5−15x4+85x3−225x2+274x−121

已知f(1.5)>0,f(2.4)<0 且方程f(x)=0 在区间[1.5,2.4] 有且只有一个根，请用二分法求出该根。

【输入】

(无)

【输出】

该方程在区间[1.5,2.4]中的根。要求四舍五入到小数点后6位。

【输入样例】

(无）

【输出样例】

（无）

程序没这么长，能学到这的应都看得懂的，归并算法没用。

只用了主函数和一个f函数

#include<bits/stdc++.h>
using namespace std;
#define MAX 100009
int n,m;
long long sum=0;
long long t=1,a[MAX],b[MAX];
void merge(long long a[],int s,int m,int e)
{
	int pb=0,p1=s,p2=m+1;
	while(p1<=m&&p2<=e)
	{
		if(a[p2]<a[p1])
		{
		b[pb++]=a[p2++];
		}
		else
		b[pb++]=a[p1++];
	}
	while(p1<=m)
	b[pb++]=a[p1++];
	while(p2<=e)
	b[pb++]=a[p2++];
	for(int i=0;i<=e-s;i++)
	a[s+i]=b[i];
 } 
 void mergesort(long long a[],int s,int e)
 {
 	if(s<e)
 	{
 		int m=s+(e-s)/2;
 		mergesort(a,s,m);
 		mergesort(a,m+1,e);
 		merge(a,s,m,e);
	 }
 }
long long numsearch(long long a[],int size,long long t)
 {
 	int l=0,r=size-1;
	 long long lastpos=-1;
 	while(l<=r)
 	{
 		int mid=l+(r-l)/2;
 		if(a[mid]>t)
 		r=mid-1;
 		else
 		{
 			lastpos=mid;
 			l=mid+1;
		 }
	 }
	 return lastpos;
 	 }
 	 long long p,q,x;
double  EPS=1e-7;
double f(double x){
 		return x*x*x*x*x-15*x*x*x*x+85*x*x*x-225*x*x+274*x-121; 
	 }
 int main()
 {
 	double root,x1=1.5,x2=2.4,y;
	 root=x1+(x2-x1)/2;
	 y=f(root);
	 while(fabs(y)>EPS) 
	 {
	 	if(y>0)
	 	x1=root;
	 	else
	 	x2=root;
	 	root=x1+(x2-x1)/2;
	 	y=f(root);
	 }
	 printf("%.6lf",root);
	 	 return 0;
 }