包含容斥原理【题目在挑战程序设计竞赛】
题目大意是这样的,给定a1,a2,...am,求从1到n的整数中至少能整除a中一个元素的数有几个?
输入:
n =100,m=2;
a={2,3};
输出:
67
输入:
n=100,m=3;
a={2,3,7}
输出:
72
因为这个学期学的离散,包含容斥原理印象还是挺深的,对这个也挺感兴趣,不知道有多少人在看挑战程序设计竞赛这本书,前面的费马小定理,线性同余方程组啥的。。。。因为没有学习数论,根本看不懂。。。唉。。。。。。有没有个大神可以教教我数论额。。。。
#include <iostream> #include <algorithm> #include <cstdio> #include <cstring> #include <string> #include <cmath> #define maxn 500 using namespace std; int n,m; int x[maxn]; int gcd(int a,int b) { if (b>a) { int temp=a; a=b; b=temp; } if (b==0)return a; return gcd(b,a%b); } int Rong_Chi() { int res=0; for (int i=1;i<1<<m;i++)/*这个运用的是位运算,把所有子集的情况表示成集合,每个二进制中的1代表列出的子集,假设m为3 ,则正好有7中情况需要计算*/ { int num=0; for (int j=i;j!=0;j>>=1) { if (j&1)num++; } int value=1; for (int j=0;j<m;j++) { if (i>>j&1) { value=value/(gcd(value,x[j]))*x[j];/*这是求最小公倍数,假设当前是算3个子集的交集, 先算1与第一个的最小公倍数,再算这个最小公倍数与第二个,继续下去,这样算完之后就是这三个子集的最小公倍数*/ } if (value>n)break; } if (num%2==0)/*根据包含容斥原理,偶数的是需要减去的,奇数是需要加上的*/ { res-=n/value; } else { res+=n/value; } } return res; } int main() { while(cin>>n>>m) { for (int i=0;i<m;i++) { scanf("%d",&x[i]); } int now=Rong_Chi(); cout<<now<<endl; } return 0; }