HDU 4548 美素数 在线打表加数状数组

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4548

解题报告：一开始本想先打个素数表，然后每次输入L 跟R 的时候都进行暴力判断，但这题测试数据太多，有1W ，这样就达到10^10级别了，很明显不能过，但是，我后来一想，数状数组不就是用来解决区间问题的吗，顿时焕然大吾。我的做法就是先打个1到1000000的素数表，然后暴力判断从1到1000000每个数是不是美素数，同时把结果存到数状数组里面去。这样就可以实现存完了就可以实现快速查找了。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>

const int maxn = 1000000;

int prim[maxn+5],tree[maxn+5];

void dabiao()
{
    memset(prim,0,sizeof(prim));
    prim[1] = 1;
    for(int i = 4;i <= maxn;i += 2)
    prim[i] = 1;
    int M = sqrt(maxn);
    for(int i = 2;i <= M;++i)
    if(prim[i] == 0)
    for(int j = i*i;j<=maxn;j+=i)
    prim[j] = 1;
}


int judge(int n)
{
    int sum = 0;
    while(n)
    {
        sum += n%10;
        n /= 10;
    }
    return sum;
}

void add(int n)
{
    if(n > maxn)
    return ;
    for(int i = n;i <= maxn;i += (i&-i))
    tree[i]++;
}
void dabiao2()
{
    memset(tree,0,sizeof(tree));
    for(int i = 1;i <= maxn;++i)
    if(!prim[i] && !prim[judge(i)])
    add(i);
}

int sum(int L,int R)
{
    int tot1 = 0,tot2 = 0;
    for(int i = L-1;i >= 1; i -= (i&-i))
    tot1 += tree[i];
    for(int i = R;i >= 1;i -= (i&-i))
    tot2 += tree[i];
    return tot2 - tot1;
}

int main()
{
    dabiao();
    dabiao2();
    int T,L,R,l = 1;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d %d",&L,&R);
        printf("Case #%d: %d
",l++,sum(L,R));
    }
    return 0;
}