HDU2098分拆素数和【打素数表+二分】

大素数表+二分

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

const int maxn = 10005;
int prm[maxn / 3 + 1], is[maxn / 32 + 1];
int n;
int k;
void get() {
    int N = maxn - 5;
    k = 0;
    for(int i = 2; i <= N; i++) {
        if(( is[i / 32] & ( 1 << ( i % 32 ) ) ) == 0) {
            prm[k++] = i;
        }
        for(int j = 0; j < k && prm[j] * i <= N; j++) {
            int num = prm[j] * i;
            is[num / 32] |= ( 1 << ( num % 32) );
            if(i % prm[j] == 0) {
                break;
            }
        }
    }
}
bool check(int num) {
    int low = 0; int high = k - 1;
    while(low <= high) {
        int mid = ( low + high ) >> 1;
        if(prm[mid] >= num) {
            high = mid - 1;
        } else {
            low = mid + 1;
        }
    }
    if(num == prm[high + 1]) return true;
    return false;
}

int main() {
    while(scanf("%d",&n) && n) {
        get();
        int m = n / 2;
        if(!(n & 1)) m--;
        int ans = 0;
        for(int i = 0; prm[i] <= m; i++) {
            if(check(n - prm[i])) {
                ans++;
//                printf("%d %d
", prm[i], n - prm[i]);
            }
        }
        printf("%d
", ans);
    }
    return 0;
}