CF 113C

求区间 [l,r]   满足 素数 = a^2+b^2    a b为任意整数

打素数表，如何证明  此时的素数满足 %4==1 

#include<iostream>
#include<stdio.h>
#include<cstring>
#include<cstdlib>
#include<bitset>
using namespace std;
const int maxn=3e8+10;

bitset<maxn>prime;
int main()
{
    int l,r,i,j;
    while(scanf("%d%d",&l,&r)>0)
    {
        prime.set();
        for(i=3;(j=i*i)<=r;i=i+2)//i*i<=r
        {
            if(prime[i]==false)continue;
            for(;j<=r;j=j+i+i) // ！！！ j=j+i+i;
                prime[j]=false;
        }
        int ans=0;
        if(l<=2 && 2<=r) ans++;
        for(i=5;i<=r;i=i+4)
            if(l<=i && prime[i]==true) ans++;
        printf("%d
",ans);
    }
    return 0;
}