扩欧,求不定方程解的个数
我们已经知道通解x0 = x +- (b/gcd(a,b))*t,那只要知道有多少个t使x在题目给定的范围中就行了
但还有y
怎么办?求交集!
分别二分确定在x取值范围内合法的t的范围和在y取值范围内合法的t的范围
然后交一下
值得一提的是,对于同一个t,由exgcd求出的x,y而找到的一组通解为
x + (b/gcd(a,b)) * t,y - (a/gcd(a,b)) * t
或
x - (b/gcd(a,b)) * t,y + (a/gcd(a,b)) * t
当a,b同号时,应注意t是互为相反数的
对此把某一组解的范围 *= -1就好了
还有一种解法是划定范围+枚举,时间上不如二分但代码更简洁
注意特判最后几组a,b == 0的数据……
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using<