双线程输出素数（双核CPU的利用率可达100%）

/* * Author: Leng_que * Date: 2009年11月14日 * E-mail: leng_que@yahoo.com.cn * Description: 一个很基础的多线程小程序，用于输出素数。 */ #include <math.h> #include <stdio.h> #include <windows.h> CRITICAL_SECTION g_cs={0}; //功能：判定一个正整数是否为素数 bool isPrime(long num) { if ( num<=1 ) { return false; } else if ( num==2 ) { return true; } else { long i = (long)sqrt(num)+1; for ( long n=2; n<i; n++ ) { if ( num%n == 0 ) { return false; } } } return true; } //输出[3,1000000]范围内的所有素数 #define N 1000000 //线程一：搜寻并输出3、7、11、15、19、……序列中的素数 DWORD WINAPI T1( LPVOID lpParameter ) { for ( int n=3; n<=N; n+=4 ) { if ( isPrime(n) ) { EnterCriticalSection(&g_cs); printf("%d ",n); LeaveCriticalSection(&g_cs); } } return 0; } //线程二：搜寻并输出5、9、13、17、……序列中的素数 DWORD WINAPI T2( LPVOID lpParameter ) { for ( int n=5; n<=N; n+=4 ) { if ( isPrime(n) ) { EnterCriticalSection(&g_cs); printf("%d ",n); LeaveCriticalSection(&g_cs); } } return 0; } #undef N //主函数 int main(void) { HANDLE handle[2]={0}; InitializeCriticalSection(&g_cs); handle[0] = CreateThread(NULL, 0, T1, NULL, 0, NULL); handle[1] = CreateThread(NULL, 0, T2, NULL, 0, NULL); WaitForMultipleObjects(2, handle, TRUE, INFINITE); DeleteCriticalSection(&g_cs); return 0; }

备注：这段代码的第一个版本被我删了，后来才传了这个更改后的版本，因为我后来发现原来那个有一些问题，主要是其中的算法在数学角度上的Bug，偶数除了2以外的绝对不可能为素数了，而之前那个版本里我写的算法就犯了这个愚蠢的错误。

 

小后记：我偶尔还是会看看自己博客以前的文章和代码的，随着时间的推移，我常常会发现以前的代码会有这样或那样的不够好的地方，甚至Bug，对于这些问题，我基本都会更正过来的，这也许也就是进步吧。