fork download
copy 
  1. #include <iostream>  
  2. #include <vector>  
  3. #include <cmath>  
  4. #include <mpi.h>  
  5.  
  6. // Function to check if a number is prime  
  7. bool isPrime(int num) {
  8. 	if (num <= 1) return false;
  9. 	for (int i = 2; i <= sqrt(num); i++) {
  10. 		if (num % i == 0) return false;
  11. 	}
  12. 	return true;
  13. }
  14.  
  15. // Function to find the largest gap between consecutive prime numbers  
  16. int findLargestGap(int start, int end) {
  17. 	int largestGap = 0;
  18. 	int prevPrime = -1;
  19. 	for (int i = start; i <= end; i++) {
  20. 		if (isPrime(i)) {
  21. 			if (prevPrime != -1) {
  22. 				largestGap = std::max(largestGap, i - prevPrime);
  23. 			}
  24. 			prevPrime = i;
  25. 		}
  26. 	}
  27. 	return largestGap;
  28. }
  29.  
  30. int main(int argc, char** argv) {
  31. 	MPI_Init(&argc, &argv);
  32.  
  33. 	int rank, size;
  34. 	MPI_Comm_rank(MPI_COMM_WORLD, &rank);
  35. 	MPI_Comm_size(MPI_COMM_WORLD, &size);
  36.  
  37. 	int n;
  38. 	if (rank == 0) {
  39. 		n = 100000;
  40. 	}
  41.  
  42. 	MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);
  43.  
  44. 	int chunkSize = n / size;
  45. 	int start = rank * chunkSize + 2; // Start from 2, the first prime number  
  46. 	int end = (rank == size - 1) ? n : (start + chunkSize - 1);
  47.  
  48. 	double startTime = MPI_Wtime();
  49. 	int largestGap = findLargestGap(start, end);
  50. 	double endTime = MPI_Wtime();
  51.  
  52. 	int globalLargestGap;
  53. 	MPI_Reduce(&largestGap, &globalLargestGap, 1, MPI_INT, MPI_MAX, 0, MPI_COMM_WORLD);
  54.  
  55. 	double timeTaken = endTime - startTime;
  56. 	double totalTimeTaken;
  57. 	MPI_Reduce(&timeTaken, &totalTimeTaken, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_WORLD);
  58.  
  59. 	if (rank == 0) {
  60. 		std::cout << "Largest gap between consecutive prime numbers less than " << n << ": " << globalLargestGap << std::endl;
  61. 		std::cout << "Time taken: " << totalTimeTaken << " seconds" << std::endl;
  62. 	}
  63.  
  64. 	MPI_Finalize();
  65. 	return 0;
  66. }
  67.  
Success #stdin #stdout #stderr 0.27s 40780KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
Standard output is empty
stderr
copy 
Error: unexpected '/' in "/"
Execution halted
https://ideone.com/AZq3X3
language:
R (R 3.5.2)
created:
16 hours ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!