191029NOIP模拟测试第一题array（扩欧）

这是我写的第一篇博客，有写的不好的地方请大家见谅并提出来，谢谢！

扩欧好题，顶！

今天模拟测试第一题爆零，没看出是一道扩欧的题，下面贴出题目：

由于n个数相对独立，所以这道题可以先对于每个数转化为解方程ax+by=c,求a、b绝对值之和的最小值，再加和得出答案。

下面贴代码：

#include<bits/stdc++.h>
using namespace std;
#define l long long

void exgcd(int a,int b,int m,l &x,l &y)//扩欧
{
	if(!b){x=m/a;y=0;}
	else
	{
		exgcd(b,a%b,m,x,y);
		swap(x,y);
		y-=a/b*x;
	}
}

int n,a,b,j;
long long x,y,p;

int main()
{
	scanf("%d%d%d",&n,&a,&b);
	int k=__gcd(a,b);
	a/=k;
	b/=k;
	if(a<b) swap(a,b);//a<b也可以，但后面就要判x<0（不判y）的情况
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&j);
		if(j%k)
		{
			printf("-1");
			return 0;
		}
		exgcd(a,b,j/k,x,y);
		if(y<0)
		{
			x-=b*((-y)/a+1);
			y+=a*((-y)/a+1);
		}
		x+=b*(y/a);
		y-=a*(y/a);//a<b的话，这两排也要改
		p+=min(abs(x)+abs(y),abs(x+b)+abs(y-a));
	}
	printf("%lld",p);
	return 0;
}

先用exgcd求出一组解(x’,y’)，那么解集是(x’+kb,y’-ka)。显然x只会取最小的正值或最大
的负值（x，y尽量接近于0）。