蓝桥杯官网题目:2.包子凑数

链接:题目点这里

首先要知道一个数学定理裴蜀定理,还有完全背包的基本运用,这里只介绍前者

也可以看一下我的个人理解,我是第一次听说这个定理,理解可能有误差。

  • 假设gcd(a,b)=d,gcd是最大公约数的意思。即a,b的最大公约数是d
  • ax+by=m(x,y是任意整数,可正可负)
  • 对于所有的m,一定会被d整除,即m%d=0

可以尝试画出ax+by=z的三维立体图,很显然是一个空间平面。

z是一个未知数,它的取值有无数个,如果在三维坐标系中看,那么是所有的z(z可以被gcd(a,b)整除)。

换句话说,ax+by表示的数的集合是{实数中所有的可以被gcd(a,b)整除的数}。

  • 以下考虑a,b互质

  • a,b如果是互质的,那么
    g c d ( a , b ) = 1 gcd(a,b)=1 gcd(a,b)=1

  • 那么ax+by可以构成所有的整数:

    {ax+by | x ∈ x\in x∈N, y ∈ y\in y∈N, a x + b y ∈ ax+by\in ax+by∈N}

  • 当x,y都是正整数的时候,包括0。ax+by不能表示的最小数是:
    ( a − 1 ) ∗ ( b − 1 ) − 1 (a-1)*(b-1)-1 (a−1)∗(b−1)−1

  • 也就是说,ax+by可以表示大于(a-1)*(b-1)-1的所有正整数

回到题目

看懂了上面的数学基础应该这个题就比较清晰了。

那我就直接上代码了,不懂的评论区留言

cpp 复制代码
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1e4+1;  //只需要比ax+by的最小值大1就可以了。
const int M = 110;
long long int dp[N];  //dp[i]=j表示取i种包子的方案数是j
//dp[i]=dp[i-a[1]]+dp[i-a[2]]+dp[i-a[3]]+....+dp[i-a[n]]; 
//即,取i种包子的方案数等于取[i-a[1]]包子的方案数+[i-a[2]]包子的方案数+...+[i-a[n]]包子的方案数
int a[M];
int sum = 0;

int gcd(int a, int b)  //辗转相除法
{
    if (a < b)swap(a, b);  //大的数是被除数
    int r = a % b;   //余数
    if (r == 0)return b;
    else
    {
        a = b;
        b = r;
    }
    return gcd(a, b);
}

int main()
{
    int n;
    cin >> n;
    int k;
    for (int i = 1; i <= n; i++)  //找最大公约数
    {
        cin >> a[i];
        if (i == 1)k = a[i];
        else
            k = gcd(k, a[i]);
    }

    if (k != 1)cout << "INF";   //如果最大公约数不是1,那么就会有无穷个数取不到。

    else  //说明最大公约数是1,那么ax+by的值是所有1的倍数,ax+by此时表示整数集(所有整数)
    {
        sort(a + 1, a + n + 1);
        dp[0] = 1;  //取0种包子的方案数是1,即不取,这个必须要有,是很重要的边界初始化
        for (int i = 1; i <= N; i++)  //枚举包子的个数
        {
            for (int j = 1; j <= n; j++)  //然后更新当前包子个数的方案数
            {
                if (i - a[j] < 0)
                    break;
                dp[i] += dp[i - a[j]];
            }
            if (dp[i] == 0)  //退出内嵌for循环判断i个包子的方案数是否是0,如果是,说明这个数不能被构成
                sum++;
        }
        cout << sum << endl;
    }
    return 0;
}

对裴蜀定理有兴趣的可以关注我这篇博客,我会从cf和leetcode等网站更新相关内容,将会以链接形式帖在本篇下面

相关推荐
惊鸿.Jh16 分钟前
【滑动窗口】3254. 长度为 K 的子数组的能量值 I
数据结构·算法·leetcode
明灯L16 分钟前
《函数基础与内存机制深度剖析:从 return 语句到各类经典编程题详解》
经验分享·python·算法·链表·经典例题
碳基学AI22 分钟前
哈尔滨工业大学DeepSeek公开课:探索大模型原理、技术与应用从GPT到DeepSeek|附视频与讲义免费下载方法
大数据·人工智能·python·gpt·算法·语言模型·集成学习
补三补四26 分钟前
机器学习-聚类分析算法
人工智能·深度学习·算法·机器学习
独好紫罗兰36 分钟前
洛谷题单3-P5718 【深基4.例2】找最小值-python-流程图重构
开发语言·python·算法
你今天刷题了吗?37 分钟前
蓝桥杯练习:翻转二叉树
蓝桥杯
正脉科工 CAE仿真1 小时前
基于ANSYS 概率设计和APDL编程的结构可靠性设计分析
人工智能·python·算法
Dovis(誓平步青云)1 小时前
【数据结构】排序算法(中篇)·处理大数据的精妙
c语言·数据结构·算法·排序算法·学习方法
2401_872945091 小时前
【补题】Xi‘an Invitational 2023 E. Merge the Rectangles
算法
暮雨哀尘2 小时前
微信小程序开发:开发实践
开发语言·算法·微信小程序·小程序·notepad++·性能·技术选型