using System;
using static System.Console;
using static System.Math;
public class main{

static (double,double) plainmc(
Func<vector,double> f,vector a,vector b,int N, Random RND = null
	){
        int dim=a.size; double V=1; for(int i=0;i<dim;i++)V*=b[i]-a[i];
        double sum=0,sum2=0;
	var x=new vector(dim);
	if(RND == null) RND = new Random();
        for(int i=0;i<N;i++){
                for(int k=0;k<dim;k++)x[k]=a[k]+RND.NextDouble()*(b[k]-a[k]);
                double fx=f(x); sum+=fx; sum2+=fx*fx;
                }
        double mean=sum/N, sigma=Sqrt(sum2/N-mean*mean);
        var result=(mean*V,sigma*V/Sqrt(N));
        return result;
}

public static void Main(string[] args){
	int N=1000;
	if(args.Length>0)N=(int)double.Parse(args[0]);
	double R=1;
	Func<vector,double> f = x =>{
		if(x.norm()<=R)return 1;
		else return 0;
		};
	vector a = new vector(0.0,0.0);
	vector b = new vector(1.0,1.0);
	(double q, double e) = plainmc(f,a,b,N);
	double exact = PI/4;
	WriteLine($"{N} {q} {e} {Abs(q-exact)}");
	}
}//class main