UVA10341解方程(二分)
题意:
      给你一个方程 F[x] = pe^-x + qsin(x) + rcos(x) + stan(x) + tx^2 + u = 0(0<=p,r<=20,-20<=q,s,t<=0),给你pqrstu然后问你在定义域内的解是多少0<=x<=1。
思路:
      在定义域内,前面5个都是减函数,所以当F[0] >= 0 && F[1] <= 0时才有唯一解,否则无解,如果有解的话既然是单调的,那么我们就可以二分去求解,二分的时候一开始根据绝对值的大小判断一直得不出答案,最后按照函数小于零的时候就往左移才过。
#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);
    }
}
UVA10341解方程(二分)的更多相关文章
- UVA 10341 Solve It 解方程 二分查找+精度
		
题意:给出一个式子以及里面的常量,求出范围为[0,1]的解,精度要求为小数点后4为. 二分暴力查找即可. e^(-n)可以用math.h里面的exp(-n)表示. 代码:(uva该题我老是出现Subm ...
 - hdu2199Can you solve this equation?(解方程+二分)
		
Can you solve this equation? Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K ( ...
 - hdu2899Strange fuction(解方程+二分)
		
Strange fuction Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)T ...
 - vijos P1915 解方程 加强版
		
背景 B酱为NOIP 2014出了一道有趣的题目, 可是在NOIP现场, B酱发现数据规模给错了, 他很伤心, 哭得很可怜..... 为了安慰可怜的B酱, vijos刻意挂出来了真实的题目! 描述 已 ...
 - HDU 4793 Collision --解方程
		
题意: 给一个圆盘,圆心为(0,0),半径为Rm, 然后给一个圆形区域,圆心同此圆盘,半径为R(R>Rm),一枚硬币(圆形),圆心为(x,y),半径为r,一定在圆形区域外面,速度向量为(vx,v ...
 - codevs3732==洛谷 解方程P2312 解方程
		
P2312 解方程 195通过 1.6K提交 题目提供者该用户不存在 标签数论(数学相关)高精2014NOIp提高组 难度提高+/省选- 提交该题 讨论 题解 记录 题目描述 已知多项式方程: a ...
 - [NOIP2014]解方程
		
3732 解方程 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题解 题目描述 Description 输入描述 Input Descrip ...
 - bzoj 3751: [NOIP2014]解方程 同余系枚举
		
3.解方程(equation.cpp/c/pas)[问题描述]已知多项式方程:a ! + a ! x + a ! x ! + ⋯ + a ! x ! = 0求这个方程在[1, m]内的整数解(n 和 ...
 - Ural 1046  Geometrical Dreams(解方程+计算几何)
		
题目链接:http://acm.timus.ru/problem.aspx?space=1&num=1046 参考博客:http://hi.baidu.com/cloudygoose/item ...
 
随机推荐
- 从JVM底层原理分析数值交换那些事
			
基础数据类型交换 这个话题,需要从最最基础的一道题目说起,看题目:以下代码a和b的值会交换么: public static void main(String[] args) { int a = 1, ...
 - Hi3559AV100 NNIE RFCN开发:V4L2->VDEC->VPSS->NNIE->VGS->VO系统整体动态调试实现
			
下面随笔将给出Hi3559AV100 NNIE RFCN开发:V4L2->VDEC->VPSS->NNIE->VGS->VO系统整体动态调试实现,最终的效果是:USB摄像 ...
 - Prometheus自定义指标
			
1. 自定义指标 为了注册自定义指标,请将MeterRegistry注入到组件中,例如: public class Dictionary { private final List<String ...
 - git的回滚与撤销【reset and revert】
			
git的工作流程-- 3个区域 工作区:我们可以看到的文件内容 在操作 git add 之前的!! 缓存区:是不可见的 已经git add操作,还没git commit -m "" ...
 - P4847 银河英雄传说V2 题解(Splay)
			
题目链接 P4847 银河英雄传说V2 解题思路 我天哪!!!\(splay\)在\(rotate\)的时候先\(upd(y)\)再\(upd(x)\)!!以后不能再因为这个\(WA\)一晚上了!!! ...
 - String 类的内存 解析
			
关于String类的内存解析 Person类的内存解析
 - Day2:Windows常用快捷键与基本的Dos命令
			
Windows常用快捷键 必须掌握: Ctrl+C:复制 Ctrl+V:粘贴 Ctrl+Z:撤销 Ctrl+S:保存 Win键+R:运行(run) alt+F4:关闭窗口/页面 Ctrl+A:全选 C ...
 - 攻防世界 reverse BABYRE
			
BABYRE XCTF 4th-WHCTF-2017 int __cdecl main(int argc, const char **argv, const char **envp) { char ...
 - 攻防世界 reverse 进阶 -gametime
			
19.gametime csaw-ctf-2016-quals 这是一个小游戏,挺有意思的 's'-->' ' 'x'-->'x' 'm'-->'m' 观察流程,发现检验函 ...
 - Ingress-nginx工作原理和实践
			
本文记录/分享 目前项目的 K8s 部署结构和请求追踪改造方案 这个图算是一个通用的前后端分离的 k8s 部署结构: Nginx Ingress 负责暴露服务(nginx前端静态资源服务), 根据十二 ...