容斥+dfs求1~n质数个数

n<=100000000

  #include<stdio.h>
 #include<string.h>
 #include<math.h>

int flag[10006];
 int num[10000],total,now,n;

void solve(int index,int Mul,int K){
       int i,t;
    if(K==0){
     now+=n/Mul;
     return;
    }
    for(i=index;i<total-K+1;i++){
      Mul*=num[i];
      t=now;
      if(Mul<=n){
               solve(i+1,Mul,K-1);
      }
      if(t==now)//优化剪枝(1return;
            Mul/=num[i];
    }
 }

int Judge(int x){
  int t,i;
  t=sqrt(x*1.0);
  for(i=2;i<=t;i++)
   if(n%i==0)
    return 0;
  return 1;
 }

int main(){
     int i,j,t,Count;
  memset(flag,0,sizeof(flag));
  for(i=2;i<=10000;i++){
   if(flag[i])
    continue;
   for(j=2;j*i<=10000;j++)
    flag[i*j]=1;
  }
  while(scanf("%d",&n)==1){
   Count=0;
   if(n<=10000){
    for(i=2;i<n;i++)
     if(!flag[i])
      Count++;
   }
   else{
    t=sqrt(1.0*n);
    total=0;
    for(i=2;i<=t;i++)
     if(!flag[i])
      num[total++]=i;
    for(i=1;i<=total;i++){
     now=0;
     solve(0,1,i);
     if(now==0)//优化(2)
      break;
     if(i&1)
        Count+=now;
     else
        Count-=now;
    }
    Count-=total;
    Count=n-Count-1;
    if(Judge(n))
     Count--;
   }
   printf("%d\n",Count);
  }
  return 0;
 }
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值