[CF1043F] Make It One

最新推荐文章于 2022-09-19 17:20:02 发布

B丶atty

最新推荐文章于 2022-09-19 17:20:02 发布

阅读量328

点赞数

分类专栏： CF版刷DP 文章标签：动态规划算法

本文链接：https://blog.csdn.net/qq_42527676/article/details/106082837

版权

CF版刷DP 专栏收录该内容

9 篇文章 0 订阅

订阅专栏

Make It One

题意：
题意很简单，给你一个长度为 $n$ 的数组 $a$ ，让你求最小的子集，使得他们的 $g c d = 1$ ，输出最小子集的长度，没有则输出-1。

分析：
这类题还是很经典的，我们先观察它的数据范围是 $\times 10^5$ ，所以我们素数最多前7个就可以了。因为前7个素数相乘就已经超过了那个范围。
我们设 $d p [i] [j]$ 表示取 $i$ 个数， $g c d$ 为 $j$ 的集合最小数。
然后考虑状态转移方程。
$C_{cnt_j}^{i} - \sum_{j|k}dp[i][k]$

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int maxn = 3e5 + 5;
const int MOD = 1e9 + 7;

ll fac[maxn], inv[maxn];

ll qpow(ll a, ll b, ll mod = MOD) {
    ll res = 1;
    while (b) {
        if (b & 1) res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}

ll C(int n, int m) {
    if (n < m) return 0;
    return fac[n] * inv[m] % MOD * inv[n - m] % MOD;
}

void init() {
    fac[0] = 1;
    for (int i = 1; i <= maxn - 5; i++) fac[i] = fac[i - 1] * i % MOD;
    inv[maxn - 5] = qpow(fac[maxn - 5], MOD - 2);
    for (int i = maxn - 6; i >= 0; i--) inv[i] = inv[i + 1] * (i + 1) % MOD;
}

int n, cnt[maxn], dp[10][maxn];

int main() {
    init();
    scanf("%d", &n);
    for (int i = 1, x; i <= n; i++) {
        scanf("%d", &x);
        cnt[x]++;
        if (x == 1) return puts("1"), 0;
    }
    for (int i = 1; i <= maxn - 5; i++) {
        for (int j = i + i; j <= maxn - 5; j += i) {
            cnt[i] += cnt[j]; // 倍数
        }
    }
    for (int i = 2; i < 10; i++) {
        for (int j = maxn - 5; j; j--) {
            dp[i][j] = C(cnt[j], i);
            for (int k = j + j; k <= maxn - 5; k += j) {
                dp[i][j] = (dp[i][j] - dp[i][k] + MOD) % MOD;
            }
        }
        if (dp[i][1] > 0) return printf("%d\n", i), 0;
    }
    puts("-1");
    return 0;
}