c_jun.cpp

c++开发的c均值算法

#include<iostream.h>
#include<math.h>
#include<stdlib.h>
typedef struct {
    float x,y;
    }element;
typedef struct{
    element e;
    int type;
    }dat;

void Type_request(element w[],dat &d,int c){
    int k,number=-1;
    float index=65535,ii;
    for(k=0;k<c;k++){
        ii=sqrt((w[k].x-d.e.x)*(w[k].x-d.e.x)+(w[k].y-d.e.y)*(w[k].y-d.e.y));
        if(ii<index){index=ii;   number=k;}
        }//for
    d.type=number;
    }//type_re

int compare(element w[],element result[],int c){
     for(int t=0;t<c;t++)
        if(w[t].x!=result[t].x||w[t].y!=result[t].y)return 0;
     return 1;//equal
    }//compare
    

void pick(dat d[],element result[],int c,int n){
    element w[c],sum;
    int i,j,num,times=1;
    for(i=0;i<c;i++){ result[i]=w[i]=d[i].e;  d[i].type=i;}
    for(i=c;i<n;i++){ Type_request(w,d[i],c);}

    for(int j=0;j<n;j++)
        cout<<d[j].e.x<<"---"<<d[j].e.y<<" --->"<<d[j].type<<endl;
    cout<<endl;
    cout<<"the "<<times<<"th time:"<<endl;
    for(j=0;j<c;j++)cout<<"Type "<<j<<":----------x="<<w[j].x<<"  Y="<<w[j].y<<endl;

  do{  cout<<"the "<<++times<<"th time:"<<endl;
    for(j=0;j<c;j++){
        num=0;sum.x=0;sum.y=0;
        for(i=0;i<n;i++)
            if(j==d[i].type){ sum.x+=d[i].e.x;  sum.y+=d[i].e.y; num++;  }
        w[j].x=sum.x/num; w[j].y=sum.y/num;
        cout<<"Type "<<j<<":----------x="<<w[j].x<<"  Y="<<w[j].y<<endl;
       }// for(j..
    if(compare(w,result,c))break;
    for(i=0;i<c;i++){result[i].x=w[i].x;  result[i].y=w[i].y;}
    for(i=0;i<n;i++){Type_request(w,d[i],c);}
   }while(1);
   
   cout<<"The result is:"<<endl;  //print the relult
    for(int m=0;m<n;m++)
   cout<<d[m].e.x<<"---"<<d[m].e.y<<"--->"<<d[m].type<<endl;
}//pick


void main()
{   int c=3;
    int n=6;
    element result[2];
    dat d[6]={0};


    d[0].e.x=1;d[0].e.y=1;
    d[1].e.x=3;d[1].e.y=4.5;
    d[2].e.x=7;d[2].e.y=7;
    d[3].e.x=7;d[3].e.y=6;
    d[4].e.x=3;d[4].e.y=3.2;
    d[5].e.x=8;d[5].e.y=7.5;
    cout<<"the date program is:"<<endl;
    for(int j=0;j<n;j++)
        cout<<d[j].e.x<<"---"<<d[j].e.y<<endl;
    cout<<"(((((((((((((((((((((((((((((((("<<endl;
    pick(d,result,c,n);
    }//main