数论(素数以及计算方程)

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

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
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<

二:

#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; iMAXL)
            {
                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
#include
#include
#include
#include
#include
#include
#include
#include
#include
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"<

你可能感兴趣的:(芝士就是力量)