实验二：证明卡布列克数6174的MPI程序设计

一、实验题目

用 MPI 技术证明卡布列克数6174猜想。

二、实验目的

熟悉 MPI 编程，加深对其编程的理解。

三、实验环境

Ubuntu22 下运行MPI。

使用实验一已经配置好的环境。
实验一：Linux下基于MPI的hello程序设计

四、原理介绍

6174猜想：一个任意的四位正整数（全相同的除外，如1111）。将数字重新组合成一个最大的数和最小的数相减，重复这个过程，最多七步，必得6174。

例如：
输入：2000
输出：
第1步：2000-0002=1998
第2步：9981-1899=8082
第3步：8820-0288=8532
第4步：8532-2358=6174
分步解题：

1. 遍历所有的四位数，判断。
2. 分解数字的千位，百位，十位，个位，放置在4个变量中，排序，将最大的放在最前面，最小的放在最后面。
3. 利用4个变量的值，构建四位数，最大值及最小值，并计算差值
4. 如差值不是6174，则重复步骤2.如果差值为6174，则程序打印输出最终结果，程序结束。
   用各自节点划分计算的数据范围，分担计算任务。
   例如，如果有4个节点机器，就将1000_{9999分成四个小范围，分给4个节点计算验证。如果有3个节点机器，就将1000}9999分成三个小范围，分给3个节点计算验证。

五、程序代码实现及测试

代码

#include <mpi.h>
#include <stdio.h>
#include <stdlib.h>

int compare(const void *a, const void *b) {
    int aa = *((int*) a);
    int bb = *((int*) b);
    if(aa == bb) return 0;
    else if(aa < bb) return -1;
    else return 1;
}

int kaprekar(int num) {
    int asc_num, desc_num;
    while (num != 6174) {
        int arr[4] = { num / 1000, (num / 100) % 10, (num / 10) % 10, num % 10 };
        qsort(arr, 4, sizeof(int), compare);
        asc_num = arr[0] * 1000 + arr[1] * 100 + arr[2] * 10 + arr[3];
        desc_num = arr[3] * 1000 + arr[2] * 100 + arr[1] * 10 + arr[0];
        num = desc_num - asc_num;
    }
    return num;
}

int main(int argc, char** argv) {
    int rank, size, name_len;
    char Processorname[20];
    int start, end;
    int count = 0;

    // 初始化MPI环境
    MPI_Init(&argc, &argv);

    // 获取总的处理器数量
    MPI_Comm_size(MPI_COMM_WORLD, &size);

    // 获取当前处理器的编号
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);

    // 获取当前处理器的名字
    MPI_Get_processor_name(Processorname, &name_len);

    // 计算每个节点分配的数据范围
    count = (9999 - 1000 + 1) / size; // 获取每个节点平均处理数量
    start = 1000 + rank * count; // 计算起始数值
    end = start + count - 1; // 计算结束数值

    if (rank == (size - 1)) {
        // 最后一个节点处理尾部剩余的数据范围
        end = 9999;
    }

    int num = start;

    while (num++ <= end) {
        if(num == 1111 || num % 1111 == 0) {
            continue;
        }
        int result = kaprekar(num);
        if (result == 6174) {
            printf("%d on %s\n",num,Processorname);
        }
    }



    MPI_Finalize();
    return 0;
}

运行结果：

mpicc test1.c -o test1
mpiexec -f /home/student/mpd_hosts ./test1