马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- #include <time.h>
- #include <math.h>
- #include <stdio.h>
- #include <stdlib.h>
- #include <iostream.h>
- #define pi 3.1415926
- #define N 50000
- void main()
- {
- double t[100][100];
- int s[1000];
- int z[7]={0,0,0,0,0,0,0};
- long a[100000];
- int m=0,n;
- double x;
- double b[1000];
- double DET=0,MTBT;
- for (int j=0;j<N;j++)
- {
- for (int i=0;i<7;i++)
- {
- time_t T;
- time(&T);
- srand(T);
- double b1,c1;
- double a1=(rand()/(RAND_MAX*1.0));
- switch (i)
- {
- case 0:
- t[i][j]=-2500*log(a1);
- break;
- case 1:
- b1=(rand()/(RAND_MAX*1.0));
- c1=(sqrt((-2)*log(a1)))*sin(2*pi*b1);
- break;
- case 2:
- t[i][j]=-2500*log(a1);
- break;
- case 3:
- t[i][j]=-2500*log(a1);
- break;
- case 4:
- b1=(rand()/(RAND_MAX*1.0));
- c1=(sqrt((-2)*log(a1)))*sin(2*pi*b1);
- break;
- case 5:
- t[i][j]=-2500*log(a1);
- break;
- case 6:
- t[i][j]=-2500*log(a1);
- break;
- }
- for (int k=0;k<7;k++)
- {
- for (int l=k;l<=k;l++)
- {
- z[l]=1;
- s[l]=(z[1]&&z[2]||z[3])&&z[6]||((z[4]||z[5])&&z[7]);
- if(s[l]==1)
- {
- a[m]=l;
- b[m]=t[l][j];
- m++;
- }
- }
- }
- for (k=0;k<m-2;k++)
- if (b[k]<b[k+1])
- {
- x=b[k];
- b[k]=b[k+1];
- b[k+1]=x;
- n=a[k];
- a[k]=a[k+1];
- a[k+1]=n;
- }
- }
- }
- DET=+b[0];
- MTBT=DET/N;
- cout<<MTBT<<endl;
- }
复制代码 帮忙看一下,出错的地方啊... |