acdream 1234-Two Cylinders（辛普森积分）

最新推荐文章于 2018-05-27 18:10:06 发布

原创最新推荐文章于 2018-05-27 18:10:06 发布 · 401 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#数论

数论专栏收录该内容

131 篇文章

订阅专栏

本篇博客介绍了一道数学题目，该题目要求计算两个不同半径的无限长圆柱体在轴线垂直相交时的相交体积。对于相同半径的情况给出了解析解，而对于不同半径的情况，则采用辛普森数值积分法进行求解。

Two Cylinders

Special Judge Time Limit: 2000/1000MS (Java/Others) Memory Limit: 128000/64000KB (Java/Others)

Submit Statistic Next Problem

Problem Description

In this problem your task is very simple.

Consider two infinite cylinders in three-dimensional space, of radii R1 and R2 respectively, located in such a way that their axes intersect and are perpendicular.

Your task is to find the volume of their intersection.

Input

Input file contains two real numbers R 1 and R 2 (1 <= R 1,R 2 <= 100).

Output

Output the volume of the intersection of the cylinders. Your answer must be accurate up to 10 -4.

Sample Input

1 1

Sample Output

5.3333

Source

题意：给你两个半径不同的圆柱体，计算垂直相交的体积

题解：当半径相同时，大一高数中学过，相交的体积为：16/3*R^3

当半径不相同时：

这里我们采用辛普森积分来进行积分（否则会越积越复杂）

#include <math.h>
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
using namespace std;
const double eps = 1e-10;
const double PI = acos(-1.0);
double r1,r2;
double f(double x)
{
	return sqrt((r1*r1-x*x)*(r2*r2-x*x)); 
}
double simpson(double a,double b)
{
	double c=a+(b-a)/2;
	return (f(a)+4*f(c)+f(b))*(b-a)/6; 
}
double asr(double a,double b) ///迭代求解数值积分  
{  
    double mid=(a+b)/2.0;  
    double res=simpson(a,b);  
    if(fabs(res-simpson(a,mid)-simpson(mid,b))<eps)  
    {  
        return res;  
    }  
    else return asr(a,mid)+asr(mid,b);  
}
int main(void)
{
	while(scanf("%lf%lf",&r1,&r2)!=EOF)
	{
		if(r1==r2)
			printf("%.4lf\n",16.0/3*r1*r1*r1);
		else
		{
			if(r1>r2)
				swap(r1,r2);
			printf("%.4lf\n",8.0*asr(0,r1));
		}
	}
	return 0;
}