#include <stdio.h>
#include <stdlib.h>
#include <mpi.h>
 
#define SIZE 8 
 
int main(int argc, char *argv[]) {
    int rank, size;
    int numbers[SIZE][3];
    int lsum = 0, gsum = 0; //local_sum and global_sum
    int count_per_process;
    int *recvbuf; 
 
    // MPI environment
    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
 
    // Process 0 reads the input numbers
    if (rank == 0) {
        printf("Enter the numbers:\n"); //inputs 3 integers
        for (int i = 0; i < SIZE; i++) {
            scanf("%d %d %d", &numbers[i][0], &numbers[i][1], &numbers[i][2]);
        }
    }
 
    MPI_Bcast(numbers, SIZE * 3, MPI_INT, 0, MPI_COMM_WORLD);
 
    count_per_process = SIZE / size;
 
    // Each process calculates its local sum
    for (int i = rank * count_per_process; i < (rank + 1) * count_per_process; i++) {
        for (int j = 0; j < 3; j++) {
            lsum += numbers[i][j];
        }
        printf("Process %d partial sum of set %d: %d\n", rank, i, lsum);
    }
 
    // Reduce local sums to global sum at process 0
    MPI_Reduce(&lsum, &gsum, 1, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD);
 
    //prints the global sum
    if (rank == 0) {
        printf("Global sum: %d\n", gsum);
    }
 
    MPI_Finalize();
    return 0;
}
 

I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxzdGRsaWIuaD4KI2luY2x1ZGUgPG1waS5oPgoKI2RlZmluZSBTSVpFIDggCgppbnQgbWFpbihpbnQgYXJnYywgY2hhciAqYXJndltdKSB7CiAgICBpbnQgcmFuaywgc2l6ZTsKICAgIGludCBudW1iZXJzW1NJWkVdWzNdOwogICAgaW50IGxzdW0gPSAwLCBnc3VtID0gMDsgLy9sb2NhbF9zdW0gYW5kIGdsb2JhbF9zdW0KICAgIGludCBjb3VudF9wZXJfcHJvY2VzczsKICAgIGludCAqcmVjdmJ1ZjsgCgogICAgLy8gTVBJIGVudmlyb25tZW50CiAgICBNUElfSW5pdCgmYXJnYywgJmFyZ3YpOwogICAgTVBJX0NvbW1fcmFuayhNUElfQ09NTV9XT1JMRCwgJnJhbmspOwogICAgTVBJX0NvbW1fc2l6ZShNUElfQ09NTV9XT1JMRCwgJnNpemUpOwoKICAgIC8vIFByb2Nlc3MgMCByZWFkcyB0aGUgaW5wdXQgbnVtYmVycwogICAgaWYgKHJhbmsgPT0gMCkgewogICAgICAgIHByaW50ZigiRW50ZXIgdGhlIG51bWJlcnM6XG4iKTsgLy9pbnB1dHMgMyBpbnRlZ2VycwogICAgICAgIGZvciAoaW50IGkgPSAwOyBpIDwgU0laRTsgaSsrKSB7CiAgICAgICAgICAgIHNjYW5mKCIlZCAlZCAlZCIsICZudW1iZXJzW2ldWzBdLCAmbnVtYmVyc1tpXVsxXSwgJm51bWJlcnNbaV1bMl0pOwogICAgICAgIH0KICAgIH0KCiAgICBNUElfQmNhc3QobnVtYmVycywgU0laRSAqIDMsIE1QSV9JTlQsIDAsIE1QSV9DT01NX1dPUkxEKTsKCiAgICBjb3VudF9wZXJfcHJvY2VzcyA9IFNJWkUgLyBzaXplOwoKICAgIC8vIEVhY2ggcHJvY2VzcyBjYWxjdWxhdGVzIGl0cyBsb2NhbCBzdW0KICAgIGZvciAoaW50IGkgPSByYW5rICogY291bnRfcGVyX3Byb2Nlc3M7IGkgPCAocmFuayArIDEpICogY291bnRfcGVyX3Byb2Nlc3M7IGkrKykgewogICAgICAgIGZvciAoaW50IGogPSAwOyBqIDwgMzsgaisrKSB7CiAgICAgICAgICAgIGxzdW0gKz0gbnVtYmVyc1tpXVtqXTsKICAgICAgICB9CiAgICAgICAgcHJpbnRmKCJQcm9jZXNzICVkIHBhcnRpYWwgc3VtIG9mIHNldCAlZDogJWRcbiIsIHJhbmssIGksIGxzdW0pOwogICAgfQoKICAgIC8vIFJlZHVjZSBsb2NhbCBzdW1zIHRvIGdsb2JhbCBzdW0gYXQgcHJvY2VzcyAwCiAgICBNUElfUmVkdWNlKCZsc3VtLCAmZ3N1bSwgMSwgTVBJX0lOVCwgTVBJX1NVTSwgMCwgTVBJX0NPTU1fV09STEQpOwoKICAgIC8vcHJpbnRzIHRoZSBnbG9iYWwgc3VtCiAgICBpZiAocmFuayA9PSAwKSB7CiAgICAgICAgcHJpbnRmKCJHbG9iYWwgc3VtOiAlZFxuIiwgZ3N1bSk7CiAgICB9CgogICAgTVBJX0ZpbmFsaXplKCk7CiAgICByZXR1cm4gMDsKfQo=