一、题目描述
二、题目简析
首先,要理解一个定理——裴蜀定理:
若任意整数
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[k−A[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;
}
完

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

被折叠的 条评论
为什么被折叠?



