CodeForces114E——Double Happiness(素数二次筛选)

Double Happiness

On the math lesson a teacher asked each pupil to come up with his own lucky numbers. As a fan of number theory Peter chose prime numbers. Bob was more original. He said that number t is his lucky number, if it can be represented as:
t = a2 + b2, 
where a, b are arbitrary positive integers.
Now, the boys decided to find out how many days of the interval [l, r] (l ≤ r) are suitable for pair programming. They decided that the day i (l ≤ i ≤ r) is suitable for pair programming if and only if the number i is lucky for Peter and lucky for Bob at the same time. Help the boys to find the number of such days.
Input
The first line of the input contains integer numbers l, r (1 ≤ l, r ≤ 3*10^8).
Output
In the only line print the number of days on the segment [l, r], which are lucky for Peter and Bob at the same time.
Sample test(s)
input
3 5
output
1
input
6 66
output
7

题目大意：

　　　　给定区间[L,R]，求区间内满足条件的数的个数：

　　　　条件1）z是素数

　　　　条件2）z=x^2+y^2 x和y为任意的正整数

解题思路：

　　　　其实就是求满足费马定理的数的个数。

　　　　费马定理：一个奇素数z可以表示成z=x^2+y^2的形式,当且仅当z可以表示成4*t+1的时候。（偶素数2=1^2+1^2）

　　　　LR的范围是1到3*10^8用普通的筛选法求素数表，时间空间都会超。使用两次筛选来求素数。

　　　　3*10^8的平方根小于17500，用17500以内的素数可以筛出范围内的所有素数。在通过判断是否满足z=t%4+1来累加。

　　　　又由于z的范围过大，但对于唯一确定的z来说，t也唯一确定，故可以用t作为数组下标即(z-1)/4。数组大小就会小4倍左右。

　　　　具体细节看代码备注。

Code：

 #include <algorithm>
 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <string>
 #include <bitset>
 #define MAXN 17500
 using namespace std;
 bitset<MAXN>ISP; //类似于bool型，可以存01
 bitset</+>result;
 int prime[MAXN];
 int is_prime()//先筛选1-17500的素数
 {
     int p=;
     ISP[]=ISP[]=;
     for (int i=; i<=MAXN; i++)
         if (ISP[i]==)
         {
             prime[p++]=i;//用prime素组将素数存起来留到后面筛选用。。
             for (int j=i+i; j<=MAXN; j+=i)
                 ISP[j]=;
         }
     return p-;
 }
 int cnt(int L,int R)
 {
     int p=is_prime();
     for (int j=; j<=p; j++)
     {
         int x=L/prime[j];
         if (x*prime[j]<L) x++;
         if (x==) x++;
         for (int k=x*prime[j]; k<=R; k+=prime[j]) //通过素数表再来筛选[L,R]中的素数
                if (k%==) //标记符合条件2且不符合条件1的数
                 result[(k-)/]=;
     }
     int cnt=;
     for (int i=L;i<=R;i+=)
         if (!result[(i-)/]) cnt++;
     return cnt;
 }
 int main()
 {
     is_prime();
     int L,R;
     int ok=;
     cin>>L>>R;
     if (L<=&&R>=) ok=; //2作为偶素数且符合条件 单独判断
     while (L%!=||L==) //将区间边界缩小到符合 %4==1
         L++;
     while (R%!=)
         R--;
     cout<<cnt(L,R)+ok;
     return ;
 }