数据挖掘原理与算法 Agnes算法

最新推荐文章于 2023-11-20 16:52:27 发布

Ayakanoinu

最新推荐文章于 2023-11-20 16:52:27 发布

阅读量919

点赞数 1

CC 4.0 BY-SA版权

分类专栏：数据挖掘原理与算法

本文链接：https://blog.csdn.net/weixin_43849505/article/details/100879782

数据挖掘原理与算法专栏收录该内容

4 篇文章

订阅专栏

本文深入探讨了一种基于暴力计算的聚类算法实现，详细解释了如何通过比较所有可能的节点组合来寻找最小距离，从而完成数据点的聚类过程。尽管此方法未进行时间优化，但在小数据集上能有效展示聚类效果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

暴力解决方法，没考虑优化时间，仅仅是暴力实现，数据量很大的话会超时。

#include<bits/stdc++.h>
using namespace std;
int n,k;
struct Node{
	int x,y;
	int id;
};
struct Node nodes[1005];
vector<struct Node> cu[1005];
double dis(int x,int y)
{
	return sqrt((nodes[x].x-nodes[y].x)*(nodes[x].x-nodes[y].x)+(nodes[x].y-nodes[y].y)*(nodes[x].y-nodes[y].y));
}
bool cmp(vector<struct Node> x,vector<struct Node> y)
{
	return x.size()<y.size();
}
bool cmp1(vector<struct Node> x,vector<struct Node> y)
{
	return x.size()>y.size();
}
bool cmp2(struct Node x,struct Node y)
{
	return x.id<
	y.id;
}
void Agnes()
{
	int N=n;
	while(N>k)
	{
		double mindis=9999;
		int tar1=-1,tar2=-1;
		for(int i=0;i<n;i++)
		{
			if(!cu[i].empty())
			for(int j=i+1;j<n;j++)
			{
				if(!cu[j].empty())
				for(int m=0;m<cu[i].size();m++)
				{
					for(int n=0;n<cu[j].size();n++)
					{
						double tempdis=dis(cu[i][m].id,cu[j][n].id);
						if(mindis>tempdis)
						{
							mindis=tempdis;
							tar1=i;
							tar2=j;
						}
					}
				}
			}
		}
		for(int i=0;i<cu[tar2].size();i++)
		{
			struct Node temp=cu[tar2][i];
			cu[tar1].push_back(temp);	
		}	
		cu[tar2].clear();
		sort(cu,cu+n,cmp);
		int cnt=0;
		for(int i=0;i<n;i++)
			if(!cu[i].empty())
				cnt++;
		N=cnt;
	}
	
}
int main()
{
	cin>>n>>k;
	for(int i=0;i<n;i++)
	{
		cin>>nodes[i].x>>nodes[i].y;
		nodes[i].id=i;
		cu[i].push_back(nodes[i]);
	}
	Agnes();
	sort(cu,cu+n,cmp1);
	for(int i=0;i<n;i++)
		if(!cu[i].empty())
		{
			sort(cu[i].begin(),cu[i].end(),cmp2);
			for(int j=0;j<cu[i].size();j++)
				if(j==0)
					cout<<"{"<<cu[i][j].id+1<<" ";
				else if(j==cu[i].size()-1)
					cout<<cu[i][j].id+1<<"}"<<endl;
				else
					cout<<cu[i][j].id+1<<" ";	
			if(cu[i].size()==1)
				cout<<"}"<<endl;
		}	
	return 0;
}