数论
gcd与lcm
#include<algorithm>
using namespace std;
int a=12,b=20;
int main(){
cout<<__gcd(a,b)<<endl;//求gcd 位于algorithm库 注意加两个下划线
cout<<a/__gcd(a,b)*b<<endl;//求lcm 注意先除再乘防止溢出
}
扩展欧几里得
#include<iostream>
using namespace std;
void exgcd(int a,int b,int &d,int &x,int &y){
if(!b) {d=a;x=1;y=0;}
else{exgcd(b,a%b,d,y,x);y-=x*(a/b);}
}
int main(){
int x,y,d,a,b;
a=6,b=15;
exgcd(a,b,d,x,y);
cout<<x<<' '<<y;
}
素数
//素数判断 s1单个 s2筛表
int n;
bool vis[10000000],flag;
bool s1(int n){
int m=sqrt(n+0.5);
for(int i=2;i<=m;i++)if(!vis[i]){
if(!(n%i))return 0;
for(int j=i*i;j<=n;j+=i)vis[j]=1;
}
return 1;
}
void s2(int n){
int m=sqrt(n+0.5);
memset(vis,0,sizeof(vis));
for(int i=2;i<=m;i++)if(!vis[i])
for(int j=i*i;j<=n;j+=i) vis[j]=1;
}
int main(){
cin>>n;
cout << s1(n);
s2(n);for(int i=1;i<=n;i++)if(!vis[i])cout<<i<<' ';
}
素数
埃氏筛法
int m=sqrt(n+0.5);
memset(vis,0,sizeod(vis);
for(int i=2;i<=m;i++)if(!vis[i])
for(int j=i*i;j<=n;j+=i)vis[j]=1;
快速幂
int pow(int a,int b,int p)
{
int ans=1;
//a%=k;
for(;b;b>>=1,a=a*a,a%=p)
{
if(b&1) ans=ans*a;
ans%=p;
}
return ans;
}
int main()
{
int a,b,p;
cin>>a>>b>>p;
cout<<pow(a,b,p);
return 0;
}
排列组合
排列
STL 中有排列算法。头文件为algorithm
next_permutation(序列第一项的地址, 序列最后一项的地址+1):产生下一排列。
prev_permutation(序列第一项的地址, 序列最后一项的地址+1):产生上一排列。
组合
int C[1001][1001];//根据实际需要开数组,必要时采用高精度类型
memset(C,0,sizeof(C));
for (int i=0;i<=n;i++)
{
C[i][0]=1;
for (int j=0;j<=i;j++)
C[i][j]=C[i-1][j-1] + C[i-1][j];
}