0和4组超强究极完整无误最新豪哥亲测完整版

1.#include <iostream>#include <string>#include <ctime>#include <iomanip>using namespace std;void ProgramInit(); // 显示本程序信息void DeckMakeDriver(); // 完成洗牌和删除牌的功能int GetRandInt(int min, int max); // 生成介于min 和max 之间的随机数void PlayCard(); // 玩牌对抗class CCard{private:char szFace[3]; // 代表1张扑克牌char chSuit; // 扑克牌的花色（黑、红、梅、方）int nNumber; // * * 牌面的大小（本成员数据用于对抗游戏）public:CCard(); // 构造函数void ShowCard(); // 显示1张牌void CardMake(int); // 根据参数生成1张牌int GetCardNum(); // * * 返回牌面的大小（本成员函数用于对抗游戏）};// 构造函数定义，初始化一些数据CCard::CCard(){int i;for(i=1; i<=3;) // 初始化扑克牌的牌面大小，初始为空{szFace[i] = ' ';i++;}chSuit = ' '; // 初始化牌面花色，初始为空}// * * 返回牌面的大小（本成员函数用于对抗游戏）int CCard::GetCardNum(){return nNumber;}// 显示牌面void CCard::ShowCard(){cout << setw(3) << szFace << setw(2) << chSuit << " ";}// 根据参数生成1张牌void CCard::CardMake(int num){int i = 0;char j;int face_num = num % 13; // 每种花色13张牌，face_num 的值的范围是：0~12 nNumber = face_num;switch( face_num ){case 0:strcpy(szFace, "A"); break; // 当牌面为非数字情况时case 9:strcpy(szFace, "10"); break;case 10:strcpy(szFace, "J"); break;case 11:strcpy(szFace, "Q"); break;case 12:strcpy(szFace, "K"); break;default: // 牌面为数字时，解释：当给出的参数是0时，49对应的ASCII 值为1，同理，参数为8时，值为9j = char(face_num + 49); // j 值的范围是：0~9，因为ASCII码表中只有0~9 这十个数字，所以把这范围的数字专门用这条语句来生成if ( i<3 ){szFace[i] = ' '; i++; // 使szFace[0] 值为' 'szFace[i] = j; i++; // 使szFace[1] 值为jszFace[i] = NULL; // 使szFace[0] 值为空break;}}// 每13张牌赋1种花色，结果就是A 到K，每张牌有4种花色，正好组成一套完整的扑克牌if ( num<=12 ) chSuit = 6; // ASCII 码中6代表的是：黑桃字符if ( num>12 && num<=25 ) chSuit = 3; // ASCII 码中3代表的是：红桃字符if ( num>25 && num<=38 ) chSuit = 5; // ASCII 码中5代表的是：梅花字符if ( num>38 && num<=51 ) chSuit = 4; // ASCII 码中4代表的是：方块字符}class CDeck{protected:CCard CardArray[52]; // 扑克牌类的对象数组，表示52张扑克牌int nLastDelt; // 临时表示最后一张牌的位数。

