永成科技C++笔试题
最后幾個題有點難度,在這里說一下:
永成科技C++筆試題
2013-11-19
1.將1億以內的質數存到一個超級大的數組中,用算法如何實現?
使用"篩法"求解1億以內的質數的程序的思路:
先動態分配1億個bit(總計12500000字節),用字節中的每一位代表每一個整數,首先將代表奇數的那些bit位置1,也就是代表偶數(合數)的位,接著再進一步從這些奇數位中篩掉合數.
篩掉合數的方法是,先從100000000(1億)的開方10000范圍內的質數i(3,5,7,11,13,17,19,23,29)開始,去找它在1億內的奇數倍數i*i,i*i+2i,i*i+4i,...,這里
沒有i*i+i是因為它可以寫成(i+1)i是2的倍數,已經被過濾掉,將代表這些合數的bit位置0,
那么最后剩下的bit為1的那些bit,就是代表質數的,統計出它們的個數就可以了.
這樣做的原理是,基于如下的事實:
(1)任意連續的6個數中,就只會測試2個而已,以6n, 6n+1, 6n+2, 6n+3, 6n+4, 6n+5為例,只需測試6n+1和6n+5, 工作量減少到1/3
(2)判斷一個數i是否是素數的話,只需要測試2->sqrt(i)之間的質數就可以了
理由如下:
按素數的定義,也就是只有 1 與本身可以整除,所以可以用 2→i-1 去除 i,如果都除不盡,i 就是素數。觀點對,但卻與上一點一樣的笨拙。當 i>2 時,有哪一個數可以被 i-1 除盡的?沒有!為什么?如果 i 不是質數,那么 i=a×b,此地 a 與 b 既不是 i 又不是 1;正因為 a>1,a 至少為 2,因此 b 最多也是 i/2 而已,去除 i 的數用不著是 2→i-1,而用 2→i/2 就可以了。不但如此,因為 i=a×b,a 與 b 不能大于 sqrt(i),為什么呢?如果 a>sqrt(i),b>sqrt(i),于是 a×b > sqrt(i)*sqrt(i) = i,因此就都不能整除i了。如果i不是質數,它的因子最大就是 sqrt(i);換言之,用 2→sqrt(i)去檢驗就行了
但是,用 2→sqrt(i) 去檢驗也是浪費。就像前面一樣,2 除不盡,2 的倍數也除不盡;同理,3 除不盡,3 的倍數也除不盡……最理想的方法就是用質數去除i。
如果只檢查 6n+1 和 6n+5 ?不難發現,它們的距離是4、2、4、2……所以,可以先定義一個變量 gab=4,然后 gab=6-gab;
(3)不能用開方而應該用平方
比較是不能用 sqrt(i),因為它不精確。舉個例子,i=121,在數學上,sqrt(i) 自然是 11,但計算機里的結果可能是 10.9999999,于是去除的數就是 2、3、5、7,而不含 11,因此 121 就變成質數了。解決這個問題的方法很簡單,不要用開方,用平方即可。例如,如果 p*p<=i,則就用 p 去除 i。而且它的效率比開方高。
為此需要先將2,3,5,7,11,13這樣的質數先定位到32bit長度的整數內,這個整數的四字節中的每個字節是10101010,將這個質數放到32bit中的唯一一個bit位上面.
最后計算的結果是5761455個素數,而且程序用時9.375秒
下面是一個老外提供的實現代碼,但是我們還有比這個更高效的處理:
//Platform: Ubuntu 12.04.3 64bit
//gcc -std=c99 -g primer_demo.c -o primer_demo
#include <stdio.h>
#include <memory.h>
#include <stdlib.h>
#include <time.h>int count(unsigned int a) //統計每個整數中的非0位個數,也就是素數的個數(素數沒被篩掉,相對應位為1)
{int sum = 0;for (unsigned int x=a; x; x>>=1) //x不斷作右移運算{if (x & 1) //x與1作與運算sum++;}return sum;
}void sieve(unsigned int* p) //篩選法求1億以內素數
{
// for(int i=2;i<=10000;i++)
// for (int i=3; i<=10000; i+=2) //只篩選奇數顯然快于原算法for(int istep=4,i=3; i < 10000 ;i+=(istep^0x6)) //進一步優化,%6余1和5時才可能是素數,即只檢查交替相隔2和4的數{if (p[i/32] & (1<<i%32)) //第i個整數對應的32位整數序位為i/32,對應無符號整數的第i%32位,該語句判斷p[i]是否為1{ //上面的1<<i%32中,%的優先級高,所以先算出i在無符號整數中在哪一位,然后將1左移這么多位
// for (int j=i*i; j<100000000; j+=i)for(int j=i*i ; j<100000000; j+=i*2) //改進版,每次跳躍檢查,j+i已經被2篩除了,每次應檢查j+2*i{p[j/32] &= ~(1<<j%32); //第j個整數置為合數,即是將相應位置0}}}
}int main()
{clock_t start = clock();unsigned int* p = (unsigned int*) malloc(12500000); //向堆空間申請一億位(12500000字節*8)if (!p){printf("Not enough memory.\n");return 1;}
// memset(p,255,12500000); //將一億位全置為1(無符號整數為255)memset(p,170,12500000); //將一億位中的奇數位索引全置為1(索引從0開始,無符號整數為170)sieve(p); //對應非素數位全清0int num = 0; // 整數1對應位置為1,整數2對應位置為0,正好計數的時候抵消了,所以NUM初始化為0。for (int i=0; i<12500000/4; i++){num += count(p[i]);}free(p);printf("within 100,000,000 primers total counter:%d, total time was %7.3f in second\n",num,(double)(clock()-start)/CLOCKS_PER_SEC);return 0;
}
2.二叉樹求值?
要求使用遞歸和循環這兩種方法實現.
3.設計簡單的文本編輯器,寫出架構設計及思路?
總結
以上是生活随笔為你收集整理的永成科技C++笔试题的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 中体骏彩C++面试题
- 下一篇: webkit入门准备