当前位置: 首页 > news >正文

题解:P12037 [USTCPC 2025] 数学分析

Solution

1. 记号

\[I_i = \int^1_0\frac{x^i}{x^2 + 1}\,dx \]

则易得

\[\begin{align*} I_0 &= \frac{\pi}{4} \\ I_1 &= \frac{\ln 2}{2} \end{align*} \]

\(I_i\) 这个式子不好直接积分,我们想办法化成我们会的样子。

2. 拆出来

\[\begin{align*} x^i &= x^{i - 2} \cdot x^2 \\ &= x^{i - 2}\cdot(x^2 + 1) - 1 \\ &= x^{i - 2}(x^2 + 1) - x^{i - 2} \end{align*} \]

3. 代回去

将这个式子带入上式,得

\[I_i = \int^1_0\frac{x^{i - 2}(x^2 + 1) - x^{i - 2}}{x^2 + 1}\,dx \]

发现了什么?中间可以拆开。把这个式子拆成两个式子,得

\[\begin{align*} I_i &= \int^1_0\frac{x^{i - 2}(x^2 + 1)}{x^2 + 1}\,dx - \int^1_0\frac{x^{i - 2}}{x^2 + 1}\,dx \\ &= \int^1_0x^{i - 2} - \int^1_0\frac{x^{i - 2}}{x^2 + 1} \\ \end{align*} \]

您猜怎么着,咱这表一查,前边儿那个式子的公式就知道了(想必学过积分的同学们这个公式都滚瓜烂熟了):

\[\int x^\alpha\,dx = \frac{x^{\alpha + 1}}{\alpha + 1} + C \]

您还记得咱最开始 \(I_i\) 的定义嘛?原式减号后边儿这个式子 \(\displaystyle\int^1_0\frac{x^{i - 2}}{x^2 + 1}\) 不就是咱的 \(I_{i - 2}\) 嘛!

所以一切就简单了:

\[\begin{align*} \int^1_0x^{i - 2} - \int^1_0\frac{x^{i - 2}}{x^2 + 1} &= \left.\frac{x^{i-1}}{i-1}\right|_{0}^{1} - I_{i - 2}\\ &= \frac{1}{i - 1} - I_{i - 2} \end{align*} \]

您瞅瞅,这递推式子不就出来了嘛!

别忘了咱现在这式子不是最终的答案,把这玩意儿算出来以后得乘上它对应的系数 \(a_i\)

啥?您问这 C++ 里边儿 \(\displaystyle\frac \pi 4\)\(\ln 2\) 咋算?这 cmath 头文件里边儿有俩函数,叫做 acos(x)log(x),这 acos(x) 就是 \(\arccos x\),咱都知道 \(\arccos(-1) = \pi\),所以这 C++ 里边儿您用 acos(-1),它就是 \(\pi\) 了。这 log(x) 就更直接了,就是 \(\log_e x\),即 \(\ln x\)

Code

#include <bits/stdc++.h>
using namespace std;typedef long double ld;
typedef long long ll;
const ld PI = acosl(-1.0), ln2 = logl(2.0);ll n, a[(int)1e5 + 5];
ld f[(int)1e5 + 5], ans;
int main() {cin >> n;for (int i = 0; i <= n; ++i) cin >> a[i];f[0] = 0.25 * PI;f[1] = 0.5 * ln2;for (int i = 2; i <= n; ++i) {f[i] = (1.0 / (i - 1.0)) - f[i - 2];ans += f[i] * a[i];}ans += f[0] * a[0] + f[1] * a[1];printf("%.12Lf", ans);	return 0;
}
http://www.gsyq.cn/news/24318.html

相关文章:

  • 8.16
  • 一文读懂隔离见证
  • 12131
  • 12131
  • FFmpeg开发笔记(八十四)使用国产的librestreaming实现RTMP直播
  • 2025 年护栏板厂家企业品牌推荐排行榜,波形,高速,镀锌,二波,三波,喷塑,国标,绳索,公路护栏板,护栏板立柱公司推荐
  • 2025 年船用锅炉厂家企业品牌推荐排行榜,基于市场口碑,评选值得信赖的船用锅炉公司推荐
  • 2025 年耐火砖厂家企业品牌推荐排行榜,绝热,轻质,莫来石,保温,莫来石轻质,氧化铝泡沫,氧化铝空心球,抗渗碳,高温轻质莫来石,高温耐火砖公司推荐
  • 2025 年启闭机厂家企业品牌推荐排行榜,四川启闭机,四川卷扬启闭机,四川螺杆启闭机,固定卷扬启闭机,手电两用螺杆启闭机,电装启闭机公司推荐
  • 2025 年铸铁闸门厂家企业品牌推荐排行榜,四川铸铁闸门,镶铜铸铁闸门,渠道铸铁闸门,圆形铸铁闸门,方形铸铁闸门公司推荐
  • Java继承的实现
  • 「学习笔记」Linux 学习笔记(持续更新中)
  • java入门代码示例
  • 极端的利己主义与虚伪的利他主义
  • WinRAR-7.13-Final-x64-烈火汉化版
  • 2025年羽绒服厂家推荐排行榜,轻薄保暖羽绒服,时尚羽绒服,户外运动羽绒服公司推荐!
  • Windows端口查看查找exe进程名
  • 12. 模型RAG评测 - Rainbow
  • 2025年卫衣厂家推荐排行榜,春秋冬款卫衣,男女同款卫衣,运动休闲卫衣公司推荐!
  • 2025年10月运动木地板厂家推荐:对比评测榜单助力场馆精准选材
  • 2025年羽绒服厂家推荐排行榜,保暖轻盈时尚羽绒服,冬季必备高品质羽绒服公司推荐!
  • 2025年10月美白精华推荐榜:OLAY水光小白瓶领衔对比评测排行
  • 2025年10月智能体公司推荐:行胜领衔五强榜对比评测
  • 完整教程:UniApp manifest.json 配置文件完全解析
  • 2025年10月北京金融街附近豪华酒店综合对比与排行榜单
  • 2025年10月金融街附近豪华酒店对比排行与评测分析
  • 大模型基础补全计划(五)---seq2seq实例与测试(编码器、解码器架构)
  • 2025年10月连锁酒店排行榜前十名深度对比与投资价值分析
  • 2025年10月亲子旅游景区综合对比与排行榜解析
  • 2025年角接触轴承厂家推荐排行榜,高精度/高承载/高精密球/机床主轴/汽车/定制/可替代进口/高转速/高刚性轴承公司推荐!