看了两个多小时,大概意思我觉得莫比乌斯是一个利用一个数的能被整除的一些数,通过容斥来计算出这个数的一些信息,很高级。这个题就当做是模板吧。然后把两种形式贴出来:

第一种:

第二种:

然后是容斥的函数:

然后有一些性质:

再耍两道题巩固一下。

#include <cstdio>
#include <cstring>
#include <vector>
#include <cmath>
#include <stack>
#include <cstdlib>
#include <queue>
#include <map>
#include <iostream>
#include <algorithm>
#include <bits/stdc++.h>using namespace std;
typedef long long LL;
const int MAXN = 1000000;
bool check[MAXN+10];
int prime[MAXN+10];
int mu[MAXN+10];
void Moblus()
{memset(check,false,sizeof(check));mu[1] = 1;int tot = 0;for(int i = 2; i <= MAXN; i++){if( !check[i] ){prime[tot++] = i;mu[i] = -1;}for(int j = 0; j < tot; j++){if(i * prime[j] > MAXN) break;check[i * prime[j]] = true;if( i % prime[j] == 0){mu[i * prime[j]] = 0;break;}else{mu[i * prime[j]] = -mu[i];}}}
}
int main()
{Moblus();LL T,ncas=1;scanf ("%I64d",&T);while (T--){LL a,b,c,d,k;scanf ("%I64d%I64d%I64d%I64d%I64d",&a,&b,&c,&d,&k);if (b>d) swap(b,d);if (k==0){printf ("Case %I64d: 0\n",ncas++);continue;}b/=k;d/=k;LL ans1=0;for (LL i=1;i<=b;i++){ans1+=1ll*mu[i]*(b/i)*(d/i);}LL ans2=0;for (LL i=1;i<=b;i++){ans2+=1ll*mu[i]*(b/i)*(b/i);}printf ("Case %I64d: %I64d\n",ncas++,ans1-ans2/2);}return 0;
}

BZOJ 2301

MDZZ......只能说有毒。。。

#include <cstdio>
#include <cstring>
#include <vector>
#include <cmath>
#include <stack>
#include <cstdlib>
#include <queue>
#include <map>
#include <iostream>
#include <algorithm>
#include <bits/stdc++.h>using namespace std;
typedef long long LL;
const long long MAXN = 1000000;
bool check[MAXN+10];
long long prime[MAXN+10];
long long mu[MAXN+10],sum[MAXN+10];
void Moblus()
{memset(check,false,sizeof(check));mu[1] = 1;long long tot = 0;for(long long i = 2; i <= MAXN; i++){if( !check[i] ){prime[tot++] = i;mu[i] = -1;}for(long long j = 0; j < tot; j++){if(i * prime[j] > MAXN) break;check[i * prime[j]] = true;if( i % prime[j] == 0){mu[i * prime[j]] = 0;break;}else{mu[i * prime[j]] = -mu[i];}}}sum[0]=0;for (long long i=1;i<=MAXN;i++)sum[i]=mu[i]+sum[i-1];
}
LL solve(long long t1,long long t2)
{if (t1>t2) swap(t1,t2);LL ans=0;long long j;for (long long i=1;i<=t1;i=j+1){j=min(t1/(t1/i),t2/(t2/i));ans+=1ll*(sum[j]-sum[i-1])*(t1/i)*(t2/i);}return ans;
}
int main()
{Moblus();LL T;int ncas=1;scanf ("%I64d",&T);while (T--){
//        printf ("Case %d: ",ncas++);
        LL a,b,c,d,k;scanf ("%I64d%I64d%I64d%I64d%I64d",&a,&b,&c,&d,&k);b/=k,d/=k;a=ceil(a*1.0/k);c=ceil(c*1.0/k);printf ("%I64d\n",solve(b,d)-solve(a-1,d)-solve(b,c-1)+solve(a-1,c-1));}return 0;
}

转载于:https://www.cnblogs.com/nj-czy/p/5900161.html

