2020暑期牛客多校训练营第九场（E）Groundhog Chasing Death（质因数分解，费马小定理）

Groundhog Chasing Death

原题请看这里

题目描述:

众所周知，“土拨鼠追逐死亡”的意思是“$GCD$”，而“$GCD$”代表“最大公约数”。
因此，您需要计算$\prod _{i=a}^b\prod _{j=c}^d\gcd (x^i，y^j)$并对$998244353$取模

输入描述：

一行包含六个整数$a，b，c，d，x，y$ $(0⩽a,b,c,d⩽3\times 10^6,0<x,y⩽10^9,a⩽b,c⩽d.)$

输出描述：

一行包含$\prod _{i=a}^b\prod _{j=c}^d\gcd (x^i，y^j)$ $mod$ $998244353$

样例：

样例输入1：

1 2 1 2 8 4

样例输出1：

2048

样例输入2：

1 2 3 4 120 180

样例输出2：

235140177

思路：

首先看看数据范围（正所谓看数据猜算法 ），发现$O(N^2)$都过不了，只有$O(NlogN)，O(N)$能过。
我们令
$x^i=p_1^{a1}\ast p_2^{a2}\ast p_3^{a3}\ast ......p_n^{an}$
$y^j=q_1^{b1}\ast q_2^{b2}\ast q_3^{b3}\ast ......q_m^{bm}$
我们发现，$x^i$和$y^j$的最大公因数一开始很快的增长，后来便不变了。假设$i<j$，因为当$x,y$中指数较小的一个和另一个数两两匹配完之后就只剩下指数较大的那个数了，也就是$y^{j-i}$。
然后我们发现，前$i-1$个$gcd$的值依次递增，所以可以用等差数列来求，后$j-i+1$个$y$和$x$的最大公因数相同，可以直接计算，如果最后$x>1$，说明$x$和$y$可能还有质因数，再做一遍即可。
时间复杂度$O(N)$

$AC$ $Code$:

#include
#define ll long long
using namespace std;
const int mod=998244353;
ll a,b,c,d,x,y,z,ans=1,n,u,v,md=mod-1;
ll ksm(ll a,ll b){
	ll ret=1;
	while(b){
		if(b&1)ret=ret*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ret;
}
ll get(ll x,ll y){
	ll k=0,p;
	for(int i=1;i<=x;i++){
		p=u*i/v;
		p=min(p,y);
		k=(k+(p+1ll)*p/2%md*v)%md;
		k=(k+1ll*i*(y-p)%md*u)%md;
	}
	return k;
}
int main(){
	scanf("%lld%lld%lld%lld%lld%lld",&a,&b,&c,&d,&x,&y);
	n=max(x,y);
	for(int i=2;i*i<=n;i++){
		while(x%i==0) u++,x/=i;
		while(y%i==0) v++,y/=i;
		if(u&&v){
			z=(2ll*md+get(b,d)+get(a-1,c-1)-get(a-1,d)-get(b,c-1))%md;
			ans=1ll*ans*ksm(i,z)%mod;
		}
		u=v=0;
	}
	if(x^0&&x==y){
		u=v=1;
		z=(2ll*md+get(b,d)+get(a-1,c-1)-get(a-1,d)-get(b,c-1))%md;
		ans=1ll*ans*ksm(x,z)%mod;
	}
	printf("%lld",ans);
}