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+;

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