SOJ 3137: Simple Computing _容斥原理

Time Limit: 2000 MS Memory Limit: 65536 K

Description

Given n integers x1 x2 ... xn, you should count how many intergers from 1 to m that can be divided by at least one of them.

Input

The first line is an integer c which shows the number of cases. Each case contains 2 lines. The first line has 2 integers n and m. The second line contains n integers. (0 < n < 11, 0 < m < 2^31, 0 < xi < 100)

Output

For each case print the answer in a single line.

Sample Input

21 1022 203 4

Sample Output

510

Source

The Star of ACMClub 2007 @ HIT

//题意，给出n个数，求1到m中有多少个数(假设这个数是i)满足n个数中至少有一个数能够整除i；

//变形的求1~n中与m互质的数的个数问题！这里将这n个数直接当做是n个素因子，可以产生多种组合，不断筛选出能够整除1~m中得数的个数，容斥！注意的是，假如n中有两个数3和6，而m=18，对于3*6|18，计算值是1，但是12显然也能够被3、6、3*6整除，所以需要求出这k个数的最小公倍数（需要用long long，没注意错了几次）。

#include<stdio.h>
#include<string.h>
using namespace std;
#define LL long long
#define maxn 15

int a[maxn];

int gcd(int a,int b)
{
int c;
while(b)
{
c=b;
b=a%b;
a=c;
}
return a;
}

int make_ans(int n,int m)
{
int ans=0,i,j,flag;
LL tmp;
for(i=1;i<(1<<n);i++)
{
tmp=1,flag=0;
for(j=0;j<n&&tmp<=m;j++)
if(i&(1<<j))
flag++,tmp=(LL)(a[j]/gcd(tmp,a[j]))*tmp;
if(flag&1)
ans+=m/tmp;
else
ans-=m/tmp;
}
return ans;
}

int main()
{
int t,n,m,i;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&m);
for(i=0;i<n;i++)
scanf("%d",&a[i]);
printf("%d\n",make_ans(n,m));
}
return 0;
}