标签:BZOJ4805 题意 int long 杜教 while 欧拉 getchar
题意
Sol
杜教筛板子题。。
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int MAXN = 1e7 + 10;
inline int read() {
char c = getchar(); int x = 0, f = 1;
while (c < '0' || c > '9') {if (c == '-')f = -1; c = getchar();}
while (c >= '0' && c <= '9') {x = x * 10 + c - '0'; c = getchar();}
return x * f;
}
int N, prime[MAXN], tot, vis[MAXN], Lim;
LL phi[MAXN];
map<LL, LL> mp;
void Get(int N) {
vis[1] = 1; phi[1] = 1;
for(int i = 2; i <= N; i++) {
if(!vis[i]) prime[++tot] = i, phi[i] = i - 1;
for(int j = 1; j <= tot && (i * prime[j] <= N); j++) {
vis[i * prime[j]] = 1;
if(!(i % prime[j])) {
phi[i * prime[j]] = 1ll * phi[i] * prime[j]; break;
} else phi[i * prime[j]] = 1ll * phi[i] * (prime[j] - 1);
}
}
for(int i = 2; i <= N; i++) phi[i] += phi[i - 1];
}
LL sieve(LL x) {
if(x <= Lim) return phi[x];
if(mp.find(x) != mp.end()) return mp[x];
LL ans = 1ll * x * (x + 1) / 2;
for(LL i = 2, nxt; i <= x; i = nxt + 1) {
nxt = min(x, x / (x / i));
ans -= 1ll * sieve(x / i) * (nxt - i + 1);
}
return mp[x] = ans;
}
signed main() {
N = read();
Get(Lim = 1e6);
cout << sieve(N);
return 0;
}
/*
1000000000
*/
标签:BZOJ4805,题意,int,long,杜教,while,欧拉,getchar 来源: https://blog.51cto.com/u_15239936/2868684
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。