即52public:CDeck() {} // 构造函数void MakeDeck(); // 生成1副扑克牌void DealOne(); // 删除1张牌void ShowDeck(); // 显示这1副牌void DeckShuff(int); // 洗牌};// 产生1副牌（52）张void CDeck::MakeDeck(){nLastDelt = 51;while ( nLastDelt > -1 ){CardArray[nLastDelt].CardMake(nLastDelt); // 调用CCard 类成员函数CardMake()，然后把参数带进去，参数从51一直到0，同时依次赋值，赋的值从12 到1，共4次12 到1nLastDelt--;}} // 数组CardArray[52] 最终产生的值，第1-4位为不同花色的“K”，然后4位是“Q”，最后4位是“A”，即是从大到小放进去的，共52张牌// 依次显示牌面void CDeck::ShowDeck(){int nIndex = nLastDelt + 1; // 这时的nLastDelt 的值应当是0 （根据上一个函数）int newline = 0; // 排列扑克用的变量，每newline 张牌换一行cout << endl;while ( nIndex<=51 ){if ( newline%13==0 )cout << endl; // 13 张牌显示在一行CardArray[nIndex].ShowCard(); // 依次显示牌面newline++;nIndex++;}}// 洗牌次数void CDeck::DeckShuff(int times){int split; // split 是分开成两部分的位置，如上部分、下部分，即随机给两方分别发一副牌cout << "\n\t\t洗牌后……\n";for(int i=0; i<=times; i++){split = GetRandInt(20, 35); // 生成介于20 和35 之间的随机数，赋值给split，假如是31CDeck topdeck; // 1副扑克牌类的对象topdeck，代表上部分的牌CDeck bottomdeck; // 1副扑克牌类的对象bottomdeck，代表下部分的牌int bottomdeckindex = 1; // 记录上部分发了多少张牌int topdeckindex = 1; // 记录下部分发了多少张牌// 给两部分发牌，两个for语句使得第1部分获得了第30张（20到35的随机数）之前的所有牌，剩下的牌全给第2部分for(int j=0; j<=split-1; j++) // 给上部分的牌赋值{topdeck.CardArray[topdeckindex] = this->CardArray[j];topdeckindex++; // 最后一次记录的值是第30张}for(int k=split; k<52; k++) // 给下部分的牌赋值（剩下的所有牌赋给它）{bottomdeck.CardArray[bottomdeckindex] = this->CardArray[k];bottomdeckindex++; // 最后一次记录的值是第51-30=21 张}int deckoutindex = 0;int numcardstomovetop;int numcardstomovebottom;// 把bottomdeckindex 和topdeckindex 的牌的数量对调topdeckindex = split; // topdeckindex 的值为31张bottomdeckindex = 52 - split; // bottomdeckindex 值就是52-31=21 张while( deckoutindex <=51 ) // 每循环一次：数组this->CardArray[] 得到随机抽的两部分的2~7张牌，这样通过这两部分的牌的运做（第1部分抽来一点，第2部分抽过来一点）// 就使得每一次循环，this->CardArray[] 得到了随机的扑克牌排列次数，一直使this->CardArray[51] 数组的牌放满{numcardstomovetop = GetRandInt(2, 7); // 从上部分31张中，抽取的张数是2~7 的随机数numcardstomovebottom = GetRandInt(2, 7); // 从下部分21张中。

