六爻蓍草的占卜方法C++随机数实现

封神电影里周文王的占卜方法 ，占卜一次的结果需要很长时间，C++生成随机数算出六爻结果

http://www.360doc.com/content/12/0121/07/1071374691_1071374691.shtml

#include <random>
#include <iostream>
#include <time.h>
#include <stdlib.h>
#include <TH1I.h>
using namespace std; 
#define N 10000000
int course(int A)
{
  
  int A1=ceil(rand()%N/(float)N*(float)A);
  //cout<<A1<<endl;a
  int A2=A-A1;
  int a1=0;
  int a2=0;
  if((A1-1)%4==0)
  {
    a1=4;
  }
  else a1=(A1-1)%4;
  if(A2%4==0)
  {
    a2=4;
  }
  else a2=A2%4;
  
  int result=a1+a2+1;
  if(A1==0)  result=0;
  
  return result;
  
  //cout<<A1<<" "<<A2<<" "<<a1<<" "<<a2<<endl;
}
int check(int yao,int X)
{
  //srand(time(NULL));
  if(yao==0)  
  {
    yao=course(X);
    check(yao,X);
  }
  if(yao>0)  yao=yao;
  return yao;
  //check(yao,X);
  
}

void sixyao6()
{
 // srand((unsigned int)time(NULL));
  int yao1=0;
  int yao2=0;
  int yao3=0;
  int X1=50-1;
  int X2=X1-check(yao1,X1);
  int X3=X2-check(yao2,X2);
  int X4=X3-check(yao3,X3);
  
  cout<<X4/4<<endl;
  //return X4/4;
}



void sixyao()
{
  
  TH1I *h=new TH1I("h","h",20,0,10);
  int a=0;
  int X=50-1;
  srand(time(NULL));
  for(int i=0;i<1000000;i++)
  {      
    //srand(time(0));
    
    int yao1=0;
    int yao2=0;
    int yao3=0;
    int X1=0;
    int X2=0;
    int X3=0;
    int X4=0;
    X1=50-1;
    X2=X1-check(yao1,X1);
    X3=X2-check(yao2,X2);
    X4=X3-check(yao3,X3);
 //   cout<<X4/4<<endl;
    
    //cout<<rand()<<endl;
    if(X4/4%2==0)	h->Fill(1);
    else	h->Fill(2);
    h->Fill(X4/4);
  }
  h->Draw();
}

srand(time(NULL));设置随机数种子，使用当前时间进行随机数发生器的初始化。

循环10000000次，将算出来的结果fill到直方图中，结果是6、7、8、9这四个数字，奇数偶数比例各占50%，阴爻阳爻的比例相当。

使用时运行6次sixyao6()，可以得到6个结果，然后怎么分析卦象我也不会了~~~~~~