1006: Redraiment猜想

时间限制: 1 Sec 内存限制: 10 MB

题目描述

redraiment在家极度无聊，于是找了张纸开始统计素数的个数。 设函数f(n)返回从1->n之间素数的个数。 redraiment发现: f(1) = 0 f(10) = 4 f(100) = 25 ... 满足g(m) = 17 * m2 / 3 - 22 * m / 3 + 5 / 3 其中m为n的位数。 他很激动，是不是自己发现了素数分布的规律了！ 请你设计一个程序，求出1->n范围内素数的个数，来验证redraiment是不是正确的，也许还可以得诺贝尔奖呢。^_^

输入

输入包括多组数据。 每组数据仅有一个整数n (1≤n≤100000000)。 输入以０结束

输出

对于每组数据输入，输出一行，为1->n(包括n)之间的素数的个数。

样例输入

110651000

样例输出

041825

提示

报告参见http://acm.zjgsu.edu.cn/Report/1006/1006.html

来源

Redraiment

值得记住的题，因为对不同地方评测机有侥幸心理，花了很多时间试着用欧拉筛进行分块打表的操作后，发现自己对于米勒-拉宾并不是很会用（当然个人比较懒，觉得知道原理就好了），后来看了别人的做法，发现考察的是容斥原理（1e8早该想到的.....）

#include <bits/stdc++.h>
#define cl(x) memset(x,0,sizeof(x))
#define _cl(a,b) memset(a,0,(b+1)<<2)
#define At(x,y) memset(x,y,sizeof x)
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define drep(i,a,b) for(int i=a;i>=b;i--)
#define ll long long
#define pb push_back
#define vct vector
#define pq priority_queue
#define all(x) x.begin(),x.end()
#define int long long
#define fi first
#define se second
//#pragma GCC optimize(2)
using namespace std;
inline int min(int a,int b) {return a<b?a:b;}
inline int max(int a,int b) {return a>b?a:b;}
//map<string,map<string,int> > mp;
//pq <ll,vct<ll>,less<ll> > now,last;
//vct<vct<int> > vt(n+1,vct<int>(m+1,0))
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3f;
const ll mod=1e9+7;
const int MOD=998244353;
const int N=1e4+10;

int n,m;
bool p[N];
int pr[N];
int cnt=0;
int ans,sqn;

void prime()
{
    memset(p,true,sizeof(p));
    p[0]=false;
    p[1]=false;
    for(int i=2;i<=N;i++)
    {
        if(p[i])pr[++cnt]=i;
        for(int j=1;j<=cnt&&pr[j]*i<=N;j++)
        {
            p[pr[j]*i]=false;
            if(i%pr[j]==0)break;
        }
    }
}


void dfs(int c,int s,int idx){
    for(int i=idx;pr[i]<=sqn;i++){
        int tmp=s*pr[i];
        if(tmp>n) return ;
        dfs(c+1,tmp,i+1);
        if(c&1) ans-=n/tmp;
        else ans+=n/tmp;
        if(c==1) ans++; 
    }
}

void solve()
{
    prime();
    while(cin>>n&&n){
        ans=n;
        sqn=sqrt(n);
        dfs(1,1,1);
        cout<<ans-1<<"\n";
    }
    return ;
}

signed main()
{
//    srand(time(NULL));
//    freopen("2022_XinShengSai_1.in", "r", stdin);
//    freopen("2022_XinShengSai_1.out", "w", stdout);
//    ios::sync_with_stdio(false); 
//    cin.tie(0); cout.tie(0);
    int t;
//    cin>>t;
    t=1;
    while(t--)
    solve();
    
    return 0;
}
/*

*/