当前位置：首页 > news >正文

【dp + 裴蜀定理】P8646 [蓝桥杯 2017 省 AB] 包子凑数题解

news 2025/7/26 9:38:45

P8646 [蓝桥杯 2017 省 AB] 包子凑数题解

题目传送门

P8646 [蓝桥杯 2017 省 AB] 包子凑数

一、题目描述

小明发现包子铺有N种蒸笼，每种能放A_i个包子（无限供应）。问有多少个正整数X无法被这些蒸笼数量的组合表示出来。若无限多个则输出INF。

二、题目分析

这是一个典型的数论+动态规划问题。需要解决两个关键点：

判断无法表示的数字是否有无限多个
有限情况下统计具体无法表示的数字个数

三、问题思考

算法分析

本题需要结合数论中的裴蜀定理和动态规划来解决。

前置知识：裴蜀定理

对于任意整数a,b，存在整数x,y使得ax+by=gcd(a,b)。推广到多个数：

当所有数的gcd=1时，只有有限多个数无法表示
当gcd>1时，所有不被gcd整除的数都无法表示（无限多个）

四、动态规划思路

a. 状态表示

f[j]表示数字j能否被表示（true/false）

b. 初始化

f[0] = true（0个包子总是可以表示）

c. 状态转移

对于每个蒸笼数量a[i]，从a[i]开始更新：

if(f[j - a[i]]) f[j] = true;

d. 最终结果

统计所有f[j]==false的j的个数

五、代码实现

#include <bits/stdc++.h>
using namespace std;
const int N = 110, M = 1e4 + 10; // N是蒸笼种类上限，M是DP数组大小
int n;
int a[N];

// 计算最大公约数
int gcd(int a, int b) {
    return b == 0 ? a : gcd(b, a % b);
}

bool f[M]; // DP数组

int main() {
    ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    cin >> n;
    for (int i = 1; i <= n; i ++) cin >> a[i];
    
    // 计算所有数的gcd
    int g = a[1];
    for (int i = 2; i <= n; i ++) {
        g = gcd(g, a[i]);
    }
    
    // 根据gcd判断是否有无限解
    if (g != 1) {
        cout << "INF";
        return 0;
    }
    
    // DP过程
    f[0] = true; // 初始化
    for (int i = 1; i <= n; i ++) { // 遍历每种蒸笼
        for (int j = a[i]; j < M; j ++) { // 更新DP数组
            if(f[j - a[i]]) f[j] = true;
        }
    }

    // 统计结果
    int res = 0;
    for (int i = 1; i < M; i ++)
        if (!f[i]) res ++;
    cout << res;
    return 0;
}