关于“钻井分布”请帮我看看有什么错误

dasheng201 · 发表于 2005-8-27 23:45:23

//刚编的关于“钻井分布”问题二的程序 //改过N次，结果总是4口井，分别为：第2，第4，第5，第10口井； //而正确答案是第1,6,7,8,9,11口井； //高手帮忙看看程序有什么错误需要改正的 #include<math.h> #include<iostream.h> const float demand=0.05,pi=3.14;//demand表示小圆的半径 double px[12],py[12],ax[12],by[12], x[]={0.50,1.41,3.00,3.37,3.40,4.72,4.72,5.43,7.57,8.38,8.98,9.50}, y[]={2.00,3.50,1.50,3.51,5.50,2.00,6.24,4.10,2.01,4.50,3.41,0.80}; int n,s=0,a1,num,nearx[12],neary[12],w[12],well[12]; void ChangeXY(int i,float a,float x1,float y1) //对旧井的坐标作转换，存于px[],py[]中 { px=x1*cos(a*pi/180)-y1*sin(a*pi/180); py=x1*sin(a*pi/180)+y1*cos(a*pi/180); }

dasheng201 · 发表于 2005-8-27 23:46:03

<

>void main() { double dex,dey,a,x3,y3; int i; cout<<"the old well:"; cin>>n; s=0; for(a=0;a<=90;a+=0.1) //角度的旋转 { for(i=0;i<n;i++) { //对角度为a的情况下，对旧井点做坐标变换，存于px,py中 nearx=int(x); neary=int(y); ax=x-nearx; by=y-neary; ChangeXY(i,a,ax,by); } for(dex=0;dex<=1.0;dex+=0.01) //用小圆在x,y范围为1的坐标内搜索，每次步长为0.01 for(dey=0;dey<=1.0;dey+=0.01) { for(i=0;i<n;i++) w=0; //给w[]清零 num=0; //给num赋初值为零 for(i=0;i<n;i++) { if(((px-dex)*(px-dex)+(py-dey)*(py-dey))<=(demand*demand)) //如果上面的表达式成立，则旧井点(px,py)在圆内 { num++; w=1;//给num加1，标志w是进入圆内的点 } } if(num>s)//s存储可用的最多的旧井点数 { s=num; for(i=0;i<n;i++) well=w;//用于说明哪些些旧井可用 } } }//for cout<<"可以用的井的总数为："<<s<<endl; for(i=0;i<n;i++) if(well) cout<<"可以用的井的序列号是："<<i+1<<" "; cout<<endl; for(i=0;i<n;i++) cout<<well<<","<<endl;//看看well[]中那些存储的是1 cout<<"此时的角度为："<<a<<endl; }

		自动登录	找回密码
密码			注-册-帐-号