HOJ 1876经理的烦恼

http://acm.hit.edu.cn/hoj/problem/view?id=1867

哎~这道题re了很多次。后来找到原因竟然是 #define maxn 1000010+5   的错误，因为这里的C++提交不支持 宏定义有+的运算改为    #define maxn 1000010就过了

解题思路：树状数组C[i]存的是A[i-lowbit[i]+1]~A[i]中的素数的个数。之前进行一下筛素数就行了。

经理的烦恼

My Tags  (Edit)

Source : HCPC 2005 SpringTime limit : 2 secMemory limit : 32 M

Submitted : 2841, Accepted : 652

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 40 1 11 4 100 11 31 1 1120 3 01 1 200 3 31 1 200 0 0

样例输出

CASE #1:02CASE #2:01

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<vector>
#include<cmath>
#include<stdlib.h>
#include<iomanip>
#include<list>
#include<deque>
#include<map>
#include <stdio.h>
#include <queue>

#define maxn 1000010
#define ull unsigned long long
#define ll long long
#define reP(i,n) for(i=1;i<=n;i++)
#define rep(i,n) for(i=0;i<n;i++)
#define cle(a) memset(a,0,sizeof(a))
#define mod 90001
#define PI 3.141592657
#define INF 1<<30
const ull inf = 1LL << 61;
const double eps=1e-5;

using namespace std;

bool cmp(int a,int b){
    return a>b;
}
int a[maxn];
int c[maxn];
int n,C,m;
bool prime[maxn*10];
void setprime()
{
memset(prime,1,sizeof prime);
prime[1]=0;
prime[0]=0;
for(int i=2;i<=maxn*10;i++)
{
if(prime[i])
{
for(int j=i+i;j<=maxn*10;j+=i)
prime[j]=0;
}
}
}
int lowbit(int i)
{
return i&(-i);
}
void add(int i,int d)
{
while(i<=C)
{
c[i]+=d;
i+=lowbit(i);
}
}
int sum(int i)
{
int ret=0;
while(i>0)
{
ret+=c[i];
i-=lowbit(i);
}
return ret;
}
int t;
int k,x,y;
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int cas=0;
    setprime();
    while(~scanf("%d%d%d",&C,&n,&m),m+C+n)
{
cas++;
   // cle(c);cle(a);
    if(prime[m])t=1;
    else t=0;
    printf("CASE #%d:\n",cas);
    for(int i=1;i<=C;i++)
{
a[i]=m;
}
for(int i=1;i<=C;i++)
{
c[i]=t*lowbit(i);
}
for(int j=1;j<=n;j++)
{
scanf("%d%d%d",&k,&x,&y);
if(k==0)
{
a[x]+=y;
if(!prime[a[x]]&&prime[a[x]-y])
{
add(x,-1);
}
else if(prime[a[x]]&&!prime[a[x]-y])
add(x,1);
}
else if(k==1)
{
printf("%d\n",sum(y)-sum(x-1));
}
}
printf("\n");
}
    return 0;
}