【51nod1415】切方块（数论毒瘤）（gcd）

传送门

---

题解：

刚拿到的时候还以为是什么类欧几里得神仙题。

这道题要抄代码倒是很好抄，推导很麻烦。

首先当$gcd$不为$1$的时候，设$g=gcd(a,b,c)$，我们可以轻易转化为$g^2$个小的平面来算。

设$F(a,b,c)$表示当$gcd(a,b,c)=1$的时候的答案，先给一个结论，只想抄代码的可以拿了直接滚：

$$2\cdot F(a,b,c)=a\cdot b+b\cdot c+c\cdot a-gcd(a,b)\cdot gcd(a,c)\cdot gcd(b,c)$$

我们发现一个问题，在三者$gcd=1$的情况下，两两的$gcd$不一定为$1$。

接下来以$a=15,b=10,c=6$为例子进行探讨。（官方题解就是这个例子）

开始官方题解的神仙推导：

首先拿到主视图（用几何画板硬生生画出来的）：

底为$a$高为$b$，首先看图中紫色的部分，这部分中所有被染到的格子（只染到顶点不算）在截面中都会被切掉一次。

显然染掉的格子总数为$S_1=\frac{a\cdot b+a+b-gcd(a,b)}{2}$。

---

简单证明一下：

显然我们可以根据对称性来考虑，则$S_1=\frac{a\cdot b+k_1}{2}$，其中$k_1$为对角线穿过的格子总数。

显然$k_1$可以考虑穿过了几次横轴和纵轴，一共$a+b$次，有几次是横轴和纵轴一起穿过的，一共$gcd(a,b)$次，则$k_1=a+b-gcd(a,b)$

---

然后，我们把主视图向后推$c$个单位，把所有$c$条截面上的交线画出来：

除了这里紫色区域的斜边外，所有剩下线段在主视图中交上了的格子全部都会在平移后对答案产生贡献需要计算。

假设线段从上到下依次为第$1$条到第$c$条，其中第$i$条线段交的格子数量为$k_i$。

则答案为$F(a,b,c)=S_1+\sum _{i=2}^c{k}_i$

$2\cdot F(a,b,c)=ab-k_1+2\cdot \sum _{i=1}^c{k}_i$

仿照计算$k_1$的时候，对于所有$k_i$，我们计算他横跨了多少行，多少列，跨过多少个顶点。

对于$i$，分别设以上三个量为$r_i,c_i,p_i$

为了方便，以下全部都求$2\cdot F(a,b,c)$

顶点我们规定必须严格在紫色区域内。

表达式为$$2\cdot F(a,b,c)=ab-k_1+2c+2\cdot (\sum _{i=1}^c{r}_i+\sum _{i=1}^c{c}_i-\sum _{i=1}^c{p}_i)$$

则现在要求$\sum _{i=1}^c{r}_i,\sum _{i=1}^c{c}_i,\sum _{i=1}^c{p}_i$

现在算一下$\sum c_i$，$\sum r_i$同理。

这个地方我觉得官方题解在口胡，我说得稍微详细一点。

直接考虑每条线的横坐标坐标可以得到$c_i=\lceil \frac{(c-i+1)\cdot a}{c}\rceil$，则：

$$\begin{array}{rlrl}\sum _{i=1}^c{c}_i& =& & \sum _{i=1}^c\lceil \frac{i\cdot a}{c}\rceil \\ & =& & \sum _{i=1}^c\frac{i\cdot a}{c}+\sum _{i=1}^c\frac{c-(i\cdot a)\%c}{c}\\ & =& & \sum _{i=1}^c\frac{i\cdot a}{c}+\sum _{i=1}^c\frac{i\cdot a\%c}{c}\end{array}$$

两部分乘$2$算一下：

$$\begin{array}{rlrl}2\cdot \sum _{i=1}^c\frac{i\cdot a}{c}& =& & (c+1)\cdot a\\ 2\cdot \sum _{i=1}^c\frac{i\cdot a\%c}{c}& =& & gcd(a,c)\cdot (\frac{c}{gcd(a,c)}-1)=c-gcd(a,c)\end{array}$$

有：${\sum }_{i=1}^c{c}_i=a\cdot c+a+c-gcd(a,c)$

同理${\sum }_{i=1}^c{r}_i=b\cdot c+b+c-gcd(b,c)$

然后考虑算${\sum }_{i=1}^c{p}_i$。

算一下紫色区域内部三线交点的个数，还是考虑计算全图，加上对角线上的$gcd(a,b)-1$个点之后除以$2$。

发现一共有$gcd(a,b)\cdot gcd(a,c)-1$条线上有$gcd(b,c)$个点，其中有$gcd(a.c)-1$条线上有$1$个点不在内部。

则$$\begin{array}{rlrl}2\cdot \sum _{i=1}^c{p}_i& =& & (gcd(a,b)\cdot gcd(a,c)-1)\cdot gcd(b,c)-(gcd(a,c)-1)+gcd(a,b)-1\\ & =& & gcd(a,b)\cdot gcd(a,c)\cdot gcd(b,c)+gcd(a,b)-gcd(b,c)-gcd(a,c)\end{array}$$

把上面一大坨东西合并起来

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

cs int mod=1e9+7,inv2=mod+1>>1;

ll a,b,c;

signed main(){
	int T;std::cin>>T;
	while(T--){
		std::cin>>a>>b>>c;
		ll g=std::__gcd(a,std::__gcd(b,c));
		a/=g,b/=g,c/=g;
		ll ab=std::__gcd(a,b),bc=std::__gcd(b,c),ac=std::__gcd(a,c);
		a%=mod,b%=mod,c%=mod;
		ab%=mod,bc%=mod,ac%=mod;
		g%=mod,g=g*g%mod*inv2%mod;
		std::cout<<(a*b+b*c+a*c-ab*bc%mod*ac%mod+mod)%mod*g%mod<<"\n";
	}
	return 0;
}