裴蜀定理是个啥
- 谨以此文纪念一位 为此定理 eat shit 的 bro
原来裴蜀是外国人
给定一个长为n的序列 a1,a2,a3……,an
,求另一个包含 n 个元素的待定整数序列 X,求 S = a1 * X1 + a2 * X2 + a3 * X3+……+ an * Xn,使得 S > 0 且 S 尽可能的小。
##代码
#include<iostream>
#include<cmath>
using namespace std;
long long n,ans;
int gcd(int x,int y){
return y?gcd(y, x%y):x;
}
int main(){
cin>>n;
cin>>ans;
for(int i=2;i<=n;i++){
int a;
cin>>a;
ans=gcd(ans,abs(a));
}
cout<<ans;
return 0;
}
这个等式的答案实际上是要求所有数的最大公约数。
等式简化一下就是 ax+by=c ,那么c就是a,b的最大公约数,比如a=38,b=57
最大公约数是19 * -1 + 38 * 1 = 19。你想想:若一个数c为两个数的最大公约数,那么这个a乘以一个整数就可以表达ax+by,比如 9 和 3 的最大公约数为 3 ,那么 3 * z就可以表达任意 9 * x + 3 * y 。
这就是我朴实无华的证明。
- 最后让我们再次纪念那位bro
哦对了,刚开始的题解因为 gcd 函数的原因被卡了了一个点
见 大佬的 gcd blog 里关于能否取0的问题