hoj 1867 哈尔滨工业大学 oj 经理的烦恼 简单的树状数组

经理的烦恼

My Tags	(Edit)   Cancel - Seperate tags with commas.

	Source : HCPC 2005 Spring
	Time limit : 2 sec		Memory limit : 32 M

Submitted : 1961, Accepted : 455

Jerry是一家公司销售部门的经理。这家公司有很多连锁店，编号为1，2，3，... Jerry每天必须关注每家连锁店的商品数量及其变化，一项很乏味的工作。在连锁店比较少的时候，Jerry喜欢计算编号在[i,j]区间内的连锁店中商品数量为素数的有多少家，但是现在连锁店的数量急剧增长，计算量很大，Jerry很难得出结果。 

输入格式
题目有多组输入。每组输入第一行有三个整数：C 连锁店的数量 N 指令的条数 M 每家连锁店初始的商品数量
接下来有N行，每行有一条指令。指令的格式为：
0 x y 连锁店x的商品数量变化值为y，y > 0商品数量增加， y < 0减少
1 i j 输出编号在[i,j]区间内的连锁店中商品数量为素数的有多少家
1 <= i, x, j < 1000000 连锁店中的商品数量a满足 0 <= a < 10000000，C = N = M = 0标志输入结束

输出格式
对于每组输入，输出它的序号。对于一组输入中的1指令输出要求的整数。每组输出后打印一行空行。

样例输入

100000 4 4
0 1 1
1 4 10
0 11 3
1 1 11

20 3 0
1 1 20
0 3 3
1 1 20

0 0 0

样例输出

CASE #1:
0
2

CASE #2:
0
1

简单的树状数组 的应用 对素数进行判断 原先是素数 后来变为不是素数 -1   原先不是素数后来变为素数 +1  注意0 1这两种不是素数的特殊情况

post code

/*This Code is Submitted by 1115332213 for Problem 1867 at 2012-09-05 23:18:52*/
#include<stdio.h>
#include<math.h>
#include<string.h>
int a[1000010],c,b[1000010];
int  isprime(int  a)   //素数的判断
{    
     int i;
     int n=int(sqrt( (double)(a) ));
     for(i=2; i<n+1; i++ )
     {
       if(a%i==0){return 0;}         
     }
     return 1;
}

int lowbit(int t)
{
  return t&(t^(t-1));    
}

int insert(int pos ,int num)
{
    while( pos<=c )
    {
      a[pos]+=num;
      pos+=lowbit(pos);             
    }    
}

int query( int pos)
{
    int sum=0;
    while( pos>0 )    
    {
      sum+=a[pos];
      pos-=lowbit(pos);       
    }
    return sum;
    
}


int main()
{
    int n,m,i,ji=0,judge,pos,num,beg,end ,flag1,flag2;
    while(scanf("%d %d %d",&c,&n,&m))
    {
        if(c==0&&n==0&&m==0)break;
        for( i=1; i<1000010; i++)
           b[i]=m;
        ji++;
        printf("CASE #%d:\n",ji);   
        memset( a, 0, sizeof(a) );
        if( isprime(m)==1 &&m!=0&&m!=1 ){
                                 for( i=1; i<=c; i++ )
                                 {
                                  insert(i,1);     
                                 }     
                                       }     
        for( i=1; i<=n; i++ )
        {
             scanf("%d",&judge);
             if(judge==0){
                           flag1=flag2=0;
                           scanf("%d %d",&pos,&num);
                           if( isprime(b[pos] )==1&&b[pos]!=0&&b[pos]!=1 )flag1=1;

                           b[pos]+=num;
                           if(isprime( b[pos] )==1&&b[pos]!=0&&b[pos]!=1) flag2=1; 
      
                           if( flag1==1 && flag2==0 )insert( pos,-1 );    //-1
                           if( flag1==0 && flag2==1 )insert( pos, 1 );    //+1
                         }
             else {
                     scanf("%d %d",&beg,&end);
                     printf("%d\n",query(end)-query(beg-1));
                  
                  
                  }
             
             
             
        }
        printf("\n");
                    
                    
    }
         
}