hdu3792

http://acm.hdu.edu.cn/showproblem.php?pid=3792
求n以内的孪生素数的对数，孪生素数: p和p+2都为素数.
这题暴力能过，用树状数组处理一下更快

#include<stdio.h>
#include<string.h>
#include <iostream>
#include <cmath>
using namespace std;

const int maxn=1e5;
int vis[maxn+10];
int prime[maxn],c;
int a[maxn+10];
void isprime()
{
    int n=100000;
    memset(vis,0,sizeof(vis));
    for(int i=2;i*i<=n;i++)
    if(!vis[i])
    {
        for(int j=i*i;j<=n;j+=i)
        vis[j]=1;
    }
}
int main()
{
    int  i,j,n,max;
    isprime();
    n=0;
    a[1]=0;a[0]=0;a[2]=0;a[3]=0;
    for(i=4;i<=100000;i++)
    {

         if(!vis[i]&&!vis[i-2])
         n++;
         a[i]=n;//不这么写也会超时
    }
//    while(i<=100000)
//    {
//        a[i++]=n;
//    }
    while(cin>>n)
    {
        if(n<0) break;
        cout<<a[n]<<endl;
    }
    return 0;
}

树状数组

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

const int maxn=100005;
int c[maxn];
int vis[maxn];
int lowbit(int x)
{
    return x&-x;
}
void add(int x,int num)
{
    while(x<=maxn)
    {
        c[x]+=num;
        x+=lowbit(x);
    }
}
int sum(int x)
{
    int ret=0;
    while(x>0)
    {
        ret+=c[x];
        x-=lowbit(x);
    }
    return ret;
}

int flag[100000],p[10000];
void init()
{
    int i,j,k,tot=0;
    memset(flag,0,sizeof(flag));
    for(i=2;i<=100000;i++)
    {
        if(!flag[i])
        {
            p[++tot]=i;
            for(j=2*i;j<=100000;j+=i)
            {
                flag[j]=1;
            }
        }
    }
    for(i=2;i<=tot;i++)
    {
        if(p[i]-p[i-1]==2)
        {
            add(p[i],1);
        }
    }
}
int main()
{
    init();
    int n;
    while(scanf("%d",&n),n>0)
    {
        if(n==0) printf("0\n");
        else
        printf("%d\n",sum(n));
    }
    return 0;
}