c语言实现k-means算法

// pree1.cpp : 定义控制台应用程序的入口点。
//


#include "stdafx.h"


#include <math.h>
#include <malloc.h>
#include <cstdlib>


#define DIM 784//定义数据的最大维数
#define NUM 10000//定义数据的最大数量
#define CLASS 10//定义最大聚类数


int initialize(double *a,int *sum,int *num,int i,int j);//初始化
void oldtonew(double *a,double *oa,int i,int j);//将新值赋给旧值
double  distance(int *a,double *b,int i);//求距离，a，b代表两个点，i代表几位数据。
int stable(double *a,double *oa,int i,int j);//判断是否稳定，a代表现在数据，oa代表之前数据，i代表数据量，j代表维数
int function(int *num,double *oa,double *a,int *sum,int b,int c,int d);
int minfind(double *l,int k);//从l中找到最小的,并返回其位置
void display(double *b,int *sum,int i,int j);


int *num = (int*)malloc(sizeof(int)*(NUM*DIM));//用来读取数据
double *a = (double *)malloc(sizeof(double)*(CLASS*DIM));//用来保存聚类中心
double *oa = (double *)malloc(sizeof(double)*(CLASS*DIM));//用来辅助进行更新
int *sum= (int  *) malloc (sizeof(int)*CLASS);//用来记录每个类的个数
char *o = (char *)malloc(sizeof(char)*(NUM*DIM));//用来将逗号吃掉




int main()
{
FILE *fp;
int i,j,k;
if((fp = fopen("ClusterSamples.csv","r"))==NULL)
    {
        printf("打开文件错误！\n");
    }
else
{
for(i = 0; i < NUM;i++)
        {
            for(j = 0;j < DIM;j++)
{
                fscanf(fp,"%d",&(*(num+i*DIM+j)));
                fscanf(fp,"%c",&(*(o+i*DIM+j)));
/* printf("%d  ",*(num+i*DIM+j));*/
}