数论去死!
数论什么都不会的蒟蒻QAQ
放两个模板就跑~
同余方程(exgcd)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
void exgcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;
y=0;
return;
}
else
{
exgcd(b,a%b,x,y);
int t=x;
x=y;
y=t-a/b*x;
}
}
int main()
{
int a,b,x,y;
scanf("%d%d",&a,&b);
exgcd(a,b,x,y);
printf("%d\n",(x+b)%b);
return 0;
}
矩阵乘法:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int a[300][300],b[300][300],c[300][300];
int main()
{
int ai,aj;
scanf("%d%d",&ai,&aj);
for(int i=1;i<=ai;i++)
for(int j=1;j<=aj;j++)
cin>>a[i][j];
int bi,bj;
scanf("%d%d",&bi,&bj);
for(int i=1;i<=bi;i++)
for(int j=1;j<=bj;j++)
cin>>b[i][j];
for(int i=1;i<=ai;i++)
{
for(int j=1;j<=aj;j++)
{
for(int k=1;k<=bj;k++)
c[i][k]+=a[i][j]*b[j][k];
}
}
for(int i=1;i<=ai;i++)
{
for(int j=1;j<=bj;j++)
cout<<c[i][j]<<" ";
cout<<endl;
}
return 0;
}
组合数:
#include<iostream>
#include<cstdio>
#include<cstring>
#ifdef unix
#define LL "%lld"
#else
#define LL "%I64d"
#endif
using namespace std;
const int MOD=100003;
int C[1000+10][1000+10];
int main()
{
int n, m;
scanf("%d%d",&n,&m);
n%=MOD,m%=MOD;
for(int i=0;i<=n;i++) C[i][0]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
C[i][j]=(C[i-1][j-1]+C[i-1][j])%MOD;
printf(LL,C[n][m]%MOD);
return 0;
}
/*
234 2
2 5
3 6
*/
欧拉筛法:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n;
int prime[10000];
bool bo[10000];
int cnt=0;
void get()
{
for(int i=2;i<=n;i++)
{
if(!bo[i])
prime[++cnt]=i;
for(int j=1;j<=cnt;j++)
{
if(i*prime[j]>n) break;
bo[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
}
int main()
{
scanf("%d",&n);
get();
for(int i=1;i<=cnt;i++) cout<<prime[i]<<' ';
puts("");
return 0;
}