HDU 1695 BZOJ 2301 莫比乌斯反演相关推荐

  1. HDU 5468 Puzzled Elena 莫比乌斯反演

    题意: 给出一棵树,每个点上有权值.然后求每棵子树中与根节点互质( \(gcd(a, b) = 1\) )的节点个数. 分析: 对于一颗子树来说,设根节点的权值为\(u\), \(count_i\)表 ...

  2. BZOJ 3309 莫比乌斯反演

    题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=3309 题意:定义f(n)为n所含质因子的最大幂指数,求 $Ans=\sum _{i=1} ...

  3. HYSBZ - 2301 莫比乌斯反演

    链接 题解:直接用公式算,用容斥来减掉重复计算的部分 但是我犯了一个非常sb的错误,直接把abcd除k了,这样算a-1的时候就错了,然后举的例子刚好还没问题= = ,结果wa了好几发 //#pragm ...

  4. BZOJ 2820 莫比乌斯反演

    思路: $\Sigma_{i=1}^n\Sigma_{j=1}^mgcd(i,j)==p(p是素数)$ $\Sigma_{p是素数}^{p<=n}\Sigma_{i=1}^{\lfloor \f ...

  5. HDU 6750 Function(莫比乌斯反演)(2020百度之星初赛1)

    Function 推式子 S(n)=∑i=1n∑d∣id[gcd(d,id)==1]=∑d=1nd∑d∣i[gcd(d,id)==1]=∑d=1nd∑i=1nd[gcd(d,i)==1]=∑d=1nd ...

  6. 莫比乌斯反演(bzoj 2301: [HAOI2011]Problem b)

    2301: [HAOI2011]Problem b Time Limit: 50 Sec  Memory Limit: 256 MB Submit: 5322  Solved: 2459 [Submi ...

  7. BZOJ 2301 Problem b(莫比乌斯反演+分块优化)

    题目链接:http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=37166 题意:对于给出的n个询问,每次求有多少个数对(x,y),满 ...

  8. BZOJ 2820 YY的GCD 莫比乌斯反演

    2820: YY的GCD Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y) ...

  9. BZOJ 5330 Luogu P4607 [SDOI2018]反回文串 (莫比乌斯反演、Pollard Rho算法)

    题目链接 (BZOJ) https://www.lydsy.com/JudgeOnline/problem.php?id=5330 (Luogu) https://www.luogu.org/prob ...

最新文章

  1. 解密 Dubbo 三大中心的部署架构
  2. 重磅!2022QS世界大学排名公布!
  3. 通过RMAN克隆11g数据库(基于active database)
  4. 备份----硬盘对拷
  5. Spring.Net 能为我们做点什么
  6. OpenGL绘制一个三角形
  7. 专业对不对口重要吗_应届生求职,专业对口到底重不重要?
  8. Redis高频面试笔记:基础+缓存雪崩+哨兵+集群+Reids场景设计
  9. ReplaceAll()和Replace的实战用法
  10. OCIEnvCreate 失败,返回代码为 -1的解决方法
  11. 光栅衍射C语言程序,基于matcom的光栅衍射仿真程序
  12. malloc分配的内存空间是连续的吗
  13. mac下安装nginx及相关配置
  14. 改变你一生的五句话 (转)
  15. JavaWeb笔记 黑马程序员课程
  16. 使用Nexus在Windows中搭建Maven私服
  17. 【Unity3D】Tank大战
  18. 局部差异二进制 LDB - Local Difference Binary
  19. [生存志] 第50节 七穆掌郑国
  20. 【SpringBoot-3】Lombok使用详解

热门文章

  1. 372. Super Pow
  2. codevs 1373 射命丸文
  3. Netty中的那些坑
  4. scanf从文件中读入,printf写入到文件
  5. 1.6.3 Uploading Data with Solr Cell using Apache Tika
  6. 安卓笔记--intent传值不更新问题
  7. TabHost选项卡的实现(一):使用TabActivity实现
  8. Char.IsDigit与Char.IsNumber的区别[转]
  9. 使用阿里云ACM简化你的Spring Cloud微服务环境配置管理 1
  10. 数据结构——Java Stack 类