from random import *
N= int(input("Taille de chaque échantillon : N = "))
P= int(input("pourcentage d’individus ayant le caractère considéré dans la population étudiée (%) : P = "))
M= int(input("Nombre d'échantillons : "))
print(" ")
print("simulation informatique des échantillons :")
print(" ")
X=[0 for i in range(M)]
MX=0
MXX=0
for j in range(M):
  X[j]=0
  E=[0 for i in range(N)]
  for i in range(N):
    roue=randrange(100)+1
    if roue<=P:
      X[j]=X[j]+1
      E[i]=1
  MX=MX+X[j]
  MXX=MXX+X[j]*X[j]
  print("échantillon numéro ",j+1," ",E)
  print(" nombre d'individus ayant le caractère considéré : X =",X[j])
MX=MX/M
MXX=MXX/M
print(" ")
print("analyse des résultats : ")
print(" ")
print("moyenne observée des nombres d'individus X : ",MX," N P =",N*P/100)
SXX=MXX-MX*MX
print("variance observée des nombres d'individus X : ",SXX," N P (1-P) = ",N*P*(1-P/100)/100)
print(" ")
print("fréquences des différentes valeurs de X : ")
print(" ")
FA=[0 for i in range(N+1)]
print("valeur de X, fréquence d'apparition")
for i in range(N+1):
  for j in range(M):
    if X[j]==i:
      FA[i]=FA[i]+1
  FA[i]=FA[i]/M
  print(" ",i," ",FA[i])
print(" ")
T=0
for i in range(N+1):
  T=T+FA[i]
print("total : ",T)