gudt 我是好人4

众所周知，我是好人！所以不会出太难的题，题意很简单

给你n个数，问你1000000000（含1e9）以内有多少个正整数不是这n个数任意一个的倍数

最后友情提供解题代码（我真是太好人了）

void solve(int p[], int n)

{

int ans = 0;

for (int i =1; i <= 1e9; i++)

{

int fl = 0;

for (int j = 0; j < n; j++)

{

if (i % p[j] == 0)

{

fl = 1;

break;

}

}

if (fl == 0)ans++;

}

printf("%d\n", ans);

}

Input

第1行是一个整数T，表示共T组数据。 接下来是T组数据，每组数据第1行是正整数n（n<=50），接下来是n个正整数（小于等于1000），任意两数用1个空格隔开，最前数前面与最后数后面无空格

Output

输出T行，对应T组数据。（T<=10） 每行输出这样的正整数有多少个

Sample Input

342 3 5 71213854 101 143 282 538 922 946 286 681 977 892 656 907

Sample Output

228571428500000000968701719

HINT

提示：数据是随机生成的，尔等尽可随意乱搞

 原理 是 容斥 ，所以用dfs 和 最小公倍数来算，不过在这之前先把给的数中 倍数之间的数给划掉，还有出现 数字1的情况 。

倍数若大于1e9 就剪枝。

#include <iostream>#include <cstring>#include <cstdio>#include <algorithm>using namespace std;typedef long long int LL;const int maxn = 1e9;int vis[55];int num[55];int a[55];int m = 0;LL ans = 0;LL p;int n;LL gcd(LL x, LL y) {return y ? gcd(y, x%y) : x;}void dfs(LL i, LL q, LL k) {if (q > maxn)return;for (; i < m; i++){p = a[i] / gcd(a[i], q)*q;ans += k * (maxn / p);dfs(i + 1, p, -k);}}int main() {int t;scanf("%d", &t);while (t--) {m = 0;memset(vis, 0, sizeof(vis));bool one = false;scanf("%d", &n);for (int i = 0; i < n; i++) {scanf("%d", &num[i]);if (num[i] == 0) {vis[num[i]] = true;}if (num[i] == 1){one = true;}}if (one) {printf("0\n");continue;}sort(num, num + n);for (int i = 0; i < n; i++){for (int j = i + 1; j < n; j++) {if (vis[j] == 0&& vis[i] == 0 && num[j] % num[i] == 0)vis[j] = 1;}if (vis[i]==0)a[m++] = num[i];}ans = 0;dfs(0, 1, 1);printf("%lld\n", maxn - ans);}return 0;}