类欧几里得算法

来源

类欧几里德算法由洪华敦大佬在 2016 年 $\text{NOI}$ 冬令营营员交流中提出。截止 2020 年 4 月 29 日，百度百科和维基百科都尚未收录该词条 $\text{Orz}$。

前置知识

取整函数

$\lfloor \frac{a}{b}\rfloor$ 为 $\frac{a}{b}$ 向下取整，$\lceil \frac{a}{b}\rceil$ 为 $\frac{a}{b}$ 向上取整。
基与下面的类欧几里得算法的推导过程，需要掌握一些取整函数的技巧。
$$\lfloor \frac{a}{b}\rfloor ⩾c\iff a⩾bc\lceil \frac{a}{b}\rceil ⩽c\iff a⩽bc\lfloor \frac{a}{b}\rfloor <c\iff a<bc\lceil \frac{a}{b}\rceil >c\iff a>bc$$
$$\lfloor \frac{a}{b}\rfloor =\lceil \frac{a+b+1}{b}\rceil \lceil \frac{a}{b}\rceil =\lfloor \frac{a+b-1}{b}\rfloor$$

证明要略

当 $a\mathrm{m}\mathrm{o}\mathrm{d}b=0$ 时，上述式子显然都成立；当 $a\mathrm{m}\mathrm{o}\mathrm{d}b\ne 0$ 时，可以令 $a=kb+m$ ，其中 $k,m\in N^{\ast }$，$1⩽m<b$，可以轻易得到证明。

引入

类欧几里得算法主要用来计算形如 $\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor$，$\sum _{i=0}^{n}i\lfloor \frac{ai+b}{c}\rfloor$，$\sum _{i=0}^n{\lfloor \frac{ai+b}{c}\rfloor }^2$的表达式。
记 $f(a,b,c,n)=\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor$，$g(a,b,c,n)=\sum _{i=0}^{n}i\lfloor \frac{ai+b}{c}\rfloor$，$h(a,b,c,n)=\sum _{i=0}^n{\lfloor \frac{ai+b}{c}\rfloor }^2$。
其结果由递归函数得到。
由于推导过程类似欧几里得算法，时间复杂度同欧几里得算法一样也是 $\log n$ ，故称类欧几里得算法。

推导

推导 $f(a,b,c,n)$

$f(a,b,c,n)$ 的形式最为简单，也最常出现；其几何意义为直角坐标系下，直线 $y=\frac{ax+b}{c}$ 与 $x,y$ 两轴包围区域内的整点个数。
当 $a⩾c$ 或 $b⩾c$时，意味着可以将 $a,b$ 对 $c$ 取模简化问题。
当 $a⩾c$ 或 $b⩾c$时
$\begin{array}{rl}& f(a,b,c,n)=\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =\sum _{i=0}^n\lfloor \frac{\left(\lfloor \frac{a}{c}\rfloor c+a\mathrm{m}\mathrm{o}\mathrm{d}c\right)i+\left(\lfloor \frac{b}{c}\rfloor c+b\mathrm{m}\mathrm{o}\mathrm{d}c\right)}{c}\rfloor \\ & =\sum _{i=0}^n\left[\lfloor \frac{a}{c}\rfloor i+\lfloor \frac{b}{c}\rfloor +\lfloor \frac{(a\mathrm{m}\mathrm{o}\mathrm{d}c)i+(b\mathrm{m}\mathrm{o}\mathrm{d}c)}{c}\rfloor \right]\\ & =\sum _{i=0}^n\left(\lfloor \frac{a}{c}\rfloor i+\lfloor \frac{b}{c}\rfloor \right)+\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 \\ & =\frac{n(n+1)}{2}\lfloor \frac{a}{c}\rfloor +(n+1)\lfloor \frac{b}{c}\rfloor +f(a\mathrm{m}\mathrm{o}\mathrm{d}c,b\mathrm{m}\mathrm{o}\mathrm{d}c,c,n)\end{array}$
那么问题就转化为 $a<c$ 且 $b<c$ 的情况。我们观察式子，发现只有 $i$ 这一个变量，因此就只能从 $i$ 下手。在推求和式时有一个常见的技巧，就是条件与贡献的放缩与转化。具体地说，在原式 $f(a,b,c,n)=\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor$ 中， $0⩽i⩽n$ 是条件，而 $\lfloor \frac{ai+b}{c}\rfloor$ 是对总和的贡献。要加快一个求和式的计算，所有的方法都可以都可以归约为 贡献合并计算。如果发现式子的贡献难以合并，就可以将贡献与条件作转化得到另一个形式的求和式。
请阁下看好接下来的神仙操作！ $\text{Orz}$
$f(a,b,c,n)=\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{ai+b}{c}\rfloor }1$
当 $a<c$ 且 $b<c$ 时
$\begin{array}{rl}& f(a,b,c,n)=\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{ai+b}{c}\rfloor }1\\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\left[\lfloor \frac{ai+b}{c}\rfloor ⩾j\right]\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\sum _{i=0}^n\left[\lfloor \frac{ai+b}{c}\rfloor ⩾j\right]\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\sum _{i=0}^n\left[ai+b⩾cj\right]\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\sum _{i=0}^n\left[i⩾\lceil \frac{cj-b}{a}\rceil \right]\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\sum _{i=\lceil \frac{cj-b}{a}\rceil }^{n}1\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\left(n-\lceil \frac{cj-b}{a}\rceil +1\right)\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\left(n-\lfloor \frac{cj+a-b-1}{a}\rfloor +1\right)\\ & =n\lfloor \frac{an+b}{c}\rfloor -\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\left(\lfloor \frac{cj+a-b-1}{a}\rfloor -1\right)\\ & =n\lfloor \frac{an+b}{c}\rfloor -\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\lfloor \frac{cj-b-1}{a}\rfloor \\ & =n\lfloor \frac{an+b}{c}\rfloor -\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\lfloor \frac{c(j+1)-b-1}{a}\rfloor \\ & =n\lfloor \frac{an+b}{c}\rfloor -\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\lfloor \frac{cj+c-b-1}{a}\rfloor \\ & =n\lfloor \frac{an+b}{c}\rfloor -f\left(c,c-b-1,a,\lfloor \frac{an+b}{c}\rfloor -1\right)\end{array}$
边界条件为，当 $a=0$ 时，$f(a,b,c,n)=(n+1)\lfloor \frac{b}{c}\rfloor$。

