UVa 11105 (筛法) Semi-prime H-numbers

题意：

你现在来到了一个所有的数都模4余1的世界，也就是除了这种数没有其他的数了。

然而素数的定义依然没变，如果一个数不能写成两个非1数字的乘积，则它是素数。

比如，在这里5就变成了最小的素数。

两个素数相乘得到一个半素数，比如5×5 = 25就是最小的半素数。

求1~h之间有多少个半素数。

分析：

虽然是要求[1, h]之间半素数的个数，但向往常筛普通素数一样先把所有的4k+1的素数筛出来。

然后二重循环枚举半素数，最后统计区间内[1, h]的半素数个数SUM_h

 #include <cstdio>
 #include <cstring>
 #include <cmath>
 using namespace std;

 const int maxn =  + ;
 bool vis[maxn + ];
 int prime[], cnt, sum[maxn + ];

 void Init()
 {
     int m = sqrt(maxn + 0.5);
     for(int i = ; i <= m; i += ) if(!vis[i])
         for(int j = i * i; j <= maxn; j += i)
             vis[j] = true;//H-素数筛选

     for(int i = ; i <= maxn; i += ) if(!vis[i]) prime[cnt++] = i;

     memset(vis, false, sizeof(vis));
     for(int i = ; i < cnt; i++)
     {//筛选H-半素数
         for(int j = i; j < cnt; j++)
         {
             long long k = (long long)prime[i] * (long long) prime[j];
             if(k > (long long)maxn) break;
             vis[k] = true;
         }
     }

     for(int i = ; i <= maxn; i++) sum[i] = sum[i-] + vis[i];
 }

 int main()
 {
     freopen("in.txt", "r", stdin);
     Init();
     int n;
     while(scanf("%d", &n) ==  && n) printf("%d %d\n", n, sum[n]);

     return ;
 }

代码君