双色球中奖率分析（python)

我的双色球分析程序

 

因为平时有买双色球，无奈屡买不中。故使用python完成了一个双色球分析程序,来提高中奖概率。 下面就此程序，进行简单说明。     程序介绍： 1、根据以往中奖纪录，进行分析，分析一段时间各个号码出现的次数，时间段分为三个，分别为一周（3期）、一月（13期）、一年（153期）。 2、根据第一条分析结果，计算每个号出现的概率。 3、给每个时间段（一周、月、年）统计出的号码出现概率，附一个权重，即用来表示，对最终分析结果的影响力。具体如下      （1）统计一年中（153期）出现最多的号码。        (2)  统计一月中（13期）出现最多的号码。        (3)  统计一周中（3期）出现最少的号码。        原因：在双色球的摇奖过程中，因为球的质量，形状问题，可能会导致某些号码总会比别的号码出现几率大。       故长期呢，分析出现次数多的，短期内，分析出现次数少的，作为将来出现概率大的号码。      备注：双色球，2年换一次球，一共4套球。 4、根据第三条计算方法，计算出每个号码的总权重（即综合各个时间段的权重），每个时间段的权重比率，采用黄金分割比率，为（4.236,1.618,2.618）作为(年、月、日）的比率。 5、根据以上统计结果，计算出将来出现几率最大的号码。       主要功能： 1、计算将来中奖概率最大的号码。 2、随机产生一注。 3、自选号码，计算自选号码的中奖概率。      代码如下： #!d:/python32/python #coding = <utf-8> #双色球计算 #Written by LiJunpin import random import re class shuangseqiu(): def __init__(self): self.quanzhong = [2.618,1.618,4.236]#权重，分别为一周统计数据，一月统计数据，一年统计数据与完全随机权重的比值，表示对最后选出的号码的影响力大小，此数据依据最伟大的黄金分割点制定。 #self.quanzhong = [1,1,1] self.qishu = [4,13,153]#一周、一月、一年的彩票期数 self.data_long = ['W','M','Y']#期数时间范围的长度表示 self.fle = '123.txt'#打开彩票文件 self.red_dict = {}#红球存储字典 self.blue_dict = {}#篮球字典 self.red_rate = {}#红球出现率字典 self.blue_rate = {}#篮球出现率字典 self.red_qz = []#红球权重 self.blue_qz = []#篮球权重 #初始化存储字典 for i in range(1,34): self.__initdict__(i,self.red_dict) self.__initdict__(i,self.red_rate) for i in range(1,17): self.__initdict__(i,self.blue_dict) self.__initdict__(i,self.blue_rate) #print(self.red_dict,self.blue_dict) def __initdict__(self,num ,color_dict): color_dict[str(num)] = {} for x in range(len(self.data_long)): color_dict[str(num)][self.data_long[x]] = 0 def __count__(self,num, color_dict, data_long):#统计一个号码出现的次数 #print(num) color_dict[str(num)][data_long] += 1 #print (num,data_long) def __countlst__(self,ball_lst,data_long): #print(ball_lst) #年数据存储 for i in range(6): self.__count__(ball_lst[i],self.red_dict, data_long) self.__count__(ball_lst[6],self.blue_dict,data_long) #print(self.red_dict) #print(self.blue_dict) def __readdata__(self): #从文件读取彩票中奖纪录 ball_file = open(self.fle, 'r') ball_lst = ball_file.readline() for i in range(1,self.qishu[2]): ball_lst = ball_file.readline().split() #print(ball_lst) #red = ball_list[5:11] #blue = ball_list[11:] #print(ball_lst) ball_lst = ball_lst[4:] if i <= self.qishu[0]: self.__countlst__(ball_lst,self.data_long[0]) if i <= self.qishu[1]: self.__countlst__(ball_lst,self.data_long[1]) self.__countlst__(ball_lst,self.data_long[2]) ball_file.close() #print(self.red_dict) #print('------------------------') #print(self.blue_dict) def __rateone__(self,qishu,data_long):#根据统计的一段时间内(以data_long为依据)的出现次数计算1-33,1-16的出现几率 #计算总出现的次数 redall = qishu * 6 blueall = qishu * 1 #计算红球出现率 for i in range(1,34): self.red_rate[str(i)][data_long] = self.red_dict[str(i)][data_long]/redall #计算篮球出现率 for i in range(1,17): self.blue_rate[str(i)][data_long] = self.blue_dict[str(i)][data_long]/blueall def __rate__(self): for i in range(len(self.data_long)): self.__rateone__(self.qishu[i],self.data_long[i]) #print(self.red_rate) def __quanzhong__(self,num,color_rate):#计算num号码对应的权重 value = (1/len(color_rate)- color_rate[str(num)][self.data_long[0]])* self.quanzhong[0]/sum(self.quanzhong)#一周出现概率高，那么后面，出现概率就会降 value += (color_rate[str(num)][self.data_long[1]] - 1/len(color_rate)) * self.quanzhong[1]/sum(self.quanzhong) value += (color_rate[str(num)][self.data_long[2]] - 1/len(color_rate)) * self.quanzhong[2]/sum(self.quanzhong) return value def make_quanzhong(self):#生成整个权重的列表 for i in range(1,34): self.red_qz.append([self.__quanzhong__(i,self.red_rate),i]) for i in range(1,17): self.blue_qz.append([self.__quanzhong__(i,self.blue_rate),i]) #print(self.red_qz) #print('--------------------------------------') #print(self.blue_qz) def init_data(self):#初始化读入数据、生成概率和权重数据结构 self.__readdata__() self.__rate__() self.make_quanzhong() def _str_format(self,s):#格式化输入的字符串。形成数字列表 a = re.compile('\d+') b = re.findall(a,s) return b def _find_qz(self,color_qz,num):#查找数字num的出现概率 for i in range(0,34): if color_qz[i][1] == num: return color_qz[i][0] def _suiji_gl(self):#计算随机概率 sjgl = 1 for x in range(28,34): sjgl *= x sjgl = 1/sjgl return sjgl * 1/16 def gailv(self,s):#计算概率 lst = self._str_format(s) #print(lst) lst = [int(x) for x in lst] #print(lst) gl = 1 if not lst or len(lst)>7 or len(lst)<7 or max(lst)>33 or min(lst)<1: return None else: for x in lst[:7]: gl *= (1+self._find_qz(self.red_qz,x)) return gl * self._suiji_gl() def caipiao_random(self):#随机红篮球列表 redball = sorted(random.sample(range(1,34),6)) blueball = random.sample(range(1,17),1) return redball + blueball def print_random(self):#打印随机结果 random_ball = self.caipiao_random() print('--------------------------------------------------') print('红球：',random_ball[:6],'篮球：',random_ball[6:7]) print('--------------------------------------------------') def print_best_number(self):#输出中奖概率最高的号码 print('根据规则选出的最可能中奖的号码如下:') print('--------------------------------------------------') print('红球：',sorted(list(x[1] for x in sorted(self.red_qz)[-6:])),'篮球：',sorted(self.blue_qz)[-1][1]) print('--------------------------------------------------') def print_gl(self): print('输入格式为-> \"(1,9,13,21,28,32)(14)\"') c = input('请输入:') if not self.gailv(c): print('输入有误!') return print('--------------------------------------------------') print(' 随机概率为：%0.10f%%'%(self._suiji_gl()*100)) print(' -------------------------') print(' 自选概率为：%0.10f%%'%(self.gailv(c)*100)) print('--------------------------------------------------') def print_all(shuangseqiu_obj): print(''' |----------功能列表-------------| |--1:打印中奖率最高代码---------| |--2:打印随机号码---------------| |--3:自助选取号码，打印中奖概率-| |--q:退出-----------------------| ''') while True: a = input('功能选择:') if a == '1': shuangseqiu_obj.print_best_number() elif a == '2': shuangseqiu_obj.print_random() elif a == '3': shuangseqiu_obj.print_gl(); elif a == 'q': break; else: print('输入有误') print('谢谢使用') if __name__ == '__main__': b = shuangseqiu() b.init_data() print_all(b) 程序中使用双色球历史中奖号码进行分析，文件为123.txt，文件格式如下： 序列    年份    期号    日期        红1    红2    红3    红4    红5    红6    蓝球 1310    2012    2012025    2012/3/1    4    12    19    21    25    28    13     1308    2012    2012023    2012/2/28    5    9    15    23    24    33    2     1307    2012    2012022    2012/2/26    4    8    12    24    26    27    4     1306    2012    2012021    2012/2/23    1    15    16    18    22    30    3     1305    2012    2012020    2012/2/21    3    8    12    18    23    29    11     1304    2012    2012019    2012/2/19    7    14    18    20    22    30    16     1303    2012    2012018    2012/2/16    3    5    6    22    26    32    15     1302    2012    2012017    2012/2/14    6    9    14    19    25    28    10     1301    2012    2012016    2012/2/12    2    5    12    17    22    25    8     1300    2012    2012015    2012/2/9    1    3    6    10    21    23    15     1299    2012    2012014    2012/2/7    1    2    5    16    28    30    12     1298    2012    2012013    2012/2/5    6    8    24    29    30    32    13     1297    2012    2012012    2012/2/2    15    17    18    20    23    27    1     1296    2012    2012011    2012/1/31    4    14    15    16    20    26    5     1295    2012    2012010    2012/1/29    1    3    13    19    25    26    10     1294    2012    2012009    2012/1/19    4    16    24    26    27    33    11     1293    2012    2012008    2012/1/17    1    12    20    23    24    29    8     1292    2012    2012007    2012/1/15    10    17    19    27    28    32    4     1291    2012    2012006    2012/1/12    2    22    25    29    32    33    8     1290    2012    2012005    2012/1/10    7    9    18    27    31    33    6     1289    2012    2012004    2012/1/8    1    5    10    11    21    23    16     1288    2012    2012003    2012/1/5    3    6    8    24    29    31    9     1287    2012    2012002    2012/1/3    2    3    7    9    10    32    13     1286    2012    2012001    2012/1/1    1    4    5    9    15    17    6     1285    2011    2011153    2011/12/29    5    8    9    10    20    25    13     1284    2011    2011152    2011/12/27    4    10    11    12    21    26    13     1283    2011    2011151    2011/12/25    7    11    16    19    31    33    10     1282    2011    2011150    2011/12/22    8    10    12    15    22    27    13     1281    2011    2011149    2011/12/20    4    5    6    7    23    31    16     1280    2011    2011148    2011/12/18    5    14    22    23    25    26    14     1279    2011    2011147    2011/12/15    4    8    12    17    18    30    10     1278    2011    2011146    2011/12/13    11    23    26    28    32    33    10     1277    2011    2011145    2011/12/11    2    4    14    15    26    30    4     1276    2011    2011144    2011/12/8    1    2    9    10    16    24    3