POJ 3292 数论 类似素数筛

POJ 3292 数论 类似素数筛

题意：一个定义的类型数 ：是A-number 是模 4 余 1的数

A-prime h是只能由 1* h（(即只能是1*A-number得到h ) 不能由其他的A-number数相乘得到）得到的数  (和素数没关系)

A-semis-primes 数是由且只能由 两个 A-prime数相乘得到

求一个数中存在多少 A-semis-primes

类似素数筛的做法

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
using namespace std;
const int Max=1000010;
int vis[1000010];
int sum[1000010];
int main()
{
    int n,t;
    memset(vis,0,sizeof(vis));
    for(int i=5;i<=1001;i+=4)
    {
        if(!vis[i])
        for(int j=5;(t=j*i)<=Max;j+=4)
        {
            if(!vis[j])                    <span id="transmark"></span>
                vis[t]=1;
            else                       //存在多个A-prime两两相乘的情况
                vis[t]=-1;
        }
    }
    sum[0]=0;
    for(int i=1;i<1000010;i++)
        sum[i]=sum[i-1]+(vis[i]==1);
    while(~scanf("%d",&n),n)
        printf("%d %d\n",n,sum[n]);
    return 0;
}