组合数学-容斥原理-求指定区间内与n互素的数的个数
2018-01-31 16:56
381 查看
求指定区间内与n互素的数的个数
给出整数n和r。求区间[1,r]中与n互素的数的个数。
去解决它的逆问题,求不与n互素的数的个数。
考虑n的所有素因子pi(i=1···k)
在[1,r]中有多少数能被pi整除呢?它就是
然而,如果我们单纯将所有结果,会得到错误答案。有些数可能被统计了多次(被好几素因子整除,如6,在计算2时,3时就重复了)
所以我们要用容斥原理来解决。我们可以用2^k的算法求出所有的pi的组合,然后计算每种组合的pi乘积,通过容斥原理来对结果进行加减处理。
关于此问题的最终实现
比如n的素因子的个数num=6,msk=000001,表示只选择了n的p[0]这个素因子
msk=011001,表示选择了n的p[0],p[3],p[4],这三个素因子。
容斥定理可以用位运算和dfs来实现
给出整数n和r。求区间[1,r]中与n互素的数的个数。
去解决它的逆问题,求不与n互素的数的个数。
考虑n的所有素因子pi(i=1···k)
在[1,r]中有多少数能被pi整除呢?它就是
然而,如果我们单纯将所有结果,会得到错误答案。有些数可能被统计了多次(被好几素因子整除,如6,在计算2时,3时就重复了)
所以我们要用容斥原理来解决。我们可以用2^k的算法求出所有的pi的组合,然后计算每种组合的pi乘积,通过容斥原理来对结果进行加减处理。
关于此问题的最终实现
int prim[maxn]; int solve(int n,int r) { int num=0; for(int i=2;i<=sqrt(n);i++) { if(n%i==0) { prim[num++]=i; while(n%i==0) n/=i; } } if(n>1) prim[num++]=n; //对n分解素因子,num是n的素因子的个数 int sum=0; for(int msk=1;msk<(1<<num);i++) {//(1<<sum)是一个长为num的二进制,msk通过从1到(1<<num)循环可以枚举这个二进制的每一种排列情况 //msk<(1<<num)是把000··的这种排列去掉了,因为得选素因子
//二进制排列中1表示选n该位置的的素因子,0表示不选
int mult=1; bits=0; for(int i=0;i<num;i++) { if((1<<i)&msk)//表示的是msk二进制的第i位数字 { bits++; mult*p[i]; } } int cur=r/mult; if(bits&1) sum+=cur; else sum-=cur; } return r-sum; }
比如n的素因子的个数num=6,msk=000001,表示只选择了n的p[0]这个素因子
msk=011001,表示选择了n的p[0],p[3],p[4],这三个素因子。
容斥定理可以用位运算和dfs来实现
相关文章推荐
- 求指定区间内与n互素的数的个数 容斥原理
- HDU 4675GCD of Sequence(大组合数取模 容斥原理 组合数学)
- 51nod 1667 概率好题 组合数学+容斥原理
- 区间DP,组合数学(Expression,HDU 5396)
- 组合数学——容斥原理和错位排列
- POJ 1091 跳蚤 [容斥原理]【组合数学】
- 区间互素(筛法欧拉函数模板+容斥原理)(1695)
- UVA.11806 Cheerleaders (组合数学 容斥原理 二进制枚举)
- HDU 1695 GCD (数论-整数和素数,组合数学-容斥原理)
- HDU 1796 How many integers can you find(组合数学-容斥原理)
- HDU5396 - Expression - 区间dp、组合数学
- 错排问题 组合数学+容斥原理
- 2014ACM/ICPC亚洲区西安站 F题 color (组合数学,容斥原理)
- HDU 1695 GCD (数论-整数和素数,组合数学-容斥原理)
- 组合数学-容斥原理专辑
- 组合数学-容斥原理
- UVALive 6469 Deranged Exams 组合数学+容斥原理
- BZOJ1853 [Scoi2010]幸运数字 [容斥原理]【组合数学】
- UVa 11806 - Cheerleaders (组合数学 容斥原理)
- 51NOD 1806 wangyurzee的树(容斥原理 + 组合数学)