LightOJ 1197 Help Hanzo(区间素数筛法)
Posted fy1999
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了LightOJ 1197 Help Hanzo(区间素数筛法)相关的知识,希望对你有一定的参考价值。
#include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #include<cmath> using namespace std; #define maxn 50000 int vis[maxn], isprime[5200], num[100005], k; void prime()//只需要把[1,sqrt(2^31)]之间的素数筛选出来就ok了。 { long long i, j; for (k=0,i=2; i<maxn; i++) if (vis[i] == 0) { isprime[k ++] = i; for (j=i*i; j<maxn; j+=i) vis[j] = 1; } //printf ("%lld ", k); } int main () { int t, l = 1; prime (); scanf ("%d", &t); while (t --) { int a, b, ans = 0; int n; scanf ("%d %d", &a, &b); n = b - a;//所求区间最大可达下标 memset (num, 0, sizeof(num)); for (int i=0; isprime[i]<=(int)sqrt(b)&& i<k; i++) { int j = 0; if (a % isprime[i] != 0 )//第一个需要筛掉的数(j+a) % isprime[i] == 0 j = j - a % isprime[i] + isprime[i]; if (a <= isprime[i])//(j+a) / isprime[i] == 1,则(j+a)是素数,要向下推一个 j += isprime[i]; for ( ; j<=n; j+=isprime[i]) { num[j] = 1; } } for (int i=0; i<=n; i++)//计算素数的数目 if (!num[i]) ans ++; if (a == 1)//对这种情况特殊处理 ans --; printf ("Case %d: %d ", l ++, ans); } return 0; }
以上是关于LightOJ 1197 Help Hanzo(区间素数筛法)的主要内容,如果未能解决你的问题,请参考以下文章
LightOJ 1197 Help Hanzo(区间素数筛法)
B - Help Hanzo (LightOJ - 1197)
M - Help Hanzo LightOJ - 1197 (大区间求素数)