hdu 1796容斥原理

最新推荐文章于 2020-06-20 17:09:33 发布

羁绊残阳

最新推荐文章于 2020-06-20 17:09:33 发布

阅读量359

点赞数

分类专栏： ACM_组合数学

本文链接：https://blog.csdn.net/jibancanyang/article/details/52675179

版权

ACM_组合数学专栏收录该内容

9 篇文章 0 订阅

订阅专栏

分析：

就是区间不互质的个数，但是这里注意 m <script type="math/tex" id="MathJax-Element-12">m</script>个数不一定都是素数，所以用lcm扩展增加，保证每次增加是素数。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define pr(x) cout << #x << ": " << x << "  " 
#define pl(x) cout << #x << ": " << x << endl;
typedef long long ll;

struct jibancanyang
{
    int n, m, A[10];
    long long ans;
    ll gcd(ll a, ll b) {
        if (!b) return a;
        return gcd(b, a % b);
    }

    ll lcm(ll a, ll b) {
        return a / gcd(a, b) * b;
    }

    void dfs(int cur, long long p, int cnt) {
        if (p >= n) return;
        if (cur == m) {
            if (cnt & 1) ans += (n - 1) / p;
            else ans -= (n - 1) / p;
        } else {
            dfs(cur + 1, p, cnt);
            dfs(cur + 1, lcm(p, A[cur]), cnt + 1);
        }
    }

    void fun() {
        while (~scanf("%d%d", &n, &m)) {
            for (int i = 0; i < m; ++i) {
                scanf("%d", &A[i]);
                if (!A[i]) --m, --i;
            }
            ans = 0;
            dfs(0, 1, 0);
            printf("%lld\n", ans + n - 1);
        }
    }

}ac;

int main()
{
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
#endif
    ac.fun();
    return 0;
}