常用的是 素筛打表
代码
离线算法 素筛
#include <bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
const int MAXN = 1e7; // ***
const int MAXM = 1e6+100;
const double PI = acos(-1.0);
const double eps = 1e-8;
inline int read(){
int x=0,f=1; char ch=getchar();
while (ch<'0'||ch>'9') { if (ch=='-') f=-1; ch=getchar(); }
while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
return x*f;
}
/***********************************************/
bool su[MAXN]={1,1,0};
void dabiao(){ //埃及筛
for(int i=2;i*i<n;i++)
if(!su[i])
for(int j=i*i;j<n;j++) // 从i*i开始就行
su[j]=1;
}
int main(){
}
代码
费马定理 在线算法
代码
#include <bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
const int maxn = 100000;
const int MAXM = 1e6+100;
const double PI = acos(-1.0);
const double eps = 1e-8;
inline int read(){
int x=0,f=1; char ch=getchar();
while (ch<'0'||ch>'9') { if (ch=='-') f=-1; ch=getchar(); }
while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
return x*f;
}
/*****************************************/
int pri[maxn],f[maxn],top;
bool vis[maxn];
int l,r,mid,n;
int G(int n,int m)
{
if (!m) return n;
if (m==1) return n-n/2;
if (n<=maxn)
{
if (f[n]<=m) return 1;
if (f[(int)sqrt(n)]<=m) return f[n]-m+1;
}
return G(n,m-1)-G(n/pri[m],m-1);
}
bool pan(int x)
{
int y=sqrt(x);
return f[y]+G(x,f[y])-1>=n;
}
int main()
{
for (int i=2;i<=maxn;++i)
{
if (!vis[i]) pri[++top]=i;
for (int j=1;j<=top&&pri[j]*i<=maxn;++j)
{
vis[i*pri[j]]=1;
if (i%pri[j]==0) break;
}
}
for (int i=2;i<=maxn;++i) f[i]=f[i-1]+(!vis[i]);
int k = 0;
/**************************************/
while(~scanf("%d",&n),n) { // 前第n个素数是什么
l=1,r=300000000;
while(l<r)
{
mid=(l+r)/2;
if (pan(mid)) r=mid;
else l=mid+1;
}
printf("Case %d: %d\n",++k,l);
}
}