推导 $g(a,b,c,n)$

当 $a⩾c$ 或 $b⩾c$ 时
$\begin{array}{rl}& g(a,b,c,n)=\sum _{i=0}^{n}i\lfloor \frac{ai+b}{c}\rfloor \\ & =\sum _{i=0}^{n}i\left[\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 \right]\\ & =\sum _{i=0}^{n}i\lfloor \frac{(a\mathrm{m}\mathrm{o}\mathrm{d}c)i+(b\mathrm{m}\mathrm{o}\mathrm{d}c)}{c}\rfloor +\sum _{i=0}^n\lfloor \frac{a}{c}\rfloor i^2+\sum _{i=0}^n\lfloor \frac{b}{c}\rfloor i\\ & =\frac{n(n+1)(2n+1)}{6}\lfloor \frac{a}{c}\rfloor +\frac{n(n+1)}{2}\lfloor \frac{b}{c}\rfloor +\sum _{i=0}^{n}i\lfloor \frac{(a\mathrm{m}\mathrm{o}\mathrm{d}c)i+(b\mathrm{m}\mathrm{o}\mathrm{d}c)}{c}\rfloor \\ & =\frac{n(n+1)(2n+1)}{6}\lfloor \frac{a}{c}\rfloor +\frac{n(n+1)}{2}\lfloor \frac{b}{c}\rfloor +g(a\mathrm{m}\mathrm{o}\mathrm{d}c,b\mathrm{m}\mathrm{o}\mathrm{d}c,c,n)\end{array}$
当 $a<c$ 且 $b<c$ 时
$\begin{array}{rl}& g(a,b,c,n)=\sum _{i=0}^{n}i\lfloor \frac{ai+b}{c}\rfloor \\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{ai+b}{c}\rfloor }i\\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }i\left[\lfloor \frac{ai+b}{c}\rfloor ⩾j\right]\\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }i\left[ai+b⩾cj\right]\\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }i\left[i⩾\lceil \frac{cj-b}{a}\rceil \right]\\ & =\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }i\left[i⩾\lfloor \frac{cj+a-b-1}{a}\rfloor \right]\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\sum _{i=0}^{n}i\left[i⩾\lfloor \frac{cj+a-b-1}{a}\rfloor \right]\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\sum _{i=\lfloor \frac{cj+a-b-1}{a}\rfloor }^{n}i\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\frac{\left(n+\lfloor \frac{cj+a-b-1}{a}\rfloor \right)\left(n-\lfloor \frac{cj+a-b-1}{a}\rfloor +1\right)}{2}\\ & =\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\frac{\left(n+\lfloor \frac{cj-b-1}{a}\rfloor +1\right)\left(n-\lfloor \frac{cj-b-1}{a}\rfloor \right)}{2}\\ & =\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\frac{\left(n+\lfloor \frac{c(j+1)-b-1}{a}\rfloor +1\right)\left(n-\lfloor \frac{c(j+1)-b-1}{a}\rfloor \right)}{2}\\ & =\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\frac{\left(n+\lfloor \frac{cj+c-b-1}{a}\rfloor +1\right)\left(n-\lfloor \frac{cj+c-b-1}{a}\rfloor \right)}{2}\\ & =\frac{1}{2}\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\left[n(n+1)-\lfloor \frac{cj+c-b-1}{a}\rfloor -{\lfloor \frac{cj+c-b-1}{a}\rfloor }^2\right]\\ & =\frac{1}{2}\left[\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}n(n+1)-\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\lfloor \frac{cj+c-b-1}{a}\rfloor -\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}{\lfloor \frac{cj+c-b-1}{a}\rfloor }^2\right]\\ & =\frac{1}{2}\left[\lfloor \frac{an+b}{c}\rfloor n(n+1)-f\left(c,c-b-1,a,\lfloor \frac{an+b}{c}\rfloor -1\right)-h\left(c,c-b-1,a,\lfloor \frac{an+b}{c}\rfloor -1\right)\right]\end{array}$
边界条件为，当 $a=0$ 时，$g(a,b,c,n)=\frac{n(n+1)}{2}\lfloor \frac{b}{c}\rfloor$ 。

