洛谷 P4549 【模板】裴蜀定理

最新推荐文章于 2023-12-02 13:52:50 发布

qq_38232157

最新推荐文章于 2023-12-02 13:52:50 发布

阅读量191

点赞数

分类专栏：洛谷数论

本文链接：https://blog.csdn.net/qq_38232157/article/details/108735146

版权

洛谷同时被 2 个专栏收录

93 篇文章 0 订阅

订阅专栏

数论

49 篇文章 0 订阅

订阅专栏

裴蜀定理指出，方程ax+by+cz+…+nm=f有解当且仅当f是所有系数的最大公约数的倍数。在给定序列{an}的情况下，寻找序列{bn}使a1b1+a2b2+…+an*bn达到最小正数值，最小值恰好是序列{an}中所有项的最大公约数。

摘要由CSDN通过智能技术生成

裴(pei)蜀定理
引理：对于给定的正整数a，b，方程ax+by=c有解的充要条件为c是gcd（a，b）的整数倍
裴(pei)蜀定理：
方程ax+by+cz+…+nm=f（其中a,b,c…n,f为整数）有解的充要条件是f为gcd（a,b,c,…,n）的整数倍
定理的应用：
给定一个序列{an}，求一个整数序列{bn}使得a1b1+a2b2+…+an*bn值最小（要求最小值为正数），求这个最小值
解：根据裴蜀定理的推广，原式最小值即为gcd(a1,a2…an)

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int MaxN = 22;
int n;
int a[MaxN], x[MaxN];

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

void solve()
{
	int ans = 0;
	for(int i = 1; i <= n; ++i)
	{
		ans = gcd(ans, a[i]);
	}
	printf("%d\n", ans);
}

int main()
{
	scanf("%d", &n);
	for(int i = 1; i <= n; ++i)
	{
		scanf("%d", &a[i]);
		if(a[i] < 0)	a[i] *= -1;
	}
	solve();
	return 0;
}

/*
2
4059 -1782
*/

/*
99
*/