python写一个双色球彩票计算器

首先声明，赌博一定不是什么好事，也完全没有意义，不要指望用彩票发财。之所以写这个，其实是用来练手的，可以参考这个来预测一些其他的东西，意在抛砖引玉。

啰嗦完了，马上开始，先上伪代码

打开网址

读取内容
内容解析

根据源码得到需爬取内容

1、开奖日期：2018年8月26日

2、红球

<li class="ball_red">03</li>

<li class="ball_red">07</li>

<li class="ball_red">08</li>

<li class="ball_red">14</li>

<li class="ball_red">25</li>

<li class="ball_red">32</li>

3、篮球

<li class="ball_blue">06</li>

打开数据库连接

爬取内容写入数据库

共3个字段

1、开奖日期

2、红球，红球使用分号“;”分隔，方便调用和导出

3、篮球

在写脚本前，建议先写伪代码，伪代码格式不是固定的，随自己喜好，主要就是在思维及算法落地前，把整个轮廓理清，可以有效降低都快写完了，发现前面有错误，结果导致整个脚本全部更新一遍这种事

的发生概率

伪代码解读：

共分为两个功能块

一就是抓取彩票数据，这个使用爬虫实现，分别抓取开奖日期、红球区、篮球区，因为考虑双色球的数据量比较庞大，所以这次使用数据库进行存储，选用的是免费又好用的mysql数据库，数据库接口文件使用MySQLdb，这个我以后会单独写一个说明，当然你也可以用文档存储，或者选择别的数据库比如oracle或者nosql的mangodb

二是用来分析彩票的，即使用彩票数据进行下一期的彩票预测，本次选用的是二项分布，这个在之前的算法里面有，就不重复说明了

下面先上脚本，再对脚本进行说明

#!/usr/bin/python
# coding: UTF-8

'''
打开网址
  读取内容
内容解析
  根据源码得到需爬取内容
  1、开奖日期：2018年8月26日
  2、红球
                                <li class="ball_red">03</li>
                                <li class="ball_red">07</li>
                                <li class="ball_red">08</li>
                                <li class="ball_red">14</li>
                                <li class="ball_red">25</li>
                                <li class="ball_red">32</li>
  3、篮球
                                <li class="ball_blue">06</li>
打开数据库连接
  爬取内容写入数据库
    共3个字段
    1、开奖日期
    2、红球，红球使用分号“;”分隔，方便调用和导出
    3、篮球
create table tow_color_ball(open_date varchar(10),red_n varchar(20),blue_n varchar(2))
'''

import urllib
import urllib2
import re
import numpy as np
import operator
import MySQLdb

# 连接mysql
def conn_db():
  db = 'pythondb'
  host = 'localhost'
  iuser = 'xxx'
  passwd = 'xxxxxx'
  conn = MySQLdb.connect(db = db, host = host, user = iuser, passwd = passwd)
  cursor = conn.cursor()
  return cursor

# 处理网页获取页面源码
def get_html_values(url):
  url_open = urllib.urlopen(url)
  url_read = url_open.read()
  return url_read

# 处理源码，获取日期、红球、篮球
def manage_html(html_values):
  red_no_re = re.compile('(?<=\<li class\=\"ball_red\"\>)[0-9]+(?=\<\/li\>)')
  blue_no_re = re.compile('(?<=\<li class\=\"ball_blue\"\>)[0-9]+(?=\<\/li\>)')
  date_re = re.compile('(?<=开奖日期：)[0-9]+年[0-9]+月[0-9]+日')
  red_no_list = re.findall(red_no_re,html_values)
  red_numbers = ';'.join(red_no_list)
  blue_number = re.search(blue_no_re,html_values)
  blue_number = blue_number.group()
  date_value = re.search(date_re,html_values)
  date_value = date_value.group()
  return date_value, red_numbers, blue_number

# 可恶的日期，竟然是YYYY年MM月DD日，需要改成YYYY-MM-DD
def manage_date(date_value):
  date_value = date_value.replace('年','-').replace('月','-').replace('日','')
  return date_value

# 处理页面编号，每次编号-1，也就是说end_page要小于url中的页码
def get_page(url,end_page):
  url_num = re.search('(?<=\/)[0-9]+(?=\.)',url)
  url_num = url_num.group()
  if int(end_page) > int(url_num):
    return 'end'
  url_num_1 = int(url_num) - 1
  url = url.replace(url_num,str(url_num_1))
  return url

