from random import *

def arrivée():
        L=[[0.1,0.9],[0.9,0.1],[0.45,0.55],[0.9,0.8]]
        for k in range(len(L)):
            X=0
            for i in range(1000):
                p=random()
                q=random()
                if q<L[k][1]:
                    X+=1
                if p<L[k][0]:
                    if X>0:
                        X-=1
            print("Pour les paramètres " + str(L[k][0])+" "+ str(L[k][1]) )
            print(X)


def S2():
        for l in range(0,1000,50):
            L=[[0.1,0.9],[0.9,0.1],[0.45,0.55],[0.9,0.8]]

            for k in range(len(L)):
                cpt=0
                for j in range(2000):
                    X=l
                    for i in range(1000):
                        p=random()
                        q=random()
                        if q<L[k][1]:
                            X+=1
                        if p<L[k][0]:
                            if X>0:
                                X-=1
                        if X==0:
                            cpt+=1
                            break
                print("Pour un départ à " + str(l) )
                print(cpt/2000)
            print()
#on remarque qu'un guichet rapide correspond aux paramètres 2 car ceux ci sont capables d'absorber rapidement une forte queue: on passe de P=1 à P=0 aux alentours de 800
#Concernant les paramètres 4, on a un basculement vers x=100
#Les autres parametres ne permettent pas d'absorber le flux

def S3(x):
        L=[[0.1,0.9],[0.9,0.1],[0.45,0.55],[0.9,0.8]]
        for k in range(len(L)):
            res=[0 for k in range(1000) ]
            for j in range(2000):
                X=x
                for i in range(10000):
                    p=random()
                    q=random()
                    if q<L[k][1]:
                        X+=1
                    if p<L[k][0]:
                        if X>0:
                            X-=1
                    if X==x:
                        res[i]+=1
                        break
            S=0
            for m in range(len(res)):
                S+=m*res[m]/2000
            print(S)

#explication pour le premier paramètre: on ne revient jamais au départ d'ou un faible score
#pour le deuxieme, on revient de suite au depart d'ou aussi un faible score
#on revient pour le troisième mais plus longtemps après
#on met aussi plus de temps =>si p et q proches on a une grande espérance

#Le choix de x importe peu