1到n里面素数的个数n<=1e11

2016沈阳网络赛的最后一题,记录一下这个模板,来自某某大神,并不懂这个的原理,先放着当模板。


//
//Created by just_sort 2016/9/12 17:01
//Copyright (c) 2016 just_sort.All Rights Reserved
//

#include <set>
#include <map>
#include <queue>
#include <stack>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define ll long long
ll f[340000],g[340000],n;
void init()
{
    ll i,j,m;
    for(m=1; m*m<=n; ++m)f[m]=n/m-1;
    for(i=1; i<=m; ++i)g[i]=i-1;
    for(i=2; i<=m; ++i)
    {
        if(g[i]==g[i-1])continue;
        ll x=min(m-1,n/i/i);
        for(j=1; j<=x; ++j)
        {
            if(i*j<m)f[j]-=f[i*j]-g[i-1];
            else f[j]-=g[n/i/j]-g[i-1];
        }
        x=i*i;
        for(j=m; j>=x; --j)g[j]-=g[j/i]-g[i-1];
    }
}
int main()
{
    while(scanf("%I64d",&n)!=EOF)
    {
        init();
        printf("%I64d\n",f[1]);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值