無標題文章

<code>

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;
}

</code>

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖或详,帶你破解...
    沈念sama閱讀 216,591評論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件量承,死亡現(xiàn)場離奇詭異,居然都是意外死亡问麸,警方通過查閱死者的電腦和手機疗杉,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,448評論 3 392
  • 文/潘曉璐 我一進店門阵谚,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事梢什∥拍担” “怎么了?”我有些...
    開封第一講書人閱讀 162,823評論 0 353
  • 文/不壞的土叔 我叫張陵绳矩,是天一觀的道長。 經(jīng)常有香客問我玖翅,道長翼馆,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,204評論 1 292
  • 正文 為了忘掉前任金度,我火速辦了婚禮应媚,結果婚禮上,老公的妹妹穿的比我還像新娘猜极。我一直安慰自己中姜,他們只是感情好,可當我...
    茶點故事閱讀 67,228評論 6 388
  • 文/花漫 我一把揭開白布跟伏。 她就那樣靜靜地躺著丢胚,像睡著了一般。 火紅的嫁衣襯著肌膚如雪受扳。 梳的紋絲不亂的頭發(fā)上携龟,一...
    開封第一講書人閱讀 51,190評論 1 299
  • 那天,我揣著相機與錄音勘高,去河邊找鬼峡蟋。 笑死,一個胖子當著我的面吹牛华望,可吹牛的內容都是我干的蕊蝗。 我是一名探鬼主播,決...
    沈念sama閱讀 40,078評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼赖舟,長吁一口氣:“原來是場噩夢啊……” “哼蓬戚!你這毒婦竟也來了?” 一聲冷哼從身側響起建蹄,我...
    開封第一講書人閱讀 38,923評論 0 274
  • 序言:老撾萬榮一對情侶失蹤碌更,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后洞慎,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體痛单,經(jīng)...
    沈念sama閱讀 45,334評論 1 310
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 37,550評論 2 333
  • 正文 我和宋清朗相戀三年劲腿,在試婚紗的時候發(fā)現(xiàn)自己被綠了旭绒。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,727評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖挥吵,靈堂內的尸體忽然破棺而出重父,到底是詐尸還是另有隱情,我是刑警寧澤忽匈,帶...
    沈念sama閱讀 35,428評論 5 343
  • 正文 年R本政府宣布房午,位于F島的核電站,受9級特大地震影響丹允,放射性物質發(fā)生泄漏郭厌。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 41,022評論 3 326
  • 文/蒙蒙 一雕蔽、第九天 我趴在偏房一處隱蔽的房頂上張望折柠。 院中可真熱鬧,春花似錦批狐、人聲如沸扇售。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,672評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽承冰。三九已至,卻和暖如春髓废,著一層夾襖步出監(jiān)牢的瞬間巷懈,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,826評論 1 269
  • 我被黑心中介騙來泰國打工慌洪, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留顶燕,地道東北人。 一個月前我還...
    沈念sama閱讀 47,734評論 2 368
  • 正文 我出身青樓冈爹,卻偏偏與公主長得像涌攻,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子频伤,可洞房花燭夜當晚...
    茶點故事閱讀 44,619評論 2 354

推薦閱讀更多精彩內容

  • 祭祀WO菛僾_3930閱讀 195評論 0 1
  • 背景 一年多以前我在知乎上答了有關LeetCode的問題, 分享了一些自己做題目的經(jīng)驗恳谎。 張土汪:刷leetcod...
    土汪閱讀 12,743評論 0 33
  • 本文主要摘抄整理自Rickjin的《LDA數(shù)學八卦》 統(tǒng)計模擬中一個很重要的問題就是給定一個概率分布$p(x)$,...
    iV0id閱讀 434評論 0 1
  • 轉至元數(shù)據(jù)結尾創(chuàng)建: 董瀟偉,最新修改于: 十二月 23, 2016 轉至元數(shù)據(jù)起始第一章:isa和Class一....
    40c0490e5268閱讀 1,709評論 0 9
  • 想寫這個話題很久了,但是真的開始寫的時候岸更,又突然會覺得像是身處大雪過后的曠野鸵膏,四周都是白茫茫的一片,無從下手怎炊。 最...
    原同學閱讀 239評論 2 2