珠排序算法

#include <iostream>
#include <cstring>
using namespace std; 
  #define BEAD(i, j) beads[i * max + j] 
void beadSort(int *a, int len) 
{ 
    // 找到最大元素
    int max = a[0]; 
    for (int i = 1; i < len; i++) 
        if (a[i] > max) 
           max = a[i]; 
  
    // 分配内存 
    unsigned char beads[max*len]; 
    memset(beads, 0, sizeof(beads)); 
  
    //标记珠子
    for (int i = 0; i < len; i++) 
        for (int j = 0; j < a[i]; j++) 
            BEAD(i, j) = 1; 
  
    for (int j = 0; j < max; j++) 
    { 
        //数一数每列有多少颗珠子
        int sum = 0; 
        for (int i = 0; i < len; i++) 
        { 
            sum += BEAD(i, j); 
            BEAD(i, j) = 0; 
        } 
  
        //向上移动珠子
        for (int i = len - sum; i < len; i++) 
            BEAD(i, j) = 1; 
    } 
  
    // 使用珠子将排序后的值放入数组中
    for (int i = 0; i < len; i++) 
    { 
        int j; 
        for (j = 0; j < max && BEAD(i, j); j++); 
  
        a[i] = j; 
    } 
} 


// 主函数
int main() 
{ 
	int a[6];
	FILE* fpread;
	fpread = fopen("data.txt", "r");
	if (fpread == NULL)
	{
		printf("file is error.");
		return -1;
	}

	
		for (int j = 0; j < 6; j++)
		{
			fscanf(fpread, "%d", &a[j]);
		}
	
	fclose(fpread);
    beadSort(a, 6); 
     for (int i = 0; i < 6; i++)
        cout << a[i] << "\n";
  
    


    return 0; 
}