自适应辛普森法

1.并查集模板2023-10-31 2.求质因数模板2023-10-31 3.二分图最大匹配模板(匈牙利算法)2023-12-02 4.欧拉函数模板2023-11-30 5.ST表模板2023-11-26 6.快速幂模板2023-10-31 7.字典树模板2023-12-04 8.矩阵模板2023-12-05 9.Dijkstra单源最短路模板2023-12-09 10.最近公共祖先模板(LCA)2023-12-12 11.拓扑排序模板2023-12-12 12.区间素数筛模板2023-12-16 13.Kruskal和Prim模板2023-12-24 14.树状数组模板2023-12-26 15.二维坐标离散化模板2024-03-07 16.单点修改区间查最值-树状数组模板2024-03-11 17.KMP模板2024-03-11 18.二叉搜索树模板2024-03-30 19.DIjkstra进阶模板路径记录按权重(结点数最小等)记录2024-03-31 20.判断负环模板2024-04-28 21.Exgcd 模板2024-05-03 22.压位高精度模板2024-05-06 23.线段树模板2024-05-12 24.扫描线模板2024-05-18 25.莫队模板2024-07-17 26.带修莫队模板2024-07-17 27.SCC缩点模板2024-07-18 28.取模+组合数2024-07-24 29.FFT 高精度乘法模板2024-07-24 30.字符串自然溢出哈希/单哈希/双哈希模板2024-07-25 31.树模板2024-07-27 32.dsu on tree 模板2024-07-28 33.线段树模板重制2024-08-03 34.主席树模板2024-08-07 35.大数质因数分解模板2024-08-08 36.线段树合并模板2024-08-09 37.int128输入输出流2024-08-09 38.Meissel_Lehmer模板2024-08-09 39.浮点高精度2024-08-13

40.自适应辛普森法2024-10-02

41.unordered_map随机底数种子2024-10-07

自适应辛普森法

 double a, b, c, d, l, r;
 
const double eps = 1e-8;
double F(double x) {
    //需要积分的公式
	return (c * x + d) / (a * x + b);
}
double simpson(double l, double r) {
	double mid = (l + r) / 2;
	return (F(l) + 4 * F(mid) + F(r)) * (r - l) / 6;
}
double asr(double l, double r, double A) {
	double mid = (l + r) / 2;
	double L = simpson(l, mid), R = simpson(mid, r);
	if (fabs(L + R - A) <= 15 * eps)
		return L + R + (L + R - A) / 15.0;
	return asr(l, mid, L) + asr(mid, r, R);
}
double asr(double l, double r) {
	return asr(l, r, simpson(l, r));
}

上一篇SMU Autumn 2024 Team Round 3（The 2024 ICPC Latin America Championship）

下一篇unordered_map随机底数种子

本文作者：Ke_scholar

本文链接：https://www.cnblogs.com/Kescholar/p/18444561

posted @ 2024-10-02 12:32 Ke_scholar 阅读(13) 评论(0) 编辑收藏举报

刷新页面返回顶部

登录后才能查看或发表评论，立即登录或者逛逛博客园首页

Ke_scholar

自适应辛普森法

自适应辛普森法

公告

最新随笔

我的标签

合集 (3)

随笔分类 (242)

随笔档案 (228)

阅读排行榜

评论排行榜

推荐排行榜

最新评论

	double a, b, c, d, l, r;

	const double eps = 1e-8;
	double F(double x) {
	//需要积分的公式
	return (c * x + d) / (a * x + b);
	}
	double simpson(double l, double r) {
	double mid = (l + r) / 2;
	return (F(l) + 4 * F(mid) + F(r)) * (r - l) / 6;
	}
	double asr(double l, double r, double A) {
	double mid = (l + r) / 2;
	double L = simpson(l, mid), R = simpson(mid, r);
	if (fabs(L + R - A) <= 15 * eps)
	return L + R + (L + R - A) / 15.0;
	return asr(l, mid, L) + asr(mid, r, R);
	}
	double asr(double l, double r) {
	return asr(l, r, simpson(l, r));
	}