UVA10341解方程(二分)
生活随笔
收集整理的這篇文章主要介紹了
UVA10341解方程(二分)
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
題意:
? ? ? 給你一個方程 F[x] = pe^-x + qsin(x) + rcos(x) + stan(x) + tx^2 + u = 0(0<=p,r<=20,-20<=q,s,t<=0),給你pqrstu然后問你在定義域內(nèi)的解是多少0<=x<=1。
思路:
? ? ? 在定義域內(nèi),前面5個都是減函數(shù),所以當F[0] >= 0 && F[1] <= 0時才有唯一解,否則無解,如果有解的話既然是單調(diào)的,那么我們就可以二分去求解,二分的時候一開始根據(jù)絕對值的大小判斷一直得不出答案,最后按照函數(shù)小于零的時候就往左移才過。
#include<math.h>
#include<stdio.h>
#define eps 0.0000000001
double p ,q ,r ,s ,t ,u;
double Fun(double x)
{
? ? return p*exp(-x) + q*sin(x) + r*cos(x) + s * tan(x) + t * x * x + u;
}
double abss(double x)
{
? ? return x > 0 ? x : -x;
}
int main ()
{
? ? while(~scanf("%lf %lf %lf %lf %lf %lf" ,&p ,&q ,&r ,&s ,&t ,&u))
? ? {
? ? ? ? if(Fun(0) < -eps || Fun(1) > eps)
? ? ? ? {
? ? ? ? ? ? printf("No solution\n");
? ? ? ? ? ? continue;
? ? ? ? }
? ? ? ? double low ,up ,mid;
? ? ? ? low = 0 ,up = 1;
? ? ? ? for(int i = 1 ;i <= 50 ;i ++)
? ? ? ? {
? ? ? ? ? ? mid = low + (up - low) / 2;
? ? ? ? ? ? if(Fun(mid) < 0) up = mid;
? ? ? ? ? ? else low = mid;
? ? ? ? }
? ? ? ? printf("%.4lf\n" ,low);
? ? }
}
? ? ? 給你一個方程 F[x] = pe^-x + qsin(x) + rcos(x) + stan(x) + tx^2 + u = 0(0<=p,r<=20,-20<=q,s,t<=0),給你pqrstu然后問你在定義域內(nèi)的解是多少0<=x<=1。
思路:
? ? ? 在定義域內(nèi),前面5個都是減函數(shù),所以當F[0] >= 0 && F[1] <= 0時才有唯一解,否則無解,如果有解的話既然是單調(diào)的,那么我們就可以二分去求解,二分的時候一開始根據(jù)絕對值的大小判斷一直得不出答案,最后按照函數(shù)小于零的時候就往左移才過。
#include<math.h>
#include<stdio.h>
#define eps 0.0000000001
double p ,q ,r ,s ,t ,u;
double Fun(double x)
{
? ? return p*exp(-x) + q*sin(x) + r*cos(x) + s * tan(x) + t * x * x + u;
}
double abss(double x)
{
? ? return x > 0 ? x : -x;
}
int main ()
{
? ? while(~scanf("%lf %lf %lf %lf %lf %lf" ,&p ,&q ,&r ,&s ,&t ,&u))
? ? {
? ? ? ? if(Fun(0) < -eps || Fun(1) > eps)
? ? ? ? {
? ? ? ? ? ? printf("No solution\n");
? ? ? ? ? ? continue;
? ? ? ? }
? ? ? ? double low ,up ,mid;
? ? ? ? low = 0 ,up = 1;
? ? ? ? for(int i = 1 ;i <= 50 ;i ++)
? ? ? ? {
? ? ? ? ? ? mid = low + (up - low) / 2;
? ? ? ? ? ? if(Fun(mid) < 0) up = mid;
? ? ? ? ? ? else low = mid;
? ? ? ? }
? ? ? ? printf("%.4lf\n" ,low);
? ? }
}
總結(jié)
以上是生活随笔為你收集整理的UVA10341解方程(二分)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UVA10294项链和手镯(等价类计数问
- 下一篇: UVA10870递推关系(矩阵乘法)