2021-04-19 寻找第k小元素

1.问题
用特定的分治策略，选第k小的元素。
设 L是n个元素的集合，从 L 中选取第k小的元素，其中 1<=k<=n.这里的第k小元素是指，当 L 按从小到大排好序之后，排在第k个位置的元素。

2.解析
以S中的某个元素m作为划分标准，将S划分为两个数组S1和S2，把这个数组中<m的数字都放入S1数组中，把这个数组中>m*的都放入S2数组中。

若k<|S1|，则原问题归纳为在数组S1中找第k小的子问题。
若k=|S1|+1，则m*就是我们所要找的第k小元素。
若k>|S1|+1，则归纳为S2中找第n−|S1|−1小的子问题。

3.设计
输入: 数组S[n],正整数k
输出: S[n]中的第k小元素

1、将S划分成5个一组,共「n/5」个组（向下取整）
2、将q组中的每一组单独排序，找出中项。把这些中项放到集合M中
3、m*←Select(M,[ |M|/2])   //M为中项集合的中项
4、把A和D中的每个元素与m*比较,小的构成S1，大的构成S2;
5、S1←S1UC; S2←S2∪B;
6、

if (k=|S1|+l)
Printf m*
else if (k≤|S1| )
Select(S1, k)
else
Select(S2, k-|S1l- 1)

int Select(int* a, int Length, int k)
{
	int n = Length / 5;
	int M[MAX], Nub = 0;
	int S1[MAX], S2[MAX];
	int Nub1 = 0, Nub2 = 0;

	if (Length <= 5)
	{
		qsort(a, Length, sizeof(int), compfun);
		return a[k - 1];
	}
	for (int i = 0; i < n; i++)
	{

		qsort(a + 5 * i, 5, sizeof(int), compfun);
		M[Nub++] = a[5 * i + 2];
	}

	int M2 = Select(M, Nub, Nub / 2 + 1);

	for (int i = 0; i < n; i++)
	{
		if (a[5 * i + 2] < M2)
		{
			S1[Nub1++] = a[5 * i];
			S1[Nub1++] = a[5 * i + 1];
			S1[Nub1++] = a[5 * i + 2];

			if (a[5 * i + 3] < M2)
				S1[Nub1++] = a[5 * i + 3];
			else
				S2[Nub2++] = a[5 * i + 3];

			if (a[5 * i + 4] < M2)
				S1[Nub1++] = a[5 * i + 4];
			else
				S2[Nub2++] = a[5 * i + 4];
		}
		else if (a[5 * i + 2] == M2)
		{
			S1[Nub1++] = a[5 * i];
			S1[Nub1++] = a[5 * i + 1];
			S2[Nub2++] = a[5 * i + 3];
			S2[Nub2++] = a[5 * i + 4];
		}
		else
		{
			S2[Nub2++] = a[5 * i + 3];
			S2[Nub2++] = a[5 * i + 4];
			S2[Nub2++] = a[5 * i + 2];

			if (a[5 * i] < M2)
				S1[Nub1++] = a[5 * i];
			else
				S2[Nub2++] = a[5 * i];

			if (a[5 * i + 1] < M2)
				S1[Nub1++] = a[5 * i + 1];
			else
				S2[Nub2++] = a[5 * i + 1];
		}
	}

	for (int i = 5 * n; i < Length; i++)
	{
		if (a[i] < M2)
			S1[Nub1++] = a[i];
		else
			S2[Nub2++] = a[i];
	}

	if (k == Nub1 + 1)
		return M2;
	else if (k <= Nub1)
		return Select(S1, Nub1, k);
	else
		return Select(S2, Nub2, k - Nub1 - 1);
}

[算法复杂度推导]
O（n）；
https://github.com/zzx114514/-.git