数论（素数以及计算方程）

判断一个数是否为质数/素数   原网址：点击打开链接

bool isPrime(int num)
{
    if(num==2||num==3)return 1;
    if(num%6!=1&&num%6!=5)return 0;
    int temp=sqrt(num);
    for(int i=5;i<=temp;i+=6)
    {
        if(num%i==0||num%(i+2)==0)
        {
            return 0;
        }
    }
    return 1;
}

欧拉线性筛：

一：

#include<bits/stdc++.h>
using namespace std;
int prime[11000000+5];//素数表
int primeSize;//素数个数
int isPrime[11000000+5];//记录i是否为素数
void getlist(int listSize)
{
    memset(isPrime,1,sizeof(isPrime));
    isPrime[1]=0;
    primeSize=1;
    for(int i=2; i<=listSize; i++)
    {
        if(isPrime[i])
        {
            prime[primeSize++]=i;
        }
        for(int j=1; j<=primeSize&&i*prime[j]<=listSize; j++)
        {
            isPrime[i*prime[j]]=0;
            if(i%prime[j]==0)break;
        }
    }
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        getlist(n);
        cout<<primeSize-1<<endl;
    }
    return 0;
}

二：

#define MAXN 10005
#define MAXL 1299710

int prime[MAXL];
int check[MAXL];

void init()
{
    int tot=0;
    memset(check,0,sizeof(check));
    for(int i=2; i<MAXL; i++)
    {
        if(!check[i])
        {
            prime[tot++]=i;
        }
        for(int j=0; j<tot; j++)
        {
            if(i*prime[j]>MAXL)
            {
                break;
            }
            check[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                break;
            }
        }
    }
}

欧几里得算法

int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}

计算ax+by=c模板

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<cmath>
#include<map>
#include<stack>
#include<vector>
#include<queue>
#include<set>
using namespace std;
const int MAX_N=10010;
typedef long long ll;
typedef long long LL;
int exgcd(int a,int b,int &x,int &y)
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    int e=exgcd(b,a%b,x,y);
    int zz=x;
    x=y;
    y=zz-a%b*y;
    return e;
}
int main()
{
    int a,b,c;
    int x,y;
    while(scanf("%d%d%d",&a,&b,&c)!=EOF)
    {
      int k=exgcd(a,b,x,y);
      if(c%k)
      {
        cout<<"Impossible"<<endl;
      }
      else
      {
        k=c/k;
        x*=k;
        y*=k;
        printf("x=%d  y=%d\n",x,y);
      }
    }
    return 0;
}