今天网上冲浪的时候突然看到一道小学六年级的数学题，如上图所示，求阴影部分的面积。
我下意识就想到了微积分，这不就建立坐标系，求出交点，计算积分就行了嘛。转念一想，小学生哪里会积分，这道题一定有简单的解法。
经过一番努力，我尝试出了一共两种精确解法和一种近似解法，但都超出了小学生的理解范畴。
首先我们把图左右翻转一下，方便建立坐标系，如下图所示：
积分方法
首先联立半圆和对角线的方程：
解出交点E的坐标是
所以可以直接求下面积分，得到要求的阴影面积：
利用简单的高等数学知识，可以解出答案是下面这样的，约等于1.252。
但是小学生肯定不会积分呀，就算高中生大多数都不会！
间接计算
既然我们不方便算这个不规则图形的面积，我们可以计算其他规则图形的面积，然后加减之后间接得到答案。
图中阴影部分面积其实等价于三角形ACE+梯形BOCE-扇形BOE，而这三部分其实都非常好算。
因为E的坐标是，所以最难算的就是扇形BOE的面积，也就是要求角OBE的大小。很轻松就可以得到角度是，所以三部分运算后可以得到和上面积分相同的答案。
这个方法初中生都会做了！但是对于小学生来说难度还是有点大了。
蒙特卡洛采样方法
如果往第一象限正方形内随机抛撒若干豆子，然后看掉落在阴影部分里的豆子所占的比例，那么就可以近似估计出阴影部分的面积了，这就叫蒙特卡洛采样。
那我们可以随机生成两个随机数，作为采样点的坐标，然后判断是否在阴影部分里，重复上述过程，记录次数就行了。
我这里写了C++和CUDA两个实现，CUDA可以并行采样，速度快很多。
执行的话，C++保存为test.cpp，然后执行g++ test.cpp -o test和./test。CUDA保存为test.cu，然后执行nvcc test.cu -o test和./test。
实测CUDA相同时间内可以采样的次数大约是100倍左右，主要受到原子写操作的局限，不然可以达到1000倍左右的加速。
但是小学生会CUDA吗？显然不可能。
C++实现
#include <iostream> typedef unsigned long long int ull; int main() { srand(time(NULL)); ull res = 0; ull cnt = 1e8; for (ull i = 0; i < cnt; ++i) { double x = double(rand()) / RAND_MAX; double y = double(rand()) / RAND_MAX; if (x * x + (y - 1) * (y - 1) > 1 && x + 2 * y - 1 < 0) { res += 1; } } printf("%.5f\n", res * 16.0 / cnt); return 0; }
CUDA实现
#include <iostream> #include <chrono> #include <curand_kernel.h> typedef unsigned long long int ull; __global__ void monte_carlo_kernel(ull* res, ull cnt, const int seed) { ull i = blockIdx.x * blockDim.x + threadIdx.x; if (i >= cnt) return; curandStatePhilox4_32_10_t state; curand_init(seed, i, 0, &state); double2 rand = curand_uniform2_double(&state); double x = rand.x; double y = rand.y; if (x * x + (y - 1) * (y - 1) > 1 && x + 2 * y - 1 < 0) { atomicAdd(res, 1); } } void launch_monte_carlo(ull* res, ull cnt, cudaStream_t& stream) { ull grid_dim = cnt >> 10; monte_carlo_kernel<<<grid_dim + 1, 1024, 0, stream>>>( res, cnt, std::chrono::duration_cast<std::chrono::microseconds>( std::chrono::system_clock::now().time_since_epoch()) .count()); } int main() { ull* res; cudaMallocManaged(&res, sizeof(ull)); res[0] = 0; ull cnt = 1e10; cudaStream_t stream; cudaStreamCreate(&stream); launch_monte_carlo(res, cnt, stream); cudaDeviceSynchronize(); printf("%.5f\n", res[0] * 16.0 / cnt); cudaFree(res); return 0; }
标准方法
其实标准的小学生解法不需要建坐标系，只需要如下图所示，做4条辅助线就行。
很容易证明，辅助线将三角形ACD划分出来的5个小三角形是全等的。
那么阴影面积就等于正方形AOBD减去四分之一圆BOD，再减去最右侧的弧边三角形ADE。而弧边三角形ADE面积等于三角形ACD减去三角形BCE，再减去扇形BDE。
所以问题的关键就是求角DBE的大小，因为小三角形的边长可以轻易求出，所以最终答案很好算：
可以发现和之前高等数学得到的结果形式并不同，其实数值是一样的，所以我们得到了一个等式：
这个可能看起来不是那么直观，但确实是成立的，根据反三角函数恒等式，可以推出等号左边等于。构造一个边长3、4、5的直角三角形，显然答案就是。
总结与拓展
其实这道题小学生也没几个人能做出来，就算给我们做也得费老鼻子劲，还不一定算的对。
最近还看到了很多类似的趣图，这里分享给大家，闲暇时间可以做了玩玩。
可以在评论区说出你的答案哦。
