#include <iostream>
#include <ctime>
#include <cstdlib>

int main(){
	int NumWin = 10000;
	int t = 3;
	srand(time(NULL));
	double eta = 0;
	for(double p = 0.00; p < 1; p += 0.05){
		//theor
		eta = (1 - p)/(1 + t*p);
		//practice
		int count = 0;
		int m = 0;
		for (int i = t+1; i <= NumWin; i += t+1, m ++){
			if ((double)rand() / RAND_MAX > p){
				i -= t;
				count ++;
				if(m > t)
					m=0;
			}
		}
		std::cout << p << '\t' << eta << '\t' <<(float)count/NumWin<< std::endl;
	}
	return 0;
}