JZOJ6025. 【GDOI2019模拟2019.2.16】加农炮

题意：

数据范围：

Analysis：

不难看出我们就是要求斜率第$K$大，一般这种题都用二分解决。
斜率可以用最简分数表出，我们采用$Stern-BrocotTree$来解决这道题。
这是一棵用来表示所有既约分数的树，且该树满足二叉搜索树的性质，中序遍历是有序的序列。
其形状大约如下：

对于中序遍历这棵树得到的序列，若$\frac{a}{b},\frac{c}{d}$中间夹着恰好一个数，那么该数为：$\frac{a+c}{b+d}$，我们可以推得这棵树上每个点的分数。
我们显然就可以在这棵树上走，每走到一个节点判断有多少条直线斜率大于当前斜率即可，这个化一下不等式就可以得到一个能用类欧解决的式子。
详细来说是：
我们已知$c,d$，要求有多少个$\frac{a}{b}>\frac{c}{d}$。
也就是$a\ast d>c\ast b$，$a>\frac{c\ast d}{b}$。
我们枚举$c$，就能得到该式子了，注意式子中的部分可能超过$m$，因此要算出什么时候超过$m$来算。
这棵树的深度并不是$\log n$的，因此还要优化，我们不妨每一次二分走多少步到下一个转折点，有一个性质就是转折点不会超过$\log$个(我不会证)。
这样就做完了，然而是真的难打。
复杂度:$O(T\ast {\log }^{3}n)$

Code：

# include<cstdio>
# include<cstring>
# include<algorithm>
using namespace std;
typedef long long ll;
struct P
{
	ll x,y;
}ans,fir,sec,nx,z;
int n,m,T;
ll k;
inline ll S(ll a,ll b,ll c,ll n)
{
	if (!a) return 0;
	if (a >= c || b >= c) return (a / c) * n * (n + 1) / 2 + (b / c) * (n + 1) + S(a % c,b % c,c,n);
	ll m = (a * n + b) / c;
	return n * m - S(c,c - b - 1,a,m - 1);
}
inline ll C(ll p,ll q)
{
	ll lim = (m - 1) * q / p;
	lim = min(lim,n - 1ll);
	if (lim < 0) return 0;
	return (lim + 1) * (m - 1) - S(p,0,q,lim);
}
int main()
{
//	freopen("a.in","r",stdin);
//	freopen("a.out","w",stdout);
	freopen("cannon.in","r",stdin);
	freopen("cannon.out","w",stdout);
	scanf("%d",&T);
	while (T--)
	{
		scanf("%d%d%lld",&n,&m,&k);
		if (k <= m - 1) { printf("%d %d\n",1,k + 1); continue; }
		if (k > (ll)(m - 1) * n) { printf("%d %d\n",k - (ll)(m - 1) * n + 1,1); continue; }
		fir = (P){1,0},sec = (P){0,1},nx = ans = (P){1,1};
		while (1)
		{
			if (nx.x >= n || nx.y >= m) break;
			int l = 1,r = 1e6,ret = r;
			if (C(nx.y,nx.x) < k)
			{
				while (l <= r)
				{
					int mid = (l + r) >> 1;
					z.x = fir.x * mid + nx.x,z.y = fir.y * mid + nx.y;
					if (C(z.y,z.x) < k) ans = z,l = mid + 1;
					else r = mid - 1,ret = mid;
				}
				nx.x += fir.x * ret,nx.y += fir.y * ret;
				sec.x += fir.x * ret,sec.y += fir.y * ret;
			}else
			{
				while (l <= r)
				{
					int mid = (l + r) >> 1;
					z.x = sec.x * mid + nx.x,z.y = sec.y * mid + nx.y;
					if (C(z.y,z.x) < k) ans = z,r = mid - 1,ret = mid;
					else l = mid + 1;
				}
				nx.x += sec.x * ret,nx.y += sec.y * ret;
				fir.x += sec.x * ret,fir.y += sec.y * ret;
			}
		}
		k -= C(ans.y,ans.x);
		ans.x *= k,ans.y *= k;
		printf("%lld %lld\n",ans.x + 1,ans.y + 1);
	}
	return 0;
}