1.求所有约数
给定 n 个正整数 ai,对于每个整数 ai,请你按照从小到大的顺序输出它的所有约数。
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
void get_prime(int x)
{
vector<int> res;
for(int i=1;i <= x / i;i++)
{
if(x % i == 0)
{
res.push_back(i);
if(i != x/i) res.push_back(x/i);
}
}
sort(res.begin(),res.end());
for(auto re :res) cout<<re<<" ";
cout<<endl;
}
int main()
{
int n;
cin>>n;
while(n--)
{
int x;
cin>>x;
get_prime(x);
}
return 0;
}
2. 约数个数
公式:所有质因数个数的乘积。[(a1+1)*(a2+1)*(a3+1)*......*(ak+1)]
一个正整数N = p1^a1 * p2^a2 * ...... pk^ak
p1,p2,p3,......,pk 是正整数的所有质因子
a1,a2,a3,......,ak 是对应的质因数存在的个数 (6 的质因子是 2 和3, 所以它有1个2,1个3,也就是2^1 * 3^1)
公式证明: 一个约数 X = p1^b1*p2^b2*......pk^bk (0<= bi <= ai)
(b1,b2,....,bk) bi 只要在每一个【0,a】区间内任选一个值,都可以得到一个值,这个值就是N的约数。
6的约数 1, 2 ,3 ,6;
1=2^0 * 3^0
2= 2^1 * 3^0
3=2^0 * 3^1
6=2^1 * 3^1
每一个质因子有a+1([0,a] 闭区间)种选法,所以k个质因子有(a1+1)*(a2+1)*(a3+1)*......*(ak+1)种选法,得到的选法种数就是N的约数个数
#include<iostream>
#include<unordered_map>
#include<algorithm>
using namespace std;
typedef long long LL;
const int mod = 1e9 + 7;
int main()
{
int n;
cin>>n;
unordered_map<int,int> prime;
while(n--)
{
int a;
cin>>a;
for(int i=2; i <= a / i; i++)
{
while(a % i == 0)
{
a /= i;
prime[i] ++;
}
}
if(a > 1) prime[a] ++;
}
LL res = 1;
for(auto p : prime)
{
res = res * (p.second + 1) % mod;
cout<<p.first<<" * "<<p.second<<endl;
}
cout<<res<<endl;
return 0;
}
3.约数之和
公式:(p1^0 + p1^1 + ..... +p1^k)(p2^0 + p2^1 + ..... +p2^k)......(pk^0 + pk^1 + ..... +pk^k)
证明:把公式按乘法分配律展开 (从上述每一个括号中选一个数相乘)+()+() + .....+ ()
可以看出展开之后每个括号都是一个约数,共有约数个(上面求过的约数个数)括号。
#include<iostream>
#include<unordered_map>
#include<algorithm>
using namespace std;
const int mod = 1e9 + 7;
typedef long long LL;
int main()
{
int n;
cin>>n;
unordered_map<int,int> prime;
while(n--)
{
int x;
cin>>x;
for(int i=2;i<=x/i;i++)
{
while(x%i==0)
{
x/=i;
prime[i] ++;
}
}
if(x>1) prime[x] ++;
}
LL res = 1;
for(auto p : prime)
{
int a = p.first ,b = p.second;
LL t = 1;
while(b--) t = (t * a + 1) % mod;//求每一项p的乘积
res = res * t % mod;
}
cout<<res<<endl;
return 0;
}
p^0+p^1+p^2+.....+p^k =
while循环
第一次循环 t = p+1
第二次 t = (p+1) *p + 1
......