六爻蓍草的占卜方法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;aint 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个结果，然后怎么分析卦象我也不会了~~~~~~