关于一个计算定积分的程序的bug
yuanfeng2019/10/18数学 IP:湖南

18687.png

请问这个的值是多少?


我需要写一个程序计算这个

我怀疑我积错了

这个程序在大部分情况都可以正常工作,但有时候会出错

#include<bits stdc++.h>
using namespace std;
int main()
{       long double a,b,c,d,l,r;
        //scanf("%0.9f %0.9f %0.9f %0.9f %0.9f %0.9f",&a,&b,&c,&d,&r,&l);
        scanf ("%10Lf %10Lf %10Lf %10Lf %10Lf %10Lf",&a,&b,&c,&d,&r,&l);
  //      printf("%.10f %.10f %.10f %.10f %.10f %.10f",a,b,c,d,r,l);
		long double R=c/a*(a*r+b)+(d-b*c/a)*log(a*r+b);
		long double L=c/a*(a*l+b)+(d-b*c/a)*log(a*l+b);
        printf("%.6Lf\n",(L-R)/a);
        return 0;
}</bits>

比如输入0.1 1 1 1 -1 0的时候是对的

输入0.3 -4 -4 0 14 34时也是对的

但输入-2 3.8 2.5 -9 2 3时就错了


所以应该积错了

大佬帮我积下吧

[修改于 4年6个月前 - 2019/10/18 15:59:36]

来自:数理化 / 数学
24
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
hello1world
4年5个月前 IP:新加坡
865213

从逻辑来看本身没错,在取第三组数时(a*r+b),(a*l+b)为负值,负数再log一下是虚数(复数).

所以你的库能自动处理这情况么?或者你需要自己重新搞个算法.

不是大佬,帮不了你.

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
UIC
4年5个月前 IP:浙江
865225

高中常见的反比例函数平移模型/常数函数

如果是常数,直接秒,return

计算出x方向偏移量,将R和L相对偏移,得到k/x+b形式,如果范围内有复数,单独减掉。其余套公式就A了。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865231
引用hello1world发表于1楼的内容
从逻辑来看本身没错,在取第三组数时(a*r+b),(a*l+b)为负值,负数再log一下是虚数(复数...

是的,是这个问题,谢谢了

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865232
引用hello1world发表于1楼的内容
从逻辑来看本身没错,在取第三组数时(a*r+b),(a*l+b)为负值,负数再log一下是虚数(复数...

这个特判一下就行了吧,我去试试

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
hello1world
4年5个月前 IP:江苏
865234
引用yuanfeng发表于4楼的内容
这个特判一下就行了吧,我去试试

没明白你的意思,如果可以的话希望能把最后成功的程序帖出来看看.

第三组数据2-3区间是有积分结果的,约为3.8455274546965374


前面那个UIC回复我也没看明白,现在高中就有这样的题了么,从函数图形来看类似双曲线旋转加平移的图形.


引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
zx-16533
4年5个月前 修改于 4年5个月前 IP:广东
865235

写成c*((r - l))/a + (a*d - b*c)/a^2*log((a*r + b)/(a*l + b))

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
hello1world
4年5个月前 IP:江苏
865236

仔细想想好象也确实没什么,可以自己判断一下自己写个log负数的功能,而且对于这个函数来看可以象上面radio的那样写成(a*r + b)/(a*l + b),如果(a*r + b)/(a*l + b)为负说明不存在结果(开区间),至少对于这个积分来说是这样.


引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865239
引用hello1world发表于5楼的内容
没明白你的意思,如果可以的话希望能把最后成功的程序帖出来看看.第三组数据2-3区间是有积分结果的,约...

可是按照我积出来的式子来看确实会出现负数的对数的情况,你是怎么积的?


引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865240
引用hello1world发表于7楼的内容
仔细想想好象也确实没什么,可以自己判断一下自己写个log负数的功能,而且对于这个函数来看可以象上面r...

问题是这个积分存在结果

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865241
引用yuanfeng发表于8楼的内容
可是按照我积出来的式子来看确实会出现负数的对数的情况,你是怎么积的?

可以了,正确的代码我发下面

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865242
#include<bits stdc++.h>
using namespace std;
int main()
{       double a,b,c,d,l,r,R,L;
        scanf("%lf%lf%lf%lf%lf%lf", &a, &b, &c, &d, &l, &r);
        if(a!=0)
        {R=c*r/a+(d/a-b*c/(a*a))*log(fabs(a*r+b));
        L=c*l/a+(d/a-b*c/(a*a))*log(fabs(a*l+b));}
        else
        {R=((c*r*r)/(2*b)+(d*r)/b);
        L=((c*l*l)/(2*b)+(d*l)/b);    }
//        cout<<(a*r+b)<<" "<<(a*l+b);         printf("%.6lf\n",(r-l));         return 0; }< pre></(a*r+b)<<" "<<(a*l+b);></bits>


引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865243
引用yuanfeng发表于11楼的内容

这个是正确的代码,特判一下就行了

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
UIC
4年5个月前 IP:浙江
865252
引用hello1world发表于5楼的内容
没明白你的意思,如果可以的话希望能把最后成功的程序帖出来看看.第三组数据2-3区间是有积分结果的,约...

反比例函数(k/x)的学习从初中开始,在高中进行了平移变换。和双曲线转1/4pi是一样的

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
UIC
4年5个月前 IP:浙江
865253

此外,高中连极限也不学,就学导数了。。。积分原理类似,但不学。在各大竞赛有所涉及。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
hello1world
4年5个月前 IP:江苏
865255
引用UIC发表于14楼的内容
此外,高中连极限也不学,就学导数了。。。积分原理类似,但不学。在各大竞赛有所涉及。

这确实是和k/x一样,当初好象不那么重视这个,不大记得了,好象当初高中讲极限,和行列式(矩阵)运算但不讲导数和积分当然那时也没什么竞赛.

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
hello1world
4年5个月前 IP:江苏
865256
引用yuanfeng发表于11楼的内容

基本就这样了,不过a!=0时要判断一下a*r+b!=0,a*l+b!=0,并且a*r+b和a*l+b符号相同(都是正数或都是负数),此时存在积分结果.然后进行下面的计算.

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865265
引用hello1world发表于15楼的内容
这确实是和k/x一样,当初好象不那么重视这个,不大记得了,好象当初高中讲极限,和行列式(矩阵)运算但...

我是搞信竞的,要学这个

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
UIC
4年5个月前 IP:浙江
865280
引用yuanfeng发表于17楼的内容
我是搞信竞的,要学这个

你是哪里的?我考完初赛,准备退役。。。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865315
引用UIC发表于18楼的内容
你是哪里的?我考完初赛,准备退役。。。

湖南,话说你过了为什么不搞复赛

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
UIC
4年5个月前 IP:浙江
865326
引用yuanfeng发表于19楼的内容
湖南,话说你过了为什么不搞复赛

目前不知道分数线,应该没过。

而且我要专心搞物理

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865328
引用UIC发表于20楼的内容
目前不知道分数线,应该没过。而且我要专心搞物理

50以上就够吧。。。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
UIC
4年5个月前 修改于 4年5个月前 IP:浙江
865329

我浙江的。。这儿都贼高,平均都有68

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
yuanfeng作者
4年5个月前 IP:湖南
865353
引用UIC发表于22楼的内容
我浙江的。。这儿都贼高,平均都有68

我67.。。。。

引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

想参与大家的讨论?现在就 登录 或者 注册

所属专业
上级专业
同级专业
文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}