/****************************************************************************
*???????????????????????????????????????????????????????????????????????????*
*??CURE?????????????????????????????????????????????????????????????????????*
*???????????????????????????????????????????????????????????????????????????*
*****************************************************************************/

#include?
#include?
#include?
#include?
#include?
#include?

#define?????????SUCCESS?????????1
#define?????????FAILURE?????????0
#define?????????TRUE????????????1
#define?????????FALSE???????????0
#define?????????NUMPATTERNS?????10???//?數據點的個數
#define?????????SIZEVECTOR??????2?????//?數據集的維數
#define?????????NUMCLUSTERS?????2?????//?類的個數
#define?????????NUMPRE??????????3????//?代表點個數
#define?????????SHRINK??????????0.5???//?收縮率

//?*****?定義結構和類?*****
struct?ClustNode?{
	int??????????Member[NUMPATTERNS];??//類中的數據項
	int??????????NumMembers;
	double???????Means[SIZEVECTOR];??//均值點
	double???????Pre[NUMPRE][SIZEVECTOR];//代表點
	int??????????PreMember[NUMPRE];//作為代表點的數據項
	int??????????closet;??//最近的簇
	double???????MinDist;//與最近的簇之間的距離
};

class?System?{
private:
	double??????Pattern[NUMPATTERNS][SIZEVECTOR];
	ClustNode?*p[NUMPATTERNS];
????void????????BuildClustList（）;?//建立初始簇鏈表
	ClustNode?*?Merge（ClustNode?*ClustNode?*）;?//合并兩個簇
	int?????????MinClust（）;????//找到最近的兩個簇
	double??????dist（ClustNode*ClustNode?*）;//兩個簇之間的最短距離
	double??????EucNorm（double?*double?*）;//代表點之間的距離
	void????????ShrinkPre（ClustNode?*）;??//收縮代表點

public:
	System（）;
	~System（）;
	int?LoadPatterns（char?*fname）;??????//?處理數據集
	void?RunCure（）;?????????????????????//?聚類的過程
	void?ShowClusters（）;????????????????//?顯示聚類結果
};
System::System（）
{
}
System::~System（）
{
	int?i;
	for（i=0;i		if（p[i]==NULL）
			continue;
		else
			delete?p[i];
	}
}
int?System::LoadPatterns（char?*fname）{
	int????ijlen;
	char?buff[300];
	char?*ptr*ptr1;
	ifstream?infile（fnameios::in）;????//將fname打開作為輸入
	if（!infile）{
		cout<<“不能打開輸入文件:\n“;
		return?FAILURE;
	}
	for（i=0;?i		infile.getline（buff300）;?????//將300個字符傳入buff中，buff是指針
		len=0;
		ptr=buff;
		for?（j=0;?j			ptr1=strchr（ptr‘‘）;???//找出sptr指向的字符串中第一次出現““的位置
			*ptr1=‘\0‘;
			Pattern[i][j]=atof（ptr）;????//將ptr指向的字符串轉化為浮點數
			len=strlen（ptr）;
			ptr=ptr+len+1;
		}?/*?endfor?*/
	}?/*?endfor?*/
	infile.close（）;
	return?SUCCESS;
}
void??System::BuildClustList（）{
	int?ij;
	double?MinDistd;
	for（i=0;i		p[i]=new?ClustNode;
		p[i]->NumMembers=1;
		p[i]->Member[0]=i;
		p[i]->PreMember[0]=i;
		for（j=0;j			p[i]->Means[j]=Pattern[i][j];//均值
		}
		for（j=0;j			p[i]->Pre[0][j]=Pattern[i][j]+（SHRINK）*（p[i]->Means[j]-Pattern[i][j]）;//計算代表點
		}
	}
	for（i=0;i		p[i]->MinDist=9.9e+99;
		for（j=0;j