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

本文介绍了如何利用裴蜀定理分析整数组合问题,讨论了当gcd(a,b)=1时的有限个数无法表达,以及gcd(a,b)>1时的无限个无法表达。通过动态规划求解在给定范围内可以表示的元素数量,并统计不能表示的元素个数。

一、题目描述

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

二、题目简析

首先,要理解一个定理——裴蜀定理
若任意整数 a a a b b b,且有 m = gcd ( a , b ) m = \text{gcd}(a, b) m=gcd(a,b),对任意整数 x x x y y y a x + b y = c ax+by=c ax+by=c,则 m   ∣   c m~|~c m  c
由该定理,我们知道 a x + b y ax+by ax+by 一定是 gcd ( a , b ) \text{gcd}(a, b) gcd(a,b) 的倍数。

利用裴蜀定理,再来看本题,可以得到两个结论:

  • 1、若 gcd ( a , b ) = = 1 \text{gcd}(a, b) == 1 gcd(a,b)==1,则无法表达的数目是有限个。
  • 2、若 gcd ( a , b ) > 1 \text{gcd}(a, b) > 1 gcd(a,b)>1,则无法表达的数目是无限个。

先判断所有元素的 gcd 是否等于1,若大于,则输出 INF;若等于,则在一定范围内(如 1e5)求出所有可以表示的元素,再统计不能表示的元素个数。我们可以采用动态规划来统计:
d p [ k ] = dp[k]= dp[k]= 能否表示该数true or false

d p [ k ] = d p [ k ]   o r   d p [ k − A [ i ] ] dp[k] = dp[k]~or~dp[k - A[i]] dp[k]=dp[k] or dp[kA[i]]

dp[0] = true;
for (int i = 0; i < n; i++)
{
	for (int j = A[i]; j <= N; j++)
		dp[j] = dp[j] | dp[j - A[i]];
}

三、AC代码

#include <bits/stdc++.h>

using namespace std;

const int N = 1e5;
const int MAX = N + 3;

int A[MAX], n;
bool dp[MAX];

int gcd(int a, int b)
{
	if (b == 0)    return a;
	return gcd(b, a % b);
}

bool check(void)
{
	int res = A[0];
	for (int i = 1; i < n; i++)
		res = gcd(res, A[i]);
	return res == 1;
}

int main()
{
	#ifdef LOCAL
	freopen("test.in", "r", stdin);
	#endif
	
	cin >> n;
	for (int i = 0; i < n; i++)
		cin >> A[i];
		
	if (check())
	{
		dp[0] = true;
		for (int i = 0; i < n; i++)
		{
			for (int j = A[i]; j <= N; j++)
				dp[j] = dp[j] | dp[j - A[i]];
		}
		int ans = 0;
		for (int i = 1; i <= N; i++)
			if (!dp[i])
				ans++;
		cout << ans << endl;
	}
	else
		cout << "INF" << endl;	
	
	return 0;	
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值