HOJ 1867

树状数组+素数打表

#include <cstdio>
#include <iostream>
#include <string.h>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=10000100;
const int maxm=1000010;
int c[maxm],now[maxm],k,n,m;
bool phi[maxn];
void getphi(){
    memset(phi,1,sizeof(phi));
    int q=(int)sqrt(maxn)+1;
    phi[0]=phi[1]=0;
    for(int i=2;i<=q;i++){
        if(phi[i]){
            for(int j=i*2;j<maxn;j+=i)
                phi[j]=0;
        }
    }
}
int lowbit(int x){return x&(-x);}
int sum(int b){
    int sum=0;
    while(b>0){
        sum+=c[b];
        b-=lowbit(b);
    }
    return sum;
}
void add(int x,int val){
    while(x<=k){
        c[x]+=val;
        x+=lowbit(x);
    }
}
int main(){
    int a,x,y,game=0;
    getphi();
    while(scanf("%d%d%d",&k,&n,&m)!=EOF){
        if(k==0&&n==0&&m==0)break;
        memset(c,0,sizeof(c));
        for(int i=1;i<=k;i++){
            now[i]=m;
            if(phi[m])add(i,1);
        }
        printf("CASE #%d:
",++game);
        for(int i=0;i<n;i++){
            scanf("%d%d%d",&a,&x,&y);
            switch (a){
                case 0:{
                    bool b=phi[now[x]],c=phi[now[x]+y];
                    if(b&&!c)add(x,-1);
                    if(!b&&c)add(x,1);
                    now[x]+=y;
                    break;
                }
                case 1:{
                    printf("%d
",sum(y)-sum(x-1));
                }
            }
        }
        printf("
");
    }
    return 0;
}