https://ac.nowcoder.com/acm/contest/5674/E
考虑x和y的每个共有的质因子在答案中的幂次q,ans*=p[k]^q
r[k]为p[k]在x中的幂次,s[k]为p[k]在y中的幂次
枚举i=a->b,然后i*r[k]<j*s[k],那么gcd中关于p[k]的幂次就是i*r[k],反之就是j*s[k]
由于q可能很大,注意幂次是要%phi(mod)也就是%(mod-1)
注意特判当x=mod,y=mod的情况
#include<bits/stdc++.h>
using namespace std;
const int mod=998244353,phi=mod-1;
int a,b,c,d;
int x,y;
int p[105],r[105],s[105],tot;
long long S(long long n){
return n*(n+1)/2;
}
long long qpow(long long x,int n){
long long ans=1;
for(;n;n>>=1,x=x*x%mod)if(n&1)ans=ans*x%mod;
return ans;
}
long long f(int k){
if(!s[k])return 1;
long long q=0;
for(int i=a;i<=b;i++){
int v=i*r[k]/s[k];
if(v<c){
q+=(long long)(d-c+1)*i%phi*r[k]%phi;
q%=phi;
}
else if(v<=d){
q+=(long long)(d-v)*i%phi*r[k]%phi;
q+=(S(v)-S(c-1))%phi*s[k]%phi;
q%=phi;
}
else{
q+=(S(d)-S(c-1))%phi*(b-i+1)%phi*s[k]%phi;
q%=phi;
break;
}
}
return qpow(p[k],q);
}
int main(){
scanf("%d%d%d%d%d%d",&a,&b,&c,&d,&x,&y);
if(x==mod&&y==mod){
if(b==0||d==0)printf("1");
else printf("0");
return 0;
}
if(x>y)swap(x,y);
for(int i=2;i*i<=x;i++){
if(x%i==0){
tot++;
p[tot]=i;
while(x%i==0)r[tot]++,x/=i;
while(y%i==0)s[tot]++,y/=i;
}
}
if(x!=1){
tot++;
p[tot]=x;
r[tot]=1;
while(y%x==0)s[tot]++,y/=x;
}
long long ans=1;
for(int i=1;i<=tot;i++)ans=ans*f(i)%mod;
printf("%lld",ans);
return 0;
}