1019. 数字黑洞

给定任一个各位数字不完全相同的4位正整数，如果我们先把4个数字按非递增排序，再按非递减排序，然后用第1个数字减第2个数字，将得到一个新的数字。一直重复这样做，我们很快会停在有“数字黑洞”之称的6174，这个神奇的数字也叫Kaprekar常数。

例如，我们从6767开始，将得到

7766 - 6677 = 1089
9810 - 0189 = 9621
9621 - 1269 = 8352
8532 - 2358 = 6174
7641 - 1467 = 6174
... ...

现给定任意4位正整数，请编写程序演示到达黑洞的过程。

输入格式：

输入给出一个(0, 10000)区间内的正整数N。

输出格式：

如果N的4位数字全相等，则在一行内输出“N - N = 0000”；否则将计算的每一步在一行内输出，直到6174作为差出现，输出格式见样例。注意每个数字按4位数格式输出。

输入样例1：

6767

输出样例1：

7766 - 6677 = 1089
9810 - 0189 = 9621
9621 - 1269 = 8352
8532 - 2358 = 6174

输入样例2：

2222

输出样例2：

2222 - 2222 = 0000

算法：

1. 算法开始。
2. 声明整型变量n，整型变量数组a，整型变量i，整型变量larger、smaller和tmp。
3. 初始化i为零，larger为零，smaller为零。
4. 读取整数，存储于n。
5. tmp等于n；
6. 如果larger-smaller等于6174，则跳到十六句。
7. 如果i大于等于四，则执行第十一句。
8. 将tmp除以十的余数存储于a[i]。
9. tmp等于tmp除以十。
10. i等于i加一。回到第七句。
11. 如果a的元素都相等，则输出n - n = 0000，跳到第十六句。
12. 对a进行非增序排列，larger等于a[0]*1000+a[1]*100+a[2]*10+a[3]。
13. 对a进行非降序排列，smaller等于a[0]*1000+a[1]*100+a[2]*10+a[3]。
14. tmp等于larger减去smaller。
15. 输出larger - smaller = tmp，返回第六句。
16. 算法结束。

注意：题目未要求末次输出是否带换行符号，则是否带换行符号对正确性没有影响。

下列代码可顺利运行于Xcode和PAT。

#include <stdio.h>
#include <stdlib.h>
#define N 4
int cmp_descent(const void* v1, const void* v2);
int cmp_increase(const void* v1, const void* v2);
int main(int argc, const char * argv[]) {
    int n, i, a[N], larger=0, smaller=0, tmp;
    scanf("%d", &n);
    tmp=n;
    while(larger-smaller!=6174){
        for(i=0; i<N; i++){
            a[i]=tmp%10;
            tmp/=10;
        }
        if(a[0]==a[1]&&a[0]==a[2]&&a[0]==a[3]){
            printf("%04d - %04d = 0000", n, n);
            break;
        }
        else{
            qsort(a, N, sizeof(int), cmp_descent);
            larger=a[0]*1000+a[1]*100+a[2]*10+a[3];
            qsort(a, N, sizeof(int), cmp_increase);
            smaller=a[0]*1000+a[1]*100+a[2]*10+a[3];
            tmp=larger-smaller;
            printf("%04d - %04d = %04d\n", larger, smaller, tmp);
        }
    }
    return 0;
}
int cmp_descent(const void* v1, const void* v2){
    if(*((int*)v1)>=*((int*)v2)){
        return 0;
    }
    else{
        return 1;
    }
}
int cmp_increase(const void* v1, const void* v2){
    if(*((int*)v1)<=*((int*)v2)){
        return 0;
    }
    else{
        return 1;
    }
}