类欧几里得算法（部分）

##Preface
欧几里得算法，就是辗转相除法。
gcd(i,j)=gcd(j,i%j)

##定义
定义函数
$F(a,b,c,n)=\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor$

##推导一波
显然当$a\ge c$或者$b\ge c$时，$F(a,b,c,n)=\sum _{i=0}^n(\lfloor \frac{(a\mathrm{m}\mathrm{o}\mathrm{d}c)i+(b\mathrm{m}\mathrm{o}\mathrm{d}c)}{c}\rfloor +\lfloor \frac{a}{c}\rfloor i+\lfloor \frac{b}{c}\rfloor$
$=F(a\%c,b\%c,c,n)+\lfloor \frac{a}{c}\rfloor (n+1)n/2+\lfloor \frac{b}{c}\rfloor (n+1)$

若当a,b均小于c怎么办？
据大佬说转换成几何意义就是一条直线与x轴、y轴以及$x=n$围成的直角梯形中的整点个数

枚举纵线，
原式可化为
设$m=\lfloor \frac{an+b}{c}\rfloor$就是上面式子i=n时的数

枚举每个数
$F(a,b,c,n)=\sum _{i=0}^n\sum _{j=1}^m[\lfloor \frac{ai+b}{c}\rfloor \ge j]$

尽量化成j=0
$=\sum _{i=0}^n\sum _{j=0}^{m-1}[\lfloor \frac{ai+b}{c}\rfloor \ge j+1]$

大于等于下整除和不整除是一样的
$=\sum _{i=0}^n\sum _{j=0}^{m-1}[(\frac{ai+b}{c})\ge j+1]$

变换一下
$=\sum _{i=0}^n\sum _{j=0}^{m-1}[(ai\ge cj+c-b]$

$=\sum _{i=0}^n\sum _{j=0}^{m-1}[(i\ge \frac{(cj+c-b)}{a}]$
注意这里没有取整

大于等于，因为左边是整数
分子减1，符号变成大于
$=\sum _{i=0}^n\sum _{j=0}^{m-1}[(i>\frac{(cj+c-b-1)}{a}]$

交换主体
$=\sum _{j=0}^{m-1}\sum _{i=0}^n[(i>\frac{(cj+c-b-1)}{a}]$

然后里面的sigma可以撤掉
$=\sum _{j=0}^{m-1}n-\frac{(cj+c-b-1)}{a}$
$=nm-\sum _{j=0}^{m-1}\frac{(cj+c-b-1)}{a}$
所以$F(a,b,c,n)=nm-F(c,c-b-1,a,m-1)$

如果只看a,c二维，变换是和欧几里得算法一样的
复杂度也和欧几里得算法一样

可以递归处理
边界条件是a=0，式子很明显。

g和h本蒟蒻并不会推，可以看看大佬的BLOG

##模版

LL fd(LL a,LL b,LL c,LL n)
{
	if(a==0) return((b/c)*(n+1));
	if(a>=c||b>=c) return fd(a%c,b%c,c,n)+(a/c)*n*(n+1)/2+(b/c)*(n+1);
	LL m=(a*n+b)/c;
	LL v=fd(c,c-b-1,a,m-1);
	return n*m-v;
}