第５話 素数探索マルチスレッド解答例
コード例
#include<iostream>
#include<math.h>
using namespace std;
using namespace System;
using namespace System::Threading;
const int N=10000000;
int cn1,cn2,cn3,cn4;
void f1();
void f2();
void f3();
void f4();
char h(int a);
int main(){
DateTime^ hj=DateTime::Now;
cn1=1;
Thread^ a=gcnew Thread(gcnew ThreadStart(f1));
Thread^ b=gcnew Thread(gcnew ThreadStart(f2));
Thread^ c=gcnew Thread(gcnew ThreadStart(f3));
Thread^ d=gcnew Thread(gcnew ThreadStart(f4));
a->Start();
b->Start();
c->Start();
d->Start();
a->Join();
b->Join();
c->Join();
d->Join();
DateTime^ ow=DateTime::Now;
TimeSpan sa=ow->Subtract(*hj);
cout<<"素数の総数は"<<cn1+cn2+cn3+cn4<<"です。"<<endl;
cout<<"計算時間は"<<sa.TotalSeconds<<"秒です。"<<endl;
}
void f1(){
int i;
//cout<<"2"<<" ";
for(i=3;i<=N;i+=8){
if(h(i)){
//cout<<i<<" ";
cn1++;
}
}
//cout<<endl;
}
void f2(){
int i;
//cout<<"2"<<" ";
for(i=5;i<=N;i+=8){
if(h(i)){
//cout<<i<<" ";
cn2++;
}
}
//cout<<endl;
}
void f3(){
int i;
//cout<<"2"<<" ";
for(i=7;i<=N;i+=8){
if(h(i)){
//cout<<i<<" ";
cn3++;
}
}
//cout<<endl;
}
void f4(){
int i;
//cout<<"2"<<" ";
for(i=9;i<=N;i+=8){
if(h(i)){
//cout<<i<<" ";
cn4++;
}
}
//cout<<endl;
}
char h(int a){
int i;
for(i=3;i<=(int)sqrt((double)a);i+=2){
if(a%i==0)return(0);
}
return(1);
}
シングルスレッド版
に比べると１２／２８です。
２倍以上速くなりました。
それでは皆さん、第２話の出来の悪かった普遍版魔方陣自動生成マルチスレッド版を.NETによって改良しましょう。
第４話へ 第６話へ
