文件读入读出
假设题目名为“add”，那么文件夹名为“add”，c++程序名为“add.cpp”，读入文件名为“add.in”，输出文件名为“add.out”。四个的拼写均不可有误，包括大小写差异。千万不要调试后就忘记修改文件读入读出了。
#include<cstdio>
int main()
{
freopen("add.in","r",stdin);//read
freopen("add.out","w",stdout);//write
}数论算法
1、线性筛（素数）：暂时代码是转载的，以后有机会会更新，看不懂请跳过
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
const int maxn=200005;
int prime[maxn];
bool not_prime[maxn];
int main()
{
int n,cnt=0;
scanf("%d",&n);
memset(not_prime,0,sizeof(not_prime));
not_prime[1]=true;
for(inti=2;i<=n;i++)
{
if(!not_prime[i])
prime[++cnt]=i;
for(intj=1;j<=cnt;j++)
{
if(prime[j]*i>n) break;
not_prime[prime[j]*i]=true;
if(i%prime[j]==0) break;
}
}
for(inti=1;i<=cnt;i++)
printf("%d ",prime[i]);
return 0;
}2、埃氏筛（素数+读写外挂）暂时代码是转载的，以后有机会会更新，看不懂请跳过
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
inline int read()
{
char c;
c=getchar();
for(;c>'9'||c<'0';c=getchar());
int sum=0;
for(;c<='9'&&c>='0';c=getchar())sum=sum*10+c-48;
return sum;
}
int n,m;
bool f[10010000];
int main()
{
n=read(),m=read();
memset(f,true,sizeof(f));
f[0]=f[1]=false;
for(inti=2;i<=n;i++)
if(f[i])
for(intj=2;i*j<=n;f[i*j]=false,j++);
for(inti=1;i<=m;i++)
{
int x;
x=read();
if(f[x])printf("Yes\n");
elseprintf("No\n");
}
return 0;
}3、拓展欧几里得算法（求逆元公式） 暂时代码是转载的，以后有机会会更新，看不懂请跳过
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
using namespace std;
int x,y;
int exgcd(int a,int b)
{
if(!b)
{
x=1;
y=0;
return a;
}
else
{
int t;
intd=exgcd(b,a%b);
t=x;
x=y;
y=t-a/b*x;
return d;
}
}
int main()
{
int a,b;
scanf("%d%d",&a,&b);
exgcd(a,b);
// cout<<__gcd(a,b)<<endl;
cout<<x<<" "<<y<<endl;
return 0;
}4、GCD&LCM（求最小公倍数与最大公约数） 暂时代码是转载的，以后有机会会更新，看不懂请跳过
#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
int gcd(int a,int b)
{
if(!b) returna;
else returngcd(b,a%b);
}
int lcm(int a,int b)
{
returna/gcd(a,b)*b;
}
int main()
{
int a,b;
scanf("%d%d",&a,&b);
cout<<gcd(a,b)<<" "<<lcm(a,b)<<endl;
return 0;
}5、分解质因数 暂时代码是转载的，以后有机会会更新，看不懂请跳过
#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
int main()
{
long long n;
scanf("%lld",&n);
for(long longi=2;i<=n;i++)
{
while(n!=i)
{
if(n%i==0)
{
printf("%lld*",i);
n=n/i;
}
elsebreak;
}
}
printf("%lld",n);
return 0;
}6、大数翻倍法 暂时代码是转载的，以后有机会会更新，看不懂请跳过
#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;
int a[233],mo[233];
int gcd(int a,int b)
{
if(!b) returna;
else returngcd(b,a%b);
}
int lcm(int a,int b)
{
returna/gcd(a,b)*b;
}
int main()
{
int n;
scanf("%d",&n);
for(inti=1;i<=n;i++)
scanf("%d%d",&a[i],&mo[i]);
intans=0,nowmo=1;
for(inti=1;i<=n;i++)
{
intother=a[i],othermo=mo[i];
if(othermo>nowmo)
{
swap(ans,other);
swap(nowmo,othermo);
}
while(ans%othermo!=other)
ans+=nowmo;
nowmo=lcm(nowmo,othermo);
}
printf("%d",ans);
}7、快速幂
模板题：luogu1226：快速幂
8、位运算
| 位运算 |
示例 |
功能 |
| x >> 1 |
(101101->10110) |
去掉最后一位 |
| x << 1 |
(101101->1011010) |
在最后加一个0 |
| x << 1+1 |
(101101->1011011) |
在最后加一个1 |
| x | 1 |
(101100->101101) |
把最后一位变成1 |
| x | 1-1 |
(101101->101100) |
把最后一位变成0 |
| x ^ 1 |
(101101->101100) |
最后一位取反 |
| x | (1 << (k-1)) |
(101001->101101,k=3) |
把右数第k位变成1 |
| x & !(1 << (k-1)) |
(101101->101001,k=3) |
把右数第k位变成0 |
| x ^ (1 << (k-1)) |
(101001->101101,k=3) |
右数第k位取反 |
| x & 7 |
(1101101->101) |
取末三位 |
| x & (1<< k -1) |
(1101101->1101,k=5) |
取末k位 |
| x >> (k-1) & 1 |
(1101101->1,k=4) |
取右数第k位 |
| x | (1 << k-1) |
(101001->101111,k=4) |
把末k位变成1 |
| x ^ (1 << k-1) |
(101001->100110,k=4) |
末k位取反 |
| x & (x+1) |
(100101111->100100000) |
把右边连续的1变成0 |
| x | (x+1) |
(100101111->100111111) |
把右起第一个0变成1 |
| x | (x-1) |
(11011000->11011111) |
把右边连续的0变成1 |
| (x ^ (x+1)) >> 1 |
(100101111->1111) |
取右边连续的1 |
| x & (x ^ (x-1)) |
(100101000->1000) |
去掉右起第一个1的左边 |