22222.cpp

是一个圆型工件的检测代码,但是有点错误。不知道如何改正。

#include<stdio.h>
#include<math.h>
void main()
{double z[36][2]={{0}},a,b,c,d,e,x,y;
int i,j,k,l,w=0;
for(i=0;i<36;i++)
scanf("%lf",&z[i][j]);
for(i=0;i<35;i++)
for(k=i+1;k<36;k++){
	if(z[i][0]-z[k][0]!=0){
		a=(z[i][0]*z[i][0]+z[i][1]*z[i][1]-z[k][0]*z[k][0]-z[k][1]*z[k][1])/(2*(z[i][0]-z[k][0]));
		b=(z[k][1]-z[i][1])/(z[i][0]-z[k][0]);
		c=b*b+1;
		d=a*b-z[i][0]*b-z[i][1];
		e=z[i][0]*z[i][0]+z[i][1]*z[i][1]+a*a-102.01-2*a*z[i][0];
		if(d*d-c*e==0){
			y=-d/c;x=a+b*y;
			for(l=0;l<36;l++)
			{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
}
if(w=36)
printf("合格");
}
}
} if(d*d-c*e>0)
{
y=(sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
y=(-sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
 else
{
if(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4==0)
{b=( z[i][1]+ z[k][1])/2;
a= z[i][0];
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else   if(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4>0)
{b=(z[i][1]+z[k][1])/2;
a= z[i][0]+sqrt(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
a=z[i][0]-sqrt(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
}
if(z[i][0]-z[k][0]!=0){
		a=(z[i][0]*z[i][0]+z[i][1]*z[i][1]-z[k][0]*z[k][0]-z[k][1]*z[k][1])/(2*(z[i][0]-z[k][0]));
		b=(z[k][1]-z[i][1])/(z[i][0]-z[k][0]);
		c=b*b+1;
		d=a*b-z[i][0]*b-z[i][1];
		e=z[i][0]*z[i][0]+z[i][1]*z[i][1]+a*a-98.01-2*a*z[i][0];
		if(d*d-c*e==0){
			y=-d/c;x=a+b*y;
			for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
}
if(w=36)
printf("合格");
}
}
else  if(d*d-c*e>0)
{
y=(sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
y=(-sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else
{
if(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4==0)
{b=(z[i][1]+z[k][1])/2;
a= z[i][0];
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else   if(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4>0)
{b=(z[i][1]+z[k][1])/2;
a= z[i][0]+sqrt(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
a=z[i][0]-sqrt(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
}
if(z[i][0]-z[k][0]!=0){
		a=(z[i][0]*z[i][0]+z[i][1]*z[i][1]-z[k][0]*z[k][0]-z[k][1]*z[k][1]+4)/(2*(z[i][0]-z[k][0]));
		b=(z[k][1]-z[i][1])/(z[i][0]-z[k][0]);
		c=b*b+1;
		d=a*b-z[i][0]*b-z[i][1];
		e=z[i][0]*z[i][0]+z[i][1]*z[i][1]+a*a-98.01-2*a*z[i][0];
		if(d*d-c*e==0){
			y=-d/c;x=a+b*y;
			for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
}
if(w=36)
printf("合格");
}
}
else  if(d*d-c*e>0)
{
y=(sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
y=(-sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else
{
if(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4==0)
{b=(z[i][1]+ z[k][1])/2;
a= z[i][0];
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else   if(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4>0)
{b=(z[i][1]+ z[k][1])/2;
a= z[i][0]+sqrt(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
a=z[i][0]-sqrt(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
}
if(z[i][0]-z[k][0]!=0){
		a=(z[i][0]*z[i][0]+z[i][1]*z[i][1]-z[k][0]*z[k][0]-z[k][1]*z[k][1]-4)/(2*(z[i][0]-z[k][0]));
		b=(z[k][1]-z[i][1])/(z[i][0]-z[k][0]);
		c=b*b+1;
		d=a*b-z[i][0]*b-z[i][1];
		e=z[i][0]*z[i][0]+z[i][1]*z[i][1]+a*a-102.01-2*a*z[i][0];
		if(d*d-c*e==0){
			y=-d/c;x=a+b*y;
			for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
}
if(w=36)
printf("合格");
}
}
else  if(d*d-c*e>0)
{
y=(sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
y=(-sqrt(d*d-c*e)-d)/c;
x=a+b*y;
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else
{
if(102.01-(z[i][1]-z[k][1]*z[i][1]-z[k][1])/4==0)
{b=( z[i][1]+ z[k][1])/2;
a= z[i][0];
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
else   if(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4>0)
{b=(z[i][1]+ z[k][1])/2;
a= z[i][0]+sqrt(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
a=z[i][0]-sqrt(102.01-(z[i][1]-z[k][1]* z[i][1]-z[k][1])/4);
for(l=0;l<36;l++)
{if(((z[l][0]-x)*(z[l][0]-x)-(z[l][1]-y)*(z[1][0]-y))<=102.01
&&((z[l][0]-x)* (z[l][0]-x)- (z[l][1]-y)*(z[1][0]-y))>=98.01)
w++;
else  break;
if(w=36)
printf("合格");
}
}
}
}