抽取的张数是2~7 的随机数for(int j=0; j<=numcardstomovebottom; j++) //从0 循环到21{if ( bottomdeckindex > 0 ) // 下部分有牌的情况下{// 将下部分的牌抽numcardstomovebottom 张给数组this->CardArray[]this->CardArray[deckoutindex] = bottomdeck.CardArray[bottomdeckindex]; // 从bottomdeck.CardArray[21] 赋给CardArray[0]开始，一直赋到bottomdeck.CardArray[0]deckoutindex++;bottomdeckindex--;}for(int l=0; l<=numcardstomovetop; l++){if ( (topdeckindex>0) && (deckoutindex<=52) ){// 将上部分的牌抽numcardstomovetop 张给数组this->CardArray[]this->CardArray[deckoutindex] = topdeck.CardArray[topdeckindex]; // 从topdeck.CardArray[31] 赋给CardArray[22]开始，一直赋到topdeck.CardArray[0]deckoutindex++;topdeckindex--;}}}}}// 这样连续250次循环以后，通过两部分的牌它抽过来一点它抽过来一点给this->CardArray[]，使得this->CardArray[51] 已经是排列的比较随机了this->nLastDelt = -1;}// 删除最后1张牌void CDeck::DealOne(){if ( nLastDelt != 51 ){nLastDelt++;cout << "\n 已删除：";CardArray[nLastDelt].ShowCard();}elsecout << "这副牌已全部被删完\n";}class CSuperDeck: public CDeck{private:int nLastCard; // 这副扑克牌中剩余的牌数public:// 构造函数，初始化剩余的牌数为52CSuperDeck(){nLastCard = 52;}int GetLastCard() // 返回剩余牌的数目{return nLastCard;}void DeleteDeckCard(int nNum); // 删除序号为nNum 的牌void ShowOneCard(int nNum); // 显示序号为nNum 的牌// 返回1张牌的牌面int GetOneCard(int nNum){return CardArray[nNum-1].GetCardNum();}};// 显示序号为nNum 的牌void CSuperDeck::ShowOneCard(int nNum){CardArray[nNum-1].ShowCard();}// 删除序号为nNum 的牌void CSuperDeck::DeleteDeckCard(int nNum){for(int i=nNum-1; i<nLastCard-1; i++)CardArray[i] = CardArray[i+1];nLastCard--;}int main(){srand( (unsigned)time(NULL) ); // 初始化随机数种子ProgramInit(); // 显示程序信息DeckMakeDriver(); // 主要步骤PlayCard(); // 对抗游戏cout << "\n\n\t\t谢谢使用本游戏，再见！\n\n";system("pause");return 0;}// 程序初始提示信息void ProgramInit(){cout << "\n\n\t\t欢迎进入扑克牌游戏！\n\n";cout << "\t\t本游戏完成洗牌和删除牌的功能\n\n";system("pause");}// 主要步骤void DeckMakeDriver(){CDeck deck; // 建立1副扑克牌对象deck.MakeDeck(); // 生成1副扑克牌deck.ShowDeck(); // 在屏幕上显示生成的牌cout << endl << endl;system("pause");deck.DeckShuff(250); // 洗牌250 次deck.ShowDeck(); // 显示洗完的牌cout << endl << endl;system("pause");string strChoice;cout << "是否删除1张牌（Y/N）？";cin >> strChoice;while( strChoice=="Y" || strChoice=="y" ){deck.DealOne(); // 删除1张牌cout << "\n\n\t余牌显示：";deck.ShowDeck(); // 余牌显示cout << endl << endl << "是否继续删除（Y/N）？";cin >> strChoice;}}// 取得介于min 和max 之间的随机数int GetRandInt(int min, int max){int numToReturn;numToReturn = rand();numToReturn = numToReturn % (max-min+1) + min;return numToReturn;}// 对抗游戏void PlayCard(){CSuperDeck deck; // 建立1副超级扑克牌对象deck.MakeDeck(); // 生成1副牌deck.ShowDeck(); // 在屏幕上显示生成的牌cout << endl << endl;system("pause");deck.DeckShuff(250); // 洗250次牌deck.ShowDeck(); // 显示洗完的牌cout << endl << endl;system("pause");string strChoice = "y";cout << "对抗游戏\n";system("pause");do // 开始对抗游戏{int nPlayer, nComputer; // 双方牌的序号int nPlayerNum, nComputerNum; // 双方牌的点数nPlayer = GetRandInt(1, deck.GetLastCard()); // 随机抽取1张牌cout << "您抽取的牌是：";deck.ShowOneCard(nPlayer); // 删除这张牌nPlayerNum = deck.GetOneCard(nPlayer); // 取得这张牌的点数deck.DeleteDeckCard(nPlayer);nComputer = GetRandInt(1, deck.GetLastCard()); // 随机抽取1张牌cout << "计算机抽取的牌是：";deck.ShowOneCard(nComputer); // 显示这张牌nComputerNum = deck.GetOneCard(nComputer); // 取得这张牌的点数deck.DeleteDeckCard(nComputer); // 删除这张牌if ( nComputerNum > nPlayerNum ) // 比较这两张牌的点数cout << "\n对不起，您输了，下次努力吧！\n\n";else if ( nComputerNum == nPlayerNum)cout << "\n平局！\n\n";elsecout << "\n恭喜，你赢了！\n\n";cout << "是否继续对抗（Y/N）？";cin >> strChoice;}while ( (strChoice=="Y" || strChoice=="y") && (deck.GetLastCard()!=0) );if ( deck.GetLastCard() == 0 ) // 扑克牌的数量为0cout << "\n\n\t\t牌已发完，欢迎下次继续！\n";cout << "\n\t\t对抗结束！\n";}2#include<cstdio>#include<cstdlib>#include<iostream>#include<cmath>#include<string>using namespace std;//定义CComplex命名空间开始namespace CComplex{class CComplex{private:double Real,Image;public:CComplex(double real=0,double image=0) //构造函数{Real=real;Image=image;}friend istream &operator >>(istream &is,CComplex &com); //重载输入friend ostream &operator <<(ostream &os,CComplex &com); //重载输出CComplex operator+(CComplex &com);//加法重载CComplex operator-(CComplex &com);//减法重载int operator ==(CComplex &com);};istream &operator>>(istream &is,CComplex &com)//重载输入{cout<<"输入复数:";is>>com.Real>>com.Image;return is;}ostream &operator <<(ostream &os,CComplex &com)//重载输出{if(fabs(com.Image)<1e-5)os<<com.Real;else if(com.Image>0)os<<com.Real<<"+"<<com.Image<<"*i";//虚部为正elseos<<com.Real<<com.Image<<"*i";return os;}CComplex CComplex::operator+(CComplex &com)//加法重载{CComplex sum;sum.Real=Real+com.Real;sum.Image=Image+com.Image;return sum;}CComplex CComplex::operator-(CComplex &com)//减法重载{CComplex sub;sub.Real=Real-com.Real;sub.Image=Image-com.Image;return sub;}int CComplex::operator==(CComplex &com){if(Real==com.Real &&Image==com.Image)return 1;elsereturn 0;}void Test(void)//{ cout <<"共10道题，作10以内的加减运算，满分100分:\n"; double real1,real2,image1,image2;CComplex answer,temp;int score=0;char op;for(int i=0;i<=9;i++){real1=rand()%10;image1=rand()%10;real2=rand()%10;image2=rand()%10;CComplex a(real1,image1),b(real2,image2);op=rand()%2;switch(op){case 0:answer=a+b;cout<<a<<"加"<<b<<"等于";break;case 1:answer=a-b;cout<<a<<"减"<<b<<"等于";break;}cin>>temp;if(answer==temp){score+=10;}else{cout<<"此题做错了\n";}}cout<<"您的最后得分是:"<<score<<endl;cout<<"按任意键继续\n";cout.flush();cin.get();cin.get();}void Add()//实现加法的函数{CComplex num1,num2,sum,Zero(0,0);cout<<"加法计算\n"<<"最少输入两个数，并且以0结束\n";cout<<"第1个数，";cin>>num1;cout<<"第2个数,";cin>>num2;sum=num1+num2;cout<<"第3个数:";cin>>num1;int i=4;while(!(num1==Zero)){sum=sum+num1;cout<<"第"<<i<<"个数:";i++;cin >>num1;}cout<<"结果是:"<<sum<<endl;cout<<"按任意键继续\n";cout.flush();cin.get();cin.get();}void Sub()//实现减法的函数{CComplex num1,num2,sub,Zero(0,0);cout<<"最少输入两个数，并且以0结束\n";cout<<"第1个数:";cin>>num1;cout<<"第2个数:";cin>>num2;sub=num1-num2;cout<<"第3个数:";cin>>num1;int i=4;while(!(num1==Zero)){sub=sub-num1;cout<<"第"<<i<<"个数:";i++;cin>>num1;}cout<<"减法运算的结果是:"<<sub<<endl;cout<<"A按任意键继续\n";cin.get();cin.get();}}//定义CComplex命名空间结束using namespace CComplex;int main(void)//进入程序后显示菜单{char strChoice[20];do{system("cls");//执行系统命令：cls为清屏cout<<"\t一个简单的复数计算器程序，实现以下功能，按对应的键(1~9)进入\n\n\n";cout<<"\t==========MENU==========\n";cout<<"\t1:多操作数加法，以0结束\n";cout<<"\t2:多操作数减法，以0结束\n";cout<<"\t3:测试10以内的加减乘法运算，1次测试10道题\n";cout<<"\t0:退出计算器程序\n\n";cout<<"\t请输入您的选择:";cin>>strChoice;if(strcmp(strChoice,"1")==0)Add();else if(strcmp(strChoice,"2")==0)Sub();else if(strcmp(strChoice,"3")==0)Test();else if(strcmp(strChoice,"0")==0){cout<<"\n\n\t欢迎下次继续使用复数计算器！\n\n";break;}else{cout<<"\n\t输入错误，请按任意键后重新输入！\n";cin.get();cin.get();}}while((strcmp(strChoice,"0")));return 0;}3.#include<iostream>#include<fstream>#include<cstdlib>#include<cassert>#include<cmath>using namespace std;class Rand{protected:long num; //平均分布的种子float r; //平均分布值public:Rand();float ran2(long&); //返回0~1之间的平均分布值float NewRand();float GetRand();};Rand::Rand(){num=-1.0;ran2(num);r=ran2(num);}float Rand::GetRand() {return r;}float Rand::NewRand() {r=ran2(num);return r;}float Rand::ran2(long&idum) //利用已知公式计算0~1之间平均分布的值{const long IM1=2147483563;const long IM2=2147483399;const double AM=(1.0/IM1);const long IMM1=(IM1-1);const long IA1=40014;const long IA2=40692;const long IQ1=53668;const long IQ2=52774;const long IR1=12211;const long IR2=3791;const long NTAB=32;const long NDIV=(1+IMM1/NTAB);const double EPS=1.2e-7;const double RNMX=(1.0-EPS);int j;long k;static long idum2=123456789;static long iy=0;static long iv[NTAB];float temp;if(idum<=0){if(-(idum)<1)idum=1;elseidum=-(idum);idum2=(idum);for(j=NTAB+7;j>=0;j--){k=(idum)/IQ1;idum=IA1*(idum-k*IQ1)-k*IR1;if(idum<0)idum+=IM1;if(j<NTAB)iv[j]=idum;}iy=iv[0];}k=idum/IQ1;idum=IA1*(idum-k*IQ1)-k*IR1;if(idum<0)idum+=IM1;k=idum2/IQ2;idum2=IA2*(idum2-k*IQ2)-k*IR2;if(idum2<0)idum2+=IM2;j=iy/NDIV;iy=iv[j]-idum2;iv[j]=idum;if(iy<1)iy+=IMM1;if((temp=AM*iy)>RNMX)return RNMX;elsereturn temp;}class Gauss:public Rand{protected:float g; //高斯分布值float mean; //均值float sigma; //高斯分布参数int count;public:Gauss(float m=0,float s=1);float box_muller(float,float);float NewGauss();float GetGauss();int GetCount();void IncrementCount();};Gauss::Gauss(float m,float s):Rand() //构造函数{mean=m;sigma=s;g=box_muller(mean,sigma);count=0;}float Gauss::NewGauss(){g=box_muller(mean,sigma);count=0;return g;}float Gauss::GetGauss(){return g;}int Gauss::GetCount(){return count;}void Gauss::IncrementCount(){count++;}float Gauss::box_muller(float m,float s) //根据已知参数求出具体符合高斯分布的值{float x1,x2,w,y1;static float y2;static int use_last=0;if(use_last){y1=y2;use_last=0;}else{do{x1=2.0*NewRand()-1.0;x2=2.0*NewRand()-1.0;w=x1*x2+x2*x2;}while(w>=1.0);w=sqrt((-2.0*log(w))/w);y1=x1*w;y2=x2*w;use_last=1;}return(m+y1*s);}class Runge{protected:double*y; //位移量的指针double*yp; //一阶微分的指针double*t; //时间的指针int num_steps; //所计算的数据个数int size; //矩阵的维数double h; //计算的时间步长public:Runge(double h=0.1,double tmax=1,double y0=0,double yp0=1,double t0=0);Runge(const Runge&);~Runge();double GetY(int n)const; //返回成员函数double GetYP(int n)const;double GetT(int n)const;int GetSize()const;double GetStepSize()const;void rk(double(*func)(double,double,double)); //Runge-Kutta算法void PrintT(); //输出时间void PrintY(); //输出位移void Print(); //输出随着时间的位移值double func(double,double,double); //计算常微分方程y''=f(x,y,y') };//构造函数Runge::Runge(double delta_t,double tmax,double y0,double yp0,double t0) {num_steps=tmax/delta_t; //总数据量size=num_steps+1; //向量维数y=new double[size]; //开辟空间，存放数据assert(y!=0);yp=new double[size];assert(yp!=0);t=new double[size];assert(t!=0);h=delta_t; //步长y[0]=y0; //常微分方程解的初值yp[0]=yp0; //常微分方程一阶初值t[0]=t0; //时间初值for(int i=1;i<size;i++) //时间的离散值t[i]=t[0]+i*h;}Runge::Runge(const Runge&from) //被拷贝的构造函数{num_steps=from.num_steps;size=from.size;h=from.h;y=new double[size];yp=new double[size];t=new double[size];for(int i=0;i<size;i++){y[i]=from.y[i];yp[i]=from.yp[i];t[i]=from.t[i];}}Runge::~Runge() //析构函数{delete[]y;delete[]yp;delete[]t;}//利用四阶Runge-Kutta方法求解常微分方程y''=f(x,y,y')void Runge::rk(double(*func)(double,double,double)){double k1,k2,k3,k4;for(int n=0;n<size-1;n++){k1=h*func(t[n],y[n],yp[n]);k2=h*func(t[n]+0.5*h,y[n]+0.5*h*yp[n]+1.0*h/8*k1,yp[n]+0.5*k1);k3=h*func(t[n]+0.5*h,y[n]+0.5*h*yp[n]+1.0*h/8*k1,yp[n]+0.5*k2);k4=h*func(t[n]+h,y[n]+h*y[n]+0.5*k3,yp[n]+k3);y[n+1]=y[n]+h*(yp[n]+1.0*1/6*(k1+k2+k3));yp[n+1]=yp[n]+1.0*1/6*(k1+2*k2+2*k3+k4);}}void Runge::PrintT() //输出离散时间{for(int l=0;l<size;l++)cout<<t[l]<<endl;cout<<"\n"<<endl;void Runge::PrintY(){for(int i=0;i<size;i++)cout<<y[i]<<endl;cout<<"\n"<<endl;}void Runge::Print(){ofstream output("runge.out");if(! output){cerr<<"File could not be opened"<<endl;exit(1);}for(int i=0;i<size;i++)output<<t[i]<<" "<<y[i]<<endl;output.close();}double Runge::GetY(int nn)const{assert(nn<size);return y[nn];double Runge::GetYP(int nn)const{assert(nn<size);return yp[nn];}int Runge::GetSize()const{return size;}double Runge::GetT(int k)const{return t[k];}double Runge::GetStepSize()const{return h;}//程序中用到的常量const double pi=3.1415927;const double K_B=1.38066e-23; //玻耳兹曼常数const double T=300; //绝对温度const double D=1e-7; //布朗粒子的半径const double n=0.001; //液体的粘滞度const double step=0.01; //时间递增的步长const double tau=step;const double m=1e-8; //粒子的质量const double sig=sqrt(6*pi*K_B*T*D*n/tau); //高斯分布的取值范围double f(double,double,double); //朗之万方程的解int main(){Runge brownian(step,1000,0.0,0.0,0.0); //生成布朗运动的对象double(*f_pointer) (double,double,double)=f;brownian.rk(f); //用Runge-Kutta方法求解二次线性朗之万方程brownian.Print(); //输出温度与对应的位移//变量的初值double sum=0; //总位移量double kinetic=0; //求解朗之万方程的平均能量double theory=0; //平均能量的理论值int count=0; //计数double Einstein=0; //爱因斯坦系数的理论值（理论位移量）double dist=0; //平均位移量double diff=0; //平均能量的相对误差double diff2=0; //平均位移的相对误差//判断是否满足布朗运动的能量关系while((count<brownian. GetSize())&&(brownian.GetYP(count)<sqrt(K_B*T/m))) count++;//计算布朗运动粒子的平均能量for(int i=count;i<brownian.GetSize();i++){sum=sum+brownian.GetYP(i)*brownian.GetYP(i);kinetic=0.5*m*sum/(brownian.GetSize()-count);}theory=0.5*K_B*T; //计算平均能量的理论值for(int j=0;j<brownian.GetSize();j++) //计算实际平均位移dist=dist+brownian.GetY(j)*brownian.GetY(j);dist=dist/brownian.GetSize();//这是布朗运动的爱因斯坦公式：求出的是布朗粒子位移的散差x^2（理论上平均位移）Einstein=(double)2*K_B*T/(3*pi*n*D)*brownian.GetT(brownian.GetSize()-1);diff=(kinetic-theory)/theory; //计算平均能量实际值与理论值的相对误差if(diff<0)diff=-diff;diff*=100;diff2=(Einstein-dist)/Einstein; //计算平均位移实际值与理论值的相对误差if(diff2<0)diff2=-diff2;diff2*=100;ofstream output2("runge.out2"); //将结果输出到文件中if(! output2){cerr<<"File could not be opened"<<endl;exit(1);}cout<<"The average kinetic energy of the particle is"<<kinetic<<"\n"<<endl;cout<<"The theoretical value is 5*K*T="<<theory<<"\n"<<endl;cout<<"The average squared displacement is"<<dist<<"\n"<<endl;cout<<"The Einstein-Stokes formula predicts a value of"<<Einstein<<"\n"<<endl;cout<<"The percent difference in energy is"<<diff<<"%"<<"\n"<<endl;cout<<"The percent difference in the squared displancement is"<<diff2<<"%"<<"\n"<<endl;cout<<"brownian.GetT(brownian.GetSize()-1)="<<brownian.GetT(brownian.GetSize( )-1)<<"\n"<<endl;cout<<"count="<<count<<endl;output2.close();}//朗之万方程的解，布朗粒子在液体中符合高斯分布，所以在程序中模拟高斯分布给出各个粒子的位置double f(double t,double y,double yp){static Gauss gauss_random(0,sig); //生成一次高斯分布的种子if(gauss_random.GetCount()%4==0)//由于解方程需要用四阶Runge-Kutta方法，所以每4次都调用相同的高斯分布值gauss_random.NewGauss();gauss_random.IncrementCount();return -3*pi*D*n*yp/m+gauss_random.GetGauss()/m; //返回朗之万方程的解}3.4.4 股票交易管理#include<iostream>#include<cstdlib>#include<string>#include<fstream>#include<iomanip>using namespace std;#define MAX 5 //可显示5支股票的交易状况class Stock //股票类{friend class Customer;private:char StockName[20]; //股票名称char StockCode[6]; //股票代号char choice;int NumOfStock; //索引long int StockVolume; //股票的总数long int VolumeAvailable; //可交易的数量public:Stock(); //构造void Interface(); //股市界面void StockExchangeMarket(); //输出信息void SwitchChoice(); //选择命令执行void InputForNewStock(); //输入新股票信息};Stock s[MAX],temp; //5支股票对象class Customer //用户类{private:char CustomerName[30]; //用户名char Password[10]; //密码char ShareHoldingName[30]; //拥有的股票名称char ShareHoldingCode[5]; //拥有的股票代码char AnyChar; //存放任意字符变量long int ShareHoldingValue;//拥有的股票数量public:Customer();void Login(char *,Customer)const; //用户登录void Register(Customer *); //用户开户void StockPortfolio(Customer)const;//用户交易操作};Customer::Customer() //默认构造，清零{ strcpy(Customer::CustomerName," ");strcpy(Customer::ShareHoldingName," ");strcpy(Customer::ShareHoldingCode," ");ShareHoldingValue=0;}void Customer::Login(char *password,Customer cust) const { ifstream read_customer("customer.dat"); //读入数据文件while (!read_customer.eof()){read_customer.read((char *)(&cust),sizeof(cust));if(read_customer.fail()){cout<<"No record at all!";break;}if(strcmp(password,cust.Password)==0) //验证密码StockPortfolio(cust); //密码正确,进入交易操作}read_customer.close();}void Customer::Register(Customer *cust) //用户开户{system("cls"); //执行系统命令;清屏cin.get();cout<<"\t\t欢迎到交易市场开户!"<<endl;cout<<"~~~~~~~~~~~~~~~~~~~"<<endl;cout<<"请输入您的密码(少于8个字符):";cin>>Password;cin.get();system("cls");cout<<endl<<endl;cout<<"您的帐户已被建立,可以开始交易......"<<endl;cout<<"按任意键开始......"<<endl;//cin.getch();cin>>AnyChar;//getch()}void Customer::StockPortfolio(Customer cust)const //用户股市操作{Stock share[MAX],temp;Customer s[MAX];int i;char ch,AnyChar,ShareCode[6];long int volume;start:system("cls");cout<<endl<<endl;cout<<"\t\t\t 欢迎返回!"<<endl;cout<<"\t\t 股票交易系统"<<endl;cout<<"\t\t~~~~~~~~~~~~~~~~~~~"<<endl;cout<<"\t\t 用户："<<cust.CustomerName<<endl<<endl;cout<<"\t\t 买股票.................................[b]"<<endl;cout<<"\t\t 卖股票.................................[s]"<<endl;cout<<"\t\t 查看用户帐户...........................[v]"<<endl;cout<<"\t\t 登出...................................[q]"<<endl;cout<<endl<<"\t\t 输入您的选择:";cin>>ch;if((ch!='v')&&(ch!='b')&&(ch!='s')&&(ch!='q'))goto start;else{switch(ch){case 'b'://买股票{system("cls");cout<<"\n\n\t\t股票交易系统"<<endl;cout<<" 名称\t\t交易代码\t在册股数\t可交易股数"<<endl<<endl;ifstream dataFile("Stock_File.dat"); //输入股票市场情况dataFile.read((char *)(&share),sizeof(share)); //将股票市场资料读入share ifstream DataFileForShare(cust.CustomerName); //用户目前资料DataFileForShare.read((char *)(&s),sizeof(s)); //将用户目前资料读入内存for(i=0;i<MAX;++i)share[i].StockExchangeMarket(); //显示股票市场行情cin.get();cout<<endl<<"输入您购买的股票的交易代码:";cin.getline(ShareCode,6);cout<<"输入您购买的股数:";cin>>volume;cin.get();i=0;while((strcmp(share[i].StockCode,ShareCode)==0)||(i<MAX)){if (strcmp(share[i].StockCode,ShareCode)==0)//输入了正确的股票代码{if(share[i].VolumeAvailable>volume) //满足交易条件{s[i].ShareHoldingValue=s[i].ShareHoldingValue+volume;//用户持有的股票数增加strncpy(s[i].ShareHoldingName,share[i].StockName,30);strcpy(s[i].ShareHoldingCode,share[i].StockCode);share[i].VolumeAvailable=share[i].VolumeAvailable-volume;//股票市场上的股票数相应减少ofstream dataFile (cust.CustomerName,ios::binary);//将交易情况存储在客户文件中for(i=0;i<5;++i)dataFile.write((char *)(&s[i]),sizeof(s[i]));ofstream file("Stock_File.dat"); //将交易结果保存在股市文件中for(i=0;i<5;++i)file.write((char *)(&share[i]),sizeof(share[i]));break;}else{cout<<share[i].StockName<<"的股数不足."<<endl;cout<<"交易无效......";cin>>AnyChar;break;}}else{i++;while(i==MAX){cout<<"交易代码错误......"<<endl;cout<<"交易无效......";cin>>AnyChar;;break;}}}break;}case 's': //卖股票{system("cls");ifstream DataFileForShare(cust.CustomerName); //打开用户资料DataFileForShare.read((char *)(&s),sizeof(s));ifstream dataFile("Stock_File.dat");dataFile.read((char *)(&share),sizeof(share));cout<<"\t\t股票交易系统"<<endl;cout<<" 证券名称\t\t交易代码\t持有股数"<<endl;for(int i=0;i<MAX;i++)cout<<" "<<s[i].ShareHoldingName<<"\t\t"<<s[i].ShareHoldingCode<<"\t"<<s[i].Share HoldingValue<<endl;cin.get();。