یک برنامهٔ C++ بنویسید که با استفاده از روش مونت‌کارلو و تابع std::srand مقدار π را با تولید N نقطهٔ تصادفی تقریب بزند.

برای تولید نقاط تصادفی از std::srand(time(0)) استفاده کنید و با std::rand() / RAND_MAX اعداد یکنواخت در [0,1) بسازید. برای هر نقطه x و y در [0,1) بررسی کنید که x*x + y*y <= 1 باشد تا داخل دایره واحد باشد و نسبت نقاط داخل را ضرب در 4 کنید تا مقدار π به دست آید. برای شمارنده‌ها از نوع دادهٔ مناسب مانند unsigned long long استفاده کنید و هرچه N بزرگ‌تر باشد دقت بهتر می‌شود. اگر به HPC فکر می‌کنید، می‌توانید با تقسیم کار یا کامپایل با -O2 کارایی را افزایش دهید و به نکتهٔ thread-safetyِ rand() نیز توجه کنید.