gmoj 6844. 【2020.11.02提高组模拟】旅途和生活

题目

https://gmoj.net/senior/#main/show/6844

前言

这道题居然在题面中给出了结论（如果不给结论，该是T4难度的吧）
本来是不打算写题解的，但是由于这题给出的结论的证明很有价值，我就写下来吧。
SLS大佬太强啦%%%

题解

考场上原本是没什么思路的。但是我看到了这句话：

而上表说3,4,5,6保证a和b都是奇数。
我们可以先把a和b都写成$2^{x}y$的形式，使得它们中至少有一个变成了奇数。接下来分类讨论：

1. 当a和b都是奇数时，根据题目给出的结论，可以只存下$a^n$和$b^n$最后的63个二进制位。因此用unsigned long long搞快速幂就好了；
2. 当a和b一奇一偶时，$\mathrm{lowbit}(a^n-b^n)$显然是1。

最后把$\mathrm{lowbit}$的结果乘上$2^{nx}$就好了。

结论

当$a,b$为奇数时，有$$\mathrm{lowbit}\left(a^n-b^n\right)=\{\begin{array}{ll}\mathrm{lowbit}(a-b),& 2\nmid n\\ \mathrm{lowbit}(a-b)\cdot \mathrm{lowbit}(a+b)\cdot \mathrm{lowbit}\left(\frac{n}{2}\right),& 2\mid n\end{array}$$
由于$a,b,n$最大都只是$10^6$，因此$\mathrm{lowbit}(a^n-b^n)$是long long存得下的。

证明

先构造三个引理：

引理1：$\forall a,b\in N,\mathrm{lowbit}(ab)=\mathrm{lowbit}(a)\cdot \mathrm{lowbit}(b)$
证明：
令$a=2^{x_1}{x}_2,b=2^{y1}y2$，其中$x_2,y_2$都是奇数。
那么$$\begin{array}{rl}\mathrm{lowbit}(ab)& =\mathrm{lowbit}(2^{x_1}{x}_2\cdot 2^{y1}y2)\\ & =\mathrm{lowbit}(2^{x_1+y_1}\cdot x_2\cdot y_2)\\ & =2^{x_1+y_1}=\mathrm{lowbit}(a)\cdot \mathrm{lowbit}(b)\end{array}$$
得证。

引理2：当n是奇数时，$\mathrm{lowbit}(a^n-b^n)=\mathrm{lowbit}(a-b)$。
证明：
$$\mathrm{lowbit}(a^n-b^n)=\mathrm{lowbit}\left((a-b)\sum _{i=0}^{n-1}{a}^i\cdot b^{n-1-i}\right)$$
那么现在要证明：$$\sum _{i=0}^{n-1}{a}^i\cdot b^{n-1-i}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}2)$$
因为$a,b$都是奇数，因此$a^i,b^{n-1-i}$是奇数，$a^i\cdot b^{n-1-i}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}2)$。又因为有$n$项，因此${\sum }_{i=0}^{n-1}{a}^i\cdot b^{n-1-i}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}2)$。
得证。

引理3：当$a,b$是奇数时，$\mathrm{lowbit}(a^2+b^2)=2$。
证明：
分类讨论$a,b$在二进制下的后两位：

1. 当$a$的后两位是01，$b$的后两位是01时，$a^2,b^2$的后两位都是01，得证；
2. 当$a$的后两位是01，$b$的后两位是11时，$a^2,b^2$的后两位都是01，得证；
3. 当$a$的后两位是11，$b$的后两位是01时，$a^2,b^2$的后两位都是01，得证；
4. 当$a$的后两位是11，$b$的后两位是11时，$a^2,b^2$的后两位都是01，得证；

