急需钻井问题的解法！！

槐树林 · 发表于 2005-7-26 04:47:35

<

>各位高手好啊！
<

>最近在做99年的《钻井问题》的题目。
<

>我感觉是用0，1分布来做，可是求解好象很麻烦哦！有那位能帮助一下啊。

wlsui_wls · 发表于 2005-7-31 05:32:30

我前两天也刚做了这个题目，编程很好实现，做几个循环就行了

hui533 · 发表于 2005-8-26 05:26:06

是的，你可以用全局搜索或局部搜索

hui533 · 发表于 2005-8-26 05:27:13

你可以用全局搜索或局部搜索

dasheng201 · 发表于 2005-8-27 23:49:16

<

>wlsui_wls
<

>能否把你编的程序贴出来给我看看
<

>或者发到我的信箱<a href="mailtdasheng0077@163.com" target="_blank" >dasheng0077@163.com</A>
谢谢

dasheng201 · 发表于 2005-8-29 18:36:54

<

>问题一的程序：
<

>#include <iostream.h> #include <math.h> const float demand=0.05; int num,well[12],w[12];//num表示计数点，well[]表示有几口井可用，w[]表示标志哪口井可用 float x[12]={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[12]={2.00,3.50,1.50,3.51,5.50,2.00,6.24,4.10,2.01,4.50,3.41,0.80}; void main() // { int i,n,a[12],b[12],sum;//sum表示计数的总数 float dex,dey,ax[12],by[12];//dex,dey表示小正方形的中心坐标 sum=0; cout<<"请输入总的井数:"; cin>>n;
<

>for(dex=0;dex<=1.0;dex+=0.01) for(dey=0;dey<=1.0;dey+=0.01) { num=0; for(i=0;i<n;i++) { a=int(x);//提取出x[],y[]的整数部分 b=int(y); ax=x-a; by=y-b;//提取小数部分 } for(i=0;i<n;i++) w=0;//给w[]赋初值
 for(i=0;i<n;i++) { if(ax-dex<=demand && by-dey<=demand&& dex-ax<=demand && dey-by<=demand) { num++; w=1; } } if(num>sum)//num对应是w[],sum对应的是well[] { sum=num; for(i=0;i<n;i++) well=w; } } cout<<"可用的旧井是:"; for(i=0;i<n;i++) if(well) cout<<i+1<<" , "; cout<<endl; cout<<"可用的旧井的总数是："<<sum<<endl; }

dasheng201 · 发表于 2005-8-29 18:37:40

<

>问题二的程序：
<

>#include<math.h> #include<iostream.h> const double demand=0.05,pi=3.1415926; float 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,num,nearx[12],neary[12],w[12],well[12]; void ChangeXY(int i,double a,double x1,double y1) { px=x1*cos(a*pi/180)-y1*sin(a*pi/180); py=x1*sin(a*pi/180)+y1*cos(a*pi/180); } float sqr(double p) { return p*p; } void main() { double dex,dey,a,a1,dex1,dey1; int i; cout<<"旧井的总数是："; cin>>n; s=0; for(a=0;a<=90;a+=0.01)// { for(i=0;i<12;i++) { ChangeXY(i,a,x,y);
<> nearx=int(px); neary=int(py); if(px<0) { ax=abs(nearx)+px+1; } else ax=px-nearx; if(py<0) { by=abs(neary)+py+1; } else by=py-neary; } for(dex=0;dex<=1.0;dex+=0.01)// for(dey=0;dey<=1;dey+=0.01)// { for(i=0;i<n;i++) w=0; num=0; for(i=0;i<n;i++) { if(sqrt(sqr(ax-dex)+sqr(by-dey))<=demand) { num++; w=1; } } if(num>=s) { s=num; for(i=0;i<n;i++) well=w; a1=a; dex1=dex; dey1=dey; } } }//for(a=...) cout<<"可以利用的旧井序列号为："; for(i=0;i<n;i++) if(well) cout<<i+1<<" "; cout<<endl; cout<<"可以利用的旧井的总数为："<<s<<endl; cout<<"原点是：("<<dex1<<","<<dey1<<")"; }

14282578 · 发表于 2005-8-31 05:27:32

感谢

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

急需钻井问题的解法！！

不要灰心