BZOJ 2005 能量采集(莫比乌斯反演+分块)
2017-08-11 21:05
555 查看
思路:
跟bzoj 2190那题差不多。在对gcd(i,j)==k的贡献做统计是分别对于每个k做统计就ok了,分块这种东西不是顺手就加上了吗。。
#include<stdio.h> #include <iostream> #include<string.h> #include<math.h> #include<algorithm> #define eps 1e-8 typedef long long int lli; using namespace std; const int maxn = 1e5+10; bool isprime[maxn]; int prime[maxn]; char miu[maxn]; void moblus(){ int cnt = 0;miu[1] = 1; for(lli i = 2;i < maxn;i++){ if(!isprime[i]){ prime[cnt++] = i,miu[i] = -1;//phi[i] = i-1; } for(lli j = 0;j < cnt && i*prime[j] < maxn;j++){ lli x = prime[j]; isprime[i*x] = 1; if(i%x){ miu[i*x] = -miu[i]; //phi[i*x] = phi[i] * phi[x]; } else{ miu[i*x] = 0; //phi[i*x] = phi[i] * x; break; } } } } int ff[maxn],sum[maxn]; int main(){ lli p,q; moblus(); for(int i = 1;i <= maxn;i++){ sum[i] = sum[i-1] + miu[i]; } lli a,b; scanf("%lld%lld",&a,&b); if(a>b) swap(a,b); lli ans = 0,l,tempans = 0; for(int i = 1;i <= a;i++){ lli n = a/i,m = b/i; for(int j=1;j<=n;j=l+1){//分块加速 l= min(n/(n/j),m/(m/j)); tempans += (lli)(sum[l]-sum[j-1])*(n/j)*(m/j); } ans += (lli)((i-1)*2+1)*tempans;tempans = 0; } printf("%lld\n",ans); return 0; }
相关文章推荐
- 【bzoj2005】[Noi2010]能量采集 莫比乌斯反演
- bzoj2005: [Noi2010]能量采集 莫比乌斯反演
- 【BZOJ2005】【NOI2010】能量采集(莫比乌斯反演,容斥原理)
- BZOJ2005: [Noi2010]能量采集 莫比乌斯反演的另一种方法——nlogn筛
- BZOJ 2005: [Noi2010]能量采集 [莫比乌斯反演]
- 【NOI2010】【BZOJ2005】能量采集(莫比乌斯反演、乱搞)
- bzoj 2005 [Noi2010]能量采集 O(n)莫比乌斯反演
- 【莫比乌斯反演】BZOJ2005 [NOI2010]能量采集
- bzoj 2005 能量采集 莫比乌斯反演
- [bzoj2005][莫比乌斯反演][Noi2010]能量采集
- 【数论】【莫比乌斯反演】【线性筛】bzoj2005 [Noi2010]能量采集
- bzoj 2005: [Noi2010]能量采集【莫比乌斯反演】
- BZOJ 2005 NOI2010 能量采集 莫比乌斯反演
- bzoj 2005: [Noi2010]能量采集 筛法||欧拉||莫比乌斯
- 【BZOJ 2005】[Noi2010]能量采集 (容斥原理| 欧拉筛+ 分块)
- 【BZOJ】2005: [Noi2010]能量采集(欧拉函数+分块)
- [NOI2010][bzoj2005] 能量采集 [欧拉函数+分块前缀和优化]
- bzoj2005 [Noi2010]能量采集(Mobius反演,欧拉函数)
- BZOJ 2005 能量采集 (Mobius反演)
- 能量采集 [Bzoj 2005,Noi2010]