10341 - Solve It 三分 + 浮点型函数用法

最新推荐文章于 2023-09-26 22:37:43 发布

原创最新推荐文章于 2023-09-26 22:37:43 发布 · 382 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#include

New idea 专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一种在给定区间内寻找特定复合函数根的方法，通过迭代过程逐步缩小搜索范围，最终准确找到根点。使用了多项数学运算如指数、三角函数和基本算术操作。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <cstdlib>
#include <cmath>
#include <stack>
#define M 1e-9

using namespace std;

int p,q,r,s,t,u;

double f(double x)
{
return p * exp(-x) + q * sin(x) + r * cos(x) + s * tan(x) + t * x * x + u;
}

void Judge()
{
double x = f(0),y = f(1);
if(x > 0 && y > 0 || x < 0 && y < 0){
printf("No solution\n");
return ;
}

if(x == 0){ printf("0.0000\n"); return;}
if(y == 0){ printf("1.0000\n"); return;}

double pre = x,rear = y,mid,midd;
x = 0; y = 1;
while(y-x >= M)
{
double m = x + (y-x)/2;
double mm = m + (y-m)/2;

mid = f(m);
midd = f(mm);
if((mid * rear) > 0){
y = m;
rear = mid;
}
else{
if( (midd * rear) > 0){
y = mm;
rear = midd;
}
else{
x = mm;
pre = midd;
}

}
}

printf("%.4lf\n",x);
}

int main()
{
// freopen("in.in","r",stdin);

while(scanf("%d%d%d%d%d%d",&p,&q,&r,&s,&t,&u) != EOF)
{
Judge();
}

return 0;
}