并行计算—OpenMP—统计完数个数

<pre class="cpp" name="code">//parallel统计完数的并行算法和OpenMP并行编程。#include "stdafx.h"#include <omp.h>#include <iostream>#include <time.h>#include <math.h>#include <stdlib.h>int isPerfect(int n){    int lyy_sum=0;    for(int i = 1; i <=n/2; i++)        if(n%i == 0) lyy_sum+= i;//统计所有真因子的和。    if(lyy_sum == n) return 1;//如果与原值相等，则该数为完数。    return 0;//不是完数。}using namespace std;int _tmain(int argc, _TCHAR* argv[]){omp_set_num_threads(2);int lyy_n,lyy_num=0;clock_t t1, t2;long lyy_sumtmp[2];cin>>lyy_n;//并行t1=clock();    #pragma omp parallel{long i;long id=omp_get_thread_num();long lyy_num1=0;for(i = id; i < lyy_n; i=i+2)        if(isPerfect(i)) lyy_num1++;lyy_sumtmp[id]=lyy_num1;}for(long i=0;i<2;i++)lyy_num+=lyy_sumtmp[i];t2=clock();cout<<lyy_num-1<<endl;cout<<"parallel time="<<(t2-t1)<<endl;//串行lyy_num=0;t1=clock();    for(int i = 1; i < lyy_n; i++)        if(isPerfect(i)) lyy_num++;t2=clock();cout<<lyy_num<<endl;cout<<"serial time="<<(t2-t1)<<endl;system("pause");return 0;}

运行结果：

分析：

    加速比计算：8253/4134=1.996