poj-2689素数筛选

这道题目我就不贴题目了，就是素数筛选题。
这道题目差点让我崩溃，自己太菜的原因，提交了很多次，总是runtime error ，后来才知道，素数筛选的范围这能是2^16，然而这道题给的范围超出了int的上界，所以，必须要用另外一种方法来筛选，看了人家的博客，感觉方法超好。这道题有个突破点，就是区间在100000之内，所以，我们找到它们的合子的位置即可。
先筛出50000以内的素数，然后再在给定的区间内以刚刚的素数为基底，筛出区间内的素数。筛的过程中算出区间左值是最小素数的几倍，然后然后不是整数倍，就往前加1（否则算出的数会不在这个区间里，没意义），然后这个整数倍乘以素数已经大于右边了，那就可以直接退出了。

这里写代码
#include <iostream>
#include <cmath>
  #include<cstring>
using namespace std;
const int N=50000;
const int Max=999999999;//定义为无穷大
int r[1000000],a[N+100],b[N+100],z;
int main()
{
    int l,u,i,j;
    for (i=2;i<=N;i++)//素数筛选
        if (!a[i])
        {
                  b[++z]=i;
                  for (j=i+i;j<=N;j+=i)
                    a[j]=1;
        }
    while (cin>>l>>u)
    {
          memset(r,0,sizeof(r));
          int dis,max=-1,min=Max,maxc,minc;
          for (i=1;i<=z;i++)//关键算法在于这里
          {
              int s,t;
              s=(l-1)/b[i]+1;
              t=u/b[i];
              for (j=s;j<=t;j++)
                  if (j>1)
                    r[j*b[i]-l]=1;
          }
          int k=-1;
          for (i=0;i<=u-l;i++)//这里就是遍历找相邻最近和最远的，其实这里理解容易些
              if (!r[i])
              {
                        if (k!=-1)
                        {
                                  dis=i-k;
                                  if (dis>max)
                                  {
                                               max=dis;
                                                maxc=i+l;
                                  }
                                  if (dis<min)
                                  {
                                               min=dis;
                                               minc=i+l;
                                  }
                        }
                        if (i+l!=1) k=i;
              }
          if (max<0)
          {
                     cout<<"There are no adjacent primes."<<endl;
          }
          else
          {
              cout<<minc-min<<','<<minc<<" are closest, ";
              cout<<maxc-max<<','<<maxc<<" are most distant."<<endl;
          }
    }
    return 0;
}

希望各位大牛们多多指教！