数论
1.定义
1.1.约数:y是x的约数(x%y==0)
1.2.质数:除了1和本身没有其他约数(x%y!=0)
1.3.合数:质数的相反(不包括1)
1.4.gcd(最大公约数):
__gcd(x,y); //c++函数:
__gcd(x,y) == 1;//互质
1.5.Lcm(最小公倍数):
a*b/__gcd(a,b);
1.6.唯一分解定理(质因分解):
x=p1a1*p2a2*p3a3…pna3
1.7.取模(防止数过大):
( 1 + 2 + … + n ) % 10
( a + b ) % 10 = ( a % 10 + b % 10 ) % 10; //加法取模
( a * b ) % 10= ( ( a % 10 ) * ( b % 10 ) ) % 10; //乘法取模
2.算法模型
2.1质因分解
所有约数有一个特点从sqrt(x)对称
sqrt(x)~x只有一个质因子
#include <bits/stdc++.h>
using namespace std;
vector<int> v;
int x;
int main(){
cin >> x;
for(int i=2;i*i<=x;i++){
if(x%i==0){
while(x%i==0){
x/=i;
v.push_back(i);
}
}
}
if(x!=1) v.push_back(x);//如果有剩下的一定就是sqrt(x)~x的质因子
return 0;
}
3.所有约数
所有约数有一个特点从sqrt(x)对称
#include <bits/stdc++.h>
using namespace std;
vector<int> v;
int x;
int main(){
cin >> x;
for(int i=2;i*i<=x;i++){
if(x%i==0){
v.push_back(i);
v.push_back(x/i);
}
}
return 0;
}
4.素数筛(求1~n之间所有的素数)
从1~n里的素数没有一个是他的约数
#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 5;
vector<int> p;
vector<int> d[maxn];
int bk[maxn];
int main()
{
int n; cin >> n;
// 方法一:O(nloglogn)
for (int i = 2 ; i <= n ; i++){
if (bk[i]) continue;
p.push_back(i);
for (int j = i + i ; j <= n ; j += i) bk[j] = 1;
}
memset (bk , 0 , sizeof bk);
// 用方法一来求1~n中每个数的质因子
for (int i = 2 ; i <= n ; i++){
if (bk[i]) continue;
p.push_back(i);
d[i].push_back(i);
for (int j = i + i ; j <= n ; j += i) {
bk[j] = 1;
d[j].push_back(i);
}
}
memset (bk , 0 , sizeof bk);
// 方法二:O(nlogn)
for (int i = 2 ; i <= n ; i++){
if (!bk[i]) p.push_back(i);
for (int j = i + i ; j <= n ; j += i) bk[j] = 1;
}
memset(bk , 0 , sizeof bk);
for (int i = 1 ; i <= n ; i++)
d[i].push_back(1);
// 方法二用来求1~n中每个数的约数.
for (int i = 2 ; i <= n ; i++){
if (!bk[i]) p.push_back(i);
d[i].push_back(i);
for (int j = i + i ; j <= n ; j += i) {
bk[j] = 1;
d[j].push_back(i);
}
}
// d[j] j
for (int i = 2 ; i <= n ; i++){
cout << i << " ";
for (auto g : d[i]) cout << g << " ";
cout << endl;
}
return 0;
}
组合数学
1.容斥定理
1.1.从1~n被5整除的数有多少
[1,n] n/5
1.2.[1~N]不能被2 3整除
n-n/2-n/3+n/(2*3)
2.抽屉原理
K倍区间
n个格子放n+1个球至少一个格子有2个球
| ai - a j |是100的倍数
( ai - aj ) % 100 = 0
ai%100=aj%100
bi=bj
n>100 Yes
n<100 O(n2)