这周继上周学习内容,主要练习了扩展欧几里得定理,终于线下上课忽然有些忙碌,但没想到在这一周的最后一天又线上了,相比于线下来说,线上上课自我学习和练习的时间会多一些,但线下效率会高。
首先来说说扩展欧几里得定理的实现:
int gcd(int a, int b, int& x, int& y) {
if (b == 0) {
x = 1;
y = 0;
return a;
}
int x1, y1;
int d = gcd(b, a % b, x1, y1);
x = y1;
y = x1 - y1 * (a / b);
return d;
}
对于两个不全为0的整数a、b,必存在一组解x,y,使得ax+by==gcd(a,b);对于扩展欧几里得算法,可以应用到解不定方程、求解线性同余方程、求模的逆元。
有一个让我看了很久的题目:荒岛野人,P2421 [NOI2002] 荒岛野人 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)
问至少有多少个山洞,才能维持岛上的和平。看了一些题解,这个比较容易理解:以有野人的居住的山洞的最大编号作为初始的山洞数,开始搜索:如果此时的山洞数符合条件就直接输出,否则山洞数加一。假设几年后有两个野人处在了同一个山洞中,设i,j为野人编号,m为此时判断的山洞数,x为经过了多少年。那么就可以得出c[i]+x*p[i]≡c[j]+x*p[j](mod m),变形得(p[i]-p[j])*x≡c[j]-c[i](mod m),即线性同余方程;用扩展欧几里得算法解出x,gcd(p[i]-p[j],m)可能为负,需要带上绝对值。
//拓展欧几里得算法
inline bool judge(int m)
{
int a,b,C,x,y,d;
for(int i=1;i<=n;i++)
{
for(int j=i+1;j<=n;j++)
{
a=p[i]-p[j];b=m;C=c[j]-c[i];
if(C%gcd(a,b)!=0) continue;//无解时
exgcd(a,b,d,x,y);
x=((x*C/d)%(abs(b/d))+abs(b/d))%(abs(b/d));//绝对值
if(x<=min(l[i],l[j])) return false;//若不符题意则直接返回false
}
}
return true;//当所有野人都枚举完后再返回true