分治法解决凸包问题（用C语言递归调用实现）

#include<stdio.h>
#include<math.h>
static float z[100][2];
shangtubao(float a[100][2])
{ int i,j,p,q,m,n,t;
   float b[100][2]={0},c[100][2]={0},d[100][2]={0},e[100][2]={0};
   for(i=0;i<100;i++)
{  if(a[i][0]==0&&a[i][1]==0) 
  break;
    else
  for(j=0;j<2;j++)
e[i][j]=a[i][j];
     }
     n=0;
for(i=0;i<100;i++)
 { if(a[i][0]!=0||a[i][1]!=0)
 n++;
 else 
  break;
 }
p=0;
 m=n-1;
 j=1;
 b[0][0]=a[0][0];
 b[0][1]=a[0][1];
 for(i=1;i<n-1;i++)
 { if((a[0][0]*a[m][1]+a[i][0]*a[0][1]+a[m][0]*a[i][1]-a[i][0]*a[m][1]-a[m][0]*a[0][1]-a[0][0]*a[i][1])>0)
  {  b[j][0]=a[i][0];
 b[j][1]=a[i][1];
 if((a[0][0]*a[m][1]+a[i][0]*a[0][1]+a[m][0]*a[i][1]-a[i][0]*a[m][1]-a[m][0]*a[0][1]-a[0][0]*a[i][1])>p)

 {p=a[0][0]*a[m][1]+a[i][0]*a[0][1]+a[m][0]*a[i][1]-a[i][0]*a[m][1]-a[m][0]*a[0][1]-a[0][0]*a[i][1];
   q=j;
 }
  j++;
  }
  }
  b[j][0]=a[m][0];
  b[j][1]=a[m][1];
  j++;
 if(j>2)
 {  for(i=0;i<100;i++)
    {  if(z[i][0]==0&&z[i][1]==0)
     {  z[i][0]=b[q][0];
         z[i][1]=b[q][1];
    break;
     }
    else  continue;
    }
  for(i=0;i<=q;i++) 
    { c[i][0]=b[i][0];
  c[i][1]=b[i][1];
    }
    for(i=0;i<100;i++)
{  if(a[i][0]==0&&a[i][1]==0) 
  break;
    else
       for(t=0;t<2;t++)
   a[i][t]=c[i][t];
}
    shangtubao(a);