模板之数论

数论去死！

数论什么都不会的蒟蒻QAQ

放两个模板就跑~

同余方程（exgcd）

#include
#include
#include
#include
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
#include
#include
#include
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<" ";
        cout<return 0;
}

组合数：

#include
#include
#include
#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
#include
#include
#include
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<' ';
    puts("");
    return 0;
}