24h購物| | PChome| 登入
2012-10-12 13:00:47| 人氣1,417| 回應1 | 上一篇 | 下一篇

[ZJ][sieve] a505. B. T-primes

推薦 0 收藏 0 轉貼0 訂閱站台

內容 :

We know that prime numbers are positive integers that have exactly two distinct positive divisors. Similarly, we'll call a positive integer t Т-prime, if t has exactly three distinct positive divisors.
You are given an array of n positive integers. For each of them determine whether it is Т-prime or not.

輸入說明 :

The first line contains a single positive integer, n (1 ≤ n ≤ 106), showing how many numbers are in the array. The next line contains nspace-separated integers xi (1 ≤ xi ≤ 1014).

輸出說明 :

Print n lines: the i-th line should contain "YES" (without the quotes), if number xi is Т-prime, and "NO" (without the quotes), if it isn't.

範例輸入 :

3
4 5 6

範例輸出 :

YES
NO
NO

提示 :

 n,xi範圍與原題不同!

 測試資料較大!請使用較快的輸出入!

2012/10/11 加強測資! 

出處 :

CodeForce 230B T-primes (管理:eddy841021)


/**********************************************************************************/
/*  Problem: a505 "B. T-primes" from CodeForce 230B T-primes                      */
/*  Language: C (687 Bytes)                                                       */
/*  Result: AC(350ms, 1.5MB) judge by this@ZeroJudge                              */
/*  Author: morris1028 at 2012-10-11 23:31:01                                     */
/**********************************************************************************/


#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define LL long long
#define maxL (10000000>>5)+1
#define GET(x) (mark[x>>5]>>(x&31)&1)
#define SET(x) (mark[x>>5] |= 1<<(x&31))
int mark[maxL];
void sieve() {
    register int i, j, k;
    SET(1);
    int n = 10000000;
    for(i = 2; i <= n; i++) {
        if(!GET(i)) {
            for(k = n/i, j = i*k; k >= i; k--, j -= i)
                SET(j);
        }
    }
}
main() {
    sieve();
    int t;
    LL n, m;
    scanf("%d", &t);
    while(t--) {
        scanf("%lld", &n), m = sqrt(n);
        if(m*m != n || GET(m))
            puts("NO");
        else    puts("YES");
    }
    return 0;
}

台長: Morris
人氣(1,417) | 回應(1)| 推薦 (0)| 收藏 (0)| 轉寄
全站分類: 不分類 | 個人分類: ZeroJudge |
此分類下一篇:[ZJ][插頭DP] a228. 就少一個插座用很不方便
此分類上一篇:[ZJ] a473. 外援計畫

卡在時限的人......
好複雜@@
可以稍微指點一下嗎?
2012-12-25 23:52:47
版主回應
首先你要先了解基本的 sieve 法去找質數,接著利用位元運算讓時間降下來,位元運算比較難一點,由於一般都是 0/1 去標記陣列,那麼多餘的 bit 就很浪費,因此將 int 32 個表示 0~31 接著是 32~63 ...,空間少了速度就快了。
2012-12-26 09:43:43
是 (若未登入"個人新聞台帳號"則看不到回覆唷!)
* 請輸入識別碼:
請輸入圖片中算式的結果(可能為0) 
(有*為必填)
TOP
詳全文