试除法求约数:
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;
//时间复杂度:O(n * sqrt(a))
vector<int> get_divisors(int n)
{
vector<int> res;
for(int i = 1; i <= n / i; ++i)
if(n % i == 0)
{
res.push_back(i);
if(i != n / i)
res.push_back(n / i);
}
sort(res.begin(), res.end());
return res;
}
int main()
{
int n;
scanf("%d", &n);
while(n--)
{
int x;
cin >> x;
auto res = get_divisors(x);
for(auto t : res)
cout << t << ' ';
cout << endl;
}
return 0;
}
求约数个数:
一个数
N =
则约数个数ans =
推导:
N的任何一个约数d =
每一项的不同,那么每一项的d就不同,每一个
都有0 ~
种方式
由乘法原理得到由ans个约数
#include <iostream>
#include <algorithm>
#include <unordered_map>
using namespace std;
typedef long long LL;
const int mod = 1e9 + 7;
//int范围内约数最多是1500多个
int main()
{
int n;
cin >> n;
unordered_map<int, int> primes;
while(n--)
{
int x;
cin >> x;
for(int i = 2; i <= x / i; ++i)
{
while(x % i == 0)
{
x /= i;
primes[i] ++;
}
}
if(x > 1)
primes[x]++;
}
LL res = 1;
for(auto ans : primes)
res = res * (ans.second + 1) % mod;
cout << res << endl;
return 0;
}
约数之和:
若: N =
则约数之和
公式:
while (b -- ) t = (t * a + 1) % mod;
t = t∗p+1
t =
t =
t =
……
t =
#include <iostream>
#include <unordered_map>
using namespace std;
typedef long long LL;
const int mod = 1e9 + 7;
int main()
{
int n;
cin >> n;
unordered_map<int, int> primes;
while(n--)
{
int x;
cin >> x;
for(int i = 2; i <= x / i; ++ i)
{
while(x % i == 0)
{
x /= i;
primes[i] ++;
}
}
if(x > 1)
primes[x] ++;
}
LL res = 1;
for(auto p : primes)
{
LL a = p.first, b = p.second;
LL t = 1;
while(b--)
t = (t * a + 1) % mod;
res = res * t % mod;
}
cout << res << endl;
return 0;
}
最大公约数:
求两个正整数 a 和 b 的 最大公约数 d
则有 gcd(a,b) = gcd(b,a%b)
证明:
设a%b = a - k*b 其中k = a/b(向下取整)
若d是(a,b)的公约数 则知 d|a 且 d|b 则易知 d|a-k*b 故d也是(b,a%b) 的公约数
若d是(b,a%b)的公约数 则知 d|b 且 d|a-k*b 则 d|a-k*b+k*b = d|a 故而d|b 故而 d也是(a,b)的公约数
因此(a,b)的公约数集合和(b,a%b)的公约数集合相同 所以他们的最大公约数也相同
#include <iostream>
using namespace std;
int gcd(int a, int b)
{
return b ? gcd(b, a % b) : a;
}
int main()
{
int n, a, b;
cin >> n;
while(n--)
{
cin >> a >> b;
cout << gcd(a, b) << endl;
}
return 0;
}