当前位置: 首页 > 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.hskmm.com/?act=detail&tid=34237

相关文章:

  • 题解:P10514 考试
  • 华为昇腾笔记之Mindspeed-LLM 中 MoE 实现机制与重写逻辑总览
  • 实时时序上下文推荐系统获KDD最佳论文奖
  • 题解:CF1010A Fly
  • 2025年精密磨床/CNC机械加工厂家推荐排行榜,覆盖铣床/车床/磨削/多轴/复合加工,专业非标定制服务首选!
  • 题解:CF1914F Programming Competition
  • 独立开发者找蓝海:新词引流实战
  • 使用云服务器搭建飞牛Frp 内网穿透服务
  • Luogu P14255 列车(train) 题解 [ 蓝 ] [ 线段树 ] [ 二维平面转化 ]
  • 使用VS2022和Unity时可能出现的问题总结
  • 2025年喷雾机器人,取件机器人,工业机器人厂家权威推荐榜单:智能高效与稳定性能的行业首选!
  • 2025年给汤机厂家推荐排行榜,高效节能/智能控制/稳定耐用的优质品牌选择!
  • 理想完美主义者的宣战:当一人面对整个时代的“合理”谎言
  • Java中的this关键字的用法
  • 网络安全威胁狩猎:主动防御的终极指南
  • C#在二合一平板电脑关于旋转模式相关设置
  • 2026 中考游记
  • MinIO 介绍(3)--MinIO 客户端 mc 管理员功能
  • 8.16
  • 2025-10-19
  • 一文读懂隔离见证
  • 12131
  • 关于火柴盒的记忆
  • PWN手的成长之路-19-int_overflow
  • FFmpeg开发笔记(八十四)使用国产的librestreaming实现RTMP直播
  • 2025 年闪测仪厂家企业品牌推荐排行榜,一键式闪测仪,卧式闪测仪,影像闪测仪,立式闪测仪,2D3D 混合式闪测仪,高精度闪测仪,大量程闪测仪,复合式闪测仪公司推荐
  • 2025 年耐火砖厂家企业品牌推荐排行榜,绝热,轻质,莫来石,保温,莫来石轻质,氧化铝泡沫,氧化铝空心球,抗渗碳,高温轻质莫来石,高温耐火砖公司推荐
  • 2025 年护栏板厂家企业品牌推荐排行榜,波形,高速,镀锌,二波,三波,喷塑,国标,绳索,公路护栏板,护栏板立柱公司推荐
  • 2025 年船用锅炉厂家企业品牌推荐排行榜,基于市场口碑,评选值得信赖的船用锅炉公司推荐
  • 2025 年反应釜厂家企业品牌推荐排行榜,实验室,高压,加氢,不锈钢,试验室,氢化,聚合,高温,钛材反应釜公司推荐