# 查看库中是否已存在开奖日期，防止重复写入
def check_open_date(open_date):
  conn = conn_db()
  check_sql = 'select 1 from tow_color_ball where open_date = %r' %open_date
  conn.execute(check_sql)
  excur = conn.fetchall()
  conn.close()
  #如过未查到excur == ()
  return excur

# 写入数据库
def write_db(date_value, red_numbers, blue_number):
  conn = conn_db()
  in_sql = "insert into tow_color_ball(open_date,red_n,blue_n) values(%r,%r,%r)" %(date_value, red_numbers, blue_number)
  conn.execute(in_sql)
  conn.execute('commit')
  conn.close()

# 彩票主程序，用来爬取彩票号码
def ball_main(url,end_page):
  while True:
    html_values = get_html_values(url)
    date_value, red_numbers, blue_number = manage_html(html_values)
    date_value = manage_date(date_value)
    data_check = check_open_date(date_value)
    if data_check == ():
      write_db(date_value, red_numbers, blue_number)
    url = get_page(url,end_page)
    if url == 'end':
      print 'url_page已到达end_page,获取完成'
      return 0

# 二项分布算法
class binomial_class(object):
  def __init__(self,case_count,real_count,p):
    self.case_count = case_count
    self.real_count = real_count
    self.p = p

  def multiply_fun(self,xlist):
    n = 1
    for x in xlist:
      n *= x
    return n

  def fact_fun(self,n):
    if n == 0:
      return 1
    n += 1
    fact_list = [i for i in range(1,n)]
    fact_num = self.multiply_fun(fact_list)
    return fact_num

  def c_n_x(self):
    fact_n = self.fact_fun(self.case_count)
    fact_x = self.fact_fun(self.real_count)
    fact_n_x = self.fact_fun(self.case_count - self.real_count)
    c_n_x_num = float(fact_n) / (fact_x * fact_n_x)
    return c_n_x_num

  def binomial_fun(self):
    c_n_k_num = self.c_n_x()
    pi = (self.p ** self.real_count) * ((1 - self.p) ** (self.case_count - self.real_count))
    binomial_num = c_n_k_num * pi
    return binomial_num

# 从库里获取彩票信息
def get_ball_infomation(start_dt,end_dt):
  conn = conn_db()
  sql = "select red_n,blue_n from tow_color_ball where date_format(open_date,'%%Y-%%m-%%d') >= %r and date_format(open_date,'%%Y-%%m-%%d') <= %r" %(start_dt,end_dt)
  conn.execute(sql)
  excur = conn.fetchall()
  conn.close()
  case_array = np.array(excur)
  row_count = case_array.shape[0]
  col_count = case_array.shape[1]
  red_ball_array = case_array[:,0]
  blue_ball_array = case_array[:,1]
  return red_ball_array,blue_ball_array,row_count,col_count

# 统计每个号码球的出现次数，这个应该在数据库里做，先放这，以后改
def every_ball_count(ball_array):
  ball_list = []
  for ball_char in ball_array:
    ball_list += ball_char.split(';')
  ball_count = {}
  for ball_num in ball_list:
    if ball_num in ball_count:
      ball_count[ball_num] += 1
    else:
      ball_count[ball_num] = ball_count.get(ball_num,0) + 1
  return ball_count

# 数据分析主函数，样本量必须大于等于7，否则不进行处理
def analysis_main(start_dt,end_dt):
  red_ball_array,blue_ball_array,row_count,col_count = get_ball_infomation(start_dt,end_dt)
  if row_count < 7:
    print '样本量不足以支持分析'
    return 1
  red_count_dict = every_ball_count(red_ball_array)
  blue_count_dict = every_ball_count(blue_ball_array)
  for red_case in red_count_dict:
    red_rate = binomial_class((red_count_dict[red_case] + 1),(row_count + 1),0.5)
    red_count_dict[red_case] = red_rate.binomial_fun()
  for blue_case in blue_count_dict:
    blue_rate = binomial_class((blue_count_dict[blue_case] + 1),(row_count + 1),0.5)
    blue_count_dict[blue_case] = blue_rate.binomial_fun()
  sorted_red_count = sorted(red_count_dict.iteritems(),key=operator.itemgetter(1),reverse=True)
  sorted_blue_count = sorted(blue_count_dict.iteritems(),key=operator.itemgetter(1),reverse=True)
  print sorted_blue_count[0]
  print '选择红球是:'
  n = 1
  for key,value in sorted_red_count:
    if n == 7:
      break
    print '%s,%s' %(key,str(value))
    n += 1
  print '选择蓝球是'
  print '%s,%s' %(sorted_blue_count[0][0],str(sorted_blue_count[0][1]))

