ioo.cpp

公路路线设计中用模式法计算插入圆心的坐标的小程序

#include<iostream.h>
#include<math.h>
double rr(double A,double R)
{
double T,t1;
t1=pow(A,4)/24/pow(R,3);
T=R+t1;
return T;
}

double tt(double A,double R)
{
double T;
T=pow(A,2)/2/R;
return T;
}

void main ()
{double A1,R1,T1,A2,R2,T2,A3,R3,T3,tt1,tt2,p,p1,d,d4,x1,x2,y1,y2,x,y,x3,y3,d1,d2,w1,w2,ddd,dd;
cout<<"输入A1,R1"<<endl;
cin>>A1>>R1;
 T1=rr(A1,R1);
cout<<"输入A2,R2"<<endl; 
cin>>A2>>R2;
 T2=rr(A2,R2);
cout<<"输入A3,R3"<<endl;
cin>>A3>>R3;
 T3=rr(A3,R3);

 tt1=tt(A1,R1);

 tt2=tt(A2,R2);
 
p=pow((T1+T2),2)+pow((tt1+tt2),2);
p1=pow(p,0.5);
cout<<"输入x1y1x2y2x3y3"<<endl;
cin>>x1>>y1>>x2>>y2>>x3>>y3;

dd=pow((x1-x2),2)+pow((y1-y2),2);
for(x=x1;x>x1-pow(A3,2)/R3;x=x-1){
	for(y=y1+R2;y<y1+T3;y=y+1){
d=pow((x-x3),2)+pow((y-y3),2);
d4=pow(d,0.5);
d1=pow((x-x1),2)+pow((y-y1),2);
d2=pow((x-x2),2)+pow((y-y2),2);
w1=d1-pow(T3,2);
w2=d2-pow(T3,2);
ddd=pow(w2,0.5)-pow(w1,0.5);
    if(fabs(p1-d4)<1&&fabs(ddd-pow(dd,0.5))<1)
	cout<<"x="<<x<<"   y="<<y<<endl;
	
	}

}
}