首先來看一個最簡單的C語言實現質因數分解的列子：
#include <stdio.h>void main( ){ int data, i = 2; scanf("%d", &data); while(data > 1) { if(data % i == 0) { printf("%d ", i); data /= i; } else i++; }}原理&&方法
把一個合數分解為若干個質因數的乘積的形式，即求質因數的過程叫做分解質因數，分解質因數只針對合數
求一個數分解質因數，要從最小的質數除起，一直除到結果為質數為止。分解質因數的算式的叫短除法，和除法的性質差不多，還可以用來求多個個數的公因式：
以24為例：
2 -- 24
2 -- 12
2 -- 6
3 （3是質數，結束）
得出 24 = 2 × 2 × 2 × 3 = 2^3 * 3
代碼
可先用素數篩選法，篩選出符合條件的質因數，然后for循環遍歷即可，通過一道題目來show一下這部分代碼
題目1
題目描述：
求正整數N(N>1)的質因數的個數。
相同的質因數需要重復計算。如120=2*2*2*3*5，共有5個質因數。
輸入：
可能有多組測試數據，每組測試數據的輸入是一個正整數N，(1<N<10^9)。
輸出：
對于每組數據，輸出N的質因數的個數。
樣例輸入：
120
樣例輸出：
5
提示：
注意：1不是N的質因數；若N為質數，N是N的質因數。
ac代碼
#include <stdio.h> int main() { int n, count, i; while (scanf("%d", &n) != EOF) { count = 0; for (i = 2; i * i <= n; i ++) { if(n % i == 0) { while (n % i == 0) { count ++; n /= i; } } } if (n > 1) { count ++; } printf("%d/n", count); } return 0; } 深入理解
我所謂的深入理解，就是通過4星的題目來靈活運用分解質因數的方法，題目如下
題目2
題目描述：
給定n，a求最大的k，使n！可以被a^k整除但不能被a^(k+1)整除。
輸入：
兩個整數n(2<=n<=1000)，a(2<=a<=1000)
輸出：
一個整數.
樣例輸入：
6 10
樣例輸出：
1
思路
a^k和n!都可能非常大，甚至超過long long int的表示范圍，所以也就不能直接用取余操作判斷它們之間是否存在整除關系，因此我們需要換一種思路，從分解質因數入手，假設兩個數a和b：
a = p1^e1 * p2^e2 * ... * pn^en, b = p1^d1 * p2^d2 * ... * pn^dn
, 則b除以a可以表示為：
b / a = (p1^d1 * p2^d2 * ... * pn^dn) / (p1^e1 * p2^e2 * ... * pn^en)
若b能被a整除，則 b / a必為整數，且兩個素數必護質，則我們可以得出如下規律：
若a存在質因數px，則b必也存在該質因數，且該素因數在b中對應的冪指數必不小于在a中的冪指數
另b = n!, a^k = p1^ke1 * p2^ke2 * ... * pn^ken，因此我們需要確定最大的非負整數k即可。要求得該k，我們只需要依次測試a中每一個素因數，確定b中該素因數是a中該素因數的冪指數的多少倍即可，所有倍數中最小的那個即為我們要求得的k
分析到這里，剩下的工作似乎只是對a和n!分解質因數，但是將n!計算出來再分解質因數，這樣n!數值太大。考慮n!中含有素因數p的個數，即確定素因數p對應的冪指數。我們知道n!包含了從1到n區間所有整數的乘積， 這些乘積中每一個p的倍數（包括其本身）都對n!貢獻至少一個p因子，且我們知道在1到n中p的倍數共有n/p個。同理，計算p^2，p^3,...即可
代碼
#include <stdio.h> #include <stdlib.h> #include <string.h> #define N 1001 int prime[N], size; /** * 素數篩選法進行預處理 */ void initProcess() { int i, j; for (prime[0] = prime[1] = 0, i = 2; i < N; i ++) { prime[i] = 1; } size = 0; for (i = 2; i < N; i ++) { if (prime[i]) { size ++; for (j = 2 * i; j < N; j += i) { prime[j] = 0; } } } } int main(void) { int i, n, a, k, num, count, base, tmp, *ansbase, *ansnum; // 預處理 initProcess(); while (scanf("%d %d", &n, &a) != EOF) { ansbase = (int *)calloc(size, sizeof(int)); ansnum = (int *)calloc(size, sizeof(int)); // 將a分解質因數 for (i = 2, num = 0; i < N && a != 1; i ++) { if (prime[i] && a % i == 0) { ansbase[num] = i; ansnum[num] = 0; while (a != 1 && a % i == 0) { ansnum[num] += 1; a = a / i; } num ++; } } // 求最小的k for (i = 0, k = 0x7fffffff; i < num; i ++) { base = ansbase[i]; count = 0; while (base <= n) { count += n / base; base *= ansbase[i]; } tmp = count / ansnum[i]; if (tmp < k) k = tmp; } printf("%d/n", k); } return 0; } /************************************************************** Problem: 1104 User: wangzhengyi Language: C Result: Accepted Time:0 ms Memory:916 kb ****************************************************************/ 約數個數定理
對于一個大于1的正整數n可以分解質因數：
n = p1^a1 * p2^a2 * p3^a3 * ... * pn^an
, 則n的正約數的個數為：
(a1 + 1) * (a2 + 1) * ... *(an + 1)
.其中p1,p2,..pn都是n的質因數，a1, a2...an是p1,p2,..pn的指數
證明
n可以分解質因數：n=p1^a1 * p2^a2 * p3^a3 * … * pk^ak,
由約數定義可知p1^a1的約數有:p1^0, p1^1, p1^2......p1^a1 ，共（a1+1）個;同理p2^a2的約數有（a2+1）個......pk^ak的約數有（ak+1）個
故根據乘法原理：n的約數的個數就是
(a1+1）*（a2+1）*（a3+1）*…* (ak+1)
題目3
題目描述：
輸入n個整數,依次輸出每個數的約數的個數
輸入：
輸入的第一行為N，即數組的個數(N<=1000)
接下來的1行包括N個整數，其中每個數的范圍為(1<=Num<=1000000000)
當N=0時輸入結束。
輸出：
可能有多組輸入數據，對于每組輸入數據，
輸出N行，其中每一行對應上面的一個數的約數的個數。
樣例輸入：
5
1 3 4 6 12
樣例輸出：
1
2
3
4
6
代碼
#include <stdio.h> #include <stdlib.h> #define N 40000 typedef long long int lint; int prime[N], size; void init() { int i, j; for (prime[0] = prime[1] = 0, i = 2; i < N; i ++) { prime[i] = 1; } size = 0; for (i = 2; i < N; i ++) { if (prime[i]) { size ++; for (j = 2 * i; j < N; j += i) prime[j] = 0; } } } lint numPrime(int n) { int i, num, *ansnum, *ansprime; lint count; ansnum = (int *)malloc(sizeof(int) * (size + 1)); ansprime = (int *)malloc(sizeof(int) * (size + 1)); for (i = 2, num = 0; i < N && n != 1; i ++) { if (prime[i] && n % i == 0) { ansprime[num] = i; ansnum[num] = 0; while (n != 1 && n % i == 0) { ansnum[num] += 1; n /= i; } num ++; } } if (n != 1) { ansprime[num] = n; ansnum[num] = 1; num ++; } for (i = 0, count = 1; i < num; i ++) { count *= (ansnum[i] + 1); } free(ansnum); free(ansprime); return count; } int main(void) { int i, n, *arr; lint count; init(); while (scanf("%d", &n) != EOF && n != 0) { arr = (int *)malloc(sizeof(int) * n); for (i = 0; i < n; i ++) { scanf("%d", arr + i); } for (i = 0; i < n; i ++) { count = numPrime(arr[i]); printf("%lld/n", count); } free(arr); } return 0; } /************************************************************** Problem: 1087 User: wangzhengyi Language: C Result: Accepted Time:190 ms Memory:1068 kb ****************************************************************/ 新聞熱點
疑難解答
圖片精選
