poj 3292 Semi-prime H-numbers (打表 同余模运算)

http://poj.org/problem?id=3292

题意:H-number为mod 4=1的数

H-prim为除了1和本身外不能被任何H-number整除的数

H-semi-prime为两个H-prim的乘积

求1-h有多少个H-semi-prime;

思路:打表

#include <iostream>
#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <limits>
#include <stack>
#include <vector>
#include <map>

using namespace std;

#define N 1000020
#define INF 0xfffffff
#define PI acos (-1.0)
#define EPS 1e-8
#define met(a, b) memset (a, b, sizeof (a))

typedef long long LL;

int Hprim[N];

void Init ()
{
    for (int i=5; i<N; i+=4)
    {
        for (int j=5; j<N; j+=4)
        {
            if (i*j>=N) break;

            if (!Hprim[i] && !Hprim[j] && i*j%4==1)
                Hprim[i*j] = 1;///<span style="font-family: KaiTi_GB2312;font-size:18px; line-height: 26px;">H-semi-prime</span>
            else Hprim[i*j] = -1;
        }
    }

    int cnt = 0;
    for (int i=5; i<N; i++)
    {
        if (Hprim[i]==1) cnt++;
        Hprim[i] = cnt;
    }
}

int main ()
{
    int n;
    Init();
    while (scanf ("%d", &n), n)
    {
        printf ("%d %d\n", n, Hprim[n]);
    }
    return 0;
}


阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/w144215160044/article/details/51547215
个人分类: 数论基础
上一篇poj 2635 The Embarrassed Cryptographer(同余模运算)
下一篇poj 1845 Sumdiv(同余模公式)
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