#!/usr/bin/env python
#-*- coding:utf-8 -*-

from math import ceil
from random import randint
from decimal import Decimal

def generate_papers():
    papers = [[1,1]for i in range(25)]
    for i in range(5):
        throw_unknown(papers)
    return papers

def throw_unknown(papers):
    unknown = randint(1,50)
    if ceil(unknown / 2.0) == unknown / 2.0:
        if papers[int(unknown / 2.0) - 1][1] == 0:
            throw_unknown(papers)
        else:
            papers[int(unknown / 2.0) - 1][1] = 0
    else:
        if papers[int(ceil(unknown / 2.0)) - 1][0] == 0:
            throw_unknown(papers)
        else:
            papers[int(ceil(unknown / 2.0)) - 1][0] = 0
    return papers

def choose_paper(papers):
    chosen = papers[randint(0,24)]
    if chosen == [1,1]:
        return 1
    else:
        return 0

def recognize_order(papers):
    unlearned_papers = []
    for paper in papers:
        if paper[0] == 0 or paper[1] == 0:
            unlearned_papers.append(paper[0] + paper[1])
    unlearned_papers.sort()
    if unlearned_papers == [1,1,1,1,1]:
        return 1
    elif unlearned_papers == [0,1,1,1]:
        return 2
    elif unlearned_papers == [0,0,1]:
        return 3
    else: print "ERROR: Impossible order! %s" % unlearned_papers


tests_quantity = int(raw_input("How many tests should we make? "))
passed = not_passed = first = second = third = 0
for i in range(tests_quantity):
    test = generate_papers()
    attempt = choose_paper(test)
    order = recognize_order(test)

    if attempt == 1: passed += 1
    elif attempt == 0: not_passed += 1

    if order == 1: first += 1
    elif order == 2: second += 1
    elif order == 3: third += 1

print("""

From %d tests %d succesed and %d failed.
Success probability is %.2f.

Unlearned papers recognition:
 1. 11111 - %d times (p = %.2f)
 2.  2111 - %d times (p = %.2f)
 3.   221 - %d times (p = %.2f)

""" % (tests_quantity, passed, not_passed,
       Decimal(passed) / Decimal(tests_quantity),
       first, Decimal(first) / Decimal(tests_quantity),
       second, Decimal(second) / Decimal(tests_quantity),
       third, Decimal(third) / Decimal(tests_quantity)))