学步园

#include <iostream><br /> #include <numeric><br /> #include <vector><br /> #include <cstdlib><br /> #include <fstream><br /> #include <ctime><br /> //很好先stl后其他library<br /> #include <boost/random/linear_congruential.hpp><br /> #include <boost/random/uniform_int.hpp><br /> #include <boost/random/uniform_real.hpp><br /> #include <boost/random/variate_generator.hpp><br /> #include <boost/random/mersenne_twister.hpp></p> <p>using namespace std;<br /> int main(int argc, char* argv[])<br /> {<br /> const int N=100; // 基础群大小<br /> const int s=4; //公母比例nf/nm<br /> const int n0=8; //每头母畜每胎所产后代数<br /> const double P1=0.16,P2=0.48,P3=0.36; //基因型频率，11，12，22<br /> const int g=20; //考察世代数<br /> int S[N]; // 个体性别（1为雄，2为雌）<br /> int G[N]; //基因型（11，12，22）<br /> double u1,u2;<br /> //产生服从（0，1）之间均匀分布的随机数<br /> //随机数生成<span class='wp_keywordlink'><a href="http://www.xuebuyuan.com/category/%E7%AE%97%E6%B3%95" title="算法" target="_blank">算法</a></span><br /> boost::mt19937 rng;<br /> //一种浮点数分布算法<br /> boost::uniform_real<> one(0,1);<br /> //得到需要的浮点数生成器，它是一个functor<br /> boost::variate_generator<boost::mt19937&, boost::uniform_real<> > die(rng, one);<br /> // 确定个体性别与基因型<br /> // for (int i=0;i<N;i++)<br /> // {<br /> // u1 = die(); // simulate rolling a die<br /> // if (u1<s/(1+s))<br /> // S[i]=1;<br /> // else<br /> // S[i]=2;<br /> // }<br /> //<br /> // //个体基因型确定（随机法）<br /> // for (int i=0;i<N;i++)<br /> // {<br /> // u2=die();<br /> //<br /> // if (u2<P1)<br /> // G[i]=11;<br /> // else {/*如果代码块是多行，请加上括号*/<br /> // if (u2>=P1&&u2<P1+P2)<br /> // G[i]=12;<br /> // else<br /> // G[i]=22;<br /> // }<br /> // }<br /> //不觉得可以一次循环搞定么？<br /> int nm=0,nf=0;<br /> for (int i=0;i<N;i++)<br /> {<br /> u1 = die(); // simulate rolling a die<br /> if (u1<s/(1+s)) {<br /> S[i]=1;<br /> ++nm;<br /> }<br /> else<br /> S[i]=2;<br /> u2=die();<br /> if (u2<P1)<br /> G[i]=11;<br /> else {/*如果代码块是多行，请加上括号*/<br /> if (u2>=P1&&u2<P1+P2)<br /> G[i]=12;<br /> else<br /> G[i]=22;<br /> }<br /> }<br /> nf=N-nm;<br /> // 随机交配，结果记录到二维数组中<br /> // 产生服从U(1,nf)的随机数u,随机挑选母畜<br /> int u;<br /> int M[nm][s]; //存放配对的数组<br /> //不需要重复声明<br /> //boost::mt19937 rng; // produces randomness out of thin air<br /> //需要重新声明一个分布算法 因为field变了<br /> boost::uniform_int<> field(1,nf); // distribution that maps to 1..6<br /> //于是分布算法+随机数算法生成的生成器对象变了，我叫die2<br /> boost::variate_generator<boost::mt19937&, boost::uniform_int<> > die2(rng, field); // glues randomness with mapping<br /> //加上括号<br /> for (int i=1;i<=nm;i++) {<br /> for (int j=1;j<=s;j++)<br /> {<br /> u = die2(); // simulate rolling a die<br /> cout<<u<<endl;<br /> M[i][j]=u;<br /> }<br /> }<br /> return 0;<br /> }