Morass

Untitled

Dec 22nd, 2016
215
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.33 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<29)
  6. #define CL(A,I) (memset(A,I,sizeof(A)))
  7. #define DEB printf("DEB!\n");
  8. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  9. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  10. typedef long long ll;
  11. typedef long double ld;
  12. typedef pair<ll,ll> pll;
  13. typedef vector<int> vi;
  14. typedef pair<int,int> ii;
  15. typedef vector<ii> vii;
  16. #define IN(n) int n;scanf("%d",&n);
  17. #define FOR(i, m, n) for (int i(m); i < n; i++)
  18. #define REP(i, n) FOR(i, 0, n)
  19. #define F(n) REP(i, n)
  20. #define FF(n) REP(j, n)
  21. #define FT(m, n) FOR(k, m, n)
  22. #define aa first
  23. #define bb second
  24. void ga(int N,int *A){F(N)scanf("%d",A+i);}
  25. #define MX (10000009)
  26. #define IS(A) ((cn[A>>6]>>(A&63))&1)
  27. #define ST(A) (cn[A>>6]|=1ULL<<((A&63)))
  28. ll cn[1+(MX>>6)];
  29. int P[MX>>1]={2},L=1;
  30. void gen(){
  31.     ST(0),ST(1);
  32.     int Q(1+sqrt(MX));
  33.     for(ll i(4);i<MX;i+=2)ST(i);
  34.     for(int k(3);k<=Q;k+=2)if(!IS(k))
  35.         for(ll h(k<<1),j(k*k);j<MX;j+=h)ST(j);
  36.     for(int i(3);i<MX;i+=2)if(!IS(i))P[L++]=i;
  37. }
  38. ll sm(int N){
  39.     ll S=0,p=upper_bound(P-1,P+L,N)-P-1;
  40.     FT(1,L&&P[k]<N){
  41.         while(p>0&&P[k]+P[p]>N)--p;
  42.         S+=p-(p>=k);
  43.     }
  44.     return S>>1;
  45. }
  46. int B,E;
  47. int main(void){
  48.     gen();
  49.     while(scanf("%d%d",&B,&E),B|E)printf("%lld\n",sm(E/2*2)-sm(--B/2*2));
  50.     return 0;
  51. }
Advertisement
Add Comment
Please, Sign In to add comment