容斥原理(二进制枚举)
2018-02-06 12:00
881 查看
转载内容:
分析:我们可以先转化下:用(1,b)区间与n互质的数的个数减去(1,a-1)区间与n互质的数的个数,那么现在就转化成求(1,m)区间于n互质的数的个数,如果要求的是(1,n)区间与n互质的数的个数的话,我们直接求出n的欧拉函数值即可,可是这里是行不通的!我们不妨换一种思路:就是求出(1,m)区间与n不互质的数的个数,假设为num,那么我们的答案就是:m-num!现在的关键就是:怎样用一种最快的方法求出(1,m)区间与n不互质的数的个数?方法实现:我们先求出n的质因子(因为任何一个数都可以分解成若干个质数相乘的),如何尽快地求出n的质因子呢?我们这里又涉及两个好的算法了!第一个:用于每次只能求出一个数的质因子,适用于题目中给的n的个数不是很多,但是n又特别大的;(http://www.cnblogs.com/jiangjing/archive/2013/06/03/3115399.html)第二个:一次求出1~n的所有数的质因子,适用于题目中给的n个数比较多的,但是n不是很大的。(http://www.cnblogs.com/jiangjing/archive/2013/06/01/3112035.html)本题适用第一个算法!举一组实例吧:假设m=12,n=30.
第一步:求出n的质因子:2,3,5;
第二步:(1,m)中是n的因子的倍数当然就不互质了(2,4,6,8,10)->n/2 6个,(3,6,9,12)->n/3 4个,(5,10)->n/5 2个。
如果是粗心的同学就把它们全部加起来就是:6+4+2=12个了,那你就大错特错了,里面明显出现了重复的,我们现在要处理的就是如何去掉那些重复的了!
第三步:这里就需要用到容斥原理了,公式就是:n/2+n/3+n/5-n/(2*3)-n/(2*5)-n/(3*5)+n/(2*3*5).
第四步:我们该如何实现呢?我在网上看到有几种实现方法:dfs(深搜),队列数组,位运算三种方法都可以!上述公式有一个特点:n除以奇数个数相乘的时候是加,n除以偶数个数相乘的时候是减。我这里就写下用队列数组如何实现吧:我们可以把第一个元素设为-1然后具体看代码如何实现吧!
同种类型的题目:hdu 2841 hdu1695
牛客:
链接:https://www.nowcoder.net/acm/contest/75/G
来源:牛客网
hdu 4135(容斥原理)
题意:就是让你求(a,b)区间于n互质的数的个数.分析:我们可以先转化下:用(1,b)区间与n互质的数的个数减去(1,a-1)区间与n互质的数的个数,那么现在就转化成求(1,m)区间于n互质的数的个数,如果要求的是(1,n)区间与n互质的数的个数的话,我们直接求出n的欧拉函数值即可,可是这里是行不通的!我们不妨换一种思路:就是求出(1,m)区间与n不互质的数的个数,假设为num,那么我们的答案就是:m-num!现在的关键就是:怎样用一种最快的方法求出(1,m)区间与n不互质的数的个数?方法实现:我们先求出n的质因子(因为任何一个数都可以分解成若干个质数相乘的),如何尽快地求出n的质因子呢?我们这里又涉及两个好的算法了!第一个:用于每次只能求出一个数的质因子,适用于题目中给的n的个数不是很多,但是n又特别大的;(http://www.cnblogs.com/jiangjing/archive/2013/06/03/3115399.html)第二个:一次求出1~n的所有数的质因子,适用于题目中给的n个数比较多的,但是n不是很大的。(http://www.cnblogs.com/jiangjing/archive/2013/06/01/3112035.html)本题适用第一个算法!举一组实例吧:假设m=12,n=30.
第一步:求出n的质因子:2,3,5;
第二步:(1,m)中是n的因子的倍数当然就不互质了(2,4,6,8,10)->n/2 6个,(3,6,9,12)->n/3 4个,(5,10)->n/5 2个。
如果是粗心的同学就把它们全部加起来就是:6+4+2=12个了,那你就大错特错了,里面明显出现了重复的,我们现在要处理的就是如何去掉那些重复的了!
第三步:这里就需要用到容斥原理了,公式就是:n/2+n/3+n/5-n/(2*3)-n/(2*5)-n/(3*5)+n/(2*3*5).
第四步:我们该如何实现呢?我在网上看到有几种实现方法:dfs(深搜),队列数组,位运算三种方法都可以!上述公式有一个特点:n除以奇数个数相乘的时候是加,n除以偶数个数相乘的时候是减。我这里就写下用队列数组如何实现吧:我们可以把第一个元素设为-1然后具体看代码如何实现吧!
同种类型的题目:hdu 2841 hdu1695
牛客:
链接:https://www.nowcoder.net/acm/contest/75/G
来源:牛客网
题目描述
给出一个数n,求1到n中,有多少个数不是2 5 11 13的倍数。输入描述:
本题有多组输入 每行一个数n,1<=n<=10^18.
输出描述:
每行输出输出不是2 5 11 13的倍数的数共有多少。
#include<iostream> #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> #define maxn 105 #define LL long long #define PI 3.1415926 #define e 2.7182818 using namespace std; LL n,a[4]; void solve() { LL ans=0; for(int i=1;i<(1<<4);i++) { LL muti=1,cnt=0; for(int j=0;j<4;j++) { if(i&(1<<j)) { cnt++; muti*=a[j]; } } if(cnt&1) ans+=n/muti; else ans-=n/muti; } cout << n-ans << endl; } int main() { a[0]=2;a[1]=5;a[2]=11;a[3]=13; while(cin>>n) { solve(); } return 0; }
相关文章推荐
- 容斥原理(枚举二进制)-HDU1796(重)
- UVA.11806 Cheerleaders (组合数学 容斥原理 二进制枚举)
- 二进制枚举(容斥原理)
- UVA.10325 The Lottery (组合数学 容斥原理 二进制枚举)
- Hdu 1796 How many integers can you find【二进制枚举容斥原理入门】
- UVA 10325 The Lottery (组合数学,容斥原理,二进制枚举)
- uva 10325 The Lottery (容斥原理)dfs枚举 或二进制枚举
- 枚举子集 增量构造法 位向量法 二进制法
- [kuangbin]专题一 简单搜索 D - Fliptile(二进制枚举)
- 二进制枚举
- hdu1557权利指数(二进制枚举)
- POJ3279 Fliptile题解 二进制枚举搜索法
- hdu 4770 Lights Against Dudely(二进制枚举情况)
- 二进制集合枚举子集
- poj 1426 Find The Multiple(dfs || 二进制枚举)
- 二进制枚举子集
- [最小生成树+二进制枚举]UVa1151 - Buy or Build
- POJ 1873 The Fortified Forest 二进制枚举 + 凸包 (final水题)
- JPA学习笔记---JPA实体Bean的建立---链接上一个博文:对实体Bean中属性进行操作:保存日期类型,设置字段的长度,名字,是否为空,可以声明枚举字段;可以存放二进制数据,可以存放
- C++之暴力算法突破-二进制枚举子集