基数排序

题目：

https://acm.zzuli.edu.cn/zzuliacm/problem.php?id=1124

题意：

Description

已知数组a中有m个按升序序排列的元素，数组b中有n个降序排列的元素，编程将a与b中的所有元素按降序存入数组c中。

Input

输入有两行，第一行首先是一个正整数m，然后是m个整数；第二行首先是一个正整数n，然后是n个整数，m, n均小于等于1000000。

Output

输出合并后的m+n个整数，数据之间用空格隔开。输出占一行。

思路：

基数排序实现，分为lsd和msd实现两种，即低位优先和高位优先。基数排序需要2倍的空间，我感觉速度很可以啊。。。

#include <bits/stdc++.h>

using namespace std;

const int N = 2000000 + 10;
const int radix = 10; //基数，即几进制
int a[N], b[N];
//四个参数：a是待排序数组，排序范围a[st,en]，len是待排序范围中最大数字的位数
void lsd_radix_sort(int a[], int st, int en, int len)
{//从低位到高位依次按每位上的数字进行计数排序
    int num[radix+1];
    int p = 1;
    while(len--)
    {
        memset(num, 0, sizeof num);
        for(int i = st; i <= en; i++) num[a[i]/p%10]++;
        for(int i = 1; i <= radix-1; i++) num[i] += num[i-1];
        for(int i = en; i >= st; i--) b[num[a[i]/p%10]--] = a[i];
        for(int i = st, j = 1; i <= en; i++, j++) a[i] = b[j];
        p *= 10;
    }
}
int get_pow(int x, int n) //求x^n，只有msd_radix_sort用到
{
    int ans = 1;
    for(int i = 1; i <= n; i++) ans *= x;
    return ans;
}
void msd_radix_sort(int a[], int st, int en, int len)
{//从高位到低位依次按每位上的数字进行计数排序后按数字分桶，然后递归
    if(len == 0 || st > en) return;
    int p = get_pow(10, len-1);
    int num[radix+1];
    memset(num, 0, sizeof num);
    for(int i = st; i <= en; i++) num[a[i]/p%10]++;
    for(int i = 1; i <= radix; i++) num[i] += num[i-1];//注意这里循环到radix
    for(int i = en; i >= st; i--) b[num[a[i]/p%10]--] = a[i];
    for(int i = st, j = 1; i <= en; i++, j++) a[i] = b[j];
    for(int i = 0; i <= radix-1; i++)//此时num[i]意为第i个桶的左边界
        msd_radix_sort(a, st + num[i], st + num[i+1] - 1, len-1);
}
int main()
{
    int n, m;
    scanf("%d", &n);
    for(int i = 1; i <= n; i++) scanf("%d", &a[i]);
    scanf("%d", &m);
    for(int i = 1; i <= m; i++) scanf("%d", &a[n+i]);
    //sort(a + 1, a + 1 + n+m);
    lsd_radix_sort(a, 1, n+m, 7);//本题数字最长为7
    //msd_radix_sort(a, 1, n+m, 7);
    for(int i = n+m; i >= 1; i--) printf("%d%c", a[i], i == 1 ? '\n' : ' ');
    return 0;
}