C++中约数定理的实例详解
对于一个大于1正整数n可以分解质因数：n = p1^a1*p2^a2*......pk^ak，则n的正约数的个数就是 ：(a1+1)*(a2+1)*......*(ak+1)
其中a1、a2、a3…ak是p1、p2、p3，…pk的指数。
用这个定理求一个数的约数个数是非常快的，贴出一道训练题目：
hdu 1492 -求约数的个数
贴出代码：
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
//约数定理的 #include <iostream> #include <algorithm> #include <iterator> #include <cstdio> #include <cstdlib> #include <cmath> #include <cstring> #include <vector> #include <queue> #include <set> using namespace std; #define ll long long int main() { // freopen("s.cpp","r",stdin); ll n; while(scanf("%lld",&n) != EOF) { if(!n) break; ll sum = 1; /* x = p1^a1*p2^a2*p3^a3...pk^ak yueshu = (a1+1)*(a2+1)*...*(ak+1)*/ for(ll i = 2; i*i <= n; i++){ int cou = 0; if(n%i==0){ cou = 1; n /= i; while(n%i==0){ cou++; n /= i; } } if(cou != 0){ sum = sum*(cou+1); } } if(n != 1){ sum = sum*2; } if(sum==1 && n==1){ sum = 1; } printf("%lld\n",sum); } return 0; } |
感谢阅读，希望能帮助到大家，谢谢大家对本站的支持！
原文链接：http://blog.csdn.net/ruangongshi/article/details/72430504
