首页 > 代码库 > HDU1796 How many integers can you find【容斥定理】
HDU1796 How many integers can you find【容斥定理】
题目链接:
http://acm.hdu.edu.cn/showproblem.php?
pid=1796
题目大意:
给你一个整数N。和M个整数的集合{A1、A2、…、Am}。集合内元素为非负数(包括零),求小于N的
正整数(1~N-1)中,能被M个整数的集合中随意一个元素整除的正整数个数。
比如N = 12。M = {2,3},在1~N-1中,能被2整除的数为{2,4,6。8。10},能被3整除的数为
{3。6,9}。则所求集合为{2,3,4。6,8,9,10},共7个,则答案为7。
思路:
就是求M个集合的并集。先看上边的样例。能被2整除的数集合S1为{2,4,6,8。10},能被3整除的数
集合S2为{3,6。9}。而两者交集S12为能被LCM(2,3) = 6整除的数为{6}。
则两者并集 S = S1 + S2 - S12。
依据容斥定理得:若有M个数,则可求出1~N-1中能被不同组合的公倍数整除的个数。
1~N-1能被公倍数整除的个数为 (N-1) / LCM。然后依据奇数项加,偶数项减的原则得出答案个数。
AC代码:
#include<iostream> #include<algorithm> #include<cstdio> #include<string> using namespace std; int Gcd(int a,int b) { if(a < b) int temp = a, a = b, b = temp; if(b == 0) return a; return Gcd(b,a%b); } int Lcm(int a,int b) { return a/Gcd(a,b)*b; } int N,M; int A[220],Select[220]; int Solve() { int ans = 0; for(int i = 1; i < (1 << M); ++i) { int odd = 0; for(int j = 0; j < M; ++j) { if((1<<j) & i) { Select[++odd] = j; } } int LCM = 1; for(int j = 1; j <= odd; ++j) LCM = Lcm(LCM,A[Select[j]]); if(odd & 1) ans += N/LCM; else ans -= N/LCM; } return ans; } int main() { int d; while(~scanf("%d%d",&N,&M)) { for(int i = 0; i < M; ++i) { scanf("%d",&d); if(d) A[i] = d; else i--,M--; } N--; printf("%d\n",Solve()); } return 0; }
HDU1796 How many integers can you find【容斥定理】
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。