Vijos P1131 最小公倍数和最大公约数问题【暴力】

一元三次方程求解

描述

有形如:ax^3+bx^2+cx+d=0 这样的一个一元三次方程。给出该方程中各项的系数(a,b,c,d 均为实数),并约定该方程存在三个不同实根(根的范围在-100至100之间),且根与根之差的绝对值>=1。要求由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后2位。

格式

输入格式

输入该方程中各项的系数(a,b,c,d 均为实数),

输出格式

由小到大依次在同一行输出这三个实根(根与根之间留有空格),并精确到小数点后2位。

样例1

样例输入1

1 -5 -4 20

样例输出1

-2.00 2.00 5.00

限制

每个测试点1s

提示

提示:记函数f(x),若存在2个实数x1和x2,且x1<x2,使f(x1)*(x2)<0,则在(x1, x2)之间一定存在实数x0使得f(x0)=0。

来源

NOIP2001第一题

题目链接:https://vijos.org/p/1116

分析:又来一道暴力题,听说这是一道很经典的二分?贴个二分的代码啊,暴力肯定可以,但是但是,一定要注意精度QAQ

暴力代码:【一个三元函数的性质,不懂的自己翻高中课本】

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 int main()
 4 {
 5     double a,b,c,d;
 6     cin>>a>>b>>c>>d;
 7     double j=-100;
 8     double i=j;
 9     for(j;j<=100;j+=0.01)
10     {
11         if(a*i*i*i+b*i*i+c*i+d<0)
12         {
13             if(a*j*j*j+b*j*j+c*j+d>0)
14             {
15                 printf("%.2lf ",j);
16                 i=j;
17             }
18         }
19         if(a*i*i*i+b*i*i+c*i+d>0)
20         {
21             if(a*j*j*j+b*j*j+c*j+d<0)
22             {
23                 printf("%.2lf ",j);
24                 i=j;
25             }
26         }
27     }
28     printf("\n");
29     return 0;
30 }

直接贴了二分代码:

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <algorithm>
 4 #include <cstring>
 5 #include <iomanip>
 6 #include <cstdlib>
 7 using namespace std;
 8 
 9 float ans[5];
10 float a,b,c,d;
11 int n=0;
12 
13 float f(float x)
14 {
15     return ((a*x+b)*x+c)*x+d;
16 }
17 
18 void solve(float l,float r)
19 {
20     if(f(l)*f(r)>0&&(((r-l)<1)||n>=2))
21         return;
22     float mid=(l+r)/2;
23 
24     if(f(mid)<=1e-4 && f(mid)>=-1e-4)
25     {
26         ans[n++]=mid;
27         return;
28     }
29 
30     solve(l,mid),solve(mid,r);
31 }
32 
33 int main()
34 {
35     cin>>a>>b>>c>>d;
36     solve(-100,100);
37     printf("%.2lf %.2lf %.2lf",ans[0],ans[1],ans[2]);
38     return 0;
39 }

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏章鱼的慢慢技术路

Direct3D 11 Tutorial 7:Texture Mapping and Constant Buffers_Direct3D 11 教程7:纹理映射和常量缓冲区

在上一个教程中,我们为项目引入了照明。 现在我们将通过向我们的立方体添加纹理来构建它。 此外,我们将介绍常量缓冲区的概念,并解释如何使用缓冲区通过最小化带宽使用...

954
来自专栏Petrichor的专栏

tensorflow编程: Neural Network

对 features 和 tf.negative(features) 分别 Relu 并 concatenate 在一起。

1773
来自专栏python3

python语句-for

722
来自专栏小樱的经验随笔

51 Nod 1007 正整数分组【类01背包】

1007 正整数分组 基准时间限制:1 秒 空间限制:131072 KB 分值: 10 难度:2级算法题 将一堆正整数分为2组,要求2组的和相差最小。 例如:1...

2867
来自专栏mwangblog

整数的存储:无符号表示法

1273
来自专栏Script Boy (CN-SIMO)

Java中随机数的产生方式与原理

查阅随机数相关资料,特做整理 首先说一下java中产生随机数的几种方式 在j2se中我们可以使用Math.random()方法来产生一个随机数,这个产生的随机数...

2650
来自专栏Python中文社区

实现属于自己的TensorFlow(一) - 计算图与前向传播

前言 前段时间因为课题需要使用了一段时间TensorFlow,感觉这种框架很有意思,除了可以搭建复杂的神经网络,也可以优化其他自己需要的计算模型,所以一直想自...

4237
来自专栏我的博客

PHP取余的那些事

1、百分号取余 $val=9.45; $result=$val*100; echo intval($result); //这里输出944 echo $re...

36710
来自专栏专注数据中心高性能网络技术研发

HERD--位运算

判断一个数是否是2的次方 1 static inline int hrd_is_power_of_2(uint32_t n) 2 { 3 retur...

3579
来自专栏AI研习社

深度学习对话系统实战篇 -- 简单 chatbot 代码实现

本文的代码都可以到我的 github 中下载:https://github.com/lc222/seq2seq_chatbot 前面几篇文章我们已经介绍了 s...

8988

扫码关注云+社区