判断一个数是否为质数/素数 原网址:点击打开链接
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;
}