NOIP 2017题解（更新ing）

D1T1：小凯的疑惑

题目：求一个最大的正整数c，使得ax+by=c（其中a,b为互质的两个正整数）没有非负正整数解。

正解：（想要直接数学推导的就去找数竞大佬吧。。。下面说说考试时怎么办——“一猜想+两验证”）

①打表找规律（不急，后面有严格证明）

观察不为-1的所有元素可不完全归纳所求最大的c即ab-a-b。

下面是赛后打表验证的代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll exgcd(ll a,ll b,ll &x,ll &y) {
	if (!b) {
		x=1,y=0;
		return a;
	}
	ll d=exgcd(b,a%b,y,x);
	y-=(a/b)*x;
	return d;
}
inline bool judge(ll a,ll b,ll c) {
	ll x,y,d=exgcd(a,b,x,y);
	if (c%d) return false;
	x*=c/d,y*=c/d;
	ll tmp=b/d;
	x=(x%tmp+tmp)%tmp;
	y=(c-a*x)/b;
	return y>=0&&x>=0;
}
inline ll gcd(ll a,ll b) {
	return !b?a:gcd(b,a%b);
}
int main() {
	for (int i=1;i<=10;++i) {
		for (int j=1;j<=10;++j) {
			if (i==1||j==1||gcd(i,j)>1) {
				printf("%-3d ",-1);
				continue;
			}
			else {
				for (int t=1000;t;--t)
					if (!judge(i,j,t)&&t) {
						printf("%-3d ",t);
						break;
					}
			}
		}
		puts("");
	}
	return 0;
}

②证明c=ab-a-b原方程没有有非负正整数解：

假设有非负正整数解：

令ax+by=ab-a-b

得x=b-1-(1+y)b/a

已知a,b互质，所以a|b不成立，又因为x为正整数所以a|(1+y)

设y+1=ka（k≥1），x=b-1-kb=(1-k)b-1，因为k≥1，所以x<0，与假设矛盾。

所以假设不成立，c=ab-a-b原方程没有非负正整数解。

③证明ab-a-b是最大的c使原方程没有非负整数解：

假设它不是，设ax+by=ab-a-b+t（t为正整数）。

设t=ap+bq，其中p,q为整数（此处不要求为正）

根据exgcd算法的性质（这一点理解后面就容易了）：设d=gcd(a,b)，由于通解为x=x0+k*(b/d)，y=y0+k*(a/d)），所以一定存在一组解(x,y)其中y>0且-b<x≤0

此处ap+bq=t（a,b互质肯定有解），令q>0且-b<p≤0。

现在观察原方程ax+by=ab-a-b+t=(b-1+p)a+(q-1)b，

因为q>0且-b<p≤0，所以b-1+p≥0且q-1≥0。

所以x=b-1+p，y=q-1即为原方程的一组非负整数解，无论t多大都成立。

证毕，下面给出正解代码（自嘲而已）。

#include<cstdio>
#include<iostream>
using namespace std;
typedef long long ll;
ll a,b;
int main() {
//	freopen("math.in","r",stdin);
//	freopen("math.out","w",stdout);
	cin>>a>>b;
	cout<<a*b-a-b<<endl;
	return 0;
}

（持续更新ing）