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

برای حل، ابتدا با std::srand(time(nullptr)) تولیدکنندهٔ اعداد تصادفی را مقداردهی اولیه کنید (نیاز به دارید)، سپس با std::rand() و تقسیم بر RAND_MAX دو عدد تصادفی در بازهٔ [0,1) بسازید و تعداد نقاطی که داخل ربع‌دایرهٔ واحد قرار می‌گیرند بشمارید. نسبت نقاط داخل به کل نمونه‌ها را در 4 ضرب کنید تا تقریب π به‌دست آید. نکات: برای دقت بهتر N را بزرگ انتخاب کنید و از نوع‌های صحیح بزرگ (مانند long long) برای شمارش استفاده کنید؛ برای کاربردهای HPC می‌توانید حلقه‌ها را بعداً با OpenMP یا توزیع کار موازی کنید تا سرعت و دقت با نمونه‌های بیشتر افزایش یابد.