推导 $h(a,b,c,n)$

当 $a⩾c$ 或 $b⩾c$ 时
$\begin{array}{rl}& h(a,b,c,n)=\sum _{i=0}^n{\lfloor \frac{ai+b}{c}\rfloor }^2\\ & =\sum _{i=0}^n{\left[\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 \right]}^2\\ & =\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 }^2+\sum _{i=0}^n{\lfloor \frac{a}{c}\rfloor }^2{i}^2+\sum _{i=0}^n{\lfloor \frac{b}{c}\rfloor }^2\\ & +2\lfloor \frac{a}{c}\rfloor \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 i+2\lfloor \frac{b}{c}\rfloor \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 +2\lfloor \frac{a}{c}\rfloor \lfloor \frac{b}{c}\rfloor \sum _{i=0}^{n}i\\ & =h(a\mathrm{m}\mathrm{o}\mathrm{d}c,b\mathrm{m}\mathrm{o}\mathrm{d}c,c,n)+2\lfloor \frac{a}{c}\rfloor g(a\mathrm{m}\mathrm{o}\mathrm{d}c,b\mathrm{m}\mathrm{o}\mathrm{d}c,c,n)+2\lfloor \frac{b}{c}\rfloor f(a\mathrm{m}\mathrm{o}\mathrm{d}c,b\mathrm{m}\mathrm{o}\mathrm{d}c,c,n)\\ & +\frac{n(n+1)(2n+1)}{6}{\lfloor \frac{a}{c}\rfloor }^2+(n+1){\lfloor \frac{b}{c}\rfloor }^2+n(n+1)\lfloor \frac{a}{c}\rfloor \lfloor \frac{b}{c}\rfloor \end{array}$
当 $a<c$ 且 $b<c$ 时
$\forall x\in N^{+}$，有如下骚操作 $x^2=2\times \frac{x(x+1)}{2}-x=2\sum _{i=1}^{x}i-x$
$\begin{array}{rl}& h(a,b,c,n)=\sum _{i=0}^n{\lfloor \frac{ai+b}{c}\rfloor }^2\\ & =\sum _{i=0}^n\left(2\sum _{j=1}^{\lfloor \frac{ai+b}{c}\rfloor }j-\lfloor \frac{ai+b}{c}\rfloor \right)\\ & =2\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{ai+b}{c}\rfloor }j-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }j\left[\lfloor \frac{ai+b}{c}\rfloor ⩾j\right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }j\left[ai+b⩾cj\right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{i=0}^n\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }j\left[i⩾\lceil \frac{cj-b}{a}\rceil \right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }j\sum _{i=0}^n\left[i⩾\lceil \frac{cj-b}{a}\rceil \right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }j\sum _{i=0}^n\left[i⩾\lfloor \frac{cj+a-b-1}{a}\rfloor \right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\left[j\left(n-\lfloor \frac{cj+a-b-1}{a}\rfloor +1\right)\right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{j=1}^{\lfloor \frac{an+b}{c}\rfloor }\left[j\left(n-\lfloor \frac{cj-b-1}{a}\rfloor \right)\right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\left[(j+1)\left(n-\lfloor \frac{c(j+1)-b-1}{a}\rfloor \right)\right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\left[(j+1)\left(n-\lfloor \frac{cj+c-b-1}{a}\rfloor \right)\right]-\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =2n\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}(j+1)-2\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}j\lfloor \frac{cj+c-b-1}{a}\rfloor -2\sum _{j=0}^{\lfloor \frac{an+b}{c}\rfloor -1}\lfloor \frac{cj+c-b-1}{a}\rfloor -\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor \\ & =\lfloor \frac{an+b}{c}\rfloor \left(\lfloor \frac{an+b}{c}\rfloor +1\right)n-2g\left(c,c-b-1,a,\lfloor \frac{an+b}{c}\rfloor -1\right)\\ & -2f\left(c,c-b-1,a,\lfloor \frac{an+b}{c}\rfloor -1\right)-f(a,b,c,n)\end{array}$
边界条件为，当 $a=0$ 时，$h(a,b,c,n)=(n+1){\lfloor \frac{b}{c}\rfloor }^2$

洛谷 P5170 【模板】类欧几里得算法 $\looparrowright$

题目描述

给定 $n,,a,b,c$，分别求 $\sum _{i=0}^n\lfloor \frac{ai+b}{c}\rfloor ,\sum _{i=0}^n{\lfloor \frac{ai+b}{c}\rfloor }^2,\sum _{i=0}^{n}i\lfloor \frac{ai+b}{c}\rfloor$ ，答案对 $998244353$ 取模。多组数据。

输入格式

第一行给出数据组数 $t$ 。
接下来 $t$ 行，每行有四个整数，分别为每组数据的 $n,a,b,c$ 。

输出格式

对于每组数据，输出一行三个整数，为三个答案对 $998244353$ 取模的结果。

输入输出样例

输入 #1

2
2 1 0 2
4 3 9 6

输出 #1

1 1 2
11 27 27

说明/提示

本题采用 $\mathrm{S}\mathrm{p}\mathrm{e}\mathrm{c}\mathrm{i}\mathrm{a}\mathrm{l}\mathrm{J}\mathrm{u}\mathrm{d}\mathrm{g}\mathrm{e}$ 。
答对所有第一问可以获得测试点 $40\%$ 的分数，答对所有第二问、第三问可以分别获得另外 $30\%$ 的分数。

测试点编号	特殊性质
$1$	$n,a,b,c⩽10$
$2\sim 3$	$n,a,b,c⩽100$
$4$	$n,a,b,c⩽10^6$
$5$	$t=1$
$6\sim 10$	无

对于所有测试点，有 $1⩽t⩽10^5,0⩽n,a,b,c⩽10^9,c\ne 0$ 。

分析

在计算的时侯，因为三个函数各有交错递归，因此可以考虑三个函数打包成结构体一起整体递归，同步计算，否则有很多项会被多次计算。
类欧几里得算法时间复杂度为 $O(\log n)$ 。

代码

#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
const ll mod=998244353;
/*-----2,6的乘法逆元-----*/
const ll inv2=499122177;
const ll inv6=166374059;
/*-----用快速幂得到------*/
/*ll quick_pow(ll a,ll b)
{
	ll res=1;
	while(b)
	{
		if(b&1) res=res*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return res;
}
int main()
{
	cout<<quick_pow(2,mod-2)<<endl;
	cout<<quick_pow(6,mod-2)<<endl;
	return 0;
}*/
struct A//将答案打包成结构体
{
	ll f;
	ll g;
	ll h;
};
A cal(ll a,ll b,ll c,ll n)//按照以上推得的公式进行递归
{
	A res,temp;
	if(!a)
	{
		res.f=(n+1)*(b/c)%mod;
		res.g=(b/c)*n%mod*(n+1)%mod*inv2%mod;
		res.h=(n+1)*(b/c)%mod*(b/c)%mod;
	}
	else if(a>=c||b>=c)
	{
		temp=cal(a%c,b%c,c,n);
		res.f=((a/c)*n%mod*(n+1)%mod*inv2%mod
		+(b/c)*(n+1)%mod+temp.f)%mod;
		res.g=((a/c)*n%mod*(n+1)%mod*(2*n+1)%mod*inv6%mod
		+(b/c)*n%mod*(n+1)%mod*inv2%mod+temp.g)%mod;
		res.h=((a/c)*(a/c)%mod*n%mod*(n+1)%mod*(2*n+1)%mod*inv6%mod
		+(b/c)*(b/c)%mod*(n+1)%mod+(a/c)*(b/c)%mod*n%mod*(n+1)%mod
		+temp.h+2*(a/c)%mod*temp.g+2*(b/c)%mod*temp.f)%mod;
	}
	else
	{
		ll m=(a*n+b)/c;
		temp=cal(c,c-b-1,a,m-1);
		res.f=(n*(m%mod)%mod-temp.f)%mod;
		res.g=(n*(n+1)%mod*(m%mod)%mod-temp.f-temp.h)%mod*inv2%mod;
		res.h=(n*(m%mod)%mod*((m+1)%mod)%mod-2*temp.g-2*temp.f-res.f)%mod;
	}
	return res;
}
int main()
{
	int _;
	for(cin>>_;_;_--)
	{
		ll n,a,b,c;
		scanf("%lld%lld%lld%lld",&n,&a,&b,&c);
		A ans=cal(a,b,c,n);
		//注意负数的问题
		printf("%lld %lld %lld\n",(ans.f+mod)%mod,(ans.h+mod)%mod,(ans.g+mod)%mod);
	}
	return 0;
}