if __name__ == '__main__':
  n = ''
  while n != '' or n != '':
    input_n = raw_input('''
请选择需要进行的功能
1、爬取页面的球号
2、进行球号分析
输入quit退出
请选择: ''')
    if input_n == '':
      url = raw_input('''
请输入需要爬取的地址(此为开始地址,因此建议选择页码较大的地址)
输入: ''')
      end_page = raw_input('''
输入结束页码
(注意:如果结束页码大于输入地址的页码,则不会爬取任何页面)
输入: ''')
      ball_main(url,end_page)
    elif input_n == '':
      analysis_main('2018-08-15','2018-09-09')
    elif input_n == 'quit':
      exit(0)

脚本共分两个功能块，一个是爬虫，用来爬取双色球号码，另一个是分析，使用二项分布对已抓取数据进行概率计算

先说说爬虫，依然是先观察页面

看过之前初级爬虫的同学应该对这个很熟悉，要爬的是1的地方，观察位置2和位置3，不难看出，每期占用一个页面，那么只要利用翻页，每次页码-1即可，下面看看1位置的源代码

在源码中搜索开奖号码，找到红圈位置，这个和伪代码中写的一样，现在看爬取正则

red_no_re = re.compile('(?<=\<li class\=\"ball_red\"\>)[0-9]+(?=\<\/li\>)')

blue_no_re = re.compile('(?<=\<li class\=\"ball_blue\"\>)[0-9]+(?=\<\/li\>)')

date_re = re.compile('(?<=开奖日期：)[0-9]+年[0-9]+月[0-9]+日')
分别爬取红球，篮球，日期
下面的事就好办了，将爬好的数据存入数据库，连接mysql，包含4个参数

db = 'pythondb'
  host = 'localhost'
  iuser = 'xxx'
  passwd = 'xxxxxx'

分别是连接数据库、连接地址、用户名、密码

因为在每期开奖之后都需要爬一次新的号码，所以我没有把这个做成一次性的，用户需要输入开始爬取页面，以及结束页面编号，当爬取到结束编号页面后则停止

为防止重复爬入数据，导致数据库数据重复，因此在数据入库前需要对库数据进行验证需要入库数据是否已存在，详见check_open_date函数

最后需要注意的是，在写入完毕后应关闭数据库连接，因为mysql默认自动commit，因此我这个没有单独commit的sql，如果是oracle不要忘记加这个

接下来是分析模块，二项分布为方便管理，单独写了个类，以后也可以单独调用，这个用法我之前的文里面有说明，就不再说了。这里有几个功能说明下

1、与爬虫一样，分析也是需要写入开始日期和结束日期的，且从库中获取case必须大于7个，否则考虑样本不足，无法分析

2、将每个号码在时间段的出现频次进行统计并+1，即二项分布的成功k次数

3、获取时间段内的期数并+1，即二项分布的实验次数

4、因为每个号码只有两种可能，要么出现，要么不出现，因此事件的成功率是0.5

结果是这样，看到概率最高的是2.34863516114e-08，这是一个科学技术法，大约是0.00000002349，也就是1/5000万，因此这个概率没有任何意义，当然你也可以尝试买2块钱的，买不中别赖我啊，哈哈

细心的同学可能注意到了，算出来的最后两个红球概率是一样的，举个例子C(5,2)其实和C(5,3)的组合数是一样的，那么二项分布的概率当然也一样，这样就需要用到泊松分布，加入历史平均值的计算，因为这个需要在库中单独建表，还没有时间做这个事，因此没有加入这项计算，如果各位有兴趣可以自己尝试把这块加上，等我有时间了，我也会把这个写出来完善下的