#include<stdio.h>
#include<math.h>
//#include<malloc.h>
#include<stdlib.h>
void heapsort(double **data,int i,int length);

main()
{
        int m=0;
        FILE *fp;
        int i=0;
        if((fp=fopen("venus.asc","r"))==NULL)
        {
                printf("no venus.asc file ");
                return;
        }
        fscanf(fp,"%d",&m);
        double **p=(double**)malloc(m*sizeof(double*));
        int j=0;
        for(j=0;j<m;j++)
        {
                p[j] = (double*)malloc(3*sizeof(double));
        }

        double arr[3];
        int n=0;
        while(!feof(fp))
        {
                int ret=fscanf(fp,"%lf %lf %lf",&arr[0],&arr[1],&arr[2]);
                if(ret!=3)
                        break;
                p[n][0]=arr[0];
                p[n][1]=arr[1];
                p[n][2]=arr[2];
                n++;
        }
        fclose(fp);
         for(j=0;j<10;j++)
        {
                printf("%lf %lf %lf\n",p[j][0],p[j][1],p[j][2]);
               
        }

        int node =0;
        for(node=m/2 -1;node>=0;node--)//生成大堆根
        {
                heapsort(p, node, m);
               
               
        }
        for(node=m-1;node>0;node--)
        {
                double tmp0=0,tmp1=0,tmp2=0;
                tmp0=p[node][0];//与最后一个记录交换
                tmp1=p[node][1];
                tmp2=p[node][2];
                p[node][0]=p[0][0];
                p[node][1]=p[0][1];
                p[node][2]=p[0][2];
                p[0][0]=tmp0;
                p[0][1]=tmp1;
                p[0][2]=tmp2;
                heapsort(p, 0, node);  //将H.r[0..i]重新调整为大根堆
                                   
        }
        fp=fopen("venus1.asc","w");
        int k=0;
        double int_arr[3];
        fprintf(fp,"%d\n",m);

        for( i=1;i<=m;i++)
        {
                int_arr[0]=p[i-1][0];
                int_arr[1]=p[i-1][1];
                int_arr[2]=p[i-1][2];
                int ret=fprintf(fp,"%lf %lf %lf\n", int_arr[0],int_arr[1],int_arr[2]);
                //if(ret!=3)
                       // break;
                n++;
        }
        fclose(fp);

}

void heapsort(double **data,int i,int length)
{
        int minnode=0;
        while((2*i+1)<length)
        {
                minnode=2*i+1;
                if((2*i+2)<length)
                { if(data[2*i+1][0]<data[2*i+2][0])//比较左子树和右子树，记录最大值的
                minnode=2*i+2;
                }
                if(data[i][0]<data[minnode][0])
                {
                        double t0=0,t1=0,t2=0;
                        t0=data[i][0];
                        t1=data[i][1];
                        t2=data[i][2];
                        data[i][0]=data[minnode][0];
                         data[i][1]=data[minnode][1];
                         data[i][2]=data[minnode][2];
                        data[minnode][0]=t0;
                        data[minnode][1]=t1;
                        data[minnode][2]=t2;
                        i=minnode;
               
               
                }
                 else
       
         {
            //比较左右孩子均大则堆未破坏，不再需要调整
            break;
         }


        }
        return;

}