引理4：当$a,b$是奇数，$\mathrm{lowbit}\left(a^{2^x}-b^{2^x}\right)=\mathrm{lowbit}(a-b)\cdot \mathrm{lowbit}(a+b)\cdot 2^{x-1}$
证明：
令$t=x-1,p=a^2,q=b^2$，则要证$$\begin{array}{rl}\mathrm{lowbit}\left(p^{2^{x-1}}-q^{2^{x-1}}\right)& =\mathrm{lowbit}(p-q)\cdot 2^{x-1}\\ \mathrm{lowbit}\left(p^{2^t}-q^{2^t}\right)& =\mathrm{lowbit}(p-q)\cdot 2^t\end{array}$$
当$t=1$时，$\mathrm{lowbit}(p^2-q^2)=\mathrm{lowbit}(p-q)\cdot \mathrm{lowbit}(p+q)$。由引理3可得$\mathrm{lowbit}(p+q)=2$，因此等式成立。

在等式两边同时乘上$\mathrm{lowbit}(p^2+q^2)$，可得$\mathrm{lowbit}(p^4-q^4)=4\mathrm{lowbit}(p-q)$；
再在等式两边同时乘上$\mathrm{lowbit}(p^2+q^2)$，可得$\mathrm{lowbit}(p^8-q^8)=8\mathrm{lowbit}(p-q)$；
……
以此类推，$$\mathrm{lowbit}\left(p^{2^{x-1}}-q^{2^{x-1}}\right)=\mathrm{lowbit}(p-q)\cdot 2^{x-1}$$

得证。

---

1. 当n是奇数时，由引理2可得$\mathrm{lowbit}(a^n-b^n)=\mathrm{lowbit}(a-b)$；
2. 当n是偶数时，令$n=2^{x}y$（y是奇数），那么现在要证$$\mathrm{lowbit}\left(a^{2^{x}y}-b^{2^{x}y}\right)=\mathrm{lowbit}(a-b)\cdot \mathrm{lowbit}(a+b)\cdot 2^{x-1}$$
   由引理4得，$$\mathrm{lowbit}\left(a^{2^x}-b^{2^x}\right)=\mathrm{lowbit}(a-b)\cdot \mathrm{lowbit}(a+b)\cdot 2^{x-1}$$
   由引理2得，$$\mathrm{lowbit}\left(a^{2^{x}y}-b^{2^{x}y}\right)=\mathrm{lowbit}\left({\left(a^{2^x}\right)}^y-{\left(b^{2^x}\right)}^y\right)=\mathrm{lowbit}\left(a^{2^x}-b^{2^x}\right)$$
   因此$$\mathrm{lowbit}\left(a^n-b^n\right)=\mathrm{lowbit}(a-b)\cdot \mathrm{lowbit}(a+b)\cdot \mathrm{lowbit}\left(\frac{n}{2}\right)$$

得证。

---

CODE

#include<cstdio>
using namespace std;
typedef unsigned long long ll;
#define P 1000000007
inline void read(int &x)
{
	char ch;while(ch=getchar(),ch<'0'||ch>'9');x=ch-48;
	while(ch=getchar(),ch>='0'&&ch<='9') x=x*10+ch-48;
}
inline long long lowbit(long long x){return x&-x;}
inline int pow1(int x,int y)
{
	int s=1;
	while(y)
	{
		if(y&1) s=1LL*s*x%P;
		x=1LL*x*x%P,y>>=1;
	}
	return s;
}
inline ll pow2(ll x,int y)
{
	ll s=1;
	while(y)
	{
		if(y&1) s*=x;
		x*=x,y>>=1;
	}
	return s;
}
int main()
{
	freopen("journey.in","r",stdin);
	freopen("journey.out","w",stdout);
	int n,a,b,t,x,y;ll u,v;long long tmp;
	read(t);
	while(t--)
	{
		read(a),read(b),read(n);
		x=a&-a,y=b&-b;
		if(x>y) x=y;a/=x,b/=x;
		if((a&1)^(b&1)) printf("%d\n",pow1(x,n));
		else
		{
			u=pow2(a,n),v=pow2(b,n);
			u^=v;if(u&(1ll<<63)) u^=1ll<<63;tmp=u;
			printf("%d\n",lowbit(tmp)%P*pow1(x,n)%P);
		}
	}
	return 0;
}