無標(biāo)題文章


#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
const int N = 5e6 + 2;
bool np[N];
int prime[N], pi[N];

int getprime() {
    int cnt = 0;
    np[0] = np[1] = true;
    pi[0] = pi[1] = 0;
    for(int i = 2; i < N; ++i) {
        if(!np[i]) prime[++cnt] = i;
        pi[i] = cnt;
        for(int j = 1; j <= cnt && i * prime[j] < N; ++j) {
            np[i * prime[j]] = true;
            if(i % prime[j] == 0)   break;
        }
    }
    return cnt;
}
const int M = 7;
const int PM = 2 * 3 * 5 * 7 * 11 * 13 * 17;
int phi[PM + 1][M + 1], sz[M + 1];
void init() {
    getprime();
    sz[0] = 1;
    for(int i = 0; i <= PM; ++i)  phi[i][0] = i;
    for(int i = 1; i <= M; ++i) {
        sz[i] = prime[i] * sz[i - 1];
        for(int j = 1; j <= PM; ++j) {
            phi[j][i] = phi[j][i - 1] - phi[j / prime[i]][i - 1];
        }
    }
}
int sqrt2(LL x) {
    LL r = (LL)sqrt(x - 0.1);
    while(r * r <= x)   ++r;
    return int(r - 1);
}
int sqrt3(LL x) {
    LL r = (LL)cbrt(x - 0.1);
    while(r * r * r <= x)   ++r;
    return int(r - 1);
}
LL getphi(LL x, int s) {
    if(s == 0)  return x;
    if(s <= M)  return phi[x % sz[s]][s] + (x / sz[s]) * phi[sz[s]][s];
    if(x <= prime[s]*prime[s])   return pi[x] - s + 1;
    if(x <= prime[s]*prime[s]*prime[s] && x < N) {
        int s2x = pi[sqrt2(x)];
        LL ans = pi[x] - (s2x + s - 2) * (s2x - s + 1) / 2;
        for(int i = s + 1; i <= s2x; ++i) {
            ans += pi[x / prime[i]];
        }
        return ans;
    }
    return getphi(x, s - 1) - getphi(x / prime[s], s - 1);
}
LL getpi(LL x) {
    if(x < N)   return pi[x];
    LL ans = getphi(x, pi[sqrt3(x)]) + pi[sqrt3(x)] - 1;
    for(int i = pi[sqrt3(x)] + 1, ed = pi[sqrt2(x)]; i <= ed; ++i) {
        ans -= getpi(x / prime[i]) - i + 1;
    }
    return ans;
}
LL lehmer_pi(LL x) {
    if(x < N)   return pi[x];
    int a = (int)lehmer_pi(sqrt2(sqrt2(x)));
    int b = (int)lehmer_pi(sqrt2(x));
    int c = (int)lehmer_pi(sqrt3(x));
    LL sum = getphi(x, a) + LL(b + a - 2) * (b - a + 1) / 2;
    for (int i = a + 1; i <= b; i++) {
        LL w = x / prime[i];
        sum -= lehmer_pi(w);
        if (i > c) continue;
        LL lim = lehmer_pi(sqrt2(w));
        for (int j = i; j <= lim; j++) {
            sum -= lehmer_pi(w / prime[j]) - (j - 1);
        }
    }
    return sum;
}

int main() {
    ios::sync_with_stdio(false);
    init();
    LL n;
    while(~scanf("%lld", &n)) {
        printf("%lld\n", lehmer_pi(n));
    }
    return 0;
}
'''
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末僵腺,一起剝皮案震驚了整個濱河市玛瘸,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌版保,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,651評論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件亏镰,死亡現(xiàn)場離奇詭異摧莽,居然都是意外死亡旺遮,警方通過查閱死者的電腦和手機医寿,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,468評論 3 392
  • 文/潘曉璐 我一進(jìn)店門栏赴,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人靖秩,你說我怎么就攤上這事须眷。” “怎么了沟突?”我有些...
    開封第一講書人閱讀 162,931評論 0 353
  • 文/不壞的土叔 我叫張陵花颗,是天一觀的道長。 經(jīng)常有香客問我惠拭,道長捎稚,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,218評論 1 292
  • 正文 為了忘掉前任求橄,我火速辦了婚禮,結(jié)果婚禮上葡公,老公的妹妹穿的比我還像新娘罐农。我一直安慰自己,他們只是感情好催什,可當(dāng)我...
    茶點故事閱讀 67,234評論 6 388
  • 文/花漫 我一把揭開白布涵亏。 她就那樣靜靜地躺著,像睡著了一般蒲凶。 火紅的嫁衣襯著肌膚如雪气筋。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,198評論 1 299
  • 那天旋圆,我揣著相機與錄音宠默,去河邊找鬼。 笑死灵巧,一個胖子當(dāng)著我的面吹牛搀矫,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播刻肄,決...
    沈念sama閱讀 40,084評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼瓤球,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了敏弃?” 一聲冷哼從身側(cè)響起卦羡,我...
    開封第一講書人閱讀 38,926評論 0 274
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后绿饵,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體欠肾,經(jīng)...
    沈念sama閱讀 45,341評論 1 311
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 37,563評論 2 333
  • 正文 我和宋清朗相戀三年蝴罪,在試婚紗的時候發(fā)現(xiàn)自己被綠了董济。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,731評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡要门,死狀恐怖虏肾,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情欢搜,我是刑警寧澤封豪,帶...
    沈念sama閱讀 35,430評論 5 343
  • 正文 年R本政府宣布,位于F島的核電站炒瘟,受9級特大地震影響吹埠,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜疮装,卻給世界環(huán)境...
    茶點故事閱讀 41,036評論 3 326
  • 文/蒙蒙 一缘琅、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧廓推,春花似錦刷袍、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,676評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至专缠,卻和暖如春雷酪,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背涝婉。 一陣腳步聲響...
    開封第一講書人閱讀 32,829評論 1 269
  • 我被黑心中介騙來泰國打工哥力, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人嘁圈。 一個月前我還...
    沈念sama閱讀 47,743評論 2 368
  • 正文 我出身青樓省骂,卻偏偏與公主長得像,于是被迫代替她去往敵國和親最住。 傳聞我的和親對象是個殘疾皇子钞澳,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 44,629評論 2 354

推薦閱讀更多精彩內(nèi)容

  • include<bits/stdc++.h> using namespace std; typedef long ...
    祭祀WO菛僾_3930閱讀 299評論 0 0
  • 背景 一年多以前我在知乎上答了有關(guān)LeetCode的問題, 分享了一些自己做題目的經(jīng)驗。 張土汪:刷leetcod...
    土汪閱讀 12,743評論 0 33
  • 本文主要摘抄整理自Rickjin的《LDA數(shù)學(xué)八卦》 統(tǒng)計模擬中一個很重要的問題就是給定一個概率分布$p(x)$,...
    iV0id閱讀 434評論 0 1
  • 轉(zhuǎn)至元數(shù)據(jù)結(jié)尾創(chuàng)建: 董瀟偉涨缚,最新修改于: 十二月 23, 2016 轉(zhuǎn)至元數(shù)據(jù)起始第一章:isa和Class一....
    40c0490e5268閱讀 1,709評論 0 9
  • 多少人欣賞你的智慧轧粟,喜歡你漂浪亮策治,善良,卻沒有幾人懂你曲折的命運兰吟,你的過去通惫,現(xiàn)在和未來。
    甜蜜蜜_b52a閱讀 